The radiative transport of photons through arbitrary three-dimensional (3D) structures of dust is a challenging problem due to the anisotropic scattering of dust grains and strong coupling between different spatial regions. The radiative transfer problem in 3D is solved using Monte Carlo or Ray Tracing techniques as no full analytic solution exists for the true 3D structures. We provide the first 3D dust radiative transfer benchmark composed of a slab of dust with uniform density externally illuminated by a star. This simple 3D benchmark is explicitly formulated to provide tests of the different components of the radiative transfer problem including dust absorption, scattering, and emission. This benchmark includes models with a range of dust optical depths fully probing cases that are optically thin at all wavelengths to optically thick at most wavelengths. This benchmark includes solutions for the full dust emission including single photon (stochastic) heating as well as two simplifying approximations: One where all grains are considered in equilibrium with the radiation field and one where the emission is from a single effective grain with size-distribution-averaged properties. A total of six Monte Carlo codes and one Ray Tracing code provide solutions to this benchmark. Comparison of the results revealed that the global SEDs are consistent on average to a few percent for all but the scattered stellar flux at very high optical depths. The image results are consistent within 10%, again except for the stellar scattered flux at very high optical depths. The lack of agreement between different codes of the scattered flux at high optical depths is quantified for the first time. We provide the first 3D dust radiative transfer benchmark and validate the accuracy of this benchmark through comparisons between multiple independent codes and detailed convergence tests.

The mass of prestellar cores is an essential ingredient to understand the onset of star formation in the core. The low level of emission from cold dust may keep parts of it hidden from observation. We aim to determine the fraction of core mass in the temperature range < 8 K that can be expected for typical low- and high-mass star formation regions. We calculate the dust temperature within standard spherically symmetric prestellar cores for a grid of density powerlaws in the outer core regions, core masses, and variation in the external multi-component radiation field. The dust is assumed to be composed of amorphous silicate and carbon, and variations of its optical properties are discussed. As measure for the distribution of cores and clumps, we use core mass functions derived for various environments. In view of the high densities in very cold central regions, dust and gas temperatures are assumed to be equal. We find that the fraction of mass with temperatures < 8 K in typical low- and high-mass cores is < 20\%. It is possible to obtain higher fractions of very cold gas by placing intermediate- or high-mass cores in a typical low-mass star formation environment. We show that the mass uncertainty arising from FIR to mm modeling of very cold dust emission is smaller than the mass uncertainty due to the unknown dust opacities.

Short version: We apply a new synergetic radiative transfer method: the derived 1D density profiles are both consistent with a cut through the Herschel PACS/SPIRE and JCMT SCUBA-2 continuum maps of L1689B and with a derived local interstellar radiation field. Choosing an appropriate cut along the filament major axis, we minimize the impact of the filament emission on the modeling. For the bulk of the core (5000-20000 au) an isothermal sphere model with a temperature of around 10 K provides the best fits. We show that the power law index of the density profile, as well as the constant temperature can be derived directly from the radial surface brightness profiles. For the inner region (< 5000 au), we find a range of densities and temperatures that are consistent with the surface brightness profiles and the local interstellar radiation field. Based on our core models, we find that pixel-by-pixel single temperature spectral energy distribution fits are incapable of determining dense core properties. We conclude that, to derive physical core properties, it is important to avoid an azimuthal average of core and filament. Correspondingly, derived core masses are too high since they include some mass of the filament, and might introduce errors when determining core mass functions. The forward radiative transfer methods also avoids the loss of information owing to smearing of all maps to the coarsest spatial resolution. We find the central core region to be colder and denser than estimated in recent inverse radiative transfer modeling, possibly indicating the start of star formation in L1689B.

