Evolution of complex organic molecules in hot molecular cores: Synthetic spectra at (sub-)mm wavebands

Hot molecular cores (HMCs) are intermediate stages of high-mass star formation and are also known for their rich emission line spectra at (sub-)mm wavebands. The observed spectral feature of HMCs such as total number of emission lines and associated line intensities are also found to vary with evolutionary stages. We developed various 3D models for HMCs guided by the evolutionary scenarios proposed by recent empirical and modeling studies. We then investigated the spatio-temporal variation of temperature and molecular abundances in HMCs by consistently coupling gas-grain chemical evolution with radiative transfer calculations. We explored the effects of varying physical conditions on molecular abundances including density distribution and luminosity evolution of the central protostar(s). The time-dependent temperature structure of the hot core models provides a realistic framework for investigating the spatial variation of ice mantle evaporation as a function of evolutionary timescales. With increasing protostellar luminosity, the water ice evaporation font ($\sim$100K) expands and the spatial distribution of gas phase abundances of these COMs also spreads out. We simulated the synthetic spectra for these models at different evolutionary timescales to compare with observations. A qualitative comparison of the simulated and observed spectra suggests that these self-consistent hot core models can reproduce the notable trends in hot core spectral variation within the typical hot core timescales of 10$^{5}$ year. These models predict that the spatial distribution of various emission line maps will also expand with evolutionary time. The model predictions can be compared with high resolution observation that can probe scales of a few thousand AU in high-mass star forming regions such as from ALMA.[Abridged]Comment: accepted for publication in A&

Triggered/sequential star formation? A multi-phase ISM study around the prominent IRDC G18.93-0.03

G18.93-0.03 is a prominent dust complex within an 0.8deg long filament, with the molecular clump G18.93/m being IR dark from near IR wavelength up to 160mu. Spitzer composite images show an IR bubble spatially associated with G18.93. We use GRS 13CO and IRAM 30m H13CO+ data to disentangle the spatial structure of the region. From ATLASGAL submm data we calculate the gas mass, while we use the H13CO+ line width to estimate its virial mass. Using HERSCHEL data we produce temperature maps from fitting the SED. With the MAGPIS 20cm and SuperCOSMOS Halpha data we trace the ionized gas, and the VGPS HI survey provides information on the atomic hydrogen gas. We show that the bubble is spatially associated with G18.93, located at a kinematic near distance of 3.6kpc. With 280Msun, the most massive clump within G18.93 is G18.93/m. The virial analysis shows that it may be gravitationally bound and has neither Spitzer young stellar objects nor mid-IR point sources within. Fitting the SED reveals a temperature distribution that decreases towards its center, but heating from the ionizing source puts it above the general ISM temperature. We find that the bubble is filled by HII gas, ionized by an O8.5 star. Between the ionizing source and the IR dark clump G18.93/m we find a layered structure, from ionized to atomic to molecular hydrogen, revealing a PDR. Furthermore, we identify an additional velocity component within the bubble's 8mu emission rim at the edge of the infrared dark cloud and speculate that it might be shock induced by the expanding HII region. While the elevated temperature allows for the build-up of larger fragments, and the shock induced velocity component may lead to additional turbulent support, we do not find conclusive evidence that the massive clump G18.93/m is prone to collapse because of the expanding HII region.Comment: Accepted for publication in A&

The Earliest Phases of Star formation (EPoS): Temperature, density, and kinematic structure of the star-forming core CB 17

Context: The initial conditions for the gravitational collapse of molecular cloud cores and the subsequent birth of stars are still not well constrained. The characteristic cold temperatures (about 10 K) in such regions require observations at sub-millimetre and longer wavelengths. The Herschel Space Observatory and complementary ground-based observations presented in this paper have the unprecedented potential to reveal the structure and kinematics of a prototypical core region at the onset of stellar birth. Aims: This paper aims to determine the density, temperature, and velocity structure of the star-forming Bok globule CB 17. This isolated region is known to host (at least) two sources at different evolutionary stages: a dense core, SMM1, and a Class I protostar, IRS. Methods: We modeled the cold dust emission maps from 100 micron to 1.2 mm with both a modified blackbody technique to determine the optical depth-weighted line-of-sight temperature and column density and a ray-tracing technique to determine the core temperature and volume density structure. Furthermore, we analysed the kinematics of CB17 using the high-density gas tracer N2H+. Results: From the ray-tracing analysis, we find a temperature in the centre of SMM1 of 10.6 K, a flat density profile with radius 9500 au, and a central volume density of n(H) = 2.3x10^5 cm-3. The velocity structure of the N2H+ observations reveal global rotation with a velocity gradient of 4.3 km/s/pc. Superposed on this rotation signature we find a more complex velocity field, which may be indicative of differential motions within the dense core. Conclusions: SMM is a core in an early evolutionary stage at the verge of being bound, but the question of whether it is a starless or a protostellar core remains unanswered.Comment: published in A&

The physical and chemical structure of Sagittarius B2 VIIIa. Dust and ionized gas contributions to the full molecular line survey of 47 hot cores

