We present observations of the Trumpler 14/Carina I region carried out using the Stratospheric Terahertz Observatory 2 (STO2). The Trumpler 14/Carina I region is in the west part of the Carina Nebula Complex, which is one of the most extreme star-forming regions in the Milky Way. We observed Trumpler 14/Carina I in the 158 $\mu$m transition of [C\,{\sc ii}] with a spatial resolution of 48$''$ and a velocity resolution of 0.17 km s$^{-1}$. The observations cover a 0.25$^\circ$ by 0.28$^\circ$ area with central position {\it l} = 297.34$^\circ$, {\it b} = -0.60$^\circ$. The kinematics show that bright [C\,{\sc ii}] structures are spatially and spectrally correlated with the surfaces of CO clouds, tracing the photodissociation region and ionization front of each molecular cloud. Along 7 lines of sight that traverse Tr 14 into the dark ridge to the southwest, we find that the [C\,{\sc ii}] luminosity from the HII region is 3.7 times that from the PDR. In same los we find in the PDRs an average ratio of 1:4.1:5.6 for the mass in atomic gas:dark-CO gas: molecular gas traced by CO. Comparing multiple gas tracers including HI 21cm, [C\,{\sc ii}], CO, and radio recombination lines, we find that the HII regions of the Carina Nebula Complex are well-described as HII regions with one-side freely expanding towards us, consistent with the champagne model of ionized gas evolution. The dispersal of the GMC in this region is dominated by EUV photoevaporation; the dispersal timescale is 20-30 Myr.

Following our recent work devoted to the effect of accretion on the pre-main-sequence evolution of low-mass stars, we perform a detailed analysis of episodic excursions of low-mass protostars in the Hertzsprung-Russell (H-R) diagram triggered by strong mass accretion bursts typical of FU Orionis-type objects (FUors). These excursions reveal themselves as sharp increases in the stellar total luminosity and/or effective temperature of the protostar and can last from hundreds to a few thousands of years, depending on the burst strength and characteristics of the protostar. During the excursions, low-mass protostars occupy the same part of the H-R diagram as young intermediate-mass protostars in the quiescent phase of accretion. Moreover, the time spent by low-mass protostars in these regions is on average a factor of several longer than that spent by the intermediate-mass stars in quiescence. During the excursions, low-mass protostars pass close to the position of most known FUors in the H-R diagram, but owing to intrinsic ambiguity the model stellar evolutionary tracks are unreliable in determining the FUor properties. We find that the photospheric luminosity in the outburst state may dominate the accretion luminosity already after a few years after the onset of the outburst, meaning that the mass accretion rates of known FUors inferred from the bolometric luminosity may be systematically overestimated, especially in the fading phase.

The Stratospheric Observatory for Infrared Astronomy (SOFIA) is a Boeing 747SP aircraft modified to accommodate a 2.7 meter gyro-stabilized telescope, which is mainly focused to studying the Universe at infrared wavelengths. As part of the Strategic Director's Discretionary Time (S-DDT) program, SOFIA performs observations of relevant science cases and immediately offers science-ready data products to the astronomical community. We present the first data release of the S-DDT program on far-infrared imaging polarimetric observations of 30 Doradus using the High-resolution Airborne Wideband Camera-Plus (HAWC+) at 53, 89, 154, and 214 micron. We present the status and quality of the observations, an overview of the SOFIA data products, and examples of working with HAWC+ polarimetric data that will enhance the scientific analysis of this, and future, data sets. These observations illustrate the potential influence of magnetic fields and turbulence in a star-forming region within the Tarantula Nebula.

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop.

