Massive stars and their associated ionized (HII) regions could play a key role in the formation and evolution of filaments that host star formation. However, the properties of filaments that interact with H regions are still poorly known. To investigate the impact of HII regions on the formation of filaments, we imaged the Galactic HII region RCW 120 and its surroundings where active star formation takes place and where the role of ionization feedback on the star formation process has already been studied. We used the ArT\'eMiS camera on the APEX telescope and combined the ArT\'eMiS data at 350 and 450 microns with Herschel-SPIRE/HOBYS. We studied the dense gas distribution around RCW 120 with a resolution of 8 arcsec (0.05 pc at a distance of 1.34 kpc). Our study allows us to trace the median radial intensity profile of the dense shell of RCW 120. This profile is asymmetric, indicating a clear compression from the HII region on the inner part of the shell. The profile is observed to be similarly asymmetric on both lateral sides of the shell, indicating a homogeneous compression over the surface. On the contrary, the profile analysis of a radial filament associated with the shell, but located outside of it, reveals a symmetric profile, suggesting that the compression from the ionized region is limited to the dense shell. The mean intensity profile of the internal part of the shell is well fitted by a Plummer like profile with a deconvolved Gaussian FWHM of 0.09 pc, as observed for filaments in low-mass star-forming regions. This study suggests that compression exerted by HII regions may play a key role in the formation of filaments and may further act on their hosted star formation. ArT\'eMiS data also suggest that RCW 120 might be a 3D ring, rather than a spherical structure

Herschel observations of nearby clouds in the Gould Belt support a paradigm for low-mass star formation, starting with the generation of molecular filaments, followed by filament fragmentation, and the concentration of mass into self-gravitating prestellar cores. With the unique far-infrared and submillimeter continuum imaging capabilities of the Herschel Space observatory, the closeby (d = 139 pc) Ophiuchus cloud was mapped at five wavelengths from 70 ${\mu}m$ to 500 ${\mu}m$ with the aim of providing a complete census of dense cores in this region, including unbound starless cores, bound prestellar cores, and protostellar cores. Taking advantage of the high dynamic range and multi-wavelength nature of the Herschel data, we used the multi-scale decomposition algorithms getsources and getfilaments to identify a complete sample of dense cores and filaments in the cloud and study their properties. The densest clouds of the Ophiuchus complex, L1688 and L1689, which thus far are only indirectly described as filamentary regions owing to the spatial distribution of their young stellar objects (YSOs), are confirmed to be dominated by filamentary structures. The tight correlation observed between prestellar cores and filamentary structures in L1688 and L1689 supports the view that solar-type star formation occurs primarily in dense filaments. While the sub clouds of the complex show disparities, L1689 being less efficient than L1688 at forming stars when considering their total mass budgets, both sub clouds share almost the same prestellar core formation efficiency in dense molecular gas. We also find evidence in the Herschel data for a remarkable concentric geometrical configuration in L1688 which is dominated by up to three arc-like compression fronts and presumably created by shockwave events emanating from the Sco OB2 association, including the neighboring massive (O9V) star ${\sigma}$ Sco.

We present a detailed study of the Orion B clouds (d~400 pc), imaged with the PACS/SPIRE cameras at 70-500 $\mu$m by the Herschel Gould Belt survey (HGBS). We release new high-res. maps of column density and dust temperature. In the filamentary sub-regions NGC2023/2024, NGC2068/2071, and L1622, 1768 starless dense cores were identified, ~28-45% of which are self-gravitating prestellar cores. A total of 76 protostellar dense cores were also found. The typical lifetime of the prestellar cores was found to be $t_{\rm pre}=1.7_{-0.6}^{+0.8}$ Myr. The prestellar core mass function (CMF) peaks at ~0.5 $M_\odot$ and is consistent with a power law with log slope -1.27$\pm$0.24 at the high-mass end, compared to the Salpeter slope of -1.35. In this region, we confirm the existence of a transition in prestellar core formation efficiency (CFE) around a fiducial value A_V_bg~7 mag in background visual extinction, similar to the trend observed with Herschel in other clouds. This is not a sharp threshold, but a smooth transition between a regime with very low prestellar CFE at A_V_bg<5 and a regime with higher, roughly constant CFE at A_V_bg$\gtrsim$10. The total mass in the form of prestellar cores represents only ~20% of the dense molecular cloud gas at A_V_bg$\gtrsim$7 mag. About 60-80% of the prestellar cores are closely associated with filaments, and this fraction increases up to >90% when a more complete sample of filamentary structures is considered. Interestingly, the median separation between nearest core neighbors corresponds to the typical inner filament width of ~0.1 pc commonly observed in nearby molecular clouds. Analysis of the CMF observed as a function of background cloud column density shows that the most massive prestellar cores are spatially segregated in the highest column density areas, and suggests that both higher- and lower-mass prestellar cores may form in denser filaments.

Herschel surveys of Galactic clouds support a paradigm for low-mass star formation in which dense filaments play a crucial role. The detailed fragmentation properties of star-forming filaments remain poorly understood, however, and the validity of the filament paradigm in the high-mass regime is still unclear. To investigate the density/velocity structure of the filament in the high-mass star-forming region NGC6334, we conducted ALMA observations in the 3mm continuum and the N2H+(1-0) line at ~3arcsec resolution. The filament was detected in both tracers. We identified 26 cores at 3mm and 5 velocity-coherent fiber-like features in N2H+ within the filament. The typical length of, and velocity difference between, the fiber-like features of the NGC6334 filament are reminiscent of the properties for the fibers of the low-mass star-forming filament B211/B213. Only 2 or 3 of the 5 velocity-coherent features are well aligned with the filament and may represent genuine, fiber sub-structures. The core mass distribution has a peak at ~10Msun. They can be divided into 7 groups of cores, closely associated with ArTeMiS clumps. The projected separation between cores and the projected spacing between clumps are roughly consistent with the effective Jeans length in the filament and a physical scale of about 4 times the filament width, respectively, suggesting a bimodal filament fragmentation process. Despite being one order of magnitude denser and more massive than the B211/B213 filament, the NGC6334 filament has a similar density/velocity structure. The difference is that the cores in NGC6334 appear to be an order of magnitude denser and more massive than the cores in Taurus. This suggests that dense filaments may evolve and fragment in a similar manner in low- and high-mass star-forming regions, and that the filament paradigm may hold in the intermediate-mass (if not high-mass) star formation regime.

Understanding the processes that determine the stellar Initial Mass Function (IMF) is a critical unsolved problem, with profound implications for many areas of astrophysics. In molecular clouds, stars are formed in cores, gas condensations which are sufficiently dense that gravitational collapse converts a large fraction of their mass into a star or small clutch of stars. In nearby star-formation regions, the core mass function (CMF) is strikingly similar to the IMF, suggesting that the shape of the IMF may simply be inherited from the CMF. Here we present 1.3 mm observations, obtained with ALMA, the world's largest interferometer, of the active star-formation region W43-MM1, which may be more representative of the Galactic-disk regions where most stars form. The unprecedented resolution of these observations reveals, for the first time, a statistically robust CMF at high masses, with a slope that is markedly shallower than the IMF. This seriously challenges our understanding of the origin of the IMF.

We present far-infrared observations of Monoceros R2 (a giant molecular cloud at approximately 830 pc distance, containing several sites of active star formation), as observed at 70 {\mu}m, 160 {\mu}m, 250 {\mu}m, 350 {\mu}m, and 500 {\mu}m by the Photodetector Array Camera and Spectrometer (PACS) and Spectral and Photometric Imaging Receiver (SPIRE) instruments on the Herschel Space Observatory as part of the Herschel imaging survey of OB young stellar objects (HOBYS) Key programme. The Herschel data are complemented by SCUBA-2 data in the submillimetre range, and WISE and Spitzer data in the mid-infrared. In addition, C18O data from the IRAM 30-m Telescope are presented, and used for kinematic information. Sources were extracted from the maps with getsources, and from the fluxes measured, spectral energy distributions were constructed, allowing measurements of source mass and dust temperature. Of 177 Herschel sources robustly detected in the region (a detection with high signal-to-noise and low axis ratio at multiple wavelengths), including protostars and starless cores, 29 are found in a filamentary hub at the centre of the region (a little over 1% of the observed area). These objects are on average smaller, more massive, and more luminous than those in the surrounding regions (which together suggest that they are at a later stage of evolution), a result that cannot be explained entirely by selection effects. These results suggest a picture in which the hub may have begun star formation at a point significantly earlier than the outer regions, possibly forming as a result of feedback from earlier star formation. Furthermore, the hub may be sustaining its star formation by accreting material from the surrounding filaments.

