The majority of stars are part of gravitationally bound stellar systems, such as binaries. Observations of protobinary systems constrain the conditions that lead to stellar multiplicity and ...subsequent orbital evolution. We report high-angular resolution observations of the circumbinary disk around BHB2007 11, a young binary protostar system. The two protostars are embedded in circumstellar disks that have radii of 2 to 3 astronomical units and probably contain a few Jupiter masses. These systems are surrounded by a complex structure of filaments connecting to the larger circumbinary disk. We also observe accretion and radio jets associated with the protobinary system. The accretion is preferentially onto the lower-mass protostar, consistent with theoretical predictions.
STATCONT is a python-based tool designed to determine the continuum emission level in spectral data, in particular for sources with a line-rich spectrum. The tool inspects the intensity distribution ...of a given spectrum and automatically determines the continuum level by using different statistical approaches. The different methods included in STATCONT are tested against synthetic data. We conclude that the sigma-clipping algorithm provides the most accurate continuum level determination, together with information on the uncertainty in its determination. This uncertainty can be used to correct the final continuum emission level, resulting in the here called ‘corrected sigma-clipping method’ or c-SCM. The c-SCM has been tested against more than 750 different synthetic spectra reproducing typical conditions found towards astronomical sources. The continuum level is determined with a discrepancy of less than 1% in 50% of the cases, and less than 5% in 90% of the cases, provided at least 10% of the channels are line free. The main products of STATCONT are the continuum emission level, together with a conservative value of its uncertainty, and datacubes containing only spectral line emission, i.e., continuum-subtracted datacubes. STATCONT also includes the option to estimate the spectral index, when different files covering different frequency ranges are provided.
ABSTRACT
The evolution of molecular clouds in galactic centres is thought to differ from that in galactic discs due to a significant influence of the external gravitational potential. We present a ...set of numerical simulations of molecular clouds orbiting on the 100-pc stream of the Central Molecular Zone (the central ${\sim }500\,{{\rm pc}}$ of the Galaxy) and characterize their morphological and kinematic evolution in response to the background potential and eccentric orbital motion. We find that the clouds are shaped by strong shear and torques, by tidal and geometric deformation, and by their passage through the orbital pericentre. Within our simulations, these mechanisms control cloud sizes, aspect ratios, position angles, filamentary structure, column densities, velocity dispersions, line-of-sight velocity gradients, spin angular momenta, and kinematic complexity. By comparing these predictions to observations of clouds on the Galactic Centre ‘dust ridge’, we find that our simulations naturally reproduce a broad range of key observed morphological and kinematic features, which can be explained in terms of well-understood physical mechanisms. We argue that the accretion of gas clouds on to the central regions of galaxies, where the rotation curve turns over and the tidal field is fully compressive, is accompanied by transformative dynamical changes to the clouds, leading to collapse and star formation. This can generate an evolutionary progression of cloud collapse with a common starting point, which either marks the time of accretion on to the tidally compressive region or of the most recent pericentre passage. Together, these processes may naturally produce the synchronized starbursts observed in numerous (extra)galactic nuclei.
Context. Stars are born deeply embedded in molecular clouds. In the earliest embedded phases, protostars emit the bulk of their radiation in the far-infrared wavelength range, where $Herschel$ is ...perfectly suited to probe at high angular resolution and dynamic range. In the high-mass regime, the birthplaces of protostars are thought to be in the high-density structures known as infrared-dark clouds (IRDCs). While massive IRDCs are believed to have the right conditions to give rise to massive stars and clusters, the evolutionary sequence of this process is not well-characterized.Aims: As part of the Earliest Phases of Star formation (EPoS) $Herschel$ guaranteed time key program, we isolate the embedded structures within IRDCs and other cold, massive molecular clouds. We present the full sample of 45 high-mass regions which were mapped at PACS 70, 100, and 160 μm and SPIRE 250, 350, and 500 μm. In the present paper, we characterize a population of cores which appear in the PACS bands and place them into context with their host molecular cloud and investigate their evolutionary stage.Methods: We construct spectral energy distributions (SEDs) of 496 cores which appear in all PACS bands, 34% of which lack counterparts at 24 μm. From single-temperature modified blackbody fits of the SEDs, we derive the temperature, luminosity, and mass of each core. These properties predominantly reflect the conditions in the cold, outer regions. Taking into account optical depth effects and performing simple radiative transfer models, we explore the origin of emission at PACS wavelengths. Results: The core population has a median temperature of 20 K and has masses and luminosities that span four to five orders of magnitude. Cores with a counterpart at 24 μm are warmer and bluer on average than cores without a 24 μm counterpart. We conclude that cores bright at 24 μm are on average more advanced in their evolution, where a central protostar(s) have heated the outer bulk of the core, than 24 μm-dark cores. The 24 μm emission itself can arise in instances where our line of sight aligns with an exposed part of the warm inner core. About 10% of the total cloud mass is found in a given cloud's core population. We uncover over 300 further candidate cores which are dark until 100 μm. These are possibly starless objects, and further observations will help us determine the nature of these very cold cores.
The temperature and density structure of molecular cloud cores are the most important physical quantities that determine the course of the protostellar collapse and the properties of the stars they ...firm. The instruments of the Herschel satellite provide us for the first time with both the spectral coverage and the spatial resolution that is needed to directly measure the dust temperature structure of nearby molecular cloud cores. With the aim of better constraining the initial physical conditions in molecular cloud cores at the onset of protostellar collapse, in particular of measuring their temperature structure, we initiated the guaranteed time key project (GTKP) "The Earliest Phases of Star Formation" (EPoS) with the Herschel satellite. This paper gives an overview of the low-mass sources in the EPoS project, the Herschel and complementary ground-based observations, our analysis method, and the initial results of the survey.
Context. Because of their short evolutionary time-scales, the earliest stages of high-mass star formation prior to the existence of any embedded heating source have barely been characterized until ...today. 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 100 μm wavelengths. The fragmentation and dynamical properties are investigated by means of Plateau de Bure Interferometer and Nobeyama 45 m single-dish spectral line and continuum observations. Results. The massive gas reservoir (between ~800 and ~1600 M⊙, depending on the assumed dust properties) fragments at spatial scales of ~18 000 AU 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.
Context. 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. Aims. We study the kinematics of the molecular gas in high-mass star-forming regions. Methods. We complemented the PACS and SPIRE far-infrared data of 16 high-mass star-forming regions from the Herschel key project EPoS with N sub(2)H super(+) molecular line data from the MOPRA and Nobeyama 45m telescope. Using the full N sub(2)H super(+) hyperfine structure, we produced column density, velocity, and linewidth maps. These were correlated with PACS 70 mu m images and PACS point sources. In addition, we searched for velocity gradients. Results. For several regions, the data suggest that the linewidth on the scale of clumps is dominated by outflows or unresolved velocity gradients. IRDC 18454 and G11.11 show two velocity components along several lines of sight. We find that all regions with a diameter larger than 1 pc 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. Conclusions. 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.
Context. Isolated starless cores within molecular clouds can be used as a testbed to investigate the conditions prior to the onset of fragmentation and gravitational proto-stellar collapse. Aims. We ...aim to determine the distribution of the dust temperature and the density of the starless core B68. Methods. In the framework of the Herschel guaranteed-time key programme “The Earliest Phases of Star formation” (EPoS), we have imaged B68 between 100 and 500 μm. Ancillary data at (sub)millimetre wavelengths, spectral line maps of the 12CO (2–1), and 13CO (2–1) transitions, as well as an NIR extinction map were added to the analysis. We employed a ray-tracing algorithm to derive the 2D mid-plane dust temperature and volume density distribution without suffering from the line-of-sight averaging effects of simple SED fitting procedures. Additional 3D radiative transfer calculations were employed to investigate the connection between the external irradiation and the peculiar crescent-shaped morphology found in the FIR maps. Results. For the first time, we spatially resolve the dust temperature and density distribution of B68, convolved to a beam size of 36.″4. We find a temperature gradient dropping from (16.7-1.0+1.3) K at the edge to (8.2-0.7+2.1) K in the centre, which is about 4 K lower than the result of the simple SED fitting approach. The column density peaks at NH = (4.3-2.8+1.4) × 1022 cm-2, and the central volume density was determined to nH = (3.4-2.5+0.9) × 105 cm-3. B68 has a mass of 3.1 M⊙ of material with AK > 0.2 mag for an assumed distance of 150 pc. We detect a compact source in the southeastern trunk, which is also seen in extinction and CO. At 100 and 160 μm, we observe a crescent of enhanced emission to the south. Conclusions. The dust temperature profile of B68 agrees well with previous estimates. We find the radial density distribution from the edge of the inner plateau outward to be nH ∝ r-3.5. Such a steep profile can arise from either or both of the following: external irradiation with a significant UV contribution or the fragmentation of filamentary structures. Our 3D radiative transfer model of an externally irradiated core by an anisotropic ISRF reproduces the crescent morphology seen at 100 and 160 μm. Our CO observations show that B68 is part of a chain of globules in both space and velocity, which may indicate that it was once part of a filament that dispersed. We also resolve a new compact source in the southeastern trunk and find that it is slightly shifted in centroid velocity from B68, lending qualitative support to core collision scenarios.
Aims. We aim to investigate the chemistry of internal photon-dominated regions (PDRs) surrounding deeply embedded hypercompact (HC) and ultracompact (UC) HII regions. We search for specific tracers ...of this evolutionary stage of massive star formation that can be detected with current astronomical facilities. Methods. We modeled hot cores with embedded HC/UCHII regions (called HII region models in the article despite the fact that we do not model the HII region itself), by coupling the astrochemical code Saptarsy to a radiative transfer framework obtaining the spatio-temporal evolution of abundances as well as time-dependent synthetic spectra. In these models where we focused on the internal PDR surrounding the HII region, the gas temperature is set to the dust temperature and we do not include dynamics thus the density structure is fixed. We compared this to hot molecular core (HMC) models and studied the effect on the chemistry of the radiation field which is included in the HII region models only during the computation of abundances. In addition, we investigated the chemical evolution of the gas surrounding HII regions with models of different densities at the ionization front, different sizes of the ionized cavity and different initial abundances. Results. We obtain the time evolution of synthetic spectra for a dozen of selected species as well as ratios of their integrated intensities. We find that some molecules such as C, N2H+, CN, and HCO do not trace the inner core and so are not good tracers to distinguish the HII/PDR regions to the HMCs phase. On the contrary, C+ and O trace the internal PDRs, in the two models starting with different initial abundances, but are unfortunately currently unobservable with the current achievable spatial resolution because of the very thin internal PDR (Δ rPDR < 100 AU). The emission of these two tracers is very dependent on the size of the HII region and on the density in the PDR. In addition, we find that the abundance profiles are highly affected by the choice of the initial abundances, hence the importance to properly define them.
Context. Stars form by the gravitational collapse of cold and dense molecular cloud cores. Constraining the temperature and density structure of such cores is fundamental for understanding the ...initial conditions of star formation. We use Herschel observations of the thermal far-infrared (FIR) dust emission from nearby and isolated molecular cloud cores and combine them with ground-based submillimeter continuum data to derive observational constraints on their temperature and density structure. Aims. The aim of this study is to verify the validity of a ray-tracing inversion technique developed to derive the dust temperature and density structure of nearby and isolated starless cores directly from the dust emission maps and to test if the resulting temperature and density profiles are consistent with physical models. Methods. We have developed a ray-tracing inversion technique that can be used to derive the temperature and density structure of starless cores directly from the observed dust emission maps without the need to make assumptions about the physical conditions. Using this ray-tracing inversion technique, we derive the dust temperature and density structure of six isolated starless molecular cloud cores from dust emission maps in the wavelengths range 100 mu m-1.2 mm. We then employ self-consistent radiative transfer modeling to the density profiles derived with the ray-tracing inversion method. In this model, the interstellar radiation field (ISRF) is the only heating source. The local strength of the ISRF as well as the total extinction provided by the outer envelope are treated as semi-free parameters which we scale within defined limits. The best-fit values of both parameters are derived by comparing the self-consistently calculated temperature profiles with those derived by the ray-tracing method. Results. We confirm earlier results and show that all starless cores are significantly colder inside than outside, with central core temperatures in the range 7.5?11.9 K and envelope temperatures that are 2.4 ? 9.6 K higher. The core temperatures show a strong negative correlation with peak column density which suggests that the thermal structure of the cores is dominated by external heating from the ISRF and shielding by dusty envelopes. We find that temperature profiles derived with the ray-tracing inversion method can be well-reproduced with self-consistent radiative transfer models if the cores have geometry that is not too complex and good data coverage with spatially resolved maps at five or more wavelengths in range between 100 mu m and 1.2 mm. We also confirm results from earlier studies that found that the usually adopted canonical value of the total strength of the ISRF in the solar neighbourhood is incompatible with the most widely used dust opacity models for dense cores. However, with the data available for this study, we cannot uniquely resolve the degeneracy between dust opacity law and strength of the ISRF.