We describe PyRaTE, a new, non-local thermodynamic equilibrium (non-LTE) line radiative transfer code developed specifically for post-processing astrochemical simulations. Population densities are estimated using the escape probability method. When computing the escape probability, the optical depth is calculated towards all directions with density, molecular abundance, temperature and velocity variations all taken into account. A very easy-to-use interface, capable of importing data from simulations outputs performed with all major astrophysical codes, is also developed. The code is written in Python using an `embarrassingly parallel' strategy and can handle all geometries and projection angles. We benchmark the code by comparing our results with those from RADEX (van der Tak et al. 2007) and against analytical solutions and present case studies using hydrochemical simulations. The code is available on GitHub (https://github.com/ArisTr/PyRaTE).

The pre-main-sequence evolution of low-mass stars and brown dwarfs is studied numerically starting from the formation of a protostellar/proto-brown dwarf seed and taking into account the mass accretion onto the central object during the initial several Myr of evolution. The stellar evolution was computed using the STELLAR evolution code developed by Yorke & Bodenheimer with recent modifications by Hosokawa et al. The mass accretion rates were taken from numerical hydrodynamics models of Vorobyov & Basu computing the circumstellar disk evolution starting from the gravitational collapse of pre-stellar cloud cores of various mass and angular momentum. The resulting stellar evolution tracks were compared with the isochrones and isomasses calculated using non-accreting models. We find that mass accretion in the initial several Myr of protostellar evolution can have a strong effect on the subsequent evolution of young stars and brown dwarfs. The disagreement between accreting and non-accreting models in terms of the total stellar luminosity L_st, stellar radius R_st and effective temperature T_eff depends on the thermal efficiency of accretion, i.e., on the fraction of accretion energy absorbed by the central object. The largest mismatch is found for the cold accretion case, in which essentially all accretion energy is radiated away. The relative deviations in L_st and R_st in this case can reach 50% for 1.0-Myr-old objects and remain notable even for 10-Myr-old objects. In the hot and hybrid accretion cases, in which a constant fraction of accretion energy is absorbed, the disagreement between accreting and non-accreting models becomes less pronounced, but still remains notable for 1.0-Myr-old objects. These disagreements may lead to the wrong age estimate for objects of (sub-)solar mass when using the isochrones based on non-accreting models (abridged).

We perform two-dimensional axially symmetric radiation-hydrodynamic simulations to assess the impact of outflows and radiative force feedback from massive protostars by varying when the protostellar outflow starts, the ratio of ejection to accretion rates, and the strength of the wide angle disk wind component. The star formation efficiency, i.e. the ratio of final stellar mass to initial core mass, is dominated by radiative forces and the ratio of outflow to accretion rates. Increasing this ratio has three effects: First, the protostar grows slower with a lower luminosity at any given time, lowering radiative feedback. Second, bipolar cavities cleared by the outflow are larger, further diminishing radiative feedback on disk and core scales. Third, the higher momentum outflow sweeps up more material from the collapsing envelope, decreasing the protostar's potential mass reservoir via entrainment. The star formation efficiency varies with the ratio of ejection to accretion rates from 50% in the case of very weak outflows to as low as 20% for very strong outflows. At latitudes between the low density bipolar cavity and the high density accretion disk, wide angle disk winds remove some of the gas, which otherwise would be part of the accretion flow onto the disk; varying the strength of these wide angle disk winds, however, alters the final star formation efficiency by only +/-6%. For all cases, the opening angle of the bipolar outflow cavity remains below 20 degree during early protostellar accretion phases, increasing rapidly up to 65 degree at the onset of radiation pressure feedback.

We present coupled stellar evolution (SE) and 3D radiation-hydrodynamic (RHD) simulations of the evolution of primordial protostars, their immediate environment, and the dynamic accretion history under the influence of stellar ionizing and dissociating UV feedback. Our coupled SE-RHD calculations result in a wide diversity of final stellar masses covering 10 Msun $\lesssim M_* \lesssim$ 1000 Msun. The formation of very massive ($\gtrsim$ 250 Msun) stars is possible under weak UV feedback, whereas ordinary massive (a few x 10 Msun) stars form when UV feedback can efficiently halt the accretion. This may explain the peculiar abundance pattern of a Galactic metal-poor star recently reported by Aoki et al. (2014), possibly the observational signature of very massive precursor primordial stars. Weak UV feedback occurs in cases of variable accretion, in particular when repeated short accretion bursts temporarily exceed 0.01 Msun/yr, causing the protostar to inflate. In the bloated state, the protostar has low surface temperature and UV feedback is suppressed until the star eventually contracts, on a thermal adjustment timescale, to create an HII region. If the delay time between successive accretion bursts is sufficiently short, the protostar remains bloated for extended periods, initiating at most only short periods of UV feedback. Disk fragmentation does not necessarily reduce the final stellar mass. Quite the contrary, we find that disk fragmentation enhances episodic accretion as many fragments migrate inward and are accreted onto the star, thus allowing continued stellar mass growth under conditions of intermittent UV feedback. This trend becomes more prominent as we improve the resolution of our simulations. We argue that simulations with significantly higher resolution than reported previously are needed to derive accurate gas mass accretion rates onto primordial protostars.

The formation of SMSs is a potential pathway to seed SMBHs in the early universe. A critical issue for forming SMSs is stellar UV feedback, which may limit the stellar mass growth via accretion. In this paper we study the evolution of an accreting SMS and its UV emissivity under conditions of realistic variable accretion from a self-gravitating circumstellar disc. First we conduct a 2D hydrodynamical simulation to follow the long-term protostellar accretion until the stellar mass exceeds $10^4~M_\odot$. The disc fragments due to gravitational instability, creating a number of small clumps that rapidly migrate inward to fall onto the star. The resulting accretion history is thus highly time-dependent: short episodic accretion bursts are followed by longer, relative quiescent phases. We show that the circumstellar disc for the so-called direct collapse model is more unstable and generates greater variability over shorter timescales than normal Pop III cases. We conduct a post-process stellar evolution calculation using the obtained accretion history. Our results show that, regardless of the strong variability of the accretion rates, the stellar radius monotonically increases with almost constant effective temperature at $T_{\rm eff} \simeq 5000$ K as the stellar mass increases. The resulting UV feedback is too weak to hinder mass accretion due to the low flux of stellar UV photons, thus verifying our implicit assumption of no stellar feedback during the hydrodynamic simulations. The insensitivity of stellar evolution to variable accretion is attributed to the fact that typical timescales of variability, $\lesssim 10^3$ years, are too short to affect the stellar structure. We argue that this evolution will continue until the SMS eventually collapses to produce a massive black hole by the general relativistic instability after the stellar mass reaches $\gtrsim 10^5~M_\odot$.

While observations of large-scale structure and the cosmic microwave background (CMB) provide strong constraints on the amplitude of the primordial power spectrum (PPS) on scales larger than 10~Mpc, the amplitude of the power spectrum on sub-galactic length scales is much more poorly constrained. We study early structure formation in a cosmological model with a blue-tilted PPS. We assume that the standard scale-invariant PPS is modified at small length scales as $P(k) \sim k^{m_{\rm s}}$ with $m_{\rm s} > 1$. We run a series of cosmological hydrodynamic simulations to examine the dependence of the formation epoch and the characteristic mass of primordial stars on the tilt of the PPS. In models with $m_{\rm s} > 1$, star-forming gas clouds are formed at $z > 100$, when formation of hydrogen molecules is inefficient because the intense CMB radiation destroys chemical intermediates. Without efficient coolant, the gas clouds gravitationally contract while keeping a high temperature. The protostars formed in such "hot" clouds grow very rapidly by accretion to become extremely massive stars that may leave massive black holes with a few hundred solar-masses at $z > 100$. The shape of the PPS critically affects the properties and the formation epoch of the first generation of stars. Future experiments of the CMB polarization and the spectrum distortion may provide important information on the nature of the first stars and their formation epoch, and hence on the shape of the small-scale power spectrum.

A promising formation channel of SMBHs at redshift 6 is the so-called DC model, which posits that a massive seed BH forms through gravitational collapse of a $\sim 10^5~M_\odot$ SMS. We study the evolution of such a SMS growing by rapid mass accretion. In particular, we examine the impact of time-dependent mass accretion of repeating burst and quiescent phases that are expected to occur with a self-gravitating circumstellar disk. We show that the stellar evolution with such episodic accretion differs qualitatively from that expected with a constant accretion rate, even if the mean accretion rate is the same. Unlike the case of constant mass accretion, whereby the star expands roughly following $R_* \simeq 2.6 \times 10^3 R_\odot (M_*/100~M_\odot)^{1/2}$, the protostar can substantially contract during the quiescent phases between accretion bursts. The stellar effective temperature and ionizing photon emissivity increase accordingly as the star contracts, which can cause strong ionizing feedback and halt the mass accretion onto the star. With a fixed duration of the quiescent phase $\Delta t_{\rm q}$, such contraction occurs in early evolutionary phases, i.e. for $M_* \lesssim 10^3~M_\odot$ with $\Delta t_{\rm q} \simeq 10^3$ yr. For later epochs and larger masses but the same $\Delta t_{\rm q}$, contraction is negligible even during quiescent phases. With larger quiescent times $\Delta t_{\rm q}$, however, the star continues to contract during quiescent phases even for the higher stellar masses. We show that such behavior is well understood by comparing the interval time and the thermal relaxation time for a bloated surface layer. We conclude that the UV radiative feedback becomes effective if the quiescent phase associated by the burst accretion is longer than $\sim 10^3$ yr, which is possible in an accretion disk forming in the direct collapse model.

We perform a large set of cosmological simulations of early structure formation and follow the formation and evolution of 1540 star-forming gas clouds to derive the mass distribution of primordial stars. The star formation in our cosmological simulations is characterized by two distinct populations, the so-called Population III.1 stars and primordial stars formed under the influence of far ultraviolet (FUV) radiation (Population III.2D stars). In this work, we determine the stellar masses by using the dependences on the physical properties of star-forming cloud and/or the external photodissociating intensity from nearby primordial stars, which are derived from the results of two-dimensional radiation hydrodynamic simulations of protostellar feedback. The characteristic mass of the Pop III stars is found to be a few hundred solar masses at z ~ 25, and it gradually shifts to lower masses with decreasing redshift. At high redshifts z > 20, about half of the star-forming gas clouds are exposed to intense FUV radiation and thus give birth to massive Pop III.2D stars. However, the local FUV radiation by nearby Pop III stars becomes weaker at lower redshifts, when typical Pop III stars have smaller masses and the mean physical separation between the stars becomes large owing to cosmic expansion. Therefore, at z < 20, a large fraction of the primordial gas clouds host Pop III.1 stars. At z =< 15, the Pop III.1 stars are formed in relatively cool gas clouds due to efficient radiative cooling by H_2 and HD molecules; such stars have masses of a few x 10 Msun. Since the stellar evolution and the final fate are determined by the stellar mass, Pop III stars formed at different epochs play different roles in the early universe.

We carry out radiation hydrodynamical simulations of the formation of massive stars in the super-Eddington regime including both their radiative feedback and protostellar outflows. The calculations start from a prestellar core of dusty gas and continue until the star stops growing. The accretion ends when the remnants of the core are ejected, mostly by the force of the direct stellar radiation in the polar direction and elsewhere by the reradiated thermal infrared radiation. How long the accretion persists depends on whether the protostellar outflows are present. We set the mass outflow rate to 1% of the stellar sink particle's accretion rate. The outflows open a bipolar cavity extending to the core's outer edge, through which the thermal radiation readily escapes. The radiative flux is funneled into the polar directions while the core's collapse proceeds near the equator. The outflow thus extends the "flashlight effect", or anisotropic radiation field, found in previous studies from the few hundred AU scale of the circumstellar disk up to the 0.1 parsec scale of the core. The core's flashlight effect allows core gas to accrete on the disk for longer, in the same way that the disk's flashlight effect allows disk gas to accrete on the star for longer. Thus although the protostellar outflows remove material near the core's poles, causing slower stellar growth over the first few free-fall times, they also enable accretion to go on longer in our calculations. The outflows ultimately lead to stars of somewhat higher mass.

The ratio of mass and magnetic flux determines the relative importance of magnetic and gravitational forces in the evolution of molecular clouds and their cores. Its measurement is thus central in discriminating between different theories of core formation and evolution. Here we discuss the effect of chemical depletion on measurements of the mass-to-flux ratio using the same molecule (OH) both for Zeeman measurements of the magnetic field and the determination of the mass of the region. The uncertainties entering through the OH abundance in determining separately the magnetic field and the mass of a region have been recognized in the literature. It has been proposed however that, when comparing two regions of the same cloud, the abundance will in both cases be the same. We show that this assumption is invalid. We demonstrate that when comparing regions with different densities, the effect of OH depletion in measuring changes of the mass-to-flux ratio between different parts of the same cloud can even reverse the direction of the underlying trends (for example, the mass-to-flux ratio may appear to decrease as we move to higher density regions). The systematic errors enter primarily through the inadequate estimation of the mass of the region.

We perform a large set of radiation hydrodynamics simulations of primordial star formation in a fully cosmological context. Our statistical sample of 100 First Stars show that the first generation of stars have a wide mass distribution M_popIII = 10 ~ 1000 M_sun. We first run cosmological simulations to generate a set of primordial star-forming gas clouds. We then follow protostar formation in each gas cloud and the subsequent protostellar evolution until the gas mass accretion onto the protostar is halted by stellar radiative feedback. The accretion rates differ significantly among the primordial gas clouds which largely determine the final stellar masses. For low accretion rates the growth of a protostar is self-regulated by radiative feedback effects and the final mass is limited to several tens of solar masses. At high accretion rates the protostar's outer envelope continues to expand and the effective surface temperature remains low; such protostars do not exert strong radiative feedback and can grow in excess to one hundred solar masses. The obtained wide mass range suggests that the first stars play a variety of roles in the early universe, by triggering both core-collapse supernovae and pair-instability supernovae as well as by leaving stellar mass black holes. We find certain correlations between the final stellar mass and the physical properties of the star-forming cloud. These correlations can be used to estimate the mass of the first star from the properties of the parent cloud or of the host halo, without following the detailed protostellar evolution.

Supermassive stars (SMSs) forming via very rapid mass accretion (Mdot >~ 0.1 Msun/yr) could be precursors of supermassive black holes observed beyond redshift of about 6. Extending our previous work, we here study the evolution of primordial stars growing under such rapid mass accretion until the stellar mass reaches 10^{4 - 5} Msun. Our stellar evolution calculations show that a star becomes supermassive while passing through the "supergiant protostar" stage, whereby the star has a very bloated envelope and a contracting inner core. The stellar radius increases monotonically with the stellar mass, until =~ 100 AU for M_* >~ 10^4 Msun, after which the star begins to slowly contract. Because of the large radius the effective temperature is always less than 10^4 K during rapid accretion. The accreting material is thus almost completely transparent to the stellar radiation. Only for M_* >~ 10^5 Msun can stellar UV feedback operate and disturb the mass accretion flow. We also examine the pulsation stability of accreting SMSs, showing that the pulsation-driven mass loss does not prevent stellar mass growth. Observational signatures of bloated SMSs should be detectable with future observational facilities such as the James Webb Space Telescope. Our results predict that an inner core of the accreting SMS should suffer from the general relativistic instability soon after the stellar mass exceeds 10^5 Msun. An extremely massive black hole should form after the collapse of the inner core.

Studies of the evolution of massive protostars and the evolution of their host molecular cloud cores are commonly treated as separate problems. However, interdependencies between the two can be significant. Here, we study the simultaneous evolution of massive protostars and their host molecular cores using a multi-dimensional radiation hydrodynamics code that incorporates the effects of the thermal pressure and radiative acceleration feedback of the centrally forming protostar. The evolution of the massive protostar is computed simultaneously using the stellar evolution code STELLAR, modified to include the effects of variable accretion. The interdependencies are studied in three different collapse scenarios. For comparison, stellar evolutionary tracks at constant accretion rates and the evolution of the host cores using pre-computed stellar evolutionary tracks are computed. The resulting interdependencies of the protostellar evolution and the evolution of the environment are extremely diverse and depend on the order of events, in particular the time of circumstellar accretion disk formation with respect to the onset of the bloating phase of the star. Feedback mechanisms affect the instantaneous accretion rate and the protostar's radius, temperature and luminosity on timescales equal or smaller than 5 kyr, corresponding to the accretion timescale and Kelvin-Helmholtz contraction timescale, respectively. Nevertheless, it is possible to approximate the overall protostellar evolution in many cases by pre-computed stellar evolutionary tracks assuming appropriate constant average accretion rates.

Numerical simulations have shown that the often cited radiation pressure barrier to accretion onto massive stars can be circumvented, when the radiation field is highly anisotropic in the presence of a circumstellar accretion disk with high optical depth. Here, these studies of the so-called flashlight effect are expanded by including the opacity of the innermost dust-free but potentially optically thick gas regions around forming massive stars. In addition to frequency-dependent opacities for the dust grains, we use temperature- and density-dependent Planck- and Rosseland mean opacities for the gas. The simulations show that the innermost dust-free parts of the accretion disks are optically thick to the stellar radiation over a substantial fraction of the solid angle above and below the disk's midplane. The temperature in the shielded disk region decreases faster with radius than in a comparison simulation with a lower constant gas opacity, and the dust sublimation front is shifted to smaller radii. The shielding by the dust-free gas in the inner disk thus contributes to an enhanced flashlight effect, which ultimately results in a smaller opening angle of the radiation pressure driven outflow and in a much longer timescale of sustained feeding of the circumstellar disk by the molecular cloud core. We conclude that it is necessary to properly account for the opacity of the inner dust-free disk regions around forming massive stars in order to correctly assess the effectiveness of the flashlight effect, the opening angle of radiation pressure driven outflows, and the lifetime and morphological evolution of the accretion disk.

The first stars in the universe ionized the ambient primordial gas through various feedback processes. "Second-generation" primordial stars potentially form from this disturbed gas after its recombination. In this Letter, we study the late formation stage of such second-generation stars, where a large amount of gas accretes onto the protostar and the final stellar mass is determined when the accretion terminates. We directly compute the complex interplay between the accretion flow and stellar ultraviolet (UV) radiation, performing radiation-hydrodynamic simulations coupled with stellar evolution calculations. Because of more efficient H2 and HD cooling in the pre-stellar stage, the accretion rates onto the star are ten times lower than in the case of the formation of the first stars. The lower accretion rates and envelope density result in the occurrence of an expanding bipolar HII region at a lower protostellar mass M_* \simeq 10Msun, which blows out the circumstellar material, thereby quenching the mass supply from the envelope to the accretion disk. At the same time the disk loses mass due to photoevaporation by the growing star. In our fiducial case the stellar UV feedback terminates mass accretion onto the star at M_* \simeq 17Msun. Although the derived masses of the second-generation primordial stars are systematically lower than those of the first generation, the difference is within a factor of only a few. Our results suggest a new scenario, whereby the majority of the primordial stars are born as massive stars with tens of solar masses, regardless of their generations.

We study the effect that non-equilibrium chemistry in dynamical models of collapsing molecular cloud cores has on measurements of the magnetic field in these cores, the degree of ionization, and the mean molecular weight of ions. We find that OH and CN, usually used in Zeeman observations of the line-of-sight magnetic field, have an abundance that decreases toward the center of the core much faster than the density increases. As a result, Zeeman observations tend to sample the outer layers of the core and consistently underestimate the core magnetic field. The degree of ionization follows a complicated dependence on the number density at central densities up to 10^5 cm^{-3} for magnetic models and 10^6 cm^{-3} in non-magnetic models. At higher central densities the scaling approaches a power-law with a slope of -0.6 and a normalization which depends on the cosmic-ray ionization rate {\zeta} and the temperature T as ({\zeta}T)^1/2. The mean molecular weight of ions is systematically lower than the usually assumed value of 20 - 30, and, at high densities, approaches a value of 3 due to the asymptotic dominance of the H3+ ion. This significantly lower value implies that ambipolar diffusion operates faster.

We combine dynamical and non-equilibrium chemical modeling of evolving prestellar molecular cloud cores, and explore the evolution of molecular abundances in the contracting core. We model both magnetic cores, with varying degrees of initial magnetic support, and non-magnetic cores, with varying collapse delay times. We explore, through a parameter study, the competing effects of various model parameters in the evolving molecular abundances, including the elemental C/O ratio, the temperature, and the cosmic-ray ionization rate. We find that different models show their largest quantitative differences at the center of the core, whereas the outer layers, which evolve slower, have abundances which are severely degenerate among different dynamical models. There is a large range of possible abundance values for different models at a fixed evolutionary stage (central density), which demonstrates the large potential of chemical differentiation in prestellar cores. However, degeneracies among different models, compounded with uncertainties induced by other model parameters, make it difficult to discriminate among dynamical models. To address these difficulties, we identify abundance ratios between particular molecules, the measurement of which would have maximal potential for discrimination among the different models examined here. In particular, we find that the ratios between NH3 and CO; NH2 and CO; NH3 and HCO+ are sensitive to the evolutionary timescale, and that the ratio between HCN and OH is sensitive to the C/O ratio. Finally, we demonstrate that measurements of the central deviation (central depletion or enhancement) of abundances of certain molecules are good indicators of the dynamics of the core.

Direct collapse of supermassive stars (SMSs) is a possible pathway for generating supermassive black holes in the early universe. It is expected that an SMS could form via very rapid mass accretion with Mdot ~ 0.1 - 1 Msun/yr during the gravitational collapse of an atomic-cooling primordial gas cloud. In this paper we study how stars would evolve under such extreme rapid mass accretion, focusing on the early evolution until the stellar mass reaches 1000 Msun. To this end we numerically calculate the detailed interior structure of accreting stars with primordial element abundances. Our results show that for accretion rates higher than 0.01 Msun/yr, stellar evolution is qualitatively different from that expected at lower rates. While accreting at these high rates the star always has a radius exceeding 100 Rsun, which increases monotonically with the stellar mass. The mass-radius relation for stellar masses exceeding ~ 100 Msun follows the same track with R_* \propto M_*^0.5 in all cases with accretion rates > 0.01 Msun/yr; at a stellar mass of 1000 Msun the radius is about 7000 Rsun (~= 30 AU). With higher accretion rates the onset of hydrogen burning is shifted towards higher stellar masses. In particular, for accretion rates exceeding Mdot > 0.1 Msun/yr, there is no significant hydrogen burning even after 1000 Msun have accreted onto the protostar. Such "supergiant" protostars have effective temperatures as low as Teff ~= 5000 K throughout their evolution and because they hardly emit ionizing photons, they do not create an HII region or significantly heat their immediate surroundings. Thus, radiative feedback is unable to hinder the growth of rapidly accreting stars to masses in excess of 1000 Msun, as long as material is accreted at rates Mdot > 0.01 Msun/yr.

The first stars fundamentally transformed the early universe by emitting the first light and by producing the first heavy elements. These effects were predetermined by the mass distribution of the first stars, which is thought to have been fixed by a complex interplay of gas accretion and protostellar radiation. We performed radiation-hydrodynamics simulations that followed the growth of a primordial protostar through to the early stages as a star with thermo-nuclear burning. The circumstellar accretion disk was evaporated by ultraviolet radiation from the star when its mass was 43 times that of the Sun. Such massive primordial stars, in contrast to the often postulated extremely massive stars, may help explain the fact that there are no signatures of the pair-instability supernovae in abundance patterns of metal-poor stars in our galaxy.

We propose a simple analytical method for estimating the central volume density of prestellar molecular cloud cores from their column density profiles. Prestellar cores feature a flat central part of the column density and volume density profiles of the same size indicating the existence of a uniform density inner region. The size of this region is set by the thermal pressure force which depends only on the central volume density and temperature of the core, and can provide a direct measurement of the central volume density. Thus a simple length measurement can immediately yield a central density estimate independent of any dynamical model for the core and without the need for fitting. Using the radius at which the column density is 90% of the central value as an estimate of the size of the flat inner part of the column density profile yields an estimate of the central volume density within a factor of 2 for well resolved cores.

We have examined methanol emission from Orion KL with of the {\em Herschel}/HIFI instrument, and detected two methanol bands centered at 524 GHz and 1061 GHz. The 524 GHz methanol band (observed in HIFI band 1a) is dominated by the isolated $\Delta$J$=$0, K$=-4\rightarrow$-3, v$_t$$=$0 Q branch, and includes 25 E-type and 2 A-type transitions. The 1061 GHz methanol band (observed in HIFI band 4b) is dominated by the $\Delta$J$=$0, K$=7\rightarrow$6, v$_t$$=$0 Q branch transitions which are mostly blended. We have used the isolated E-type v$_t$$=$0 methanol transitions to explore the physical conditions in the molecular gas. With HIFI's high velocity resolution, the methanol emission contributed by different spatial components along the line of sight toward Orion KL (hot core, low velocity flow, and compact ridge) can be distinguished and studied separately. The isolated transitions detected in these bands cover a broad energy range (upper state energy ranging from 80 K to 900 K), which provides a unique probe of the thermal structure in each spatial component. The observations further show that the compact ridge is externally heated. These observations demonstrate the power of methanol lines as probes of the physical conditions in warm regions in close proximity to young stars.

We report observations of the reactive molecular ions OH$^+$, H$_2$O$^+$, and H$_3$O$^+$ towards Orion KL with Herschel/HIFI. All three $N=1-0$ fine-structure transitions of OH$^+$ at 909, 971, and 1033GHz and both fine-structure components of the doublet {\it ortho}-H$_2$O$^+$ $1_{11}-0_{00}$ transition at 1115 and 1139GHz were detected; an upper limit was obtained for H$_3$O$^+$. OH$^+$ and H$_2$O$^+$ are observed purely in absorption, showing a narrow component at the source velocity of 9 kms$^{-1}$, and a broad blueshifted absorption similar to that reported recently for HF and {\it para}-H$_{2}^{18}$O, and attributed to the low velocity outflow of Orion KL. We estimate column densities of OH$^+$ and H$_2$O$^+$ for the 9 km s$^{-1}$ component of $9 \pm 3 \times 10^{12}$cm$^{-2}$ and $7 \pm 2 \times 10^{12}$cm$^{-2}$, and those in the outflow of $1.9 \pm 0.7 \times 10^{13}$cm$^{-2}$ and $1.0 \pm 0.3 \times 10^{13}$cm$^{-2}$. Upper limits of $2.4\times 10^{12}$cm$^{-2}$ and $8.7\times 10^{12}$cm$^{-2}$ were derived for the column densities of {\it ortho} and {\it para}-H$_3$O$^+$ from transitions near 985 and 1657GHz. The column densities of the three ions are up to an order of magnitude lower than those obtained from recent observations of W31C and W49N. The comparatively low column densities may be explained by a higher gas density despite the assumption of a very high ionization rate.

We present results from numerical simulations of the cooling-core cluster A2199 produced by the two-dimensional (2-D) resistive magnetohydrodynamics (MHD) code MACH2. In our simulations we explore the effect of anisotropic thermal conduction on the energy balance of the system. The results from idealized cases in 2-D axisymmetric geometry underscore the importance of the initial plasma density in ICM simulations, especially the near-core values since the radiation cooling rate is proportional to ${n_e}^2$. Heat conduction is found to be non-effective in preventing catastrophic cooling in this cluster. In addition we performed 2-D planar MHD simulations starting from initial conditions deliberately violating both thermal balance and hydrostatic equilibrium in the ICM, to assess contributions of the convective terms in the energy balance of the system against anisotropic thermal conduction. We find that in this case work done by the pressure on the plasma can dominate the early evolution of the internal energy over anisotropic thermal conduction in the presence of subsonic flows, thereby reducing the impact of the magnetic field. Deviations from hydrostatic equilibrium near the cluster core may be associated with transient activity of a central active galactic nucleus and/or remnant dynamical activity in the ICM and warrant further study in three dimensions.

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10^-6 is placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display asymmetric profiles, which may arise from infall around the hot core.

We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence on the illuminating UV field. Using Herschel-HIFI observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. Our approach consists in determining the energetics of the region by combining the information carried by the mid-IR spectrum (extinction by classical grains, emission from very small dust particles) with that of the main gas coolant lines. In this letter, we discuss more specifically the intensity and line profile of the 158 micron (1901 GHz) [CII] line measured by HIFI and provide information on the emitting gas. We show that both the [CII] emission and the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) arise from the regions located in the transition zone between atomic and molecular gas. Using the Meudon PDR code and a simple transfer model, we find good agreement between the calculated and observed [CII] intensities. HIFI observations of NGC 7023 provide the opportunity to constrain the energetics at the surface of PDRs. Future work will include analysis of the main coolant line [OI] and use of a new PDR model that includes PAH-related species.

We report the detection of more than 48 velocity-resolved ground rotational state transitions of H2(16)O, H2(18)O, and H2(17)O - most for the first time - in both emission and absorption toward Orion KL using Herschel/HIFI. We show that a simple fit, constrained to match the known emission and absorption components along the line of sight, is in excellent agreement with the spectral profiles of all the water lines. Using the measured H2(18)O line fluxes, which are less affected by line opacity than their H2(16)O counterparts, and an escape probability method, the column densities of H2(18)O associated with each emission component are derived. We infer total water abundances of 7.4E-5, 1.0E-5, and 1.6E-5 for the plateau, hot core, and extended warm gas, respectively. In the case of the plateau, this value is consistent with previous measures of the Orion-KL water abundance as well as those of other molecular outflows. In the case of the hot core and extended warm gas, these values are somewhat higher than water abundances derived for other quiescent clouds, suggesting that these regions are likely experiencing enhanced water-ice sublimation from (and reduced freeze-out onto) grain surfaces due to the warmer dust in these sources.

Observations of HDO are an important complement for studies of water, because they give strong constraints on the formation processes -- grain surfaces versus energetic process in the gas phase, e.g. in shocks. The HIFI observations of multiple transitions of HDO in Sgr~B2(M) presented here allow the determination of the HDO abundance throughout the envelope, which has not been possible before with ground-based observations only. The abundance structure has been modeled with the spherical Monte Carlo radiative transfer code RATRAN, which also takes radiative pumping by continuum emission from dust into account. The modeling reveals that the abundance of HDO rises steeply with temperature from a low abundance ($2.5\times 10^{-11}$) in the outer envelope at temperatures below 100~K through a medium abundance ($1.5\times 10^{-9}$) in the inner envelope/outer core, at temperatures between 100 and 200~K, and finally a high abundance ($3.5\times 10^{-9}$) at temperatures above 200~K in the hot core.

To investigate the accretion and feedback processes in massive star formation, we analyze the shapes of emission lines from hot molecular cores, whose asymmetries trace infall and expansion motions. The high-mass star forming region SgrB2(M) was observed with Herschel/HIFI (HEXOS key project) in various lines of HCN and its isotopologues, complemented by APEX data. The observations are compared to spherically symmetric, centrally heated models with density power-law gradient and different velocity fields (infall or infall+expansion), using the radiative transfer code RATRAN. The HCN line profiles are asymmetric, with the emission peak shifting from blue to red with increasing J and decreasing line opacity (HCN to H$^{13}$CN). This is most evident in the HCN 12--11 line at 1062 GHz. These line shapes are reproduced by a model whose velocity field changes from infall in the outer part to expansion in the inner part. The qualitative reproduction of the HCN lines suggests that infall dominates in the colder, outer regions, but expansion dominates in the warmer, inner regions. We are thus witnessing the onset of feedback in massive star formation, starting to reverse the infall and finally disrupting the whole molecular cloud. To obtain our result, the THz lines uniquely covered by HIFI were critically important.

Using Herschel's HIFI instrument we have observed [C II] along a cut through S140 and high-J transitions of CO and HCO+ at two positions on the cut, corresponding to the externally irradiated ionization front and the embedded massive star forming core IRS1. The HIFI data were combined with available ground-based observations and modeled using the KOSMA-tau model for photon dominated regions. Here we derive the physical conditions in S140 and in particular the origin of [C II] emission around IRS1. We identify three distinct regions of [C II] emission from the cut, one close to the embedded source IRS1, one associated with the ionization front and one further into the cloud. The line emission can be understood in terms of a clumpy model of photon-dominated regions. At the position of IRS1, we identify at least two distinct components contributing to the [C II] emission, one of them a small, hot component, which can possibly be identified with the irradiated outflow walls. This is consistent with the fact that the [C II] peak at IRS1 coincides with shocked H2 emission at the edges of the outflow cavity. We note that previously available observations of IRS1 can be well reproduced by a single-component KOSMA-tau model. Thus it is HIFI's unprecedented spatial and spectral resolution, as well as its sensitivity which has allowed us to uncover an additional hot gas component in the S140 region.

The [CII] fine--structure line at 158um is an excellent tracer of the warm diffuse gas in the ISM and the interfaces between molecular clouds and their surrounding atomic and ionized envelopes. Here we present the initial results from Galactic Observations of Terahertz C+ (GOTC+), a Herschel Key Project devoted to study the [CII] fine structure emission in the galactic plane using the HIFI instrument. We use the [CII] emission together with observations of CO as a probe to understand the effects of newly--formed stars on their interstellar environment and characterize the physical and chemical state of the star-forming gas. We collected data along 16 lines--of--sight passing near star forming regions in the inner Galaxy near longitudes 330 degrees and 20 degrees. We identify fifty-eight [CII] components that are associated with high--column density molecular clouds as traced by 13CO emission. We combine [CII], 12CO, and 13CO observations to derive the physical conditions of the [CII]--emitting regions in our sample of high--column density clouds based on comparison with results from a grid of Photon Dominated Region (PDR) models. From this unbiased sample, our results suggest that most of [CII] emission originates from clouds with H2 volume densities between 10e3.5 and 10e5.5 cm^-3 and weak FUV strength (CHI_0=1-10). We find two regions where our analysis suggests high densities >10e5 cm^-3 and strong FUV fields (CHI=10e4-10e6), likely associated with massive star formation. We suggest that [CII] emission in conjunction with CO isotopes is a good tool to differentiate between regions of massive star formation (high densities/strong FUV fields) and regions that are distant from massive stars (lower densities/weaker FUV fields) along the line--of--sight

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.

We present the first high spectral resolution observations of Orion KL in the frequency ranges 1573.4 - 1702.8 GHz (band 6b) and 1788.4 - 1906.8 GHz (band 7b) obtained using the HIFI instrument on board the Herschel Space Observatory. We characterize the main emission lines found in the spectrum, which primarily arise from a range of components associated with Orion KL including the hot core, but also see widespread emission from components associated with molecular outflows traced by H2O, SO2, and OH. We find that the density of observed emission lines is significantly diminished in these bands compared to lower frequency Herschel/HIFI bands.

We present observations of twelve rotational transitions of H2O-16, H2O-18, and H2O-17 toward the massive star-forming region NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical conditions in these components. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of 200 K, whereas water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 4 10^-5 (outflow) and 10^-8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust grains. The H2O-18/H2O-17 abundance ratio is 3.2, which is consistent with the O-18/O-17 ratio determined from CO isotopologues. The ortho/para ratio in water appears to be relatively low 1.6(1) in the cold, quiescent gas, but close to the equilibrium value of three in the warmer outflow material (2.5(0.8)).

We present the first results on the diffuse transition clouds observed in [CII] line emission at 158 microns (1.9 THz) towards Galactic longitudes near 340deg (5 LOSs) and 20deg (11 LOSs) as part of the GOT C+ survey. Out of the total 146 [CII] velocity components detected by profile fitting we identify 53 as diffuse molecular clouds with associated $^{12}$CO emission but without $^{13}$CO emission and characterized by A$_V$ < 5 mag. We estimate the fraction of the [CII] emission in the diffuse HI layer in each cloud and then determine the [CII] emitted from the molecular layers in the cloud. We show that the excess [CII] intensities detected in a few clouds is indicative of a thick H$_2$ layer around the CO core. The wide range of clouds in our sample with thin to thick H$_2$ layers suggests that these are at various evolutionary states characterized by the formation of H$_2$ and CO layers from HI and C$^+$, respectively. In about 30% of the clouds the H$_2$ column densities (''dark gas'') traced by the [CII] is 50% or more than that traced by $^{12}$CO emission. On the average about 25% of the total H$_2$ in these clouds is in an H$_2$ layer which is not traced by CO. We use the HI, [CII], and $^{12}$CO intensities in each cloud along with simple chemical models to obtain constraints on the FUV fields and cosmic ray ionization rates.

We present the first results of the Herschel open time key program, Galactic Observations of Terahertz C$^+$ (GOT C+) survey of the [CII] fine-structure line at 1.9 THz (158 microns) using the HIFI instrument on Herschel. We detected 146 interstellar clouds along sixteen lines-of-sight towards the inner Galaxy. We also acquired HI and CO isotopologue data along each line-of-sight for analysis of the physical conditions in these clouds. Here we analyze 29 diffuse clouds (A$_{V}$ < 1.3 mag.) in this sample characterized by having [CII] and HI emission, but no detectable CO. We find that [CII] emission is generally stronger than expected for diffuse atomic clouds, and in a number of sources is much stronger than anticipated based on their HI column density. We show that excess [CII] emission in these clouds is best explained by the presence of a significant diffuse warm H$_2$, dark gas, component. This first [CII] 158 micron detection of warm dark gas demonstrates the value of this tracer for mapping this gas throughout the Milky Way and in galaxies.

We report a detection of the fundamental rotational transition of hydrogen fluoride in absorption towards Orion KL using Herschel/HIFI. After the removal of contaminating features associated with common molecules ("weeds"), the HF spectrum shows a P-Cygni profile, with weak redshifted emission and strong blue-shifted absorption, associated with the low-velocity molecular outflow. We derive an estimate of 2.9 x 10^13 cm^-2 for the HF column density responsible for the broad absorption component. Using our best estimate of the H2 column density within the low-velocity molecular outflow, we obtain a lower limit of ~1.6 x 10^-10 for the HF abundance relative to hydrogen nuclei, corresponding to 0.6% of the solar abundance of fluorine. This value is close to that inferred from previous ISO observations of HF J=2--1 absorption towards Sgr B2, but is in sharp contrast to the lower limit of 6 x 10^-9 derived by Neufeld et al. (2010) for cold, foreground clouds on the line of sight towards G10.6-0.4.

We present initial results from the Herschel GT key program: Herschel observations of EXtra-Ordinary Sources (HEXOS) and outline the promise and potential of spectral surveys with Herschel/HIFI. The HIFI instrument offers unprecedented sensitivity, as well as continuous spectral coverage across the gaps imposed by the atmosphere, opening up a largely unexplored wavelength regime to high-resolution spectroscopy. We show the spectrum of Orion KL between 480 and 560 GHz and from 1.06 to 1.115 THz. From these data, we confirm that HIFI separately measures the dust continuum and spectrally resolves emission lines in Orion KL. Based on this capability we demonstrate that the line contribution to the broad-band continuum in this molecule-rich source is ~20-40% below 1 THz and declines to a few percent at higher frequencies. We also tentatively identify multiple transitions of HD18O in the spectra. The first detection of this rare isotopologue in the interstellar medium suggests that HDO emission is optically thick in the Orion hot core with HDO/H2O ~ 0.02. We discuss the implications of this detection for the water D/H ratio in hot cores.

We have observed CH absorption lines ($J=3/2, N=1 \leftarrow J=1/2, N=1$) against the continuum source Sgr~B2(M) using the \textit{Herschel}/HIFI instrument. With the high spectral resolution and wide velocity coverage provided by HIFI, 31 CH absorption features with different radial velocities and line widths are detected and identified. The narrower line width and lower column density clouds show `spiral arm' cloud characteristics, while the absorption component with the broadest line width and highest column density corresponds to the gas from the Sgr~B2 envelope. The observations show that each `spiral arm' harbors multiple velocity components, indicating that the clouds are not uniform and that they have internal structure. This line-of-sight through almost the entire Galaxy offers unique possibilities to study the basic chemistry of simple molecules in diffuse clouds, as a variety of different cloud classes are sampled simultaneously. We find that the linear relationship between CH and H$_2$ column densities found at lower $A_V$ by UV observations does not continue into the range of higher visual extinction. There, the curve flattens, which probably means that CH is depleted in the denser cores of these clouds.