Modern high-resolution images obtained with space observatories display extremely strong intensity variations across images on all spatial scales. Source extraction in such images with methods based on global thresholding may bring unacceptably large numbers of spurious sources in bright areas while failing to detect sources in low-background or low-noise areas. It would be highly beneficial to subtract background and equalize the levels of small-scale fluctuations in the images before extracting sources or filaments. This paper describes getimages, a new method of background derivation and image flattening. It is based on median filtering with sliding windows that correspond to a range of spatial scales from the observational beam size up to a maximum structure width $X_{\lambda}$. The latter is a single free parameter of getimages that can be evaluated manually from the observed image $\mathcal{I}_{\lambda}$. The median filtering algorithm provides a background image $\tilde{\mathcal{B}}_{\lambda}$ for structures of all widths below $X_{\lambda}$. The same median filtering procedure applied to an image of standard deviations $\mathcal{D}_{\lambda}$ derived from a background-subtracted image $\tilde{\mathcal{S}}_{\lambda}$ results in a flattening image $\tilde{\mathcal{F}}_{\lambda}$. Finally, a flattened detection image $\mathcal{I}_{{\lambda}\mathrm{D}}{\,=\,}\tilde{\mathcal{S}}_{\lambda}{/}\tilde{\mathcal{F}}_{\lambda}$ is computed, whose standard deviations are uniform outside sources and filaments. Detecting sources in such greatly simplified images results in much cleaner extractions that are more complete and reliable. As a bonus, getimages reduces various observational and map-making artifacts and equalizes noise levels between independent tiles of mosaicked images.

Filaments are key for star formation models. As part of the study carried out by the Herschel GCC Programme, here we study the filament properties presented in GCC.VII in context with theoretical models of filament formation and evolution. A conservative sample of filaments at a distance D<500pc was extracted with the Getfilaments algorithm. Their physical structure was quantified according to two main components: the central (Gaussian) region (core component), and the power-law like region dominating the filament column density profile at larger radii (wing component). The properties and behaviour of these components relative to the total linear mass density of the filament and its environmental column density were compared with theoretical models describing the evolution of filaments under gravity-dominated conditions. The feasibility of a transition to supercritical state by accretion is dependent on the combined effect of filament intrinsic properties and environmental conditions. Reasonably self-gravitating (high Mline-core) filaments in dense environments (av\sim3mag) can become supercritical in timescales of t\sim1Myr by accreting mass at constant or decreasing width. The trend of increasing Mline-tot (Mline-core and Mline-wing), and ridge Av with background also indicates that the precursors of star-forming filaments evolve coevally with their environment. The simultaneous increase of environment and filament Av explains the association between dense environments and high Mline-core values, and argues against filaments remaining in constant single-pressure equilibrium states. The simultaneous growth of filament and background in locations with efficient mass assembly, predicted in numerical models of collapsing clouds, presents a suitable scenario for the fulfillment of the combined filament mass-environment criterium that is in quantitative agreement with Herschel observations.

Triggered star formation around HII regions could be an important process. The Galactic HII region RCW 79 is a prototypical object for triggered high-mass star formation. We take advantage of Herschel data from the surveys HOBYS, "Evolution of Interstellar Dust", and Hi-Gal to extract compact sources in this region, complemented with archival 2MASS, Spitzer, and WISE data to determine the physical parameters of the sources (e.g., envelope mass, dust temperature, and luminosity) by fitting the spectral energy distribution. We obtained a sample of 50 compact sources, 96% of which are situated in the ionization-compressed layer of cold and dense gas that is characterized by the column density PDF with a double-peaked lognormal distribution. The 50 sources have sizes of 0.1-0.4 pc with a typical value of 0.2 pc, temperatures of 11-26 K, envelope masses of 6-760 $M_\odot$, densities of 0.1-44 $\times$ $10^5$ cm$^{-3}$, and luminosities of 19-12712 $L_\odot$. The sources are classified into 16 class 0, 19 intermediate, and 15 class I objects. Their distribution follows the evolutionary tracks in the diagram of bolometric luminosity versus envelope mass (Lbol-Menv) well. A mass threshold of 140 $M_\odot$, determined from the Lbol-Menv diagram, yields 12 candidate massive dense cores that may form high-mass stars. The core formation efficiency (CFE) for the 8 massive condensations shows an increasing trend of the CFE with density. This suggests that the denser the condensation, the higher the fraction of its mass transformation into dense cores, as previously observed in other high-mass star-forming regions.

The expansion of HII regions can trigger the formation of stars. An overdensity of young stellar objects (YSOs) is observed at the edges of HII regions but the mechanisms that give rise to this phenomenon are not clearly identified. Moreover, it is difficult to establish a causal link between HII-region expansion and the star formation observed at the edges of these regions. A clear age gradient observed in the spatial distribution of young sources in the surrounding might be a strong argument in favor of triggering. We have observed the Galactic HII region RCW120 with herschel PACS and SPIRE photometers at 70, 100, 160, 250, 350 and 500$\mu$m. We produced temperature and H$_2$ column density maps and use the getsources algorithm to detect compact sources and measure their fluxes at herschel wavelengths. We have complemented these fluxes with existing infrared data. Fitting their spectral energy distributions (SEDs) with a modified blackbody model, we derived their envelope dust temperature and envelope mass. We computed their bolometric luminosities and discuss their evolutionary stages. The herschel data, with their unique sampling of the far infrared domain, have allowed us to characterize the properties of compact sources observed towards RCW120 for the first time. We have also been able to determine the envelope temperature, envelope mass and evolutionary stage of these sources. Using these properties we have shown that the density of the condensations that host star formation is a key parameter of the star-formation history, irrespective of their projected distance to the ionizing stars.

Fitting spectral distributions of total fluxes or image intensities are two standard methods for estimating the masses of starless cores and protostellar envelopes. These mass estimates, which are the main source and basis of our knowledge of the origin and evolution of self-gravitating cores and protostars, are uncertain. In this model-based study, a grid of radiative transfer models of starless cores and protostellar envelopes was computed and their total fluxes and image intensities were fitted to derive the model masses. To investigate intrinsic effects related to the physical objects, all observational complications were explicitly ignored. Known true values of the numerical models allow us to assess the qualities of the methods and fitting models, as well as the effects of nonuniform temperatures, far-infrared opacity slope, selected subsets of wavelengths, background subtraction, and angular resolutions. The method of fitting image intensities gives more accurate masses for more resolved objects than the method of fitting total fluxes. With the latter, a fitting model that assumes optically thin emission gives much better results than the one allowing substantial optical depths. Temperature excesses within the objects above the mass-averaged values skew their spectral shapes towards shorter wavelengths, leading to masses underestimated typically by factors 2-5. With a fixed opacity slope deviating from the true value by a factor of 1.2, derived masses are inaccurate within a factor of 2. The most accurate masses of the models are estimated by fitting just two or three longest-wavelength fluxes. Conventional algorithm of background subtraction is a likely source of large systematic errors. The absolute values of masses of the unresolved or poorly resolved objects in star-forming regions are uncertain to within at least a factor of 2-3.

