Search

Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures. Results: We infer density distributions n~r^{-p} of the regions with typical power-law slopes p around ~1.5. There is no obvious correlation between the power-law slopes of the density structures on larger clump scales (~1pc) and the number of fragments on smaller core scales (<0.1pc). Comparing the large-scale single-dish density profiles to those derived earlier from interferometric observations at smaller spatial scales, we find that the smaller-scale power-law slopes are steeper, typically around ~2.0. The flattening toward larger scales is consistent with the star-forming regions being embedded in larger cloud structures that do not decrease in density away from a particular core. Regarding the magnetic field, for several regions it appears aligned with filamentary structures leading toward the densest central cores. Furthermore, we find different polarization structures with some regions exhibiting central polarization holes whereas other regions show polarized emission also toward the central peak positions. Nevertheless, the polarized intensities are inversely related to the Stokes I intensities. We estimate magnetic field strengths between ~0.2 and ~4.5mG, and we find no clear correlation between magnetic field strength and the fragmentation level of the regions. Comparison of the turbulent to magnetic energies shows that they are of roughly equal importance in this sample. The mass-to-flux ratios range between ~2 and ~7, consistent with collapsing star-forming regions., Comment: Accepted for Astronomy & Astrophysics, 14 pages, 14 figures plus appendices, also download option at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in the 1.37 mm wavelength regime in both continuum and line emission, reaching $\sim$0.4" resolution (800 au at 2 kpc). Using the dense gas tracer CH$_3$CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. Specific angular momentum ($j$) radial profiles are on average $\sim10^{-3}$ km /s pc and follow $j \propto r^{1.7}$, consistent with a poorly resolved rotating and infalling envelope/disk model. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of $\sim 10-25$ $M_\odot$. Modelling the level population of CH$_3$CN lines, we present temperature maps and find median gas temperatures in the range $70-210$ K. We create Toomre $Q$ maps to study the stability of the disks and find almost all (11 of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations. In particular, disks with masses greater than $\sim10-20\%$ of the mass of their host (proto)stars are Toomre unstable, and more luminous protostellar objects tend to have disks that are more massive and hence more prone to fragmentation. Our finings show that most disks around high-mass protostars are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars., Comment: 27 pages, 12 figures, 6 appendices - accepted for publication in Astronomy and Astrophysics

We aim to characterize the physical and chemical properties of fragmented cores during the earliest evolutionary stages in the very young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953. NOEMA 1.3 mm data are used in combination with archival mid- and far-infrared observations to construct and fit the SEDs of individual fragmented cores. The radial density profiles are inferred from the 1.3 mm continuum visibility profiles and the radial temperature profiles are estimated from H2CO rotation temperature maps. Molecular column densities are derived with the line fitting tool XCLASS. The physical and chemical properties are combined by applying the physical-chemical model MUSCLE in order to constrain the chemical timescales of a few line-rich cores. The morphology and spatial correlations of the molecular emission are analyzed using the HOG method. The mid-infrared data show that both regions contain a cluster of young stellar objects. Bipolar molecular outflows are observed in the CO 2-1 transition toward the strong mm cores indicating protostellar activity. We find strong molecular emission of SO, SiO, H2CO, and CH3OH in locations which are not associated with the mm cores. These shocked knots can be either associated with the bipolar outflows or, in the case of ISOSS J23053+5953, with a colliding flow that creates a large shocked region between the mm cores. The mean chemical timescale of the cores is lower (20 000 yr) compared to that of the sources of the more evolved CORE sample (60 000 yr). With the HOG method, we find that the spatial emission of species tracing the extended emission and of shock-tracing molecules are well correlated within transitions of these groups., Comment: 34 pages, 17 figures, accepted for publication in A&A

Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk around the high-mass star forming region AFGL 2591-VLA 3 which was hypothesized to be fragmenting based on the observations that show multiple outflow directions. We use a new set of high-resolution (0.19 arcsec) IRAM/NOEMA observations at 843 micron towards VLA 3 which allow us to resolve its disk, characterize the fragmentation, and study its kinematics. In addition to the 843 micron continuum emission, our spectral setup targets warm dense gas and outflow tracers such as HCN, HC$_3$N and SO$_2$, as well as vibrationally excited HCN lines. The high resolution continuum and line emission maps reveal multiple fragments with subsolar masses within the inner 1000 AU of VLA 3. Furthermore, the velocity field of the inner disk observed at 843 micron shows a similar behavior to that of the larger scale velocity field studied in the CORE project at 1.37 mm. We present the first observational evidence for disk fragmentation towards AFGL 2591-VLA 3, a source that was thought to be a single high-mass core. While the fragments themselves are low-mass, the rotation of the disk is dominated by the protostar with a mass of 10.3$\pm 1.8~M_{\odot}$. These data also show that NOEMA Band 4 can obtain the highest currently achievable spatial resolution at (sub-)mm wavelengths in observations of strong northern sources., Comment: 16 pages, 11 figures, Accepted for publication in A&A