Mon R2, at a distance of 830 pc, is the only ultracompact HII region (UC HII) where the photon-dominated region (PDR) between the ionized gas and the molecular cloud can be resolved with Herschel. HIFI observations of the abundant compounds 13CO, C18O, o-H2-18O, HCO+, CS, CH, and NH have been used to derive the physical and chemical conditions in the PDR, in particular the water abundance. The 13CO, C18O, o-H2-18O, HCO+ and CS observations are well described assuming that the emission is coming from a dense (n=5E6 cm-3, N(H2)>1E22 cm-2) layer of molecular gas around the UC HII. Based on our o-H2-18O observations, we estimate an o-H2O abundance of ~2E-8. This is the average ortho-water abundance in the PDR. Additional H2-18O and/or water lines are required to derive the water abundance profile. A lower density envelope (n~1E5 cm-3, N(H2)=2-5E22 cm-2) is responsible for the absorption in the NH 1_1-0_2 line. The emission of the CH ground state triplet is coming from both regions with a complex and self-absorbed profile in the main component. The radiative transfer modeling shows that the 13CO and HCO+ line profiles are consistent with an expansion of the molecular gas with a velocity law, v_e =0.5 x (r/Rout)^{-1} km/s, although the expansion velocity is poorly constrained by the observations presented here.