The association of filaments with protostellar objects has made these structures a priority target in star formation studies. The datasets of the Herschel Galactic Cold Cores Key Programme allow for a statistical study of filaments with a wide range of intrinsic and environmental characteristics. Characterisation of this sample can be used to identify key physical parameters and quantify the role of environment in the formation of supercritical filaments. Filaments were extracted from fields at D<500pc with the getfilaments algorithm and characterised according to their column density profiles and intrinsic properties. Each profile was fitted with a beam-convolved Plummer-like function and quantified based on the relative contributions from the filament 'core', represented by a Gaussian, and 'wing' component, dominated by the power-law of the Plummer-like function. These parameters were examined for populations associated with different background levels. We find that filaments increase their core (Mcore) and wing (Mwing) contributions while increasing their total linear mass density (Mtot). Both components appear to be linked to the local environment, with filaments in higher backgrounds having systematically more massive Mcore and Mwing. This dependence on the environment supports an accretion-based model for filament evolution in the local neighbourhood (D<500pc). Structures located in the highest backgrounds develop the highest central Av, Mcore, and Mwing as Mtot increases with time, favoured by the local availability of material and the enhanced gravitational potential. Our results indicate that filaments acquiring a significantly massive central region with Mcore>Mcrit/2 may become supercritical and form stars. This translates into a need for filaments to become at least moderately self-gravitating in order to undergo localised star formation or become star-forming filaments.

We present a catalogue of dense cores in a $\sim 4^\circ\times2^\circ$ field of the Taurus star-forming region, inclusive of the L1495 cloud, derived from Herschel SPIRE and PACS observations in the 70 $\mu$m, 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, and 500 $\mu$m continuum bands. Estimates of mean dust temperature and total mass are derived using modified blackbody fits to the spectral energy distributions. We detect 525 starless cores of which $\sim10$-20% are gravitationally bound and therefore presumably prestellar. Our census of unbound objects is $\sim85$% complete for $M>0.015\,M_\odot$ in low density regions ($A_V\stackrel{<}{_\sim}5$ mag), while the bound (prestellar) subset is $\sim85$% complete for $M>0.1\,M_\odot$ overall. The prestellar core mass function (CMF) is consistent with lognormal form, resembling the stellar system initial mass function, as has been reported previously. All of the inferred prestellar cores lie on filamentary structures whose column densities exceed the expected threshold for filamentary collapse, in agreement with previous reports. Unlike the prestellar CMF, the unbound starless CMF is not lognormal, but instead is consistent with a power-law form below $0.3\,M_\odot$ and shows no evidence for a low-mass turnover. It resembles previously reported mass distributions for CO clumps at low masses ($M\stackrel{<}{_\sim}0.3\,M_\odot$). The volume density PDF, however, is accurately lognormal except at high densities. It is consistent with the effects of self-gravity on magnetized supersonic turbulence. The only significant deviation from lognormality is a high-density tail which can be attributed unambiguously to prestellar cores.

We present a detailed study of V899 Mon (a new member in the FUors/EXors family of young low-mass stars undergoing outburst), based on our long-term monitoring of the source starting from November 2009 to April 2015. Our optical and near-infrared photometric and spectroscopic monitoring recorded the source transitioning from its first outburst to a short duration quiescence phase ($<$ 1 year), and then returning to a second outburst. We report here the evolution of the outflows from inner region of the disk as the accretion rate evolved in various epochs. Our high resolution (R$\sim$37000) optical spectrum could resolve interesting clumpy structures in the outflow traced by various lines. Change in far-infrared flux was also detected between two outburst epochs. Based on our observations we constrained various stellar and envelope parameters of V899 Mon, as well as the kinematics of its accretion and outflow. The photometric and spectroscopic properties of this source fall between classical FUors and EXors. Our investigation of V899 Mon hints instability associated with magnetospheric accretion to be the physical cause of sudden short duration pause of outburst in 2011. It is also a good candidate to explain similar short duration pauses in outburst of some other FUors/EXors sources.

We present and discuss the results of the Herschel Gould Belt survey observations in a ~11 deg^2 area of the Aquila molecular cloud complex at d~260 pc, imaged with the SPIRE/PACS cameras from 70 to 500 micron. We identify a complete sample of starless dense cores and embedded protostars in this region, and analyze their global properties and spatial distributions. We find a total of 651 starless cores, ~60% of which are gravitationally bound prestellar cores, and they will likely form stars in the future. We also detect 58 protostellar cores. The core mass function (CMF) derived for the prestellar cores is very similar in shape to the stellar initial mass function (IMF), supporting the earlier view that there is a close physical link between the IMF and the CMF. The global shift in mass scale observed between the CMF and the IMF is consistent with a typical star formation efficiency of ~40%. By comparing the numbers of starless cores to the number of young stellar objects, we estimate that the lifetime of prestellar cores is ~1 Myr. We find a strong correlation between the spatial distribution of prestellar cores and the densest filaments. About 90% of the Herschel-identified prestellar cores are located above a background column density corresponding to A_V~7, and ~75% of them lie within filamentary structures with supercritical masses per unit length >~16 M_sun/pc. These findings support a picture wherein the cores making up the peak of the CMF (and probably responsible for the base of the IMF) result primarily from the gravitational fragmentation of marginally supercritical filaments. Given that filaments appear to dominate the mass budget of dense gas at A_V>7, our findings also suggest that the physics of prestellar core formation within filaments is responsible for a characteristic "efficiency" SFR/M_dense ~5+-2 x 10^-8 yr^-1 for the star formation process in dense gas.

The W3 GMC is a prime target for investigating the formation of high-mass stars and clusters. This second study of W3 within the HOBYS Key Program provides a comparative analysis of subfields within W3 to further constrain the processes leading to the observed structures and stellar population. Probability density functions (PDFs) and cumulative mass distributions (CMDs) were created from dust column density maps, quantified as extinction Av. The shape of the PDF, typically represented with a lognormal function at low Av "breaking" to a power-law tail at high Av, is influenced by various processes including turbulence and self-gravity. The breaks can also be identified, often more readily, in the CMDs. The PDF break from lognormal (Av(SF)= 6-10 mag) appears to shift to higher Av by stellar feedback, so that high-mass star-forming regions tend to have higher PDF breaks. A second break at Av > 50 mag traces structures formed or influenced by a dynamic process. Because such a process has been suggested to drive high-mass star formation in W3, this second break might then identify regions with potential for hosting high-mass stars/clusters. Stellar feedback appears to be a major mechanism driving the local evolution and state of regions within W3. A high initial star formation efficiency in a dense medium could result in a self-enhancing process, leading to more compression and favourable star-formation conditions (e.g., colliding flows), a richer stellar content, and massive stars. This scenario would be compatible with the "convergent constructive feedback" model introduced in our previous Herschel study.

We use PACS and SPIRE continuum data at 160 um, 250 um, 350 um, and 500 um from the Herschel Gould Belt Survey to sample seven clumps in Perseus: B1, B1-E, B5, IC348, L1448, L1455, and NGC1333. Additionally, we identify and characterize the embedded Class 0 protostars using detections of compact Herschel sources at 70 um as well as archival Spitzer catalogues and SCUBA 850 um photometric data. We identify 28 candidate Class 0 protostars, four of which are newly discovered sources not identified with Spitzer. We find that the star formation efficiency of clumps, as traced by Class 0 protostars, correlates strongly with the flatness of their respective column density distributions at high values. This correlation suggests that the fraction of high column density material in a clump reflects only its youngest protostellar population rather than its entire source population. We propose that feedback from either the formation or evolution of protostars changes the local density structure of clumps.

The Chamaeleon molecular cloud complex is one of the nearest star-forming sites encompassing three molecular clouds with a different star-formation history, from quiescent (Cha III) to actively forming stars (Cha II), and reaching the end of star-formation (Cha I). To charactize its large-scale structure, we derived column density and temperature maps using PACS and SPIRE observations from the Herschel Gould Belt Survey, and applied several tools, such as filament tracing, power-spectra, \Delta-variance, and probability distribution functions of column density (PDFs), to derive physical properties. The column density maps reveal a different morphological appearance for the three clouds, with a ridge-like structure for Cha I, a clump-dominated regime for Cha II, and an intricate filamentary network for Cha III. The filament width is measured to be around 0.12\pm0.04 pc in the three clouds, and the filaments found to be gravitationally unstable in Cha I and II, but mostly subcritical in Cha III. Faint filaments (striations) are prominent in Cha I showing a preferred alignment with the large-scale magnetic field. The PDFs of all regions show a lognormal distribution at low column densities. For higher densities, the PDF of Cha I shows a turnover indicative of an extended higher density component, culminating with a power-law tail. Cha II shows a power-law tail with a slope characteristic of gravity. The PDF of Cha III can be best fit by a single lognormal. The turbulence properties of the three regions are found to be similar, pointing towards a scenario where the clouds are impacted by large-scale processes. The magnetic field could possibly play an important role for the star-formation efficiency in the Chamaeleon clouds if proven that it can effectively channel material on Cha I, and possibly Cha II, but probably less efficiently on the quiescent Cha III cloud.