Sagittarius B2 (Sgr B2) is a giant molecular cloud complex in the central molecular zone of our Galaxy hosting several sites of high-mass star formation. The two main centers of activity are Sgr B2(M) and Sgr B2(N), which contain 27 and 20 continuum sources, respectively. Our analysis aims to be a comprehensive modeling of each core spectrum, where we take the complex interaction between molecular lines, dust attenuation, and free-free emission arising from HII regions into account. In this work, we determine the dust and, if HII regions are contained, the parameters of the free-free thermal emission of the ionized gas for each core, and derive a self-consistent description of the continuum levels of each core. Using the high sensitivity of ALMA, we characterize the physical and chemical structure of these continuum sources and gain better insight into the star formation process within the cores. We used ALMA to perform an unbiased spectral line survey of all 47 sources in ALMA band 6 with a frequency coverage from 211 GHz to 275 GHz. In order to model the free-free continuum contribution of a specific core, we fit the contained recombination lines (RRLs) to obtain the electron temperatures and the emission measures, where we use an extended XCLASS program to describe RRLs and free-free continuum simultaneously. In contrast to previous analyses, we derived the corresponding parameters here not only for each core, but also for their local surrounding envelope, and determined their physical properties. The distribution of RRLs we found in the core spectra closely fits the distribution of HII regions described in previous analyses. For the cores we determine average dust temperatures of around 236 K (Sgr B2(M)) and 225 K (Sgr B2(N)), while the electronic temperatures are located in a range between 3800 K and 23800 K

Kinematic structure of massive star-forming regions - I. Accretion along filaments

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. We 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.Comment: accepted by A&

Fragmentation and dynamical collapse of the starless high-mass star-forming region IRDC18310-4

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.Comment: 11 pages, 10 figures, accepted to Astronomy and Astrophysics, high-resolution version with all figures included can be found at http://www.mpia.de/homes/beuther/papers.htm

Gas-phase CO depletion and N2H+ abundances in starless cores

Seven isolated, nearby low-mass starless molecular cloud cores have been observed as part of the Herschel key program Earliest Phases of Star formation (EPoS). By applying a ray-tracing technique to the obtained continuum emission and complementary (sub)mm emission maps, we derive the physical structure (density, dust temperature) of these cloud cores. We present observations of the 12CO, 13CO, and C18O (2-1) and N2H+ (1-0) transitions towards the same cores. Based on the density and temperature profiles, we apply time-dependent chemical and line-radiative transfer modeling and compare the modeled to the observed molecular emission profiles. CO is frozen onto the grains in the center of all cores in our sample. The level of CO depletion increases with hydrogen density and ranges from 46% up to more than 95% in the core centers in the core centers in the three cores with the highest hydrogen density. The average hydrogen density at which 50% of CO is frozen onto the grains is 1.1+-0.4 10^5 cm^-3. At about this density, the cores typically have the highest relative abundance of N2H+. The cores with higher central densities show depletion of N2H+ at levels of 13% to 55%. The chemical ages for the individual species are on average 2+-1 10^5 yr for 13CO, 6+-3 10^4 yr for C18O, and 9+-2 10^4 yr for N2H+. Chemical modeling indirectly suggests that the gas and dust temperatures decouple in the envelopes and that the dust grains are not yet significantly coagulated. We observationally confirm chemical models of CO-freezeout and nitrogen chemistry. We find clear correlations between the hydrogen density and CO depletion and the emergence of N2H+. The chemical ages indicate a core lifetime of less than 1 Myr.Comment: 24 pages, 25 figures, Accepted for publication in Astronomy and Astrophysic

Herschel observations of extraordinary sources: Full Herschel/HIFI molecular line survey of Sagittarius B2(M)

We present a full analysis of a broadband spectral line survey of Sagittarius B2 (Main), one of the most chemically rich regions in the Galaxy located within the giant molecular cloud complex Sgr B2 in the Central Molecular Zone. Our goal is to derive the molecular abundances and temperatures of the high-mass star-forming region Sgr B2(M) and thus its physical and astrochemical conditions. Sgr B2(M) was observed using the Heterodyne Instrument for the Far-Infrared (HIFI) on board the Herschel Space Observatory in a spectral line survey from 480 to 1907 GHz at a spectral resolution of 1.1 MHz, which provides one of the largest spectral coverages ever obtained toward this high-mass star-forming region in the submillimeter with high spectral resolution and includes frequencies > 1 THz unobservable from the ground. We model the molecular emission from the submillimeter to the far-IR using the XCLASS program. For each molecule, a quantitative description was determined taking all emission and absorption features of that species across the entire spectral range into account. Additionally, we derive velocity resolved ortho / para ratios for those molecules for which ortho and para resolved molecular parameters are available. Finally, the temperature and velocity distributions are analyzed and the derived abundances are compared with those obtained for Sgr B2(N) from a similar HIFI survey. A total of 92 isotopologues were identified, arising from 49 different molecules, ranging from free ions to complex organic compounds and originating from a variety of environments from the cold envelope to hot and dense gas within the cores. Sulfur dioxide, methanol, and water are the dominant contributors. For the ortho / para ratios we find deviations from the high temperature values between 13 and 27 %. In total 14 % of all lines remain unidentified.Comment: 67 pages, 102 figures, submitted to A&