We report the first detection of the ground-state rotational transition of the methylidyne cation CH+ towards the massive star-forming region DR21 with the HIFI instrument onboard the Herschel satellite. The line profile exhibits a broad emission line, in addition to two deep and broad absorption features associated with the DR21 molecular ridge and foreground gas. These observations allow us to determine a CH+ J=1-0 line frequency of 835137 +/- 3 MHz, in good agreement with a recent experimental determination. We estimate the CH+ column density to be a few 1e13 cm^-2 in the gas seen in emission, and > 1e14 cm^-2 in the components responsible for the absorption, which is indicative of a high line of sight average abundance [CH+]/[H] > 1.2x10^-8. We show that the CH+ column densities agree well with the predictions of state-of-the-art C-shock models in dense UV-illuminated gas for the emission line, and with those of turbulent dissipation models in diffuse gas for the absorption lines.

We present Herschel/HIFI observations of the fundamental rotational transitions of ortho- and para-H$_2^{16}$O and H$_2^{18}$O in absorption towards Sagittarius~B2(M) and W31C. The ortho/para ratio in water in the foreground clouds on the line of sight towards these bright continuum sources is generally consistent with the statistical high-temperature ratio of 3, within the observational uncertainties. However, somewhat unexpectedly, we derive a low ortho/para ratio of $2.35 \pm 0.35$, corresponding to a spin temperature of $\sim$27~K, towards Sagittarius~B2(M) at velocities of the expanding molecular ring. Water molecules in this region appear to have formed with, or relaxed to, an ortho/para ratio close to the value corresponding to the local temperature of the gas and dust.

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