Far-infrared imaging surveys of Galactic star-forming regions with Herschel have shown that a substantial part of the cold interstellar medium appears as a fascinating web of omnipresent filamentary structures. This highly anisotropic ingredient of the interstellar material further complicates the difficult problem of the systematic detection and measurement of dense cores in the strongly variable but (relatively) isotropic backgrounds. Observational evidence that stars form in dense filaments creates severe problems for automated source extraction methods that must reliably distinguish sources not only from fluctuating backgrounds and noise, but also from the filamentary structures. A previous paper presented the multi-scale, multi-wavelength source extraction method getsources based on a fine spatial scale decomposition and filtering of irrelevant scales from images. In this paper, a multi-scale, multi-wavelength filament extraction method getfilaments is presented that solves this problem, substantially improving the robustness of source extraction with getsources in filamentary backgrounds. The main difference is that the filaments extracted by getfilaments are now subtracted by getsources from detection images during source extraction, greatly reducing the chances of contaminating catalogs with spurious sources. The intimate physical relationship between forming stars and filaments seen in Herschel observations demands that accurate filament extraction methods must remove the contribution of sources and that accurate source extraction methods must be able to remove underlying filamentary structures. Source extraction with getsources now provides researchers also with clean images of filaments, free of sources, noise, and isotropic backgrounds.

We present the first overview of the Herschel observations of the nearby high-mass star-forming region NGC 7538, taken as part of the Herschel imaging study of OB Young Stellar objects (HOBYS) Key Programme. These PACS and SPIRE maps cover an approximate area of one square degree at five submillimeter and far-infrared wavebands. We have identified 780 dense sources and classified 224 of those. With the intention of investigating the existence of cold massive starless or class 0-like clumps that would have the potential to form intermediate- to high-mass stars, we further isolate 13 clumps as the most likely candidates for followup studies. These 13 clumps have masses in excess of 40 M_sun and temperatures below 15 K. They range in size from 0.4 pc to 2.5 pc and have densities between 3x10^3 cm^-3 to 4x10^4 cm^-3. Spectral energy distributions are then used to characterize their energetics and evolutionary state through a luminosity-mass diagram. NGC 7538 has a highly filamentary structure, previously unseen in the dust continuum of existing submillimeter surveys. We report the most complete imaging to date of a large, evacuated ring of material in NGC 7538 which is bordered by many cool sources.

The whole of the Taurus region (a total area of 52 sq. deg.) has been observed by the Herschel SPIRE and PACS instruments at wavelengths of 70, 160, 250, 350 and 500 {\mu}m as part of the Herschel Gould Belt Survey. In this paper we present the first results from the part of the Taurus region that includes the Barnard 18 and L1536 clouds. A new source-finding routine, the Cardiff Source-finding AlgoRithm (CSAR), is introduced, which is loosely based on CLUMPFIND, but that also generates a structure tree, or dendrogram, which can be used to interpret hierarchical clump structure in a complex region. Sources were extracted from the data using the hierarchical version of CSAR and plotted on a mass-size diagram. We found a hierarchy of objects with sizes in the range 0.024-2.7 pc. Previous studies showed that gravitationally bound prestellar cores and unbound starless clumps appeared in different places on the mass-size diagram. However, it was unclear whether this was due to a lack of instrumental dynamic range or whether they were actually two distinct populations. The excellent sensitivity of Herschel shows that our sources fill the gap in the mass-size plane between starless and pre-stellar cores, and gives the first clear supporting observational evidence for the theory that unbound clumps and (gravitationally bound) prestellar cores are all part of the same population, and hence presumably part of the same evolutionary sequence (c.f. Simpson et al. 2011).

We present first results from the Herschel Gould Belt survey for the B211/L1495 region in the Taurus molecular cloud. Thanks to their high sensitivity and dynamic range, the Herschel images reveal the structure of the dense, star-forming filament B211 with unprecedented detail, along with the presence of striations perpendicular to the filament and generally oriented along the magnetic field direction as traced by optical polarization vectors. Based on the column density and dust temperature maps derived from the Herschel data, we find that the radial density profile of the B211 filament approaches a power-law behavior {\rho} {\propto} r^(-2.0{\pm}0.4) at large radii and that the temperature profile exhibits a marked drop at small radii. The observed density and temperature profiles of the B211 filament are in good agreement with a theoretical model of a cylindrical filament undergoing gravitational contraction with a polytropic equation of state: P {\propto} {\rho}^{\gamma} and T {\propto} {\rho}^({\gamma}-1), with {\gamma}=0.97{\pm}0.01<1 (i.e. not strictly isothermal). The morphology of the column density map, where some of the perpendicular striations are apparently connected to the B211 filament, further suggests that the material may be accreting along the striations onto the main filament. The typical velocities expected for the infalling material in this picture are ~0.5-1 km/s, which are consistent with the existing kinematical constraints from previous CO observations.

In order to characterise the cloud structures responsible for the formation of high-mass stars, we present Herschel observations of the DR21 environment. Maps of the column density and dust temperature unveil the structure of the DR21 ridge and several connected filaments. The ridge has column densities larger than 1e23/cm^2 over a region of 2.3 pc^2. It shows substructured column density profiles and branching into two major filaments in the north. The masses in the studied filaments range between 130 and 1400 Msun whereas the mass in the ridge is 15000 Msun. The accretion of these filaments onto the DR21 ridge, suggested by a previous molecular line study, could provide a continuous mass inflow to the ridge. In contrast to the striations seen in e.g., the Taurus region, these filaments are gravitationally unstable and form cores and protostars. These cores formed in the filaments potentially fall into the ridge. Both inflow and collisions of cores could be important to drive the observed high-mass star formation. The evolutionary gradient of star formation running from DR21 in the south to the northern branching is traced by decreasing dust temperature. This evolution and the ridge structure can be explained by two main filamentary components of the ridge that merged first in the south.

We present a multi-scale, multi-wavelength source extraction algorithm called getsources. Although it has been designed primarily for use in the far-infrared surveys of Galactic star-forming regions with Herschel, the method can be applied to many other astronomical images. Instead of the traditional approach of extracting sources in the observed images, the new method analyzes fine spatial decompositions of original images across a wide range of scales and across all wavebands. It cleans those single-scale images of noise and background, and constructs wavelength-independent single-scale detection images that preserve information in both spatial and wavelength dimensions. Sources are detected in the combined detection images by following the evolution of their segmentation masks across all spatial scales. Measurements of the source properties are done in the original background-subtracted images at each wavelength; the background is estimated by interpolation under the source footprints and overlapping sources are deblended in an iterative procedure. In addition to the main catalog of sources, various catalogs and images are produced that aid scientific exploitation of the extraction results. We illustrate the performance of getsources on Herschel images by extracting sources in sub-fields of the Aquila and Rosette star-forming regions. The source extraction code and validation images with a reference extraction catalog are freely available.