Context: The formation of high-mass star-forming regions from their parental gas cloud and the subsequent fragmentation processes lie at the heart of star formation research. Aims: We aim to study the dynamical and fragmentation properties at very early evolutionary stages of high-mass star formation. Methods: Employing the NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30m telescope, we observed two young high-mass star-forming regions, ISOSS22478 and ISOSS23053, in the 1.3mm continuum and spectral line emission at a high angular resolution (~0.8''). Results: We resolved 29 cores that are mostly located along filament-like structures. Depending on the temperature assumption, these cores follow a mass-size relation of approximately M~r^2.0, corresponding to constant mean column densities. However, with different temperature assumptions, a steeper mass-size relation up to M~r^3.0, which would be more likely to correspond to constant mean volume densities, cannot be ruled out. The correlation of the core masses with their nearest neighbor separations is consistent with thermal Jeans fragmentation. We found hardly any core separations at the spatial resolution limit, indicating that the data resolve the large-scale fragmentation well. Although the kinematics of the two regions appear very different at first sight - multiple velocity components along filaments in ISOSS22478 versus a steep velocity gradient of more than 50km/s/pc in ISOSS23053 - the findings can be explained within the framework of a dynamical cloud collapse scenario. Conclusions: While our data are consistent with a dynamical cloud collapse scenario and subsequent thermal Jeans fragmentation, the importance of additional environmental properties, such as the magnetization of the gas or external shocks triggering converging gas flows, is nonetheless not as well constrained and would require future investigation., Comment: 16 pages, 16 figures, accepted for Astronomy and Astrophysics, a higher-resolution version can also be found at https://www.mpia.de/homes/beuther/papers.html

We use sub-arcsecond resolution ($\sim$0.4$''$) observations with NOEMA at 1.37 mm to study the dust emission and molecular gas of 18 high-mass star-forming regions. We combine the derived physical and chemical properties of individual cores in these regions to estimate their ages. The temperature structure of these regions are determined by fitting H2CO and CH3CN line emission. The density profiles are inferred from the 1.37 mm continuum visibilities. The column densities of 11 different species are determined by fitting the emission lines with XCLASS. Within the 18 observed regions, we identify 22 individual cores with associated 1.37 mm continuum emission and with a radially decreasing temperature profile. We find an average temperature power-law index of q = 0.4$\pm$0.1 and an average density power-law index of p = 2.0$\pm$0.2 on scales on the order of several 1 000 au. Comparing these results with values of p derived in the literature suggest that the density profiles remain unchanged from clump to core scales. The column densities relative to N(C18O) between pairs of dense gas tracers show tight correlations. We apply the physical-chemical model MUSCLE to the derived column densities of each core and find a mean chemical age of $\sim$60 000 yrs and an age spread of 20 000-100 000 yrs. With this paper we release all data products of the CORE project available at https://www.mpia.de/core. The CORE sample reveals well constrained density and temperature power-law distributions. Furthermore, we characterize a large variety in molecular richness that can be explained by an age spread confirmed by our physical-chemical modeling. The hot molecular cores show the most emission lines, but we also find evolved cores at an evolutionary stage, in which most molecules are destroyed and thus the spectra appear line-poor again., Comment: 29 pages, 14 figures, accepted for publication in A&A