Recently discovered scattered light from molecular cloud cores in the wavelength range 3-5 {\mu}m (called "coreshine") seems to indicate the presence of grains with sizes above 0.5 {\mu}m. We aim to analyze 3.6 and 4.5 {\mu}m coreshine from molecular cloud cores to probe the largest grains in the size distribution. We analyzed dedicated deep Cycle 9 Spitzer IRAC observations in the 3.6 and 4.5 {\mu}m bands for a sample of 10 low-mass cores. We used a new modeling approach based on a combination of ratios of the two background- and foreground-subtracted surface brightnesses and observed limits of the optical depth. The dust grains were modeled as ice-coated silicate and carbonaceous spheres. We discuss the impact of local radiation fields with a spectral slope differing from what is seen in the DIRBE allsky maps. For the cores L260, ecc806, L1262, L1517A, L1512, and L1544, the model reproduces the data with maximum grain sizes around 0.9, 0.5, 0.65, 1.5, 0.6, and > 1.5 {\mu}m, respectively. The maximum coreshine intensities of L1506C, L1439, and L1498 in the individual bands require smaller maximum grain sizes than derived from the observed distribution of band ratios. Additional isotropic local radiation fields with a spectral shape differing from the DIRBE map shape do not remove this discrepancy. In the case of Rho Oph 9, we were unable to reliably disentangle the coreshine emission from background variations and the strong local PAH emission. Considering surface brightness ratios in the 3.6 and 4.5 {\mu}m bands across a molecular cloud core is an effective method of disentangling the complex interplay of structure and opacities when used in combination with observed limits of the optical depth.

We define an appropriate problem for benchmarking dust emissivity calculations in the context of radiative transfer (RT) simulations, specifically including the emission from stochastically heated dust grains. Our aim is to provide a self-contained guide for implementors of such functionality, and to offer insights in the effects of the various approximations and heuristics implemented by the participating codes to accelerate the calculations. The benchmark problem definition includes the optical and calorimetric material properties, and the grain size distributions, for a typical astronomical dust mixture with silicate, graphite and PAH components; a series of analytically defined radiation fields to which the dust population is to be exposed; and instructions for the desired output. We process this problem using six RT codes participating in this benchmark effort, and compare the results to a reference solution computed with the publicly available dust emission code DustEM. The participating codes implement different heuristics to keep the calculation time at an acceptable level. We study the effects of these mechanisms on the calculated solutions, and report on the level of (dis)agreement between the participating codes. For all but the most extreme input fields, we find agreement within 10% across the important wavelength range from 3 to 1000 micron. We conclude that the relevant modules in RT codes can and do produce fairly consistent results for the emissivity spectra of stochastically heated dust grains.