For many years feedback processes generated by OB-stars in molecular clouds, including expanding ionization fronts, stellar winds, or UV-radiation, have been proposed to trigger subsequent star formation. However, hydrodynamic models including radiation and gravity show that UV-illumination has little or no impact on the global dynamical evolution of the cloud. The Rosette molecular cloud, irradiated by the NGC2244 cluster, is a template region for triggered star-formation, and we investigated its spatial and density structure by applying a curvelet analysis, a filament-tracing algorithm (DisPerSE), and probability density functions (PDFs) on Herschel column density maps, obtained within the HOBYS key program. The analysis reveals not only the filamentary structure of the cloud but also that all known infrared clusters except one lie at junctions of filaments, as predicted by turbulence simulations. The PDFs of sub-regions in the cloud show systematic differences. The two UV-exposed regions have a double-peaked PDF we interprete as caused by shock compression. The deviations of the PDF from the log-normal shape typically associated with low- and high-mass star-forming regions at Av~3-4m and 8-10m, respectively, are found here within the very same cloud. This shows that there is no fundamental difference in the density structure of low- and high-mass star-forming regions. We conclude that star-formation in Rosette - and probably in high-mass star-forming clouds in general - is not globally triggered by the impact of UV-radiation. Moreover, star formation takes place in filaments that arose from the primordial turbulent structure built up during the formation of the cloud. Clusters form at filament mergers, but star formation can be locally induced in the direct interaction zone between an expanding HII--region and the molecular cloud.

A growing body of evidence indicates that the formation of filaments in interstellar clouds is a key component of the star formation process. In this paper, we present new Herschel PACS and SPIRE observations of the B59 and Stem regions in the Pipe Nebula complex, revealing a rich, organized network of filaments. The asymmetric column density profiles observed for several filaments, along with the bow-like edge of B59, indicates that the Pipe Nebula is being compressed from its western side, most likely by the winds from the nearby Sco OB2 association. We suggest that this compressive flow has contributed to the formation of some of the observed filamentary structures. In B59, the only region of the entire Pipe complex showing star formation activity, the same compressive flow has likely enhanced the initial column density of the clump, allowing it to become globally gravitationally unstable. Although more speculative, we propose that gravity has also been responsible for shaping the converging filamentary pattern observed in B59. While the question of the relative impact of large-scale compression and gravity remains open in B59, large-scale compression appears to be a plausible mechanism for the initial formation of filamentary structures in the rest of the complex

A new adaptive mesh refinement (AMR) version of the ZEUS-3D astrophysical magnetohydrodynamical (MHD) fluid code, AZEuS, is described. The AMR module in AZEuS has been completely adapted to the staggered mesh that characterises the ZEUS family of codes, on which scalar quantities are zone-centred and vector components are face-centred. In addition, for applications using static grids, it is necessary to use higher-order interpolations for prolongation to minimise the errors caused by waves crossing from a grid of one resolution to another. Finally, solutions to test problems in 1-, 2-, and 3-dimensions in both Cartesian and spherical coordinates are presented.

The filament IRDC G035.39--00.33 in the W48 molecular complex is one of the darkest infrared clouds observed by \textit{Spitzer}. It has been observed by the PACS (70 and 160\,$\micron$) and SPIRE (250, 350, and 500\,$\micron$) cameras of the \textit{Herschel} Space Observatory as part of the W48 molecular cloud complex in the framework of the HOBYS key programme. The observations reveal a sample of 28 compact sources (deconvolved FWHM sizes $<$0.3 pc) complete down to $\sim$$5 \msun$ in G035.39--00.33 and its surroundings. Among them, 13 compact sources are massive dense cores with masses $>$$20 \msun$. The cloud characteristics we derive from the analysis of their spectral energy distributions are masses of $20-50 \msun$, sizes of 0.1--0.2 pc, and average densities of $2-20 \times 10^{5} \cmc$, which make these massive dense cores excellent candidates to form intermediate- to high-mass stars. Most of the massive dense cores are located inside the G035.39--00.33 ridge and host IR-quiet high-mass protostars. The large number of protostars found in this filament suggests that we are witnessing a mini-burst of star formation with an efficiency of $\sim$15% and a rate density of $\sim$$40 \msun\,$yr$^{-1}\,$kpc$^{-2}$ within $\sim$8 pc$^2$, a large area covering the full ridge. Part of the extended SiO emission observed towards G035.39--00.33 is not associated with obvious protostars and may originate from low-velocity shocks within converging flows, as advocated by previous studies.

Abridged -- We present an analysis of the Aquila Rift complex which addresses the questions of the star formation rate (SFR), star formation efficiency (SFE) and typical lifetime of the Class 0 protostellar phase in two nearby cluster-forming clumps: the Serpens South and W40 protoclusters. We carried out a 1.2 mm dust continuum mapping of the Aquila Rift complex with the MAMBO bolometer array on the IRAM 30m telescope. We perform a systematic source extraction in our millimeter continuum map. Based on complementary data from the Herschel Gould Belt survey and Spitzer maps, we characterize the SEDs of the 77 mm continuum sources detected with MAMBO and estimate their evolutionary stages. Taking advantage of the comprehensive dataset available for the Serpens South region, spanning wavelengths from 2 microns to 1.2 mm, we estimate the numbers of young stellar objects (YSOs) at different evolutionary stages and find a ratio of Class 0 to Class I protostars N(0)/N(I)=0.19-0.27. This low ratio supports a scenario of relatively fast accretion at the beginning of the protostellar phase, and leads to a Class 0 lifetime of ~4-9x10^{4} yr. We also show that both the Serpens South and W40 protoclusters are characterized by large fractions of protostars and high SFRs ~20-50 Msun.Myr^{-1}pc^{-2}, in agreement with the idea that these two nearby clumps are active sites of clustered star formation currently undergoing bursts of star formation, and have the potential ability to form bound star clusters. While the formation of these two protoclusters is likely to have been initiated in a very different manner, the resulting protostellar populations are observed to be very similar. This suggests that after the onset of gravitational collapse, the detailed manner in which the collapse has been initiated does not affect much the ability of a clump to form stars.

We present the first Herschel PACS and SPIRE results of the Vela C molecular complex in the far-infrared and submillimetre regimes at 70, 160, 250, 350, and 500 um, spanning the peak of emission of cold prestellar or protostellar cores. Column density and multi-resolution analysis (MRA) differentiates the Vela C complex into five distinct sub-regions. Each sub-region displays differences in their column density and temperature probability distribution functions (PDFs), in particular, the PDFs of the `Centre-Ridge' and `South-Nest' sub-regions appear in stark contrast to each other. The Centre-Ridge displays a bimodal temperature PDF representative of hot gas surrounding the HII region RCW 36 and the cold neighbouring filaments, whilst the South-Nest is dominated by cold filamentary structure. The column density PDF of the Centre-Ridge is flatter than the South-Nest, with a high column density tail, consistent with formation through large-scale flows, and regulation by self-gravity. At small to intermediate scales MRA indicates the Centre-Ridge to be twice as concentrated as the South-Nest, whilst on larger scales, a greater portion of the gas in the South-Nest is dominated by turbulence than in the Centre-Ridge. In Vela C, high-mass stars appear to be preferentially forming in ridges, i.e., dominant high column density filaments.

We provide a first look at the results of the Herschel Gould Belt survey toward the IC5146 molecular cloud and present a preliminary analysis of the filamentary structure in this region. The column density map, derived from our 70-500 micron Herschel data, reveals a complex network of filaments, and confirms that these filaments are the main birth sites of prestellar cores. We analyze the column density profiles of 27 filaments and show that the underlying radial density profiles fall off as r^{-1.5} to r^{-2.5} at large radii. Our main result is that the filaments seem to be characterized by a narrow distribution of widths having a median value of 0.10 +- 0.03 pc, which is in stark contrast to a much broader distribution of central Jeans lengths. This characteristic width of ~0.1 pc corresponds to within a factor of ~2 to the sonic scale below which interstellar turbulence becomes subsonic in diffuse gas, supporting the argument that the filaments may form as a result of the dissipation of large-scale turbulence.