We present observations of the rotational ortho-water ground transition, the two lowest para-water transitions, and the ground transition of ionised ortho-water in the archetypal starburst galaxy M82, performed with the HIFI instrument on the Herschel Space Observatory. These observations are the first detections of the para-H2O(111-000) (1113\,GHz) and ortho-H2O+(111-000) (1115\,GHz) lines in an extragalactic source. All three water lines show different spectral line profiles, underlining the need for high spectral resolution in interpreting line formation processes. Using the line shape of the para-H2O(111-000) and ortho-H2O+(111-000) absorption profile in conjunction with high spatial resolution CO observations, we show that the (ionised) water absorption arises from a ~2000 pc^2 region within the HIFI beam located about ~50 pc east of the dynamical centre of the galaxy. This region does not coincide with any of the known line emission peaks that have been identified in other molecular tracers, with the exception of HCO. Our data suggest that water and ionised water within this region have high (up to 75%) area-covering factors of the underlying continuum. This indicates that water is not associated with small, dense cores within the ISM of M82 but arises from a more widespread diffuse gas component.

H2O+ has been observed in its ortho- and para- states toward the massive star forming core Sgr B2(M), located close to the Galactic center. The observations show absorption in all spiral arm clouds between the Sun and Sgr B2. The average o/p ratio of H2O+ in most velocity intervals is 4.8, which corresponds to a nuclear spin temperature of 21 K. The relationship of this spin temperature to the formation temperature and current physical temperature of the gas hosting H2O+ is discussed, but no firm conclusion is reached. In the velocity interval 0 to 60 km/s, an ortho/para ratio of below unity is found, but if this is due to an artifact of contamination by other species or real is not clear.