Context: Observations of scattered light in the 1-5 $\mu$m range have revealed dust grains in molecular cores with sizes larger than commonly inferred for the diffuse interstellar medium. It is currently unclear whether these grains are grown within the molecular cores or are an ubiquitous component of the interstellar medium. Aims: We investigate whether the large grains necessary for efficient scattering at 1-5 mum are associated with the abundance of water-ice within molecular clouds and cores. Methods: We combined water-ice abundance measurements for sight lines through the Lupus IV molecular cloud complex with measurements of the scattered light at 3.6 mum for the same sight lines. Results: We find that there is a similar threshold for the cores in emission in scattered light at 3.6 mum (tau_9.7=0.15pm0.05, A_K=0.4pm0.2 as water-ice (tau_9.7=0.11pm0.01, A_K=0.19pm0.04) and that the scattering efficiency increases as the relative water-ice abundance increases. The ice layer increases the average grain size, which again strongly increases the albedo. Conclusions: The higher scattering efficiency is partly due to layering of ice on the dust grains. Although the layer can be relatively thin it can enhance the scattering substantially.

Context. Using observations to deduce dust properties, grain size distribution, and physical conditions in molecular clouds is a highly degenerate problem. Aims. The coreshine phenomenon, a scattering process at 3.6 and 4.5 $\mu$m that dominates absorption, has revealed its ability to explore the densest parts of clouds. We want to use this effect to constrain the dust parameters. The goal is to investigate to what extent grain growth (at constant dust mass) inside molecular clouds is able to explain the coreshine observations. We aim to find dust models that can explain a sample of Spitzer coreshine data. We also look at the consistency with near-infrared data we obtained for a few clouds. Methods. We selected four regions with a very high occurrence of coreshine cases: Taurus-Perseus, Cepheus, Chameleon and L183/L134. We built a grid of dust models and investigated the key parameters to reproduce the general trend of surface bright- nesses and intensity ratios of both coreshine and near-infrared observations with the help of a 3D Monte-Carlo radiative transfer code. The grid parameters allow to investigate the effect of coagulation upon spherical grains up to 5 $\mu$m in size derived from the DustEm diffuse interstellar medium grains. Fluffiness (porosity or fractal degree), ices, and a handful of classical grain size distributions were also tested. We used the near- and mostly mid-infrared intensity ratios as strong discriminants between dust models. Results. The determination of the background field intensity at each wavelength is a key issue. In particular, an especially strong background field explains why we do not see coreshine in the Galactic plane at 3.6 and 4.5 $\mu$m. For starless cores, where detected, the observed 4.5 $\mu$m / 3.6 $\mu$m coreshine intensity ratio is always lower than $\sim$0.5 which is also what we find in the models for the Taurus-Perseus and L183 directions. Embedded sources can lead to higher fluxes (up to four times greater than the strongest starless core fluxes) and higher coreshine ratios (from 0.5 to 1.1 in our selected sample). Normal interstellar radiation field conditions are sufficient to find suitable grain models at all wavelengths for starless cores. The standard interstellar grains are not able to reproduce observations and, due to the multi-wavelength approach, only a few grain types meet the criteria set by the data. Porosity does not affect the flux ratios while the fractal dimension helps to explain coreshine ratios but does not seem able to reproduce near-infrared observations without a mix of other grain types. Conclusions. Combined near- and mid-infrared wavelengths confirm the potential to reveal the nature and size distribution of dust grains. Careful assessment of the environmental parameters (interstellar and background fields, embedded or nearby reddened sources) is required to validate this new diagnostic.

The recently discovered coreshine effect can aid in exploring the core properties and in probing the large grain population of the ISM. We discuss the implications of the coreshine detected from the molecular cloud core L1506C in the Taurus filament for the history of the core and the existence of a primitive ISM component of large grains becoming visible in cores. The coreshine surface brightness of L1506C is determined from IRAC Spitzer images at 3.6 micron. We perform grain growth calculations to estimate the grain size distribution in model cores similar in gas density, radius, and turbulent velocity to L1506C. Scattered light intensities at 3.6 micron are calculated for a variety of MRN and grain growth distributions to compare with the observed coreshine. For a core with the overall physical properties of L1506C, no detectable coreshine is predicted for an MRN size distribution. Extending the distribution to grain radii of about 0.65 $\mu$m allows to reproduce the observed surface brightness level in scattered light. Assuming the properties of L1506C to be preserved, models for the growth of grains in cores do not yield sufficient scattered light to account for the coreshine within the lifetime of the Taurus complex. Only increasing the core density and the turbulence amplifies the scattered light intensity to a level consistent with the observed coreshine brightness. The grains could be part of primitive omni-present large grain population becoming visible in the densest part of the ISM, could grow under the turbulent dense conditions of former cores, or in L1506C itself. In the later case, L1506C must have passed through a period of larger density and stronger turbulence. This would be consistent with the surprisingly strong depletion usually attributed to high column densities, and with the large-scale outward motion of the core envelope observed today.

Recently discovered scattered light at 3-5 $\mu$m from low-mass cores (so-called "coreshine") reveals the presence of grains around 1 $\mu$m, which is larger than the grains found in the low-density interstellar medium. But only about half of the 100+ cores investigated so far show the effect. This prompts further studies on the origin of this detection rate. From the 3D continuum radiative transfer equation, we derive the expected scattered light intensity from a core placed in an arbitrary direction seen from Earth. We use the approximation of single scattering, consider extinction up to 2nd-order Taylor approximation, and neglect spatial gradients in the dust size distribution. The impact of the directional characteristics of the scattering on the detection of scattered light from cores is calculated for a given grain size distribution, and local effects like additional radiation field components are discussed. The surface brightness profiles of a core with a 1D density profile are calculated for various Galactic locations, and the results are compared to the approximate detection limits. We find that for optically thin radiation and a constant size distribution, a simple limit for detecting scattered light from a low-mass core can be derived that holds for grains with sizes smaller than 0.5 $\mu$m. The extinction by the core prohibits detection in bright parts of the Galactic plane, especially near the Galactic center. For scattered light received from low-mass cores with grain sizes beyond 0.5 $\mu$m, the directional characteristics of the scattering favors the detection of scattered light above and below the Galactic center, and to some extent near the Galactic anti-center. We identify the local incident radiation field as the major unknown causing deviations from this simple scheme.

Context. The largest grains (0.5-1 micron) in the interstellar size distribution are efficient in scattering near- and mid-infrared radiation. These wavelengths are therefore particularly well suited to probe the still uncertain high-end of the size distribution. Aims. We investigate the change in appearance of a typical low-mass molecular core from the Ks (2.2 micron) band to the Spitzer IRAC 3.6 and 8 micron bands, and compare with model calculations, which include variations of the grain size distribution. Methods. We combine Spitzer IRAC and ground-based near-infrared observations to characterize the scattered light observed at the near- and mid-infrared wavelengths from the core L260. Using a spherical symmetric model core, we perform radiative transfer calculations to study the impact of various dust size distributions on the intensity profiles across the core. Results. The observed scattered light patterns in the Ks and 3.6 micron bands are found to be similar. By comparison with radiative transfer models the two profiles place constraints on the relative abundance of small and large (more than 0.25 micron) dust grains. The scattered light profiles are found to be inconsistent with an interstellar silicate grain distribution extending only to 0.25 micron and large grains are needed to reach the observed fluxes and the flux ratios. The shape of the Ks band surface brightness profile limits the largest grains to 1-1.5 micron. Conclusions. In addition to observing coreshine in the Spitzer IRAC channels, the combination with ground-based near-infrared observations are suited to constrain the properties of large grains in cores.

Cosmic dust is present in many astrophysical objects, and recent observations across the electromagnetic spectrum show that the dust distribution is often strongly three-dimensional (3D). Dust grains are effective in absorbing and scattering ultraviolet (UV)/optical radiation, and they re-emit the absorbed energy at infrared wavelengths. Understanding the intrinsic properties of these objects, including the dust itself, therefore requires 3D dust radiative transfer (RT) calculations. Unfortunately, the 3D dust RT problem is nonlocal and nonlinear, which makes it one of the hardest challenges in computational astrophysics. Nevertheless, significant progress has been made in the past decade, with an increasing number of codes capable of dealing with the complete 3D dust RT problem. We discuss the complexity of this problem, the two most successful solution techniques [ray-tracing (RayT) and Monte Carlo (MC)], and the state of the art in modeling observational data using 3D dust RT codes. We end with an outlook on the bright future of this field.

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 form. Nevertheless, density profiles often rely either on the simplifying assumption of isothermality or on observationally poorly constrained model temperature profiles. With the aim of better constraining the initial physical conditions in molecular cloud cores at the onset of protostellar collapse, 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, including all observations, the analysis method, and the initial results of the survey. We study the thermal dust emission of 12 previously well-characterized, isolated, nearby globules using FIR and submm continuum maps at up to eight wavelengths between 100 micron and 1.2 mm. Our sample contains both globules with starless cores and embedded protostars at different early evolutionary stages. The dust emission maps are used to extract spatially resolved SEDs, which are then fit independently with modified blackbody curves to obtain line-of-sight-averaged dust temperature and column density maps. We find that the thermal structure of all globules is dominated by external heating from the interstellar radiation field and moderate shielding by thin extended halos. All globules have warm outer envelopes (14-20 K) and colder dense interiors (8-12 K). The protostars embedded in some of the globules raise the local temperature of the dense cores only within radii out to about 5000 AU, but do not significantly affect the overall thermal balance of the globules.

(Abriged) In the framework of the Herschel GTKP "The earliest phases of star formation", we have imaged B68 between 100 and 500 um. Ancillary (sub)mm data, spectral line maps of the 12/13CO(2-1) transitions as well as a 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 LoS 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. For the first time, we spatially resolve the dust temperature and density distribution of B68. We find T_dust dropping from 16.7 K at the edge to 8.2 K in the centre, which is about 4 K lower than the result of the simple SED fitting approach. N_H peaks at 4.3x10^22 cm^-2 and n_H at 3.4x10^5 cm^-3 in the centre. B68 has a mass of 3.1 M_sun of material with A_K > 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. We find the radial density distribution from the edge of the inner plateau outward to be n_H ~ 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. 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 which dispersed. We also resolve a new compact source in the SE trunk and find that it is slightly shifted in centroid velocity from B68, lending qualitative support to core collision scenarios.

(Abridged) We present an overview of the sample of high-mass star and cluster forming regions observed as part of the Earliest Phases of Star Formation (EPoS) Herschel Guaranteed Time Key Program. A sample of 45 infrared-dark clouds (IRDCs) were mapped at PACS 70, 100, and 160 micron and SPIRE 250, 350, and 500 micron. In this paper, we characterize a population of cores which appear in the PACS bands and place them into context with their host cloud and investigate their evolutionary stage. We construct spectral energy distributions (SEDs) of 496 cores which appear in all PACS bands, 34% of which lack counterparts at 24 micron. 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. The core population has a median temperature of 20K and has masses and luminosities that span four to five orders of magnitude. Cores with a counterpart at 24 micron are warmer and bluer on average than cores without a 24 micron counterpart. We conclude that cores bright at 24 micron are on average more advanced in their evolution, where a central protostar(s) have heated the outer bulk of the core, than 24 micron-dark cores. The 24 micron 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 micron. These are candidate starless objects, and further observations will help us determine the nature of these very cold cores.

Cold molecular clouds are the birthplaces of stars and planets, where dense cores of gas collapse to form protostars. The dust mixed in these clouds is thought to be made of grains of an average size of 0.1 micrometer. We report the widespread detection of the coreshine effect as a direct sign of the existence of grown, micrometer-sized dust grains. This effect is seen in half of the cores we have analyzed in our survey, spanning all Galactic longitudes, and is dominated by changes in the internal properties and local environment of the cores, implying that the coreshine effect can be used to constrain fundamental core properties such as the three-dimensional density structure and ages and also the grain characteristics themselves.

We present status and results of AstroGrid-D, a joint effort of astrophysicists and computer scientists to employ grid technology for scientific applications. AstroGrid-D provides access to a network of distributed machines with a set of commands as well as software interfaces. It allows simple use of computer and storage facilities and to schedule or monitor compute tasks and data management. It is based on the Globus Toolkit middleware (GT4). Chapter 1 describes the context which led to the demand for advanced software solutions in Astrophysics, and we state the goals of the project. We then present characteristic astrophysical applications that have been implemented on AstroGrid-D in chapter 2. We describe simulations of different complexity, compute-intensive calculations running on multiple sites, and advanced applications for specific scientific purposes, such as a connection to robotic telescopes. We can show from these examples how grid execution improves e.g. the scientific workflow. Chapter 3 explains the software tools and services that we adapted or newly developed. Section 3.1 is focused on the administrative aspects of the infrastructure, to manage users and monitor activity. Section 3.2 characterises the central components of our architecture: The AstroGrid-D information service to collect and store metadata, a file management system, the data management system, and a job manager for automatic submission of compute tasks. We summarise the successfully established infrastructure in chapter 4, concluding with our future plans to establish AstroGrid-D as a platform of modern e-Astronomy.

This paper has been withdrawn by the authors.

We present multi-wavelengths observations and a radiative transfer model of a newly discovered massive circumstellar disk of gas and dust which is one of the largest disks known today. Seen almost edge-on, the disk is resolved in high-resolution near-infrared (NIR) images and appears as a dark lane of high opacity intersecting a bipolar reflection nebula. Based on molecular line observations we estimate the distance to the object to be 3.5 kpc. This leads to a size for the dark lane of ~10500 AU but due to shadowing effects the true disk size could be smaller. In Spitzer/IRAC 3.6 micron images the elongated shape of the bipolar reflection nebula is still preserved and the bulk of the flux seems to come from disk regions that can be detected due to the slight inclination of the disk. At longer IRAC wavelengths, the flux is mainly coming from the central regions penetrating directly through the dust lane. Interferometric observations of the dust continuum emission at millimeter wavelengths with the SMA confirm this finding as the peak of the unresolved mm-emission coincides perfectly with the peak of the Spitzer/IRAC 5.8 micron flux and the center of the dark lane seen in the NIR images. Simultaneously acquired CO data reveal a molecular outflow along the northern part of the reflection nebula which seems to be the outflow cavity. An elongated gaseous disk component is also detected and shows signs of rotation. The emission is perpendicular to the molecular outflow and thus parallel to but even more extended than the dark lane in the NIR images. Based on the dust continuum and the CO observations we estimate a disk mass of up to a few solar masses depending on the underlying assumptions. Whether the disk-like structure is an actual accretion disk or rather a larger-scale flattened envelope or pseudodisk is difficult to discriminate with the current dataset (abridged).

Aims: Our aim is to understand the evolutionary sequence of high-mass star formation from the earliest evolutionary stage of high-mass starless cores, via high-mass cores with embedded low- to intermediate-mass objects, to finally high-mass protostellar objects. Methods: Herschel far-infrared PACS and SPIRE observations are combined with existing data at longer and shorter wavelengths to characterize the spectral and physical evolution of massive star-forming regions. Results: The new Herschel images spectacularly show the evolution of the youngest and cold high-mass star-forming regions from mid-infrared shadows on the Wien-side of the spectral energy distribution (SED), via structures almost lost in the background emission around 100mum, to strong emission sources at the Rayleigh-Jeans tail. Fits of the SEDs for four exemplary regions covering evolutionary stages from high-mass starless cores to high-mass protostellar objects reveal that the youngest regions can be fitted by single-component black-bodies with temperatures on the order of 17K. More evolved regions show mid-infrared excess emission from an additional warmer component, which however barely contributes to the total luminosities for the youngest regions. Exceptionally low values of the ratio between bolometric and submm luminosity additionally support the youth of the infrared-dark sources. Conclusions: The Herschel observations reveal the spectral and physical properties of young high-mass star-forming regions in detail. The data clearly outline the evolutionary sequence in the images and SEDs. Future work on larger samples as well as incorporating full radiative transfer calculations will characterize the physical nature at the onset of massive star formation in even more depth.

We present Herschel observations of the isolated, low-mass star-forming Bok globule CB244. It contains two cold sources, a low-mass Class 0 protostar and a starless core, which is likely to be prestellar in nature, separated by 90 arcsec (~ 18000 AU). The Herschel data sample the peak of the Planck spectrum for these sources, and are therefore ideal for dust-temperature and column density modeling. With these data and a near-IR extinction map, the MIPS 70 micron mosaic, the SCUBA 850 micron map, and the IRAM 1.3 mm map, we model the dust-temperature and column density of CB244 and present the first measured dust-temperature map of an entire star-forming molecular cloud. We find that the column-averaged dust-temperature near the protostar is ~ 17.7 K, while for the starless core it is ~ 10.6K, and that the effect of external heating causes the cloud dust-temperature to rise to ~ 17 K where the hydrogen column density drops below 10^21 cm^-2. The total hydrogen mass of CB244 (assuming a distance of 200 pc) is 15 +/- 5 M_sun. The mass of the protostellar core is 1.6 +/- 0.1 M_sun and the mass of the starless core is 5 +/- 2 M_sun, indicating that ~ 45% of the mass in the globule is participating in the star-formation process.

The intermediate-mass star-forming core UYSO 1 has previously been found to exhibit intriguing features. While deeply embedded and previously only identified by means of its (sub-)millimeter emission, it drives two powerful, dynamically young, molecular outflows. Although the process of star formation has obviously started, the chemical composition is still pristine. We present Herschel PACS and SPIRE continuum data of this presumably very young region. The now complete coverage of the spectral energy peak allows us to precisely constrain the elevated temperature of 26 - 28 K for the main bulge of gas associated with UYSO1, which is located at the interface between the hot HII region Sh 2-297 and the cold dark nebula LDN 1657A. Furthermore, the data identify cooler compact far-infrared sources of just a few solar masses, hidden in this neighbouring dark cloud.

Theoretical arguments suggest that dust grains should grow in the dense cold parts of molecular clouds. Evidence of larger grains has so far been gathered in near/mid infrared extinction and millimeter observations. Interpreting the data is, however, aggravated by the complex interplay of density and dust properties (as well as temperature for thermal emission). We present new Spitzer data of L183 in bands that are sensitive and insensitive to PAHs. The visual extinction AV map derived in a former paper was fitted by a series of 3D Gaussian distributions. For different dust models, we calculate the scattered MIR radiation images of structures that agree agree with the AV map and compare them to the Spitzer data. The Spitzer data of L183 show emission in the 3.6 and 4.5 micron bands, while the 5.8 micron band shows slight absorption. The emission layer of stochastically heated particles should coincide with the layer of strongest scattering of optical interstellar radiation, which is seen as an outer surface on I band images different from the emission region seen in the Spitzer images. Moreover, PAH emission is expected to strongly increase from 4.5 to 5.8 micron, which is not seen. Hence, we interpret this emission to be MIR cloudshine. Scattered light modeling when assuming interstellar medium dust grains without growth does not reproduce flux measurable by Spitzer. In contrast, models with grains growing with density yield images with a flux and pattern comparable to the Spitzer images in the bands 3.6, 4.5, and 8.0 micron.

We investigated the nature of the hitherto unresolved elliptical infrared emission in the centre of the ~20000 AU disc silhouette in M 17. We combined high-resolution JHKsL'M' band imaging carried out with NAOS/CONICA at the VLT with [Fe II] narrow band imaging using SOFI at the NTT. The analysis is supported by Spitzer/GLIMPSE archival data and by already published SINFONI/VLT Integral Field Spectroscopy data. For the first time, we resolve the elongated central infrared emission into a point-source and a jet-like feature that extends to the northeast in the opposite direction of the recently discovered collimated H2 jet. They are both orientated almost perpendicular to the disc plane. In addition, our images reveal a curved southwestern emission nebula whose morphology resembles that of the previously detected northeastern one. Both nebulae are located at a distance of 1500 AU from the disc centre. We describe the infrared point-source in terms of a protostar that is embedded in circumstellar material producing a visual extinction of 60 <= Av <= 82. The observed Ks band magnitude is equivalent to a stellar mass range of 2.8 Msun <= Mstar <= 8 Msun adopting conversions for a main-sequence star. Altogether, we suggest that the large M 17 accretion disc is forming an intermediate to high-mass protostar. Part of the accreted material is expelled through a symmetric bipolar jet/outflow.

We investigate the morphology and the evolutionary stage of the hypercompact HII region M17-UC1 using observations at infrared wavelengths and NIR radiative transfer modeling. For the first time, this region is resolved into two emission areas separated by a dark lane reminiscent of an obscuring silhouette of a circumstellar disk. So far, the observational data as well as model calculations suggest that M17-UC1 is surrounded by a disk of cool dust. This direct direct detection of a circumstellar disk candidate around a hypercompact HII region is in agreement with the expectations of the disk accretion model for high-mass star formation.

The formation of low-mass protostars and especially of brown dwarfs currently are ``hot topics'' in cool star research. The talks contributed to this splinter session discussed how low in mass and how low in luminosity objects might exist, if these substellar objects show evidence for circum(sub)stellar disks, and how the bottom of the mass function in young clusters after the formation process looks like. In a lively open discussion, a vast majority of the speakers and the audience expressed why, given the available data, a stellar-like formation mechanism down to the lowest masses should be preferred.

At the onset of high-mass star formation, accreting protostars are deeply embedded in massive cores made of gas and dust. Their spectral energy distribution is still dominated by the cold dust and rises steeply from near-to far-infrared wavelengths. The young massive star-forming region IRDC18223-3 is a prototypical Infrared-Dark-Cloud with a compact mm continuum core that shows no protostellar emission below 8mum. However, based on outflow tracers, early star formation activity was previously inferred for this region. Here, we present recent Spitzer observations from the MIPSGAL survey that identify the central protostellar object for the first time at 24 and 70mum. Combining the mid- to far-infrared data with previous mm continuum observations and the upper limits below 8mum, one can infer physical properties of the central source. At least two components with constant gas mass M and dust temperature T are necessary: one cold component (~15K and ~576M_sun) that contains most of the mass and luminosity, and one warmer component (>=51K and >=0.01M_sun) to explain the 24mum data. The integrated luminosity of ~177L_sun can be used to constrain additional parameters of the embedded protostar from the turbulent core accretion model for massive star formation. The data of IRDC18223-3 are consistent with a massive gas core harboring a low-mass protostellar seed of still less than half a solar mass with high accretion rates of the order 10^-4M_sun/yr. In the framework of this model, the embedded protostar is destined to become a massive star at the end of its formation processes.

We present a ray-tracing technique for radiative transfer modeling of complex three-dimensional (3D) structures which include dense regions of high optical depth like in dense molecular clouds, circumstellar disks, envelopes of evolved stars, and dust tori around active galactic nuclei. The corresponding continuum radiative transfer problem is described and the numerical requirements for inverse 3D density and temperature modeling are defined. We introduce a relative intensity and transform the radiative transfer equation along the rays to solve machine precision problems and to relax strong gradients in the source term. For the optically thick regions where common ray-tracers are forced to perform small trace steps, we give two criteria for making use of a simple approximative solver crossing the optically thick region quickly. Using an example of a density structure with optical depth changes of 6 orders of magnitude and sharp temperature variations, we demonstrate the accuracy of the proposed scheme using a common 5th-order Runge-Kutta ray-tracer with adaptive step size control. In our test case, the gain in computational speed is about a factor of 870. The method is applied to calculate the temperature distribution within a massive molecular cloud core for different boundary conditions for the radiation field.

Constraints on the density and thermal 3D structure of the dense molecular cloud core Rho Oph D are derived from a detailed 3D radiative transfer modeling. Two ISOCAM images at 7 and 15 micron are fitted simultaneously by representing the dust distribution in the core with a series of 3D Gaussian density profiles. Size, total density, and position of the Gaussians are optimized by simulated annealing to obtain a 2D column density map. The projected core density has a complex elongated pattern with two peaks. We propose a new method to calculate an approximate temperature in an externally illuminated complex 3D structure from a mean optical depth. This T(tau)-method is applied to a 1.3 mm map obtained with the IRAM 30m telescope to find the approximate 3D density and temperature distribution of the core Rho Oph D. The spatial 3D distribution deviates strongly from spherical symmetry. The elongated structure is in general agreement with recent gravo-turbulent collapse calculations for molecular clouds. We discuss possible ambiguities of the background determination procedure, errors of the maps, the accuracy of the T(tau)-method, and the influence of the assumed dust particle sizes and properties.

We analyze a three-dimensional smoothed particle hydrodynamics simulation of an evolving and later collapsing pre-stellar core. Using a three-dimensional continuum radiative transfer program, we generate images at 7 micron, 15 micron, 175 micron, and 1.3 mm for different evolutionary times and viewing angles. We discuss the observability of the properties of pre-stellar cores for the different wavelengths. For examples of non-symmetric fragments, it is shown that, misleadingly, the density profiles derived from a one-dimensional analysis of the corresponding images are consistent with one-dimensional core evolution models. We conclude that one-dimensional modeling based on column density interpretation of images does not produce reliable structural information and that multidimensional modeling is required.

We present benchmark problems and solutions for the continuum radiative transfer (RT) in a 2D disk configuration. The reliability of three Monte-Carlo and two grid-based codes is tested by comparing their results for a set of well-defined cases which differ for optical depth and viewing angle. For all the configurations, the overall shape of the resulting temperature and spectral energy distribution is well reproduced. The solutions we provide can be used for the verification of other RT codes.We also point out the advantages and disadvantages of the various numerical techniques applied to solve the RT problem.

We report the detection of D2CO in a sample of starless dense cores, in which we previously measured the degree of CO depletion. The deuterium fractionation is found extremely high, [D2CO]/[H2CO] ~ 1-10 %, similar to that reported in low-mass protostars. This provides convincing evidence that D2CO is formed in the cold pre-stellar cores, and later desorbed when the gas warms up in protostars. We find that the cores with the highest CO depletions have also the largest [D2CO]/[H2CO] ratios, supporting the theoretical prediction that deuteration increases with increasing CO depletion.

We present new results on CO depletion in a sample of nearby pre-stellar cores, based on observations of the millimeter C17O and C18O lines and the 1.3 mm dust emission with the IRAM 30m telescope. In most cases, the distribution of CO is much flatter than that of the dust, whereas other tracers, like N2H+, still probe the latter. In the centre of these objects, we estimate CO to be underabundant by a factor 4-15 depending on the cores. The CO underabundance is more pronounced in the central regions and appears to decrease with increasing distance from the core centre. This underabundance is most likely due to the freezing out of CO onto the dust grains in the cold, dense parts of the cores. We find evidence for an increase of the CO depletion degree with the core density.