We summarize the first results from the Gould Belt survey, obtained toward the Aquila Rift and Polaris Flare regions during the 'science demonstration phase' of Herschel. Our 70-500 micron images taken in parallel mode with the SPIRE and PACS cameras reveal a wealth of filamentary structure, as well as numerous dense cores embedded in the filaments. Between ~ 350 and 500 prestellar cores and ~ 45-60 Class 0 protostars can be identified in the Aquila field, while ~ unbound starless cores and no protostars are observed in the Polaris field. The prestellar core mass function (CMF) derived for the Aquila region bears a strong resemblance to the stellar initial mass function (IMF), already confirming the close connection between the CMF and the IMF with much better statistics than earlier studies. Comparing and contrasting our Herschel results in Aquila and Polaris, we propose an observationally-driven scenario for core formation according to which complex networks of long, thin filaments form first within molecular clouds, and then the densest filaments fragment into a number of prestellar cores via gravitational instability.

The Rosette molecular cloud is promoted as the archetype of a triggered star-formation site. This is mainly due to its morphology, because the central OB cluster NGC 2244 has blown a circular-shaped cavity into the cloud and the expanding HII-region now interacts with the cloud. Studying the spatial distribution of the different evolutionary states of all star-forming sites in Rosette and investigating possible gradients of the dust temperature will help to test the 'triggered star-formation' scenario in Rosette. We use continuum data obtained with the PACS (70 and 160 micron) and SPIRE instruments (250, 350, 500 micron) of the Herschel telescope during the Science Demonstration Phase of HOBYS. Three-color images of Rosette impressively show how the molecular gas is heated by the radiative impact of the NGC 2244 cluster. A clear negative temperature gradient and a positive density gradient (running from the HII-region/molecular cloud interface into the cloud) are detected. Studying the spatial distribution of the most massive dense cores (size scale 0.05 to 0.3 pc), we find an age-sequence (from more evolved to younger) with increasing distance to the cluster NGC 2244. No clear gradient is found for the clump (size-scale up to 1 pc) distribution. The existence of temperature and density gradients and the observed age-sequence imply that star formation in Rosette may indeed be influenced by the radiative impact of the central NGC 2244 cluster. A more complete overview of the prestellar and protostellar population in Rosette is required to obtain a firmer result.

We present a preliminary analysis of the small-scale structure found in new 70-520 micron continuum maps of the Rosette molecular cloud (RMC), obtained with the SPIRE and PACS instruments of the Herschel Space Observatory. We find 473 clumps within the RMC using a new structure identification algorithm, with sizes up to ~1.0 pc in diameter. A comparison with recent Spitzer maps reveals that 371 clumps are "starless" (without an associated young stellar object), while 102 are "protostellar." Using the respective values of dust temperature, we determine the clumps have masses (M_C) over the range -0.75 <= log (M_C/M_sun) <= 2.50. Linear fits to the high-mass tails of the resulting clump mass spectra (CMS) have slopes that are consistent with those found for high-mass clumps identified in CO emission by other groups.

Our PACS and SPIRE images of the Aquila Rift and part of the Polaris Flare regions, taken during the science demonstration phase of Herschel discovered fascinating, omnipresent filamentary structures that appear to be physically related to compact cores. We briefly describe a new multi-scale, multi-wavelength source extraction method used to detect objects and measure their parameters in our Herschel images. All of the extracted starless cores (541 in Aquila and 302 in Polaris) appear to form in the long and very narrow filaments. With its combination of the far-IR resolution and sensitivity, Herschel directly reveals the filaments in which the dense cores are embedded; the filaments are resolved and have deconvolved widths of 35 arcsec in Aquila and 59 arcsec in Polaris (9000 AU in both regions). Our first results of observations with Herschel enable us to suggest that in general dense cores may originate in a process of fragmentation of complex networks of long, thin filaments, likely formed as a result of an interplay between gravity, interstellar turbulence, and magnetic fields. To unravel the roles of the processes, one has to obtain additional kinematic and polarization information; these follow-up observations are planned.

The Herschel OB young stellar objects survey (HOBYS) has observed the Rosette molecular cloud, providing an unprecedented view of its star formation activity. These new far-infrared data reveal a population of compact young stellar objects whose physical properties we aim to characterise. We compiled a sample of protostars and their spectral energy distributions that covers the near-infrared to submillimetre wavelength range. These were used to constrain key properties in the protostellar evolution, bolometric luminosity, and envelope mass and to build an evolutionary diagram. Several clusters are distinguished including the cloud centre, the embedded clusters in the vicinity of luminous infrared sources, and the interaction region. The analysed protostellar population in Rosette ranges from 0.1 to about 15 Msun with luminosities between 1 and 150 Lsun, which extends the evolutionary diagram from low-mass protostars into the high-mass regime. Some sources lack counterparts at near- to mid-infrared wavelengths, indicating extreme youth. The central cluster and the Phelps & Lada 7 cluster appear less evolved than the remainder of the analysed protostellar population. For the central cluster, we find indications that about 25% of the protostars classified as Class I from near- to mid-infrared data are actually candidate Class 0 objects. As a showcase for protostellar evolution, we analysed four protostars of low- to intermediate-mass in a single dense core, and they represent different evolutionary stages from Class 0 to Class I. Their mid- to far-infrared spectral slopes flatten towards the Class I stage, and the 160 to 70um flux ratio is greatest for the presumed Class 0 source. This shows that the Herschel observations characterise the earliest stages of protostellar evolution in detail.

The origin and possible universality of the stellar initial mass function (IMF) is a major issue in astrophysics. One of the main objectives of the Herschel Gould Belt Survey is to clarify the link between the prestellar core mass function (CMF) and the IMF. We present and discuss the core mass function derived from Herschel data for the large population of prestellar cores discovered with SPIRE and PACS in the Aquila Rift cloud complex at d ~ 260 pc. We detect a total of 541 starless cores in the entire ~11 deg^2 area of the field imaged at 70-500 micron with SPIRE/PACS. Most of these cores appear to be gravitationally bound, and thus prestellar in nature. Our Herschel results confirm that the shape of the prestellar CMF resembles the stellar IMF, with much higher quality statistics than earlier submillimeter continuum ground-based surveys.

As part of the science demonstration phase of the Herschel mission of the Gould Belt Key Program, the Aquila Rift molecular complex has been observed. The complete ~ 3.3deg x 3.3deg imaging with SPIRE 250/350/500 micron and PACS 70/160 micron allows a deep investigation of embedded protostellar phases, probing of the dust emission from warm inner regions at 70 and 160 micron to the bulk of the cold envelopes between 250 and 500 micron. We used a systematic detection technique operating simultaneously on all Herschel bands to build a sample of protostars. Spectral energy distributions are derived to measure luminosities and envelope masses, and to place the protostars in an M_env - L_bol evolutionary diagram. The spatial distribution of protostars indicates three star-forming sites in Aquila, with W40/Sh2-64 HII region by far the richest. Most of the detected protostars are newly discovered. For a reduced area around the Serpens South cluster, we could compare the Herschel census of protostars with Spitzer results. The Herschel protostars are younger than in Spitzer with 7 Class 0 YSOs newly revealed by Herschel. For the entire Aquila field, we find a total of ~ 45-60 Class 0 YSOs discovered by Herschel. This confirms the global statistics of several hundred Class~0 YSOs that should be found in the whole Gould Belt survey.

The Polaris Flare cloud region contains a great deal of extended emission. It is at high declination and high Galactic latitude. It was previously seen strongly in IRAS Cirrus emission at 100 microns. We have detected it with both PACS and SPIRE on Herschel. We see filamentary and low-level structure. We identify the five densest cores within this structure. We present the results of a temperature, mass and density analysis of these cores. We compare their observed masses to their virial masses, and see that in all cases the observed masses lie close to the lower end of the range of estimated virial masses. Therefore, we cannot say whether they are gravitationally bound prestellar cores. Nevertheless, these are the best candidates to be potentialprestellar cores in the Polaris cloud region.