We present high resolution HIFI spectroscopy of the nucleus of the archetypical starburst galaxy M82. Six 12CO lines, 2 13CO lines and 4 fine-structure lines are detected. Besides showing the effects of the overall velocity structure of the nuclear region, the line profiles also indicate the presence of multiple components with different optical depths, temperatures and densities in the observing beam. The data have been interpreted using a grid of PDR models. It is found that the majority of the molecular gas is in low density (n=10^3.5 cm^-3) clouds, with column densities of N_H=10^21.5 cm^-2 and a relatively low UV radiation field (GO = 10^2). The remaining gas is predominantly found in clouds with higher densities (n=10^5 cm^-3) and radiation fields (GO = 10^2.75), but somewhat lower column densities (N_H=10^21.2 cm^-2). The highest J CO lines are dominated by a small (1% relative surface filling) component, with an even higher density (n=10^6 cm^-3) and UV field (GO = 10^3.25). These results show the strength of multi-component modeling for the interpretation of the integrated properties of galaxies.

Recent observations of column densities in molecular clouds find lognormal distributions with power-law high-density tails. These results are often interpreted as indications that supersonic turbulence dominates the dynamics of the observed clouds. We calculate and present the column-density distributions of three clouds, modeled with very different techniques, none of which is dominated by supersonic turbulence. The first star-forming cloud is simulated using smoothed particle hydrodynamics (SPH); in this case gravity, opposed only by thermal-pressure forces, drives the evolution. The second cloud is magnetically subcritical with subsonic turbulence, simulated using nonideal MHD; in this case the evolution is due to gravitationally-driven ambipolar diffusion. The third cloud is isothermal, self-gravitating, and has a smooth density distribution analytically approximated with a uniform inner region and an r^-2 profile at larger radii. We show that in all three cases the column-density distributions are lognormal. Power-law tails develop only at late times (or, in the case of the smooth analytic profile, for strongly centrally concentrated configurations), when gravity dominates all opposing forces. It therefore follows that lognormal column-density distributions are generic features of diverse model clouds, and should not be interpreted as being a consequence of supersonic turbulence.

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

Mass accretion onto (proto-)stars at high accretion rates > 10^-4 M_sun/yr is expected in massive star formation. We study the evolution of massive protostars at such high rates by numerically solving the stellar structure equations. In this paper we examine the evolution via disk accretion. We consider a limiting case of "cold" disk accretion, whereby most of the stellar photosphere can radiate freely with negligible backwarming from the accretion flow, and the accreting material settles onto the star with the same specific entropy as the photosphere. We compare our results to the calculated evolution via spherically symmetric accretion, the opposite limit, whereby the material accreting onto the star contains the entropy produced in the accretion shock front. We examine how different accretion geometries affect the evolution of massive protostars. For cold disk accretion at 10^-3 M_sun/yr the radius of a protostar is initially small, about a few R_sun. After several solar masses have accreted, the protostar begins to bloat up and for M \simeq 10 M_sun the stellar radius attains its maximum of 30 - 400 R_sun. The large radius about 100 R_sun is also a feature of spherically symmetric accretion at the same accreted mass and accretion rate. Hence, expansion to a large radius is a robust feature of accreting massive protostars. At later times the protostar eventually begins to contract and reaches the Zero-Age Main-Sequence (ZAMS) for M \simeq 30 M_sun, independent of the accretion geometry. For accretion rates exceeding several 10^-3 M_sun/yr the protostar never contracts to the ZAMS. The very large radius of several 100s R_sun results in a low effective temperature and low UV luminosity of the protostar. Such bloated protostars could well explain the existence of bright high-mass protostellar objects, which lack detectable HII regions.