In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we derive a mass infall rate of about 10^{-4}~M_{sun}~yr^{-1}. The most massive cores (labeled 1, 2, and 3) are found at the center of the cluster, and these are the only ones that present a signature of protostellar activity in terms of emission from high-excitation molecular lines or a molecular outflow. We reveal an extended (size about 0.1~pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by a (previously discovered) compact (size <= 1000~au) radio jet, ejected by a YSO embedded in core 1 (named YSO-1), since the molecular outflow and the radio jet are almost parallel and have a comparable momentum rate. By means of high-excitation lines, we find a large (14~km/s over 500~au) Vlsr gradient at the position of YSO-1, oriented approximately perpendicular to the radio jet. Assuming this is an edge-on, rotating disk and fitting a Keplerian rotation pattern, we determine the YSO-1 mass to be 5.6+/-2.0~M_{sun}. The water masers (previously observed with VLBI) emerge within 100-300~au from YSO-1 and are unique tracers of the jet kinematics. Their three-dimensional (3D) velocity pattern reveals that the gas flows along, and rotates about, the jet axis. We show that the 3D maser velocities are fully consistent with the magneto-centrifugal disk-wind models predicting a cylindrical rotating jet. Under this hypothesis, we determine the jet radius to be about 16~au and the corresponding launching radius and terminal velocity to be about 2.2~au and 200~km/s, respectively., Comment: 25 pages, 16 figures

We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth., Comment: 22 pages, 13 figures, accepted for publication in A&A

The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions., Comment: 23 pages, 16 figures, 5 tables, accepted for publication in A&A

This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars.

The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales., Comment: 23 pages, 26 figures, accepted for publication in Astronomy and Astrophysics

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation., Comment: 29 pages, 17 figures, Astronomy & Astrophysics in press, for a higher-resolution version of the paper see http://www.mpia.de/homes/beuther/papers.html

We analyse simulations of turbulent, magnetised molecular cloud cores focussing on the formation of Class 0 stage protostellar discs and the physical conditions in their surroundings. We show that for a wide range of initial conditions Keplerian discs are formed in the Class 0 stage already. In particular, we show that even subsonic turbulent motions reduce the magnetic braking efficiency sufficiently in order to allow rotationally supported discs to form. We therefore suggest that already during the Class 0 stage the fraction of Keplerian discs is significantly higher than 50%, consistent with recent observational trends but significantly higher than predictions based on simulations with misaligned magnetic fields, demonstrating the importance of turbulent motions for the formation of Keplerian discs. We show that the accretion of mass and angular momentum in the surroundings of protostellar discs occurs in a highly anisotropic manner, by means of a few narrow accretion channels. The magnetic field structure in the vicinity of the discs is highly disordered, revealing field reversals up to distances of 1000 AU. These findings demonstrate that as soon as even mild turbulent motions are included, the classical disc formation scenario of a coherently rotating environment and a well-ordered magnetic field breaks down. Hence, it is highly questionable to assess the magnetic braking efficiency based on non-turbulent collapse simulation. We strongly suggest that, in addition to the global magnetic field properties, the small-scale accretion flow and detailed magnetic field structure have to be considered in order to assess the likelihood of Keplerian discs to be present., Comment: 14 pages, 6 figures, accepted for publication in MNRAS, updated to final version

We present collapse simulations of strongly magnetised, turbulent molecular cloud cores with masses ranging from 2.6 to 1000 M_sun in order to study the influence of the initial conditions on the turbulence-induced disc formation mechanism proposed recently by Seifried et al. 2012. We find that Keplerian discs are formed in all cases independently of the core mass, the strength of turbulence, or the presence of global rotation. The discs appear within a few kyr after the formation of the protostar, are 50 - 150 AU in size, and have masses between 0.05 and a few 0.1 M_sun. During the formation of the discs the mass-to-flux ratio stays well below the critical value of 10 for Keplerian disc formation. Hence, flux-loss alone cannot explain the formation of Keplerian discs. The formation of rotationally supported discs at such early phases is rather due to the disordered magnetic field structure and due to turbulent motions in the surroundings of the discs, two effects lowering the classical magnetic braking efficiency. Binary systems occurring in the discs are mainly formed via the disc capturing mechanism rather than via disc fragmentation, which is largely suppressed by the presence of magnetic fields., Comment: 13 pages, 11 figures, accepted for publication in MNRAS, updated to final version

We present collapse simulations of strongly magnetised, 100 M_sun, turbulent cloud cores. Around the protostars formed during the collapse Keplerian discs with typical sizes of up to 100 AU build up in contrast to previous simulations neglecting turbulence. Analysing the condensations in which the discs form, we show that the magnetic flux loss is not sufficient to explain the build-up of Keplerian discs. The average magnetic field is strongly inclined to the disc which might reduce the magnetic braking efficiency. However, the main reason for the reduced magnetic braking efficiency is the highly disordered magnetic field in the surroundings of the discs. Furthermore, due to the lack of a coherently rotating structure in the turbulent environment of the disc no toroidal magnetic field necessary for angular momentum extraction can build up. Simultaneously the angular momentum inflow remains high due to local shear flows created by the turbulent motions. We suggest that the "magnetic braking catastrophe" is an artefact of the idealised non-turbulent initial conditions and that turbulence provides a natural mechanism to circumvent this problem., Comment: 4 pages, 2 figures. To appear in the proceedings of 'The Labyrinth of Star Formation' (18-22 June 2012, Chania, Greece), published by Springer

Long-lived, stable jets are observed in a wide variety of systems, from protostars, through Galactic compact objects to active galactic nuclei (AGN). Magnetic fields play a central role in launching, accelerating, and collimating the jets through various media. The termination of jets in molecular clouds or the interstellar medium deposits enormous amounts of mechanical energy and momentum, and their interactions with the external medium, as well, in many cases, as the radiation processes by which they are observed, are intimately connected with the magnetic fields they carry. This review focuses on the properties and structures of magnetic fields in long-lived jets, from their launch from rotating magnetized young stars, black holes, and their accretion discs, to termination and beyond. We compare the results of theory, numerical simulations, and observations of these diverse systems and address similarities and differences between relativistic and non-relativistic jets in protostellar versus AGN systems. On the observational side, we focus primarily on jets driven by AGN because of the strong observational constraints on their magnetic field properties, and we discuss the links between the physics of these jets on all scales., Comment: Refereed review accepted by Space Science Reviews; contribution to ISSI workshop on 'Large-Scale Magnetic Fields in the Universe'. 48 pages. v2 fixes a couple of typos and updates references

We present collapse simulations of 100 M_{\sun}, turbulent cloud cores threaded by a strong magnetic field. During the initial collapse phase filaments are generated which fragment quickly and form several protostars. Around these protostars Keplerian discs with typical sizes of up to 100 AU build up in contrast to previous simulations neglecting turbulence. We examine three mechanisms potentially responsible for lowering the magnetic braking efficiency and therefore allowing for the formation of Keplerian discs. Analysing the condensations in which the discs form, we show that the build-up of Keplerian discs is neither caused by magnetic flux loss due to turbulent reconnection nor by the misalignment of the magnetic field and the angular momentum. It is rather a consequence of the turbulent surroundings of the disc which exhibit no coherent rotation structure while strong local shear flows carry large amounts of angular momentum. We suggest that the "magnetic braking catastrophe", i.e. the formation of sub-Keplerian discs only, is an artefact of the idealised non-turbulent initial conditions and that turbulence provides a natural mechanism to circumvent this problem., Comment: 6 pages, 5 figures, accepted by MNRAS Letters, updated to final version

Numerical simulations of outflows formed during the collapse of 100 M_sun cloud cores are presented. We derive a generalised criterion from MHD wind theory to analyse the launching mechanism of these outflows. The criterion is successfully applied to the whole outflow structure and cases with sub-Keplerian disc rotation. It allows us to decide whether an outflow is driven centrifugally or by the toroidal magnetic pressure. We show that quantities such as the magnetic field line inclination or the ratio of the toroidal to poloidal magnetic field alone are insufficient to determine the driving mechanism of outflows. By performing 12 runs we are able to study the influence of the initial conditions on the properties of outflows around massive protostars in detail. Our simulations reveal a strong effect of the magnetic field strength on the outflow morphology. In runs with weak fields or high rotational energies, well collimated, fast jets are observed whereas for strong fields poorly collimated, low-velocity outflows are found. We show that the occurrence of a fast jet is coupled to the existence of a Keplerian protostellar disc. Despite the very different morphologies all outflows are launched from the discs by centrifugal acceleration with the toroidal magnetic field increasingly contributing to the gas acceleration further away from the discs. The poor collimation of the outflows in runs with strong fields is a consequence of the weak hoop stresses. This in turn is caused by the slow build-up of a toroidal field due to sub-Keplerian disc rotation. The mass ejection/accretion ratios scatter around a mean of 0.3 in accordance with observational and analytical results. We suggest an evolutionary scenario for the earliest stage of massive star formation in which initially poorly collimated outflows develop which progressively get better collimated due to the generation of fast jets., Comment: 22 pages, 15 figures, accepted by MNRAS, updated to final version

We present simulations of collapsing 100 M_\sun mass cores in the context of massive star formation. The effect of variable initial rotational and magnetic energies on the formation of massive stars is studied in detail. We focus on accretion rates and on the question under which conditions massive Keplerian disks can form in the very early evolutionary stage of massive protostars. For this purpose, we perform 12 simulations with different initial conditions extending over a wide range in parameter space. The equations of magnetohydrodynamics (MHD) are solved under the assumption of ideal MHD. We find that the formation of Keplerian disks in the very early stages is suppressed for a mass-to-flux ratio normalised to the critical value \mu below 10, in agreement with a series of low-mass star formation simulations. This is caused by very efficient magnetic braking resulting in a nearly instantaneous removal of angular momentum from the disk. For weak magnetic fields, corresponding to \mu > 10, large-scale, centrifugally supported disks build up with radii exceeding 100 AU. A stability analysis reveals that the disks are supported against gravitationally induced perturbations by the magnetic field and tend to form single stars rather than multiple objects. We find protostellar accretion rates of the order of a few 10^-4 M_\sun yr^-1 which, considering the large range covered by the initial conditions, vary only by a factor of ~ 3 between the different simulations. We attribute this fact to two competing effects of magnetic fields. On the one hand, magnetic braking enhances accretion by removing angular momentum from the disk thus lowering the centrifugal support against gravity. On the other hand, the combined effect of magnetic pressure and magnetic tension counteracts gravity by exerting an outward directed force on the gas in the disk thus reducing the accretion onto the protostars., Comment: 22 pages, 17 figures, accepted for publication in MNRAS, updated to final version

We examine the case for Quark-Novae (QNe) as possible sources for the reionization and early metal enrichment of the universe. Quark-Novae are predicted to arise from the explosive collapse (and conversion) of sufficiently massive neutron stars into quark stars. A Quark-Nova (QN) can occur over a range of time scales following the supernova event. For QNe that arise days to weeks after the supernovae, we show that dual-shock that arises as the QN ejecta encounter the supernova ejecta can produce enough photons to reionize hydrogen in most of the Inter-Galactic medium (IGM) by z ~ 6. Such events can explain the large optical depth tau_e ~ 0.1 as measured by WMAP, if the clumping factor, C, of the material being ionized is smaller than 10. We suggest a way in which a normal initial mass function (IMF) for the oldest stars can be reconciled with a large optical depth as well as the mean metallicity of the early IGM post reionization. We find that QN also make a contribution to r-process element abundances for atomic numbers A > 130. We predict that the main cosmological signatures of Quark-Novae are the gamma-ray bursts that announce their birth. These will be clustered at redshifts in the range z ~ 7-8 in our model., Comment: Accepted for publication in the Astrophysical Journal. 9 journal pages, 2 figures. This version includes more discussion, 3 new appendices and extended literature

Given the importance of simulating hydromagnetic processes that impact star formation, we have earlier developed a 3D adaptive mesh approach that allows us to include hydromagnetic processes during the formation and evolution of cores, discs, and stars in observed regions of star formation. In this paper, we take the next step in this program - namely - to develop a modified version of the 3D adaptive mesh refinement (AMR) code FLASH in which the ambipolar diffusion of the magnetic field in poorly ionized molecular gas is implemented. We approach the problem using a single-fluid approximation to simplify numerical calculations. In this paper, we present a series of test cases including oblique isothermal and non--isothermal C-shocks. We also present a study of the quasi-static collapse of an initial uniform, self-gravitating, magnetized sphere that is initially supported by its magnetic field against collapse (i.e. magnetically subcritical). Applications to the collapse of a pre-stellar Bonnor-Ebert sphere are presented in a companion paper., Comment: to appear in MNRAS; 18 pages, 10 figures

We develop models for the self-similar collapse of magnetized isothermal cylinders. We find solutions for the case of a fluid with a constant toroidal flux-to-mass ratio (Gamma_phi=constant) and the case of a fluid with a constant gas to magnetic pressure ratio (beta=constant). In both cases, we find that a low magnetization results in density profiles that behave as rho ~ r^{-4} at large radii, and at high magnetization we find density profiles that behave as rho ~ r^{-2}. This density behaviour is the same as for hydrostatic filamentary structures, suggesting that density measurements alone cannot distinguish between hydrostatic and collapsing filaments--velocity measurements are required. Our solutions show that the self-similar radial velocity behaves as v_r ~ r during the collapse phase, and that unlike collapsing self-similar spheres, there is no subsequent accretion (i.e. expansion-wave) phase. We also examine the fragmentation properties of these cylinders, and find that in both cases, the presence of a toroidal field acts to strengthen the cylinder against fragmentation. Finally, the collapse time scales in our models are shorter than the fragmentation time scales. Thus, we anticipate that highly collapsed filaments can form before they are broken into pieces by gravitational fragmentation., Comment: 20 pages, 4 figures, accepted to ApJ

We have performed fully 3D simulations of the collapse of molecular cloud cores which obey the logatropic equation of state. By following the collapse of these cores from states of near hydrodynamic equilibrium, we are able to produce accretion histories which closely resemble those of observed cores. The accretion proceeds in four distinct stages: an initial period of very slow accretion; a period of vigorous accretion following the development of a central density singularity, with a mass accretion rate proportional to t^3, as predicted by self-similar models; a period of relatively stable, vigorous accretion; and finally a gradual decrease in the accretion rate once about 50% of the available mass of the molecular cloud core has been accreted. These results may explain the accretion histories of cores as they pass through the pre-protostellar, Class 0, and Class I stages., Comment: To be published in: Chemistry as a diagnostic of star formation, eds. C.L. Curry & M. Fich (NRC Press). 3 pages, 3 figures

We present the extension of previous two-dimensional simulations of the time-dependent evolution of non-relativistic outflows from the surface of Keplerian accretion disks, to three dimensions. The accretion disk itself is taken to provide a set of fixed boundary conditions for the problem. The 3-D results are consistent with the theory of steady, axisymmetric, centrifugally driven disk winds up to the Alfv\'en surface of the outflow. Beyond the Alfv\'en surface however, the jet in 3-D becomes unstable to non-axisymmetric, Kelvin-Helmholtz instabilities. We show that jets maintain their long-term stability through a self-limiting process wherein the average Alfv\'enic Mach number within the jet is maintained to order unity. This is accomplished in at least two ways. First, poloidal magnetic field is concentrated along the central axis of the jet forming a ``backbone'' in which the Alfv\'en speed is sufficiently high to reduce the average jet Alfv\'enic Mach number to unity. Second, the onset of higher order Kelvin-Helmholtz ``flute'' modes (m \ge 2) reduce the efficiency with which the jet material is accelerated, and transfer kinetic energy of the outflow into the stretched, poloidal field lines of the distorted jet. This too has the effect of increasing the Alfv\'en speed, and thus reducing the Alfv\'enic Mach number. The jet is able to survive the onset of the more destructive m=1 mode in this way. Our simulations also show that jets can acquire corkscrew, or wobbling types of geometries in this relatively stable end-state, depending on the nature of the perturbations upon them. Finally, we suggest that jets go into alternating periods of low and high activity as the disappearance of unstable modes in the sub-Alfv\'enic regime enables another cycle of acceleration to super-Alfv\'enic speeds., Comment: 57 pages, 22 figures, submitted to ApJ

We present the results of fully 3-D hydrodynamic simulations of the gravitational collapse of isolated, turbulent molecular cloud cores. Starting from initial states of hydrostatic equilibrium, we follow the collapse of both singular and nonsingular logatropic cores until the central protostar has accreted > 90% of the total available mass. We find that, in the collapse of a singular core with access to a finite mass reservoir, the mass of the central protostar increases as M_acc proportional to t^4 until it has accreted about 35% of the total available mass. For nonsingular cores of fiducial masses 1, 2.5, and 5 M_solar, we find that protostellar accretion proceeds slowly prior to the formation of a singular density profile. Immediately thereafter, the accretion rate in each case increases to about 10^{-6} M_solar/yr, for cores with central temperature T_c= 10 K and truncation pressure P_s = 1.3E5 k_B K/cm^3. It remains at that level until half the available mass has been accreted. After this point, the accretion rate falls steadily as the remaining material is accreted onto the growing protostellar core. We suggest that this general behaviour of the protostellar accretion rate may be indicative of evolution from the Class 0 to the Class I protostellar phase., Comment: 18 pages, 11 figures. Accepted for publication in the Astrophysical Journal

Young stars form on a wide range of scales, producing aggregates and clusters with various degrees of gravitational self-binding. The loose aggregates have a hierarchical structure in both space and time that resembles interstellar turbulence, suggesting that these stars form in only a few turbulent crossing times with positions that map out the previous gas distribution. Dense clusters, on the other hand, are often well mixed, as if self-gravitational motion has erased the initial fine structure. Nevertheless, some of the youngest dense clusters also show sub-clumping, so it may be that all stellar clustering is related to turbulence. Some of the densest clusters may also be triggered. The evidence for mass segregation of the stars inside clusters is reviewed, along with various explanations for this effect. Other aspects of the theory of cluster formation are reviewed as well, including many specific proposals for cluster formation mechanisms. The conditions for the formation of bound clusters are discussed. Critical star formation efficiencies can be as low as 10% if the gas removal process is slow and the stars are born at sub-virial speeds. Environmental conditions, particularly pressure, may affect the fraction and masses of clusters that end up bound. Globular clusters may form like normal open clusters but in conditions that prevailed during the formation of the halo and bulge, or in interacting and starburst galaxies today. Various theories for the formation of globular clusters are summarized., Comment: 36 pages, 5 figures, To be published in Protostars and Planets IV, eds. V. G. Mannings, A. P. Boss, and S. S. Russell, from the conference at Santa Barbara, CA, July 6-11, 1998

New wide-field images of the Rho Ophiuchus molecular cloud at 850 and 450 microns obtained with SCUBA on the James Clerk Maxwell Telescope reveal a wide variety of large-scale features that were previously unknown. Two linear features each 4' (0.2 pc) in length extend to the north of the bright emission region containing SM1 and VLA 1623. These features may correspond to the walls of a previously unidentified outflow cavity, or the boundary of a photon dominated region powered by a nearby B star. A previously unidentified source is located in the north-east region of the image. The properties of this source (diameter ~5000 AU, mass ~0.3-1 Msun) suggest that it is a pre-protostellar core. Two arcs of emission are seen in the direction of the north-west extension of the VLA 1623 outflow. The outer arc appears relatively smooth at 850 microns and is estimated to have a mass of ~0.3 Msun, while the inner arc breaks up into a number of individual clumps, some of which are previously identified protostars., Comment: 9 pages, 2 color figures; accepted for publication in ApJ Letters

Observational advances over the last decade reveal that star formation is associated with the simultaneous presence of gravitationally collapsing gas, bipolar outflow, and an accretion disk. Two theoretical views of star formation suppose that either stellar mass is determined from the outset by gravitational instability, or by the outflow which sweeps away the collapsing envelope of initially singular density distributions. Neither picture appears to explain all of the facts. This contribution examines some of the key issues facing star formation theory., Comment: 12 pages with 2 ps figures. Invited review to appear in Computational Astrophysics, Proceedings of the 12th Kingston Meeting, ed. D. A. Clarke and M. J. West (San Francisco: ASP)

Continuum emission at 220 and 355 GHz from the prototype Class 0 source VLA 1623 has been detected using the JCMT-CSO interferometer. Gaussian fits to the data place an upper limit of 70 AU on the half-width half-maximum radius of the emission, which implies an upper limit of ~175 AU for the cutoff radius of the circumstellar disk in the system. In the context of existing collapse models, this disk could be magnetically supported on the largest scales and have an age of ~6x10^4 yr, consistent with previous suggestions that Class 0 sources are quite young. The innermost region of the disk within ~6 AU is likely to be in centrifugal support, which is likely large enough to provide a drive for the outflow according to current theoretical models. Alternatively, if 175 AU corresponds to the centrifugal radius of the disk, the age of the system is ~2x10^5 yr, closer to age estimates for Class I sources., Comment: 14 pages, aastex, 2 postscript figures; accepted for publication in ApJ Letters; also available at http://www.physics.mcmaster.ca/Wilson_Preprints