[abridged] Our knowledge of the inner structure of embedded massive young stellar objects is still quite limited. We attempt here to overcome the spatial resolution limitations of conventional thermal infrared imaging. We employed mid-infrared interferometry using the MIDI instrument on the ESO/VLTI facility to investigate M8E-IR, a well-known massive young stellar object suspected of containing a circumstellar disk. Spectrally dispersed visibilities in the 8-13 micron range were obtained at seven interferometric baselines. We resolve the mid-infrared emission of M8E-IR and find typical sizes of the emission regions of the order of 30 milli-arcseconds (~45 AU). Radiative transfer simulations have been performed to interpret the data. The fitting of the spectral energy distribution, in combination with the measured visibilities, does not provide evidence for an extended circumstellar disk with sizes > 100 AU but requires the presence of an extended envelope. The data are not able to constrain the presence of a small-scale disk in addition to an envelope. In either case, the interferometry measurements indicate the existence of a strongly bloated, relatively cool central object, possibly tracing the recent accretion history of M8E-IR. In addition, we present 24.5 micron images that clearly distinguish between M8E-IR and the neighbouring ultracompact HII region and which show the cometary-shaped infrared morphology of the latter source. Our results show that IR interferometry, combined with radiative transfer modelling, can be a viable tool to reveal crucial structure information on embedded massive young stellar objects and to resolve ambiguities arising from fitting the SED.

The very inner structure of massive young stellar objects (YSOs) is difficult to trace. With conventional observational methods we identify structures still several hundreds of AU in size. However, the (proto-)stellar growth takes place at the innermost regions (<100 AU) where the actual mass transfer onto the forming high-mass star occurs. We present results from our programme toward massive YSOs at the VLTI, utilising the two-element interferometer MIDI. To date, we observed 10 well-known massive YSOs down to scales of 20 mas (typically corresponding to 20 - 40 AU for our targets) in the 8-13 micron region. We clearly resolve these objects which results in low visibilities and sizes in the order of 30-50 mas. For two objects, we show results of our modelling. We demonstrate that the MIDI data can reveal decisive structure information for massive YSOs. They are often pivotal in order to resolve ambiguities still immanent in model parameters derived from sole SED fitting.

The very inner structure of massive YSOs is difficult to trace. With conventional observational methods we often identify structures still several hundreds of AU in size. But we also need information about the innermost regions where the actual mass transfer onto the forming high-mass star occurs. An innovative way to probe these scales is to utilise mid-infrared interferometry. Here, we present first results of our MIDI GTO programme at the VLTI. We observed 10 well-known massive YSOs down to scales of 20 mas. We clearly resolve these objects which results in low visibilities and sizes in the order of 30 - 50 mas. Thus, with MIDI we can for the first time quantify the extent of the thermal emission from the warm circumstellar dust and thus calibrate existing concepts regarding the compactness of such emission in the pre-UCHII region phase. Special emphasis will be given to the BN-type object M8E-IR where our modelling is most advanced and where there is indirect evidence for a strongly bloated central star.

We investigate the pulsation-phase dependent properties of the circumstellar dust shell (CDS) of the OH/IR star OH104.9+2.4 based on radiative transfer modeling (RTM) using the code DUSTY. Our previous study concerning simultaneous modeling of the spectral energy distribution (SED) and near-infrared (NIR) visibilities (Riechers et al. 2004) has now been extended by means of a more detailed analysis of the pulsation-phase dependence of the model parameters of OH104.9+2.4. In order to investigate the temporal variation in the spatial structure of the CDS, additional NIR speckle interferometric observations in the K' band were carried out with the 6 m telescope of the Special Astrophysical Observatory (SAO). At a wavelength of 2.12 micron the diffraction-limited resolution of 74 mas was attained. Several key parameters of our previous best-fitting model had to be adjusted in order to be consistent with the newly extended amount of observational data. It was found that a simple rescaling of the bolometric flux F_bol is not sufficient to take the variability of the source into account, as the change in optical depth over a full pulsation cycle is rather high. On the other hand, the impact of a change in effective temperature T_eff on SED and visibility is rather small. However, observations, as well as models for other AGB stars, show the necessity of including a variation of T_eff with pulsation phase in the radiative transfer models. Therefore, our new best-fitting model accounts for these changes.

We present near-infrared speckle interferometry of the OH/IR star OH 104.9+2.4 in the K' band obtained with the 6m telescope of the Special Astrophysical Observatory (SAO). At a wavelength of lambda = 2.12 micron the diffraction-limited resolution of 74 mas was attained. The reconstructed visibility reveals a spherically symmetric, circumstellar dust shell (CDS) surrounding the central star. The visibility function shows that the stellar contribution to the total flux at lambda = 2.12 micron is less than ~50%, indicating a rather large optical depth of the CDS. The azimuthally averaged 1-dimensional Gaussian visibility fit yields a diameter of 47 +/- 3mas (FHWM), which corresponds to 112 +/- 13 AU for an adopted distance of D = 2.38 +/- 0.24 kpc. To determine the structure and the properties of the CDS of OH 104.9+2.4, radiative transfer calculations using the code DUSTY were performed to simultaneously model its visibility and the spectral energy distribution (SED). We found that both the ISO spectrum and the visibility of OH 104.9+2.4 can be well reproduced by a radiative transfer model with an effective temperature T_eff = 2500 +/- 500 K of the central source, a dust temperature T_in = 1000 +/- 200 K at the inner shell boundary R_in = 9.1 R_star = 25.4 AU, an optical depth tau = 6.5 +/- 0.3 at 2.2 micron, and dust grain radii ranging from a_min = 0.005 +/- 0.003 micron to a_max = 0.2 +/- 0.02 micron with a power law with index -3.5. It was found that even minor changes in a_max have a major impact on both the slope and the curvature of the visibility function, while the SED shows only minor changes. Our detailed analysis demonstrates the potential of dust shell modeling constrained by both the SED and visibilities.