The molecular gas in the DR21 massive star formation region is known to be affected by the strong UV field from the central star cluster and by a fast outflow creating a bright shock. The relative contribution of both heating mechanisms is the matter of a long debate. By better sampling the excitation ladder of various tracers we provide a quantitative distinction between the different heating mechanisms. HIFI observations of mid-J transitions of CO and HCO+ isotopes allow us to bridge the gap in excitation energies between observations from the ground, characterizing the cooler gas, and existing ISO LWS spectra, constraining the properties of the hot gas. Comparing the detailed line profiles allows to identify the physical structure of the different components. In spite of the known shock-excitation of H2 and the clearly visible strong outflow, we find that the emission of all lines up to > 2 THz can be explained by purely radiative heating of the material. However, the new Herschel/HIFI observations reveal two types of excitation conditions. We find hot and dense clumps close to the central cluster, probably dynamically affected by the outflow, and a more widespread distribution of cooler, but nevertheless dense, molecular clumps.

We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum background and allows us to derive the velocity profile of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the frequency of the transition and the conditions for its formation. In DR21, the velocity distribution of H2O+ matches that of the [CII] line at 158\mu\m and of OH cm-wave absorption, both stemming from the hot and dense clump surfaces facing the HII-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H2O+ column density of 7.2e12 cm^-2 in NGC 6334, 2.3e13 cm^-2 in DR21, and 1.1e15 cm^-2 in Sgr B2.

Although fundamental for astrophysics, the processes that produce massive stars are not well understood. Large distances, high extinction, and short timescales of critical evolutionary phases make observations of these processes challenging. Lacking good observational guidance, theoretical models have remained controversial. This review offers a basic description of the collapse of a massive molecular core and a critical discussion of the three competing concepts of massive star formation: - monolithic collapse in isolated cores - competitive accretion in a protocluster environment - stellar collisions and mergers in very dense systems We also review the observed outflows, multiplicity, and clustering properties of massive stars, the upper initial mass function and the upper mass limit. We conclude that high-mass star formation is not merely a scaled-up version of low-mass star formation with higher accretion rates, but partly a mechanism of its own, primarily owing to the role of stellar mass and radiation pressure in controlling the dynamics.

By providing sensitive sub-arcsecond images and integral field spectroscopy in the 25 - 400 micron wavelength range, a far-IR interferometer will revolutionize our understanding of planetary system formation, reveal otherwise-undetectable planets through the disk perturbations they induce, and spectroscopically probe the atmospheres of extrasolar giant planets in orbits typical of most of the planets in our solar system. The technical challenges associated with interferometry in the far-IR are greatly relaxed relative to those encountered at shorter wavelengths or when starlight nulling is required. A structurally connected far-IR interferometer with a maximum baseline length of 36 m can resolve the interesting spatial structures in nascent and developed exoplanetary systems and measure exozodiacal emission at a sensitivity level critical to TPF-I mission planning. The Space Infrared Interferometric Telescope was recommended in the Community Plan for Far-IR/Submillimeter Space Astronomy, studied as a Probe-class mission, and estimated to cost 800M dollars. The scientific communities in Europe, Japan, and Canada have also demonstrated a keen interest in far-IR interferometry through mission planning workshops and technology research, suggesting the possibility of an international collaborative effort.

We show that the optically-thick dusty envelopes surrounding young high-mass stars are subject to the photon bubble instability. The infrared radiation passing through the envelope amplifies magnetosonic disturbances, with growth rates in our local numerical radiation MHD calculations that are consistent with a linear analysis. Modes with wavelengths comparable to the gas pressure scale height grow by more than two orders of magnitude in a thousand years, reaching non-linear amplitudes within the envelope lifetime. If the magnetic pressure in the envelope exceeds the gas pressure, the instability develops into trains of propagating shocks. Radiation escapes readily through the low-density material between the shocks, enabling accretion to continue despite the Eddington limit imposed by the dust opacity. The supersonic motions arising from the photon bubble instability can help explain the large velocity dispersions of hot molecular cores, while conditions in the shocked gas are suitable for maser emission. We conclude that the photon bubble instability may play a key role in the formation of massive stars.

We continue our numerical analysis of the morphological and energetic influence of massive stars on their ambient interstellar medium for a 35 solar mass star that evolves from the main sequence through red supergiant and Wolf-Rayet phases, until it ultimately explodes as a supernova. We find that structure formation in the circumstellar gas during the early main-sequence evolution occurs as in the 60 solar mass case but is much less pronounced because of the lower mechanical wind luminosity of the star. Since on the other hand the shell-like structure of the HII region is largely preserved, effects that rely on this symmetry become more important. At the end of the stellar lifetime 1% of the energy released as Lyman continuum radiation and stellar wind has been transferred to the circumstellar gas. From this fraction 10% is kinetic energy of bulk motion, 36% is thermal energy, and the remaining 54% is ionization energy of hydrogen. The sweeping up of the slow red supergiant wind by the fast Wolf-Rayet wind produces remarkable morphological structures and emission signatures, which are compared with existing observations of the Wolf-Rayet bubble S308. Our model reproduces the correct order of magnitude of observed X-ray luminosity, the temperature of the emitting plasma as well as the limb brightening of the intensity profile. This is remarkable, because current analytical and numerical models of Wolf-Rayet bubbles fail to consistently explain these features. A key result is that almost the entire X-ray emission in this stage comes from the shell of red supergiant wind swept up by the shocked Wolf-Rayet wind rather than from the shocked Wolf-Rayet wind itself as hitherto assumed and modeled. This offers a possible solution to what is called the ``missing wind problem'' of Wolf-Rayet bubbles.

The effects of turbulence on the mixing of gases and dust in the outer Solar nebula are examined using 3-D MHD calculations in the shearing-box approximation with vertical stratification. The turbulence is driven by the magneto-rotational instability. The magnetic and hydrodynamic stresses in the turbulence correspond to an accretion time at the midplane about equal to the lifetimes of T Tauri disks, while accretion in the surface layers is thirty times faster. The mixing resulting from the turbulence is also fastest in the surface layers. The mixing rate is similar to the rate of radial exchange of orbital angular momentum, so that the Schmidt number is near unity. The vertical spreading of a trace species is well-matched by solutions of a damped wave equation when the flow is horizontally-averaged. The damped wave description can be used to inexpensively treat mixing in 1-D chemical models. However, even in calculations reaching a statistical steady state, the concentration at any given time varies substantially over horizontal planes, due to fluctuations in the rate and direction of the transport. In addition to mixing species that are formed under widely varying conditions, the turbulence intermittently forces the nebula away from local chemical equilibrium. The different transport rates in the surface layers and interior may affect estimates of the grain evolution and molecular abundances during the formation of the Solar system.

Spectroscopic studies indicate that gas in the photospheres of young O stars moves at speeds up to the sound speed. We show, using two-dimensional radiation MHD calculations and results from a local linear analysis, that the motions may be due to photon bubble instability if young O stars have magnetic fields.

Current thinking about the Moon as a destination has revitalized interest in lunar astronomical observatories. Once seen by a large scientific community as a highly enabling site, the dramatic improvement in capabilities for free-space observatories prompts reevaluation of this interest. Whereas the lunar surface offers huge performance advantages for astronomy over terrestrial sites, free-space locales such as Earth orbit or Lagrange points offer performance that is superior to what could be achieved on the Moon. While astronomy from the Moon may be cost effective once infrastructure is there, it is in many respects no longer clearly enabling compared to free-space.

The collapse of massive molecular clumps can produce high mass stars, but the evolution is not simply a scaled-up version of low mass star formation. Outflows and radiative effects strongly hinder the formation of massive stars via accretion. A necessary condition for accretion growth of a hydrostatic object up to high masses M > 20 M_sun (rather than coalescence of optically thick objects) is the formation of and accretion through a circumstellar disk. Once the central object has accreted approximately 10 M_sun it has already evolved to core hydrogen-burning; the resultant main sequence star continues to accrete material as it begins to photoevaporate its circumstellar disk (and any nearby disks) on a timescale of 100,000 years, similar to the accretion timescale. Until the disk(s) is (are) completely photoevaporated, this configuration is observable as an ultra-compact HII region (UCHII). The final mass of the central star (and any nearby neighboring systems) is determined by the interplay between radiation acceleration, UV photoevaporation, stellar winds and outflows, and the accretion through the disk. Several aspects of this evolutionary sequence have been simulated numerically, resulting in a "proof of concept". This scenario places strong constraints on the accretion rate necessary to produce high mass stars and offers an opportunity to test the accretion hypothesis.

We consider the atmospheric flow on short-period extra-solar planets through two-dimensional numerical simulations of hydrodynamics with radiation transfer. One side is always exposed to the irradiation from the host star. The other is always in shadow. The temperature of the day side is determined by the equilibrium which the planetary atmosphere establishes with stellar radiation. Part of the thermal energy deposited on the day side is advected to the night side by a current. The radiation transfer, the night-side temperature distribution and by this the spectroscopic signature of the planet are sensitive functions of the atmospheric opacity. If the atmosphere contains grains with an abundanceand size distribution comparable to that of the interstellar medium, shallow heating occurs on the day side and the night side cools well below the day side. The temperature difference decreases as the abundance of grains is reduced. A simple analytic model of the dissipation of the circulation flow and associated kinetic heating is considered. This heating effect occurs mostly near the photosphere, not deep enough to significantly affect the size of planets. We show that the surface irradiation suppresses convection near the photospheric region on the day side. In some cases convection zones appear near the surface on the night side. This structural modification may influence the response and dissipation of tidal disturbances and alter the circularization and synchronization time scales.

We present results of numerical simulations carried out with a 2D radiation hydrodynamics code in order to study the impact of massive stars on their surrounding interstellar medium. This first paper deals with the evolution of the circumstellar gas around an isolated 60 Msun star. The interaction of the photoionized HII region with the stellar wind bubble forms a variety of interesting structures like shells, clouds, fingers, and spokes. These results demonstrate that complex structures found in HII regions are not necessarily relics from the time before the gas became ionized but may result from dynamical processes during the course of the HII region evolution. We have also analyzed the transfer and deposit of the stellar wind and radiation energy into the circumstellar medium until the star explodes as a supernova. Although the total mechanical wind energy supplied by the star is negligible compared to the accumulated energy of the Lyman continuum photons, the kinetic energy imparted to the circumstellar gas over the star's lifetime is 4 times higher than for a comparable windless simulation. Furthermore, the thermal energy of warm photoionized gas is lower by some 55%. Our results document the necessity to consider both ionizing radiation and stellar winds for an appropriate description of the interaction of OB stars with their circumstellar environment.

To understand the cosmic history of element synthesis it will be important to obtain extinction-free measures of the heavy element contents of high-redshift objects and to chart two monumental events: the collapse of the first metal-free clouds to form stars, and the initial seeding of the universe with dust. The information needed to achieve these objectives is uniquely available in the far-infrared/submillimeter (FIR/SMM) spectral region. Following the Decadal Report and anticipating the development of the Single Aperture Far-IR (SAFIR) telescope and FIR/SMM interferometry, we estimate the measurement capabilities of a large-aperture, background-limited FIR/SMM observatory and an interferometer on a boom, and discuss how such instruments could be used to measure the element synthesis history of the universe.

The question "How did we get here and what will the future bring?" captures the human imagination and the attention of the National Academy of Science's Astronomy and Astrophysics Survey Commitee (AASC). Fulfillment of this "fundamental goal" requires astronomers to have sensitive, high angular and spectral resolution observations in the far-infrared/submillimeter (far-IR/sub-mm) spectral region. With half the luminosity of the universe and vital information about galaxy, star and planet formation, observations in this spectral region require capabilities similar to those currently available or planned at shorter wavelengths. In this paper we summarize the scientific motivation, some mission concepts and technology requirements for far-IR/sub-mm space interferometers that can be developed in the 2010-2020 timeframe.

We calculate numerically the collapse of slowly rotating, non-magnetic, massive molecular clumps, which conceivably could lead to the formation of massive stars. Because radiative acceleration on dust grains plays a critical role in the clump's dynamical evolution, we utilize a wavelength-dependent radiation transfer and a three component dust model: amorphous carbon particles, silicates and "dirty ice"-coated silicates. We do not spatially resolve the innermost regions of the molecular clump and assume that all material in the innermost grid cell accretes onto a single object. We introduce a semi-analytical scheme for augmenting existing evolution tracks of pre-main sequence protostars by including the effects of accretion. By considering an open outermost boundary, an arbitrary amount of material could, in principal, be accreted onto this central star. However, for the three cases considered (30, 60, and 120 solar masses originally within the computation grid), radiation acceleration limited the final masses to 31.6, 33.6, and 42.9 solar masses, respectively, for wavelength-dependent radiation transfer and to 19.1, 20.1, and 22.9 solar masses for comparison simulations with grey radiation transfer. We demonstrate that massive stars can in principle be formed via accretion through a disk. We conclude with the warning that a careful treatment of radiation transfer is a mandatory requirement for realistic simulations of the formation of massive stars.

We investigate dust dynamics and evolution during the formation of a protostellar accretion disk around intermediate mass stars via 2D numerical simulations. Using three different detailed dust models, compact spherical particles, fractal BPCA grains, and BCCA grains, we find that even during the early collapse and the first 10,000 yr of dynamical disk evolution, the initial dust size distribution is strongly modified. Close to the disk's midplane coagulation produces dust particles of sizes of several 10 micons (for compact spherical grains) up to several mm (for fluffy BCCA grains), whereas in the vicinity of the accretion shock front (located several density scale heights above the disk), large velocity differences inhibit coagulation. Dust particles larger than about 1 micron segregate from the smaller grains behind the accretion shock. Due to the combined effects of coagulation and grain segregation the infrared dust emission is modified. Throughout the accretion disk a MRN dust distribution provides a poor description of the general dust properties. Estimates of the consequences of the "freezing out" of molecules in protostellar disks should consider strongly modified grains. Physical model parameters such as the limiting sticking strength and the grains' resistivity against shattering are crucial factors determining the degree of coagulation reached. In dense regions (e.g. in the mid-plane of the disk) a steady-state is quickly attained. High above the equatorial plane coagulation equilibrium is not reached due to the much lower densities. Here, the dust size distribution is affected primarily by differential advection, rather than coagulation.

We report the discovery of three proplyd-like structures in the giant HII region NGC 3603. The emission nebulae are clearly resolved in narrow-band and broad-band HST/WFPC2 observations in the optical and broad-band VLT/ISAAC observations in the near-infrared. All three nebulae are tadpole shaped, with the bright ionization front at the head facing the central cluster and a fainter ionization front around the tail pointing away from the cluster. Typical sizes are 6,000 A.U. x 20,000 A.U. The nebulae share the overall morphology of the proplyds (``PROto PLanetarY DiskS'') in Orion, but are 20 to 30 times larger in size. Additional faint filaments located between the nebulae and the central ionizing cluster can be interpreted as bow shocks resulting from the interaction of the fast winds from the high-mass stars in the cluster with the evaporation flow from the proplyds. The striking similarity of the tadpole shaped emission nebulae in NGC 3603 to the proplyds in Orion suggests that the physical structure of both types of objects might be the same. We present 2D radiation hydrodynamical simulations of an externally illuminated star-disk-envelope system, which was still in its main accretion phase when first exposed to ionizing radiation from the central cluster. The simulations reproduce the overall morphology of the proplyds in NGC 3603 very well, but also indicate that mass-loss rates of up to 10^-5 Mo/yr are required in order to explain the size of the proplyds. (abbreviated)

We present theoretical continuum emission spectra (SED's), isophotal maps and line profiles for several models of photoevaporating disks at different orientations with respect to the observer. The hydrodynamic evolution of these models has been the topic of the two previous papers of this series. We discuss in detail the numerical scheme used for these diagnostic radiation transfer calculations. Our results are qualitatively compared to observed ultracompact Hii-regions (UCHii's). Our conclusion is that the high fraction of ``unresolved'' UCHii's from the catalogues of Wood & Churchwell (1989) and Kurtz et al. (1994) cannot be explained by disks around massive stars. In particular, the observed infrared spectra of these objects indicate dust temperatures which are about one order of magnitude lower than expected. We suggest that disks around close companions to OB stars may be necessary to resolve this inconsistency. Alternatively, strong stellar winds and radiative acceleration could remove disk material from the immediate vicinity of luminous O stars, whereas for the lower luminosity sources considered here this will not occur. We also find that line profiles tracing the evaporated material originating from the disk are not influenced significantly by the existence of stellar winds over a wide range of wind velocities (400 - 1000 km/s). We compare our results to the bright IRAS source MWC349A. Many of its properties, especially its spatial appearance in high-resolution radio maps, can be well explained by a disk surrounding a UV luminous star with a high velocity stellar wind.