We present two-dimensional (2D) radiative transfer modeling of IRC+10216 at selected moments of its evolution in 1995-2001, which correspond to three epochs of our series of 8 near-infrared speckle images (Osterbart et al. 2000, Weigelt et al. 2002). The high-resolution images obtained over the last 5.4 years revealed the dynamic evolution of the subarcsecond dusty environment of IRC+10216 and our recent time-independent 2D radiative transfer modeling reconstructed its physical properties at the single epoch of January 1997 (Men'shchikov et al. 2001). Having documented the complex changes in the innermost bipolar shell of the carbon star, we incorporate the evolutionary constraints into our new modeling to understand the physical reasons for the observed changes. The new calculations imply that during the last 50 years, we have been witnessing an episode of a steadily increasing mass loss from the central star, from Mdot ~ 10^-5 Msun/yr to the rate of Mdot ~ 3x10^-4 Msun/yr in 2001. The rapid increase of the mass loss of IRC+10216 and continuing time-dependent dust formation and destruction caused the observed displacement of the initially faint components C and D and of the bright cavity A from the star which has almost disappeared in our images in 2001. Increasing dust optical depths are causing strong backwarming that leads to higher temperatures in the dust formation zone, displacing the latter outward with a velocity v_T ~ 27 km/s due to the evaporation of the recently formed dust grains. This shift of the dust density peak in the bipolar shell mimics a rapid radial expansion, whereas the actual outflow has probably a lower speed v < v_inf ~ 15 km/s. The model predicts that the star will remain obscured until Mdot starts to drop back to lower values in the dust formation zone.

We present new near-infrared (JHK) bispectrum speckle-interferometry monitoring of the carbon star IRC+10216 obtained between 1999 and 2001 with the SAO 6m telescope. The J-, H-, and K-band resolutions are 50mas, 56mas, and 73mas, resp. The total sequence of K-band observations covers now 8 epochs from 1995 to 2001 and shows the dynamic evolution of the inner dust shell. The present observations show that the appearance of the dust shell has considerably changed compared to the epochs of 1995 to 1998. Four main components within a 0.2" radius can be identified in the K-band images. The apparent separation of the two initially brightest components A and B increased from ~191mas in 1995 to ~351mas in 2001. Simultaneously, component B has been fading and almost disappeared in 2000 whereas the initially faint components C and D became brighter (relative to peak intensity). These changes can be related to changes of the optical depth caused, e.g., by mass-loss variations or new dust condensation in the wind. Our 2D radiative transfer model suggests that the observed relative motion of A and B is not consistent with the known terminal wind velocity of 15 km/s. The apparent motion with a deprojected velocity of 19 km/s on average and of recently 27 km/s appears to be caused by adisplacement of the dust density peak due to dust evaporation in the optically thicker and hotter environment. Our monitoring, covering more than 3 pulsation periods, shows that the structural variations are not related to the stellar pulsation cycle in a simple way. This is consistent with the predictions of hydrodynamical models that enhanced dust formation takes place on a timescale of several pulsation periods. The timescale of the fading of component B can well be explained by the formation of new dust in the circumstellar envelope.

New diffraction-limited speckle images of the Red Rectangle in the wavelength range 2.1--3.3 microns with angular resolutions of 44--68 mas and previous speckle images at 0.7--2.2 microns revealed well-resolved bright bipolar outflow lobes and long X-shaped spikes originating deep inside the outflow cavities. This set of high-resolution images stimulated us to reanalyze all infrared observations of the Red Rectangle using our two-dimensional radiative transfer code. The new detailed modeling, together with estimates of the interstellar extinction in the direction of the Red Rectangle enabled us to more accurately determine one of the key parameters, the distance D=710 pc with model uncertainties of 70 pc, which is twice as far as the commonly used estimate of 330 pc. The central binary is surrounded by a compact, massive (M=1.2 Msun), very dense dusty torus with hydrogen densities reaching n_H=2.5x10^12 cm^-3 (dust-to-gas mass ratio rho_d/rho~0.01). The bright component of the spectroscopic binary HD 44179 is a post-AGB star with mass M*=0.57 Msun, luminosity L*=6000 Lsun, and effective temperature T*=7750 K. Based on the orbital elements of the binary, we identify its invisible component with a helium white dwarf with Mwd~0.35 Msun, Lwd~100 Lsun, and Twd~6x10^4 K. The hot white dwarf ionizes the low-density bipolar outflow cavities inside the dense torus, producing a small HII region observed at radio wavelengths. We propose an evolutionary scenario for the formation of the Red Rectangle nebula, in which the binary initially had 2.3 and 1.9 Msun components at a separation of 130 Rsun. The nebula was formed in the ejection of a common envelope after Roche lobe overflow by the present post-AGB star.

We present spectrophotometric ISO imaging with the LWS and the CAM-CVF of the Serpens molecular cloud core. The LWS map is centred on the far infrared and submillimetre source SMM1 and its size is 8' x 8'. The fine structure line emission in [OI] and [CII] is extended and can be successfully modelled to originate in a PDR with G_0 = 15 and n(H2) about 10^4 - 10^5 cm^-3. Extended emission is also observed in the rotational line emission of H2O and high-J CO. However, lack of sufficient angular resolution prevents us from excluding the possibility that the emssion regions of these lines are point like, which could be linked to the embedded objects SMM9 and SMM4. Toward the Class0 source SMM1, the LWS observations reveal, in addition to fine structure line emission, a rich spectrum of molecular lines. The sub-thermally excited and optically thick CO, H2O and OH lines are tracing an about 10^3 AU source with temperatures higher than 300 K and densities above 10^6 cm^-3. We show that geometry is of concern for the correct interpretation of the data and based on 2D-radiative transfer modelling of the disk/torus around SMM1, which successfully reproduces the entire observed SED and the observed line profiles of CO isotopomers, we can exclude the disk to be the source of the LWS-molecular line emission. The CAM-CVF permits us to see a region of rotational H2 emission. This H2 gas has a temperature of 10^3 K, which suggests that the heating of the gas is achieved through relatively slow shocks. Although we are not able to establish any firm conclusion regarding the detailed nature of the shock waves, our observations of the molecular line emission from SMM1 can be explainable in terms of an admixture of J-shocks and of C-shocks.

Structural and chemical changes during the AGB and post-AGB evolution are discussed with respect to two recent observational and theoretical findings. On the one hand, high-resolution infrared observations revealed details of the dynamical evolution of the fragmented, bipolar dust shell around the far-evolved carbon star IRC+10216 giving evidence for rapid changes of an already PPN-like structure during the very end of the AGB evolution. On the other hand, stellar evolution calculations considering convective overshoot have shown how thermal pulses during the post-AGB stage lead to the formation of hydrogen-deficient post-AGB stars with abundance patterns consistent with those observed for Wolf-Rayet central stars.

ISO-LWS mapping observations of the Serpens molecular cloud core are presented. The spectral range is 50 - 200 micron and the map size is 8' X 8'. These observations suffer from severe source confusion at FIR wavelengths and we employ a Maximum Likelihood Method for the spectro-spatial deconvolution. The strong and fairly isolated source SMM1 FIRS1 presented a test case, whose modelled spectral energy distribution (SED), within observational errors, is identical to the observed one. The model results for the other infrared and submillimetre sources are therefore likely to represent their correct SEDs. Simulations demonstrating the reliability and potential of the developed method support this view. For the majority of sources the peak of the SEDs is found within the spectral range of the LWS and derived temperatures are generally higher (>30K) than have been found by earlier deconvolution attempts using IRAS data. SMM sizes are found to be only a few arcsec in diameter. In addition, the SMMs are generally optically thick even at LWS wavelengths. Self-consistent radiative transfer calculations indicates that the SMMs are optically thick out to even longer wavelengths. Models were calculated for five sources, for which sufficient data were available, viz. SMM 1, 2, 3, 4 and 9. These models are optically thick out to millimetre wavelengths (wavelength of unit optical depth 900 to 1400 micron). Envelope masses for these SMMs are in the range 2-6 Msun, which is of course considerably more massive than estimates based on the optically thin assumption. The luminosities are in the range 10-70 Lsun, suggesting the formation of low-mass to intermediate mass stars, so that the existence of such massive envelopes argues for extreme youth of the SMMs in the Serpens cloud core.

We present high-resolution J-, H-, and K-band observations and the first H-K color image of the carbon star IRC +10216. The images were reconstructed from 6m telescope speckle interferograms using the bispectrum speckle interferometry method. The H and K images with resolutions between 70mas and 92mas consist of several compact components within a 0.2" radius and a fainter asymmetric nebula. The brightest four components are denoted with A to D in the order of decreasing brightness in the 1996 image. A comparison of our images from 1995, 1996, 1997, and 1998 gives - almost like a movie of five frames - insight into the dynamical evolution of the inner nebula. For instance, the separation of the two brightest components A and B increased from 191 mas in 1995 to 265 mas in 1998. At the same time, component B is fading and the components C and D become brighter. The X-shaped bipolar structure of the nebula, most prominently present in the J-band image, implies an asymmetric mass loss. Such asymmetries are often present in protoplanetary nebulae but are unexpected for AGB stars. IRC +10216 is thus likely to be very advanced in its AGB evolution, shortly before turning into a protoplanetary nebula. The cometary shapes of A in the H and J images and in the 0.79 micron and 1.06 micron HST images suggest that the core of A is not the central star, but the southern lobe of a bipolar structure. The position of the central star is probably at or near the position of component B, where the H-K color has a value of 4.2. If the star is at or near B, then the components A, C, and D are likely to be located at the inner boundary of the dust shell.

We present high-resolution J-, H-, and K-band observations of the carbon star IRC+10216. The images were reconstructed from 6 m telescope speckle interferograms using the bispectrum speckle interferometry method. The H and K images consist of several compact components within a 0.2" radius and a fainter asymmetric nebula. The brightest four components are denoted with A to D in the order of decreasing brightness. A comparison of our images gives - almost like a movie of five frames - insight to the dynamical evolution of the inner nebula. For instance, the separation of the two brightest components A and B increased by almost 40% from 191 mas in 1995 to 265 mas in 1998. At the same time, component B is fading and the components C and D become brighter. The X-shaped bipolar structure of the nebula implies an asymmetric mass-loss suggesting that IRC+10216 is very advanced in its AGB evolution, shortly before turning into a protoplanetary nebula. The cometary shape of component A suggests that the core of A is not the central star, but the southern lobe of a bipolar structure. The position of the central star is probably at or near the position of component B.