Source author record

L. Testi

L. Testi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

148works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

148 published item(s)

preprint2025arXiv

ALMAGAL V. Relations between the core populations and the parent clump physical properties

Context. The fragmentation of massive molecular clumps into smaller, potentially star-forming cores plays a key role in the processes of high-mass star formation. The ALMAGAL project offers high-resolution data to investigate these processes across various evolutionary stages in the Galactic plane. Aims. This study aims at correlating the fragmentation properties of massive clumps, obtained from ALMA observations, with their global physical parameters (e.g., mass, surface density, and temperature) and evolutionary indicators (such as luminosity-to-mass ratio and bolometric temperature) obtained from Herschel observations. It seeks to assess whether the cores evolve in number and mass in tandem with their host clumps, and to determine the possible factors influencing the formation of massive cores (M > 24M_\odot). Methods. We analyzed the masses of 6348 fragments, estimated from 1.4 mm continuum data for 1007 ALMAGAL clumps. Leveraging this unprecedentedly large data set, we evaluated statistical relationships between clump parameters, estimated over about 0.1 pc scales, and fragment properties, corresponding to scales of a few 1000 au, while accounting for potential biases related to distance and observational resolution. Our results were further compared with predictions from numerical simulations. Results. The fragmentation level correlates preferentially with clump surface density, supporting a scenario of density-driven fragmentation, whereas it does not show any clear dependence on total clump mass. Both the mass of the most massive core and the core formation efficiency show a broad range and increase on average by an order of magnitude in the intervals spanned by evolutionary indicators such as clump dust temperature and the luminosity-to-mass ratio. This suggests that core growth continues throughout the clump evolution, favoring clump-fed over core-fed theoretical scenarios.

preprint2025arXiv

FAUST. XXVIII. High-Resolution ALMA Observations of Class 0/I Disks: Structure, Optical Depths, and Temperatures

We present high-resolution (~7.5 au) ALMA observations at 1.3 and 3 mm of 16 disks around Class 0/I protostars across multiple star-forming regions and a variety of multiplicities, showing a range of disk sizes (~2-100 au) and including circumbinary disks (CBDs) in binaries with separations <100 au. The disk properties show similarities to Class II disks, including (a) low spectral index (SI) values (alpha=2.1) that increase with disk radius, (b) 3 mm disk sizes only marginally smaller than at 1.3 mm (<10%), and (c) radial intensity profiles well described by modified self-similar profiles. We also find key differences: (i) SI values increasing with radius, but exceeding 2 only at the disk edge (ii) higher brightness temperatures Tb, in some cases higher than the predicted temperatures due to irradiation, and (iii) ~10x higher luminosity at a given size compared to the Class II disks. These results confirm significant optical depth in the observed Class 0/I disks, at both 1.3 and 3 mm, helping to explain their higher luminosities, but higher temperatures are also required for the most compact (< 40 au) disks, suggesting additional viscous heating. Considering optical depth, most disk dust masses are estimated in the range 30-900 Mearth (0.01-0.3 Msun in gas), resulting in some disks reaching marginal gravitational instability. The median location of the water iceline is ~3 au, but it can extend beyond 10-20 au for the hottest disks. CBDs exhibit lower optical depths at both wavelengths and hence higher SI values (alpha=3.0), dust masses of 100 Mearth, and beta~1.5 (2 Class 0 CBDs) and beta~1 (1 Class I CBD), suggesting substantial grain growth only in the more evolved CBD. The inferred high optical depths provide a compelling explanation for the apparent scarcity of dust substructures in the younger disks at ~ 1 mm, despite mounting evidence for early planet formation (ABRIDGED).

preprint2022arXiv

A giant planet shaping the disk around the very low-mass star CIDA 1

(Abridged) Exoplanetary research has provided us with exciting discoveries of planets around very low-mass (VLM) stars (e.g., TRAPPIST-1 and Proxima Centauri). However, current theoretical models strive to explain planet formation in these conditions and do not predict the development of giant planets. Recent high-resolution observations from ALMA of the disk around CIDA 1, a VLM star in Taurus, show substructures hinting at the presence of a massive planet. We aim to reproduce the dust ring of CIDA 1, observed in the dust continuum emission in ALMA Band 7 (0.9 mm) and Band 4 (2.1 mm), along with its $^{12}$CO (J=3-2) and $^{13}$CO (J=3-2) channel maps, assuming the structures are shaped by the interaction of the disk with a massive planet. We seek to retrieve the mass and position of the putative planet. We model the protoplanetary disk with a set of hydrodynamical simulations, varying the mass and locations of the embedded planet. We compute the dust and gas emission using radiative transfer simulations, and, finally, we obtain the synthetic observations treating the images as the actual ALMA observations. Our models indicate that a planet with a minimum mass of $\sim1.4\,\text{M}_\text{Jup}$ orbiting at a distance of $\sim 9-10$ au can explain the morphology and location of the observed dust ring at Band 7 and Band 4. We can reproduce the low spectral index ($\sim 2$) observed where the dust ring is detected. Our synthetic images reproduce the morphology of the $^{12}$CO and $^{13}$CO observed channel maps where the cloud absorption allowed a detection. Applying an empirical relation between planet mass and gap width in the dust, we predict a maximum planet mass of $\sim4 - 8\,\text{M}_\text{Jup}$. Our results suggest the presence of a massive planet orbiting CIDA 1, thus challenging our understanding of planet formation around VLM stars.

preprint2022arXiv

CHEMOUT: CHEMical complexity in star-forming regions of the OUTer Galaxy. I. Organic molecules and tracers of star-formation activity

The outer Galaxy is an environment with metallicity lower than the Solar one. Because of this, the formation and survival of molecules in star-forming regions located in the inner and outer Galaxy is expected to be different. To gain understanding on how chemistry changes throughout the Milky Way, it is crucial to observe outer Galaxy star-forming regions to constrain models adapted for lower metallicity environments. In this paper we present a new observational project: chemical complexity in star-forming regions of the outer Galaxy (CHEMOUT). The goal is to unveil the chemical composition in 35 dense molecular clouds associated with star-forming regions of the outer Galaxy through observations obtained with the IRAM 30m telescope. In this first paper, we present the sample, and report the detection at 3~mm of simple organic species HCO+, H13CO+, HCN, c-C3H2, HCO, C4H, and HCS+, of the complex hydrocarbon CH3CCH, and of SiO, CCS and SO. From c-C3H2, we estimate new kinematic heliocentric and Galactocentric distances based on an updated rotation curve of the Galaxy. The detection of the molecular tracers does not seem to have a clear dependence on the Galactocentric distance. We also analyse the HCO+ line profiles. We find high velocity wings in ~71% of the targets, and their occurrence does not depend on the Galactocentric distance. Our results, confirmed by a statistical analysis, show that the presence of organic molecules and tracers of protostellar activity is ubiquitous in the low-metallicity environment of the outer Galaxy. Based on this, and on the additional evidence that small, terrestrial planets are omnipresent in the Galaxy, we support previous claims that the definition of Galactic Habitable Zone should be rediscussed in view of the ubiquitous capacity of the interstellar medium to form organic molecules.

preprint2022arXiv

Distribution of solids in the rings of the HD 163296 disk: a multiwavelength study

In this paper we analyze new observations from ALMA and VLA, at a high angular resolution corresponding to 5 - 8 au, of the protoplanetary disk around HD 163296 to determine the dust spatial distribution and grain properties. We fit the spectral energy distribution as a function of the radius at five wavelengths from 0.9 to 9\,mm, using a simple power law and a physical model based on an analytic description of radiative transfer that includes isothermal scattering. We considered eight dust populations and compared the models' performance using Bayesian evidence. Our analysis shows that the moderately high optical depth ($τ$>1) at $λ\leq$ 1.3 mm in the dust rings artificially lower the millimeter spectral index, which should therefore not be considered as a reliable direct proxy of the dust properties and especially the grain size. We find that the outer disk is composed of small grains on the order of 200 $μ$m with no significant difference between rings at 66 and 100 au and the adjacent gaps, while in the innermost 30 au, larger grains ($\geq$mm) could be present. We show that the assumptions on the dust composition have a strong impact on the derived surface densities and grain size. In particular, increasing the porosity of the grains to 80\% results in a total dust mass about five times higher with respect to grains with 25\% porosity. Finally, we find that the derived opacities as a function of frequency deviate from a simple power law and that grains with a lower porosity seem to better reproduce the observations of HD163296. While we do not find evidence of differential trapping in the rings of HD163296, our overall results are consistent with the postulated presence of giant planets affecting the dust temperature structure and surface density, and possibly originating a second-generation dust population of small grains.

preprint2022arXiv

Hot methanol in the [BHB2007] 11 protobinary system: hot corino versus shock origin? : FAUST V

Methanol is a ubiquitous species commonly found in the molecular interstellar medium. It is also a crucial seed species for the building-up of the chemical complexity in star forming regions. Thus, understanding how its abundance evolves during the star formation process and whether it enriches the emerging planetary system is of paramount importance. We used new data from the ALMA Large Program FAUST (Fifty AU STudy of the chemistry in the disk/envelope system of Solar-like protostars) to study the methanol line emission towards the [BHB2007] 11 protobinary system (sources A and B), where a complex structure of filaments connecting the two sources with a larger circumbinary disk has been previously detected. Twelve methanol lines have been detected with upper energies in the range [45-537] K along with one 13CH3OH transition. The methanol emission is compact and encompasses both protostars, separated by only 28 au and presents three velocity components, not spatially resolved by our observations, associated with three different spatial regions, with two of them close to 11B and the third one associated with 11A. A non-LTE radiative transfer analysis of the methanol lines concludes that the gas is hot and dense and highly enriched in methanol with an abundance as high as 1e-5. Using previous continuum data, we show that dust opacity can potentially completely absorb the methanol line emission from the two binary objects. Although we cannot firmly exclude other possibilities, we suggest that the detected hot methanol is resulting from the shocked gas from the incoming filaments streaming towards [BHB2007] 11 A and B, respectively. Higher spatial resolution observations are necessary to confirm this hypothesis.

preprint2022arXiv

Mass ejection and time variability in protostellar outflows: Cep E. SOLIS XVI

Protostellar jets are an important agent of star formation feedback, tightly connected with the mass-accretion process. The history of jet formation and mass-ejection provides constraints on the mass accretion history and the nature of the driving source. We want to characterize the time-variability of the mass-ejection phenomena at work in the Class 0 protostellar phase, in order to better understand the dynamics of the outflowing gas and bring more constraints on the origin of the jet chemical composition and the mass-accretion history. We have observed the emission of the CO 2-1 and SO N_J=5_4-4_3 rotational transitions with NOEMA, towards the intermediate-mass Class 0 protostellar system Cep E. The CO high-velocity jet emission reveals a central component associated with high-velocity molecular knots, also detected in SO, surrounded by a collimated layer of entrained gas. The gas layer appears to accelerate along the main axis over a length scale delta_0 ~700 au, while its diameter gradually increases up to several 1000au at 2000au from the protostar. The jet is fragmented into 18 knots of mass ~10^-3 Msun, unevenly distributed between the northern and southern lobes, with velocity variations up to 15 km/s close to the protostar, well below the jet terminal velocities. The knot interval distribution is approximately bimodal with a scale of ~50-80yr close to the protostar and ~150-200yr at larger distances >12". The mass-loss rates derived from knot masses are overall steady, with values of 2.7x10^-5 Msun/yr (8.9x10^-6 Msun/yr) in the northern (southern) lobe. The interaction of the ambient protostellar material with high-velocity knots drives the formation of a molecular layer around the jet, which accounts for the higher mass-loss rate in the north. The jet dynamics are well accounted for by a simple precession model with a period of 2000yr and a mass-ejection period of 55yr.

preprint2022arXiv

Negative and Positive Feedback from a Supernova Remnant with SHREC: A detailed Study of the Shocked Gas in IC443

Supernova remnants (SNRs) contribute to regulate the star formation efficiency and evolution of galaxies. As they expand into the interstellar medium (ISM), they transfer vast amounts of energy and momentum that displace, compress and heat the surrounding material. Despite the extensive work in galaxy evolution models, it remains to be observationally validated to what extent the molecular ISM is affected by the interaction with SNRs. We use the first results of the ESO-ARO Public Spectroscopic Survey SHREC, to investigate the shock interaction between the SNR IC443 and the nearby molecular clump G. We use high sensitivity SiO(2-1) and H$^{13}$CO$^+$(1-0) maps obtained by SHREC together with SiO(1-0) observations obtained with the 40m telescope at the Yebes Observatory. We find that the bulk of the SiO emission is arising from the ongoing shock interaction between IC443 and clump G. The shocked gas shows a well ordered kinematic structure, with velocities blue-shifted with respect to the central velocity of the SNR, similar to what observed toward other SNR-cloud interaction sites. The shock compression enhances the molecular gas density, n(H$_2$), up to $>$10$^5$ cm$^{-3}$, a factor of >10 higher than the ambient gas density and similar to values required to ignite star formation. Finally, we estimate that up to 50\% of the momentum injected by IC443 is transferred to the interacting molecular material. Therefore the molecular ISM may represent an important momentum carrier in sites of SNR-cloud interactions.

preprint2022arXiv

The protoplanetary disk population in the rho-Ophiuchi region L1688 and the time evolution of Class II YSOs

(Abridged) We present a study of the disk population in L1688, the densest and youngest region in Ophiuchus, and we compare it with other nearby regions of different age, namely Lupus, Chamaeleon I, Corona Australis, Taurus and Upper Scorpius. We select our L1688 sample using a combination of criteria (ALMA data, Gaia, optical/near-IR spectroscopy) and determine stellar and disk properties, specifically stellar mass (Mstar), average population age, mass accretion rate (Macc) and disk dust mass (Mdust). a) In L1688 the relations between Macc and Mstar, Mdust and Mstar, and Macc and Mdust have a roughly linear trend with slopes 1.8-1.9 for the first two relations and ~1 for the third, similarly to what found in the other regions. b) When ordered according to the characteristic age of each region, Macc decreases as 1/t, when corrected for the different stellar mass content; Mdust follows roughly the same trend between 0.5 and 5 Myr, but has an increase of a factor ~3 at ages of 2-3 Myr. We suggest that this could result from an earlier planet formation, followed by collisional fragmentation that temporarily replenishes the millimeter-size grain population. c) The dispersion of Macc and Mdust around the best-fitting relation with Mstar, as well as that of Macc versus Mdust are large: we find that the dispersions have continuous distributions with a log-normal shape and similar width (~0.8 dex). The amount of dust observed at ~1 Myr does not appear to be sufficient to assemble the majority of planetary systems, which suggests an earlier planetary cores formation. The dust mass traces to a large extent the disk gas mass evolution. Two properties remain puzzling: the steep dependence of Macc and Mdust on Mstar and the cause of the large dispersion in the three relations analyzed in this paper, in particular the one of the Macc versus Mdust relation.

preprint2021arXiv

ALMA chemical survey of disk-outflow sources in Taurus (ALMA-DOT) VI: Accretion shocks in the disk of DG Tau and HL Tau

Planet-forming disks are not isolated systems. Their interaction with the surrounding medium affects their mass budget and chemical content. In the context of the ALMA-DOT program, we obtained high-resolution maps of assorted lines from six disks that are still partly embedded in their natal envelope. In this work, we examine the SO and SO$_2$ emission that is detected from four sources: DG Tau, HL Tau, IRAS 04302+2247, and T Tau. The comparison with CO, HCO$^+$, and CS maps reveals that the SO and SO$_2$ emission originates at the intersection between extended streamers and the planet-forming disk. Two targets, DG Tau and HL Tau, offers clear cases of inflowing material inducing an accretion shock on the disk material. The measured rotational temperatures and radial velocities are consistent with this view. In contrast to younger Class 0 sources, these shocks are confined to the specific disk region impacted by the streamer. In HL Tau, the known accreting streamer induces a shock in the disk outskirt, and the released SO and SO$_2$ molecules spiral toward the star in a few hundreds years. These results suggest that shocks induced by late accreting material may be common in the disks of young star-forming regions with possible consequences on the chemical composition and mass content of the disk. They also highlight the importance of SO and SO$_2$ line observations to probe accretion shocks from a larger sample.

preprint2020arXiv

ALMA chemical survey of disk-outflow sources in Taurus (ALMA-DOT) III: The interplay between gas and dust in the protoplanetary disk of DG Tau

Planets form in protoplanetary disks and inherit their chemical composition. It is therefore crucial to understand the disks molecular content. We aim to characterize the distribution and abundance of molecules in the disk of DG Tau. In the context of the ALMA chemical survey of Disk-Outflow sources in Taurus (ALMA-DOT) we analyse ALMA observations of the disk of DG Tau in H2CO 3(1,2)-2(1,1), CS 5-4, and CN 2-1 at ~0.15", i.e. ~18 au at 121 pc. H2CO and CS originate from a disk ring at the edge of the 1.3mm dust continuum, with CS probing an outer disk region with respect to H2CO (peaking at ~70 and ~60 au, respectively). CN originates from an outermost disk/envelope region peaking at ~80 au. H2CO is dominated by disk emission, while CS probes also two streams of material possibly accreting onto the disk with a peak of emission where the stream connects to the disk. The ring- and disk-height- averaged column densities are ~2.4-8.6e13 cm-2 (H2CO), ~1.7-2.5e13 cm-2 (CS), and ~1.9-4.7e13 cm-2 (CN). Unsharp masking reveals a ring of enhanced dust emission at ~40 au, i.e. just outside the CO snowline (~30 au). CS and H2CO emissions are co-spatial suggesting that they are chemically linked. The observed rings of molecular emission at the edge of the 1.3mm continuum may be due to dust opacity effects and/or continnum over-subtraction in the inner disk; as well as to increased UV penetration and/or temperature inversion at the edge of the mm-dust which would cause an enhanced gas-phase formation and desorption of these molecules. Moreover, H2CO and CS originate from outside the ring of enhanced dust emission, which also coincides with a change of the linear polarization at 0.87mm. This suggests that outside the CO snowline there could be a change of the dust properties which would reflect in the increase of the intensity (and change of polarization) of continuum, and of molecular emission.

preprint2020arXiv

ALMA chemical survey of disk-outflow sources in Taurus (ALMA-DOT). IV. Thioformaldehyde (H$_2$CS) in protoplanetary disks: spatial distributions and binding energies

Aims: To trace the radial and vertical spatial distribution of H2CS, a key species of the S-bearing chemistry, in protoplanetary disks. To analyse the observed distributions in light of the H2CS binding energy, in order to discuss the role of thermal desorption in enriching the gas disk component. Methods: In the context of the ALMA chemical survey of Disk-Outflow sources in the Taurus star forming region (ALMA-DOT), we observed five Class I or early Class II sources with the o-H2CS(7_1,6-6_1,5) line on a 40 au scale. We estimated the binding energy (BEs) of H2CS using quantum mechanical calculations, for the first time, for an extended, periodic, crystalline ice. Results: We imaged H2CS in two rotating molecular rings in the HL Tau and IRAS04302+2247 disks. The outer radii are about 140 au (HL Tau), and 115 au (IRAS 04302+2247). The edge-on geometry of IRAS 04302+2247 reveals that H2CS emission peaks, at radii of 60-115 au, at z = +- 50 au from the equatorial plane. The column densities are about 10^14 cm^-2. For HL Tau, we derive, for the first time, the [H2CS]/[H] abundance in a protoplanetary disk (about 10^-14). The BEs of H2CS computed for extended crystalline ice and amorphous ices is 4258 K and 3000-4600 K, respectively, implying a thermal evaporation where dust temperature is larger than 50-80 K. Conclusions: H2CS traces the so-called warm molecular layer, a region previously sampled using CS, and H2CO. Thioformaldehyde peaks closer to the protostar than H2CO and CS, plausibly due to the relatively high-excitation level of observed 7_1,6-6_1,5 line (60 K). The H2CS BEs implies that thermal desorption dominates in thin, au-sized, inner and/or upper disk layers, indicating that the observed H2CS emitting up to radii larger than 100 au is likely injected in the gas due to non-thermal processes.

preprint2020arXiv

Angular momentum profiles of Class 0 protostellar envelopes

[abridged] Understanding how the infalling gas redistribute most of its initial angular momentum inherited from prestellar cores before reaching the stellar embryo is a key question. Disk formation has been naturally considered as a possible solution to this "angular momentum problem". However, how the initial angular momentum of protostellar cores is distributed and evolves during the main accretion phase and the beginning of disk formation has largely remained unconstrained up to now. In the framework of the IRAM CALYPSO survey, we used high dynamic range C$^{18}$O (2-1) and N$_2$H$^+$ (1-0) observations to quantify the distribution of specific angular momentum along the equatorial axis in a sample of 12 Class 0 protostellar envelopes from scales ~50 to 10000 au. The radial distributions of specific angular momentum in the CALYPSO sample suggest two distinct regimes within protostellar envelopes: the specific angular momentum decreases as $j \propto r^{1.6 \pm 0.2}$ down to ~1600 au and then tends to become relatively constant around 6 $\times$ 10$^{-4}$ km s$^{-1}$ pc down to ~50 au. The values of specific angular momentum measured in the inner Class 0 envelopes, namely that of the material directly involved in the star formation process ($<$1600 au), is on the same order of magnitude as what is inferred in small T-Tauri disks. Thus, disk formation appears to be a direct consequence of angular momentum conservation during the collapse. Our analysis reveals a dispersion of the directions of velocity gradients at envelope scales $>$1600 au, suggesting that they may not be related to rotational motions of the envelopes. We conclude that the specific angular momentum observed at these scales could find its origin in core-forming motions (infall, turbulence) or trace an imprint of the initial conditions for the formation of protostellar cores.

preprint2020arXiv

Bipolar molecular outflow of the very low-mass star Par-Lup3-4

Very low-mass stars are known to have jets and outflows, which is indicative of a scaled-down version of low-mass star formation. However, only very few outflows in very low-mass sources are well characterized. We characterize the bipolar molecular outflow of the very low-mass star Par-Lup3-4, a 0.12 M$_{\odot}$ object known to power an optical jet. We observed Par-Lup3-4 with ALMA in Bands 6 and 7, detecting both the continuum and CO molecular gas. In particular, we studied three main emission lines: CO(2-1), CO(3-2), and $^{13}$CO(3-2). Our observations reveal for the first time the base of a bipolar molecular outflow in a very low-mass star, as well as a stream of material moving perpendicular to the primary outflow of this source. The primary outflow morphology is consistent with the previously determined jet orientation and disk inclination. The outflow mass is $9.5\times10^{-7}\mathrm{M}_{\odot}$ , with an outflow rate of $4.3\times10^{-9}\mathrm{M}_{\odot}\mathrm{yr}^{-1}$ A new fitting to the spectral energy distribution suggests that Par-Lup3-4 may be a binary system. We have characterized Par-Lup3-4 in detail, and its properties are consistent with those reported in other very low-mass sources. This source provides further evidence that very low-mass sources form as a scaled-down version of low-mass stars.

preprint2020arXiv

Detectability of embedded protoplanets from hydrodynamical simulations

We predict magnitudes for young planets embedded in transition discs, still affected by extinction due to material in the disc. We focus on Jupiter-size planets at a late stage of their formation, when the planet has carved a deep gap in the gas and dust distributions and the disc starts being transparent to the planet flux in the infrared (IR). Column densities are estimated by means of three-dimensional hydrodynamical models, performed for several planet masses. Expected magnitudes are obtained by using typical extinction properties of the disc material and evolutionary models of giant planets. For the simulated cases located at $5.2$ AU in a disc with local unperturbed surface density of $127$ $\mathrm{g} \cdot \mathrm{cm}^{-2}$, a $1$ $M_J$ planet is highly extincted in J-, H- and K-bands, with predicted absolute magnitudes $\ge 50$ mag. In L- and M-bands extinction decreases, with planet magnitudes between $25$ and $35$ mag. In the N-band, due to the silicate feature on the dust opacities, the expected magnitude increases to $40$ mag. For a $2$ $M_J$ planet, the magnitudes in J-, H- and K-bands are above $22$ mag, while for L-, M- and N-bands the planet magnitudes are between $15$ and $20$ mag. For the $5$ $M_J$ planet, extinction does not play a role in any IR band, due to its ability to open deep gaps. Contrast curves are derived for the transition discs in CQ Tau, PDS70, HL Tau, TW Hya and HD163296. Planet mass upper-limits are estimated for the known gaps in the last two systems.

preprint2020arXiv

FAUST I. The hot corino at the heart of the prototypical Class I protostar L1551 IRS5

The study of hot corinos in Solar-like protostars has been so far mostly limited to the Class 0 phase, hampering our understanding of their origin and evolution. In addition, recent evidence suggests that planet formation starts already during Class I phase, which, therefore, represents a crucial step in the future planetary system chemical composition. Hence, the study of hot corinos in Class I protostars has become of paramount importance. Here we report the discovery of a hot corino towards the prototypical Class I protostar L1551 IRS5, obtained within the ALMA Large Program FAUST. We detected several lines from methanol and its isopotologues ($^{13}$CH$_{\rm 3}$OH and CH$_{\rm 2}$DOH), methyl formate and ethanol. Lines are bright toward the north component of the IRS5 binary system, and a possible second hot corino may be associated with the south component. The methanol lines non-LTE analysis constrains the gas temperature ($\sim$100 K), density ($\geq$1.5$\times$10$^{8}$ cm$^{-3}$), and emitting size ($\sim$10 au in radius). All CH$_{\rm 3}$OH and $^{13}$CH$_{\rm 3}$OH lines are optically thick, preventing a reliable measure of the deuteration. The methyl formate and ethanol relative abundances are compatible with those measured in Class 0 hot corinos. Thus, based on the present work, little chemical evolution from Class 0 to I hot corinos occurs.

preprint2020arXiv

Questioning the spatial origin of complex organic molecules in young protostars with the CALYPSO survey

Complex organic molecules (COMs) have been detected in a few Class 0 protostars but their origin is not well understood. Going beyond studies of individual objects, we want to investigate the origin of COMs in young protostars on a statistical basis. We use the CALYPSO survey performed with the IRAM PdBI to search for COMs at high angular resolution in a sample of 26 solar-type protostars, including 22 Class 0 and four Class I objects. Methanol is detected in 12 sources and tentatively in one source, which represents half of the sample. Eight sources (30%) have detections of at least three COMs. We find a strong chemical differentiation in multiple systems with five systems having one component with at least three COMs detected but the other component devoid of COM emission. The internal luminosity is found to be the source parameter impacting the most the COM chemical composition of the sources, while there is no obvious correlation between the detection of COM emission and that of a disk-like structure. A canonical hot-corino origin may explain the COM emission in four sources, an accretion-shock origin in two or possibly three sources, and an outflow origin in three sources. The CALYPSO sources with COM detections can be classified into three groups on the basis of the abundances of oxygen-bearing molecules, cyanides, and CHO-bearing molecules. These chemical groups correlate neither with the COM origin scenarii, nor with the evolutionary status of the sources if we take the ratio of envelope mass to internal luminosity as an evolutionary tracer. We find strong correlations between molecules that are a priori not related chemically (for instance methanol and methyl cyanide), implying that the existence of a correlation does not imply a chemical link. [abridged]

preprint2020arXiv

Seeds of Life in Space (SOLIS).VII. Discovery of a cold dense methanol blob toward the L1521F VeLLO system

The SOLIS (Seeds Of Life In Space) IRAM/NOEMA Large Program aims at studying a set of crucial complex organic molecules in a sample of sources, with well-known physical structure, covering the various phases of Solar-type star formation. One representative object of the transition from the prestellar core to the protostar phases has been observed toward the Very Low Luminosity Object (VeLLO) called L1521F. This type of source is important to study to make the link between prestellar cores and Class 0 sources and also to constrain the chemical evolution during the process of star formation. Two frequency windows (81.6-82.6 GHz and 96.65-97.65 GHz) were used to observe the emission from several complex organics toward the L1521F VeLLO. Only 2 transitions of methanol (A+, E2) have been detected in the narrow window centered at 96.7 GHz (with an upper limit on E1) in a very compact emission blob (~7'' corresponding to ~1000au) toward the NE of the L1521F protostar. The CS 2-1 transition is also detected within the WideX bandwidth. Consistently, with what has been found in prestellar cores, the methanol emission appears ~1000au away from the dust peak. The location of the methanol blob coincides with one of the filaments previously reported in the literature. The Tex of the gas inferred from methanol is (10$\pm$2) K, while the H2 gas density (estimated from the detected CS 2-1 emission and previous CS 5-4 ALMA obs.) is a factor >25 higher than the density in the surrounding environment (n(H2) >10$^{7}$ cm$^{-3}$). From its compactness, low excitation temperature and high gas density, we suggest that the methanol emission detected with NOEMA is either a cold and dense shock-induced blob, recently formed ($\leq$ few hundred years) by infalling gas or a cold and dense fragment that may have just been formed as a result of the intense gas dynamics found within the L1521F VeLLO system.

preprint2020arXiv

The Forgotten Quadrant Survey. $^{12}$CO and $^{13}$CO (1-0) survey of the Galactic Plane in the range 220°$<l<$240° -2.5°$<b<$0°

We present the Forgotten Quadrant Survey (FQS), an ESO large project that used the 12m antenna of the Arizona Radio Observatory to map the Galactic Plane in the range 220°$<l<$240°and -2.5°$<b<$0°, both in $^{12}$CO(1-0) and $^{13}$CO(1-0), at a spectral resolution of 0.65 km s$^{-1}$ and 0.26 km s$^{-1}$. Our dataset allows us to easily identify how the molecular dense gas is organised at different spatial scales: from the giant clouds with their denser filamentary networks, down to the clumps and cores that host the newborn stars and to obtain reliable estimates of their key physical parameters. We present the first release of the FQS data and discuss their quality. Spectra with 0.65 km s$^{-1}$ velocity channels have a noise ranging from 0.8 K to 1.3 K for $^{12}$CO(1-0) and from 0.3 K to 0.6 K for $^{13}$CO(1-0). In this paper, we used the $^{12}$CO(1-0) spectral cubes to produce a catalogue of 263 molecular clouds. This is the first selfconsistent, statistical catalogue of molecular clouds of the outer Galaxy, obtained with a subarcminute spatial resolution and therefore able to detect not only the classical giant molecular clouds, but also the small clouds and to resolve the cloud structure at the subparsec scale up to a distance of a few kpc. We found two classes of objects: structures with size above a few parsecs that are typical molecular clouds and may be self-gravitating, and subparsec structures that cannot be in gravitational equilibrium and are likely transient or confined by external pressure. We used the ratio between the Herschel H$_2$ column density and the integrated intensity of the CO lines to calculate the CO conversion factor and we found mean values of (3.3$\pm$1.4)$\times 10^{20}$ cm$^{-2}$(K km s$^{-1})^{-1}$ and (1.2$\pm$0.4)$\times 10^{21}$ cm$^{-2}$(K km s$^{-1})^{-1}$, for $^{12}$CO(1-0) and $^{13}$CO(1-0), respectively.

preprint2020arXiv

The GUAPOS project: G31.41+0.31 Unbiased ALMA sPectral Observational Survey -- I. Isomers of C$_{2}$H$_{4}$O$_{2}$

Understanding the degree of chemical complexity that can be reached in star-forming regions, together with the identification of precursors of the building blocks of life in the interstellar medium, is one of the goals of astrochemistry. Unbiased spectral surveys with large bandwidth and high spectral resolution are thus needed, to resolve line blending in chemically rich sources and identify complex organic molecules. This kind of observations has been successfully carried out, mainly towards the Galactic Center, a region that shows peculiar environmental conditions. We present an unbiased spectral survey at 3mm of one of the most chemically rich hot molecular cores located outside the Galactic Center, in the high-mass star-forming region G31.41+0.31. In this first paper, we present the survey and discuss the detection of the 3 isomers of C$_{2}$H$_{4}$O$_{2}$: methyl formate, glycolaldehyde and acetic acid. Observations were carried out with ALMA and cover the entire Band 3 from 84 to 116 GHz with an angular resolution of $1.2^{''}$x$1.2^{''}$ and a spectral resolution of $\sim0.488$ MHz. The transitions of the 3 molecules have been analyzed with the software XCLASS. All three isomers were detected and methyl formate and acetic acid abundances in G31 are the highest detected up to now, if compared to sources in literature. The size of the emission varies among the three isomers with acetic acid showing the most compact emission while methyl formate the most extended. The comparison with chemical models suggests the necessity of grain-surface routes for the formation of methyl formate in G31, while for glycolaldehyde both scenarios could be feasible. Proposed grain-surface reaction for acetic acid is not able to reproduce the observed abundance in this work, while gas-phase scenario should be further tested due to large uncertainties.

preprint2020arXiv

The Science Case for ALMA Band 2 and Band 2+3

We discuss the science drivers for ALMA Band 2 which spans the frequency range from 67 to 90 GHz. The key science in this frequency range are the study of the deuterated molecules in cold, dense, quiescent gas and the study of redshifted emission from galaxies in CO and other species. However, Band 2 has a range of other applications which are also presented. The science enabled by a single receiver system which would combine ALMA Bands 2 and 3 covering the frequency range 67 to 116 GHz, as well as the possible doubling of the IF bandwidth of ALMA to 16 GHz, are also considered.

preprint2020arXiv

Wideband 67-116 GHz receiver development for ALMA Band 2

ALMA has been operating since 2011, but has not yet been populated with the full suite of intended frequency bands. In particular, ALMA Band 2 (67-90 GHz) is the final band in the original ALMA band definition to be approved for production. We aim to produce a wideband, tuneable, sideband-separating receiver with 28 GHz of instantaneous bandwidth per polarisation operating in the sky frequency range 67-116 GHz. Our design anticipates new ALMA requirements following the recommendations in the 2030 ALMA Development Roadmap. The cryogenic cartridge is designed to be compatible with the ALMA Band 2 cartridge slot, where the coldest components -- the feedhorns, orthomode transducers, and cryogenic low noise amplifiers -- operate at a temperature of 15 K. We use multiple simulation methods and tools to optimise our designs for both the passive optics and the active components. The cryogenic cartridge interfaces with a room temperature cartridge hosting the local oscillator (LO) and the downconverter module. This warm cartridge is largely based on GaAs semiconductor technology and is optimised to match the cryogenic receiver bandwidth with the required instantaneous LO tuning range. Our collaboration has designed, fabricated, and tested multiple technical solutions for each of the components, producing a state-of-the-art receiver covering the full ALMA Band 2 & 3 atmospheric window. The receiver is suitable for deployment on ALMA in the coming years, and is capable of dual-polarisation, sideband-separating observations in intermediate frequency bands spanning 4-18 GHz, for a total of 28 GHz on-sky bandwidth per polarisation channel. We conclude that the 67-116 GHz wideband implementation for ALMA Band 2 is now feasible, and this receiver is a compelling instrumental upgrade that will enhance observational capabilities and scientific reach.

preprint2020arXiv

X-Shooter survey of disk accretion in Upper Scorpius I. Very high accretion rates at age>5 Myr

Determining the mechanisms that drive the evolution of protoplanetary disks is a necessary step to understand how planets form. Here we measured the mass accretion rate for young stellar objects at age >5 Myr, a critical test for the current models of disk evolution. We present the analysis of the spectra of 36 targets in the ~5-10 Myr old Upper Scorpius region for which disk masses were measured with ALMA. We find that the mass accretion rates in this sample of old but still survived disks are similarly high as those of the younger (<3 Myr old) star-forming regions of Lupus and Cha I, when considering the dependence on stellar and disk mass. In particular, several disks show high mass accretion rates >10^-9 Msun/yr while having low disk masses. Furthermore, the median values of the measured mass accretion rates in the disk mass ranges where our sample is complete at a level ~60-80% are compatible in these three regions. At the same time, the spread of mass accretion rates at any given disk mass is still >0.9 dex even at age>5 Myr. These results are in contrast with simple models of viscous evolution, which would predict that the values of the mass accretion rate diminish with time, and a tighter correlation with disk mass at age>5 Myr. Similarly, simple models of internal photoevaporation cannot reproduce the observed mass accretion rates, while external photoevaporation might explain the low disk masses and high accretion rates. A partial possible solution to the discrepancy with the viscous models is that the gas-to-dust ratio of the disks at >5 Myr is significantly different and higher than the canonical 100, as suggested by some dust and gas disk evolution models. The results shown here require the presence of several inter-playing processes, such as detailed dust evolution, external photoevaporation and possibly MHD winds, to explain the secular evolution of protoplanetary disks.

preprint2019arXiv

ALMA and ROSINA detections of phosphorus-bearing molecules: the interstellar thread between star-forming regions and comets

To understand how Phosphorus-bearing molecules are formed in star-forming regions, we have analysed ALMA observations of PN and PO towards the massive star-forming region AFGL 5142, combined with a new analysis of the data of the comet 67P/Churyumov-Gerasimenko taken with the ROSINA instrument onboard Rosetta. The ALMA maps show that the emission of PN and PO arises from several spots associated with low-velocity gas with narrow linewidths in the cavity walls of a bipolar outflow. PO is more abundant than PN in most of the spots, with the PO/PN ratio increasing as a function of the distance to the protostar. Our data favor a formation scenario in which shocks sputter phosphorus from the surface of dust grains, and gas-phase photochemistry induced by UV photons from the protostar allows efficient formation of the two species in the cavity walls. Our analysis of the ROSINA data has revealed that PO is the main carrier of P in the comet, with PO/PN>10. Since comets may have delivered a significant amount of prebiotic material to the early Earth, this finding suggests that PO could contribute significantly to the phosphorus reservoir during the dawn of our planet. There is evidence that PO was already in the cometary ices prior to the birth of the Sun, so the chemical budget of the comet might be inherited from the natal environment of the Solar System, which is thought to be a stellar cluster including also massive stars.

preprint2019arXiv

Demographics of disks around young very low-mass stars and brown dwarfs in Lupus

We present new 890 $μm$ continuum ALMA observations of 5 brown dwarfs (BDs) with infrared excess in Lupus I and III -- which, in combination with 4 BDs previously observed, allowed us to study the mm properties of the full known BD disk population of one star-forming region. Emission is detected in 5 out of the 9 BD disks. Dust disk mass, brightness profiles and characteristic sizes of the BD population are inferred from continuum flux and modeling of the observations. Only one source is marginally resolved, allowing for the determination of its disk characteristic size. We conduct a demographic comparison between the properties of disks around BDs and stars in Lupus. Due to the small sample size, we cannot confirm or disprove if the disk mass over stellar mass ratio drops for BDs, as suggested for Ophiuchus. Nevertheless, we find that all detected BD disks have an estimated dust mass between 0.2 and 3.2 $M_{\bigoplus}$; these results suggest that the measured solid masses in BD disks can not explain the observed exoplanet population, analogous to earlier findings on disks around more massive stars. Combined with the low estimated accretion rates, and assuming that the mm-continuum emission is a reliable proxy for the total disk mass, we derive ratios of $\dot{M}_{\mathrm{acc}} / M_{\mathrm{disk}}$ significantly lower than in disks around more massive stars. If confirmed with more accurate measurements of disk gas masses, this result could imply a qualitatively different relationship between disk masses and inward gas transport in BD disks.

preprint2019arXiv

The Hi-GAL catalogue of dusty filamentary structures in the Galactic Plane

The recent data collected by {\it Herschel} have confirmed that interstellar structures with filamentary shape are ubiquitously present in the Milky Way. Filaments are thought to be formed by several physical mechanisms acting from the large Galactic scales down to the sub-pc fractions of molecular clouds, and they might represent a possible link between star formation and the large-scale structure of the Galaxy. In order to study this potential link, a statistically significant sample of filaments spread throughout the Galaxy is required. In this work we present the first catalogue of $32,059$ candidate filaments automatically identified in the Hi-GAL survey of the entire Galactic Plane. For these objects we determined morphological (length, $l^{a}$, and geometrical shape) and physical (average column density, $N_{\rm H_{2}}$, and average temperature, $T$) properties. We identified filaments with a wide range of properties: 2$'$\,$\leq l^{a}\leq$\, 100$'$, $10^{20} \leq N_{\rm H_{2}} \leq 10^{23}$\,cm$^{-2}$ and $10 \leq T\leq$ 35\,K. We discuss their association with the Hi-GAL compact sources, finding that the most tenuous (and stable) structures do not host any major condensation and we also assign a distance to $\sim 18,400$ filaments for which we determine mass, physical size, stability conditions and Galactic distribution. When compared to the spiral arms structure, we find no significant difference between the physical properties of on-arm and inter-arm filaments. We compared our sample with previous studies, finding that our Hi-GAL filament catalogue represents a significant extension in terms of Galactic coverage and sensitivity. This catalogue represents an unique and important tool for future studies devoted to understanding the filament life-cycle.

preprint2016arXiv

A Steeper than Linear Disk Mass-Stellar Mass Scaling Relation

The disk mass is among the most important input parameter for every planet formation model to determine the number and masses of the planets that can form. We present an ALMA 887micron survey of the disk population around objects from 2 to 0.03Msun in the nearby 2Myr-old Chamaeleon I star-forming region. We detect thermal dust emission from 66 out of 93 disks, spatially resolve 34 of them, and identify two disks with large dust cavities of about 45AU in radius. Assuming isothermal and optically thin emission, we convert the 887micron flux densities into dust disk masses, hereafter Mdust. We find that the Mdust-Mstar relation is steeper than linear with power law indices 1.3-1.9, where the range reflects two extremes of the possible relation between the average dust temperature and stellar luminosity. By re-analyzing all millimeter data available for nearby regions in a self-consistent way, we show that the 1-3 Myr-old regions of Taurus, Lupus, and Chamaeleon I share the same Mdust-Mstar relation, while the 10Myr-old Upper Sco association has a steeper relation. Theoretical models of grain growth, drift, and fragmentation reproduce this trend and suggest that disks are in the fragmentation-limited regime. In this regime millimeter grains will be located closer in around lower-mass stars, a prediction that can be tested with deeper and higher spatial resolution ALMA observations.

preprint2016arXiv

Dust properties across the CO snowline in the HD 163296 disk from ALMA and VLA observations

To characterize the mechanisms of planet formation it is crucial to investigate the properties and evolution of protoplanetary disks around young stars, where the initial conditions for the growth of planets are set. Our goal is to study grain growth in the disk of the young, intermediate mass star HD163296 where dust processing has already been observed, and to look for evidence of growth by ice condensation across the CO snowline, already identified in this disk with ALMA. Under the hypothesis of optically thin emission we compare images at different wavelengths from ALMA and VLA to measure the opacity spectral index across the disk and thus the maximum grain size. We also use a Bayesian tool based on a two-layer disk model to fit the observations and constrain the dust surface density. The measurements of the opacity spectral index indicate the presence of large grains and pebbles ($\geq$1 cm) in the inner regions of the disk (inside $\sim$50 AU) and smaller grains, consistent with ISM sizes, in the outer disk (beyond 150 AU). Re-analysing ALMA Band 7 Science Verification data we find (radially) unresolved excess continuum emission centered near the location of the CO snowline at $\sim$90 AU. Our analysis suggests a grain size distribution consistent with an enhanced production of large grains at the CO snowline and consequent transport to the inner regions. Our results combined with the excess in infrared scattered light found by Garufi et al. (2014) suggests the presence of a structure at 90~AU involving the whole vertical extent of the disk. This could be evidence for small scale processing of dust at the CO snowline.

preprint2016arXiv

Evidence for a correlation between mass accretion rates onto young stars and the mass of their protoplanetary disks

A relation between the mass accretion rate onto the central young star and the mass of the surrounding protoplanetary disk has long been theoretically predicted and observationally sought. For the first time, we have accurately and homogeneously determined the photospheric parameters, mass accretion rate, and disk mass for an essentially complete sample of young stars with disks in the Lupus clouds. Our work combines the results of surveys conducted with VLT/X-Shooter and ALMA. With this dataset we are able to test a basic prediction of viscous accretion theory, the existence of a linear relation between the mass accretion rate onto the central star and the total disk mass. We find a correlation between the mass accretion rate and the disk dust mass, with a ratio that is roughly consistent with the expected viscous timescale when assuming an interstellar medium (ISM) gas-to-dust ratio. This confirms that mass accretion rates are related to the properties of the outer disk. We find no correlation between mass accretion rates and the disk mass measured by CO isotopologues emission lines, possibly owing to the small number of measured disk gas masses. This suggests that the mm-sized dust mass better traces the total disk mass and that masses derived from CO may be underestimated, at least in some cases.

preprint2016arXiv

First image of the L1157 molecular jet by the CALYPSO IRAM-PdBI survey

Fast jets are thought to be a crucial ingredient of star formation because they might extract angular momentum from the disk and thus allow mass accretion onto the star. However, it is unclear whether jets are ubiquitous, and likewise, their contribution to mass and angular momentum extraction during protostar formation remains an open question. Our aim is to investigate the ejection process in the low-mass Class 0 protostar L1157. This source is associated with a spectacular bipolar outflow, and the recent detection of high-velocity SiO suggests the occurrence of a jet. Observations of CO 2-1 and SiO 5-4 at 0.8" resolution were obtained with the IRAM Plateau de Bure Interferometer as part of the CALYPSO large program. The jet and outflow structure were fit with a precession model. We derived the column density of CO and SiO, as well as the jet mass-loss rate and mechanical luminosity. High-velocity CO and SiO emission resolve for the first time the first 200 au of the outflow-driving molecular jet. The jet is strongly asymmetric, with the blue lobe 0.65 times slower than the red lobe. This suggests that the large-scale asymmetry of the outflow is directly linked to the jet velocity and that the asymmetry in the launching mechanism has been at work for the past 1800 yr. Velocity asymmetries are common in T Tauri stars, which suggests that the jet formation mechanism from Class 0 to Class II stages might be similar. Our model simultaneously fits the inner jet and the clumpy 0.2 pc scale outflow by assuming that the jet precesses counter-clockwise on a cone inclined by 73 degree to the line of sight with an opening angle of 8 degree on a period of 1640 yr. The estimated jet mass flux and mechanical luminosity are 7.7e-7 Msun/yr, and 0.9 Lsun, indicating that the jet could extract at least 25% of the gravitational energy released by the forming star.

preprint2016arXiv

Hi-GAL, the Herschel infrared Galactic Plane Survey: photometric maps and compact source catalogues. First data release for Inner Milky Way: +68°> l > -70°

(Abridged) We present the first public release of high-quality data products (DR1) from Hi-GAL, the {\em Herschel} infrared Galactic Plane Survey. Hi-GAL is the keystone of a suite of continuum Galactic Plane surveys from the near-IR to the radio, and covers five wavebands at 70, 160, 250, 350 and 500 micron, encompassing the peak of the spectral energy distribution of cold dust for 8 < T < 50K. This first Hi-GAL data release covers the inner Milky Way in the longitude range 68° > l > -70° in a |b|<1° latitude strip. Photometric maps have been produced with the ROMAGAL pipeline, that optimally capitalizes on the excellent sensitivity and stability of the bolometer arrays of the {\em Herschel} PACS and SPIRE photometric cameras, to deliver images of exquisite quality and dynamical range, absolutely calibrated with {\em Planck} and {\em IRAS}, and recovering extended emission at all wavelengths and all spatial scales. The compact source catalogues have been generated with the CuTEx algorithm, specifically developed to optimize source detection and extraction in the extreme conditions of intense and spatially varying background that are found in the Galactic Plane in the thermal infrared. Hi-GAL DR1 images will be accessible via a dedicated web-based image cutout service. The DR1 Compact Source Catalogues are delivered as single-band photometric lists containing, in addition to source position, peak and integrated flux and source sizes, a variety of parameters useful to assess the quality and reliability of the extracted sources, caveats and hints to help this assessment are provided. Flux completeness limits in all bands are determined from extensive synthetic source experiments and depend on the specific line of sight along the Galactic Plane. Hi-GAL DR1 catalogues contain 123210, 308509, 280685, 160972 and 85460 compact sources in the five bands, respectively.

preprint2016arXiv

Hot and dense water in the inner 25 AU of SVS13-A

In the context of the ASAI (Astrochemical Surveys At IRAM) project, we carried out an unbiased spectral survey in the millimeter window towards the well known low-mass Class I source SVS13-A. The high sensitivity reached (3-12 mK) allowed us to detect at least 6 HDO broad (FWHM ~ 4-5 km/s) emission lines with upper level energies up to Eu = 837 K. A non-LTE LVG analysis implies the presence of very hot (150-260 K) and dense (> 3 10^7 cm-3) gas inside a small radius ($\sim$ 25 AU) around the star, supporting, for the first time, the occurrence of a hot corino around a Class I protostar. The temperature is higher than expected for water molecules are sublimated from the icy dust mantles (~ 100 K). Although we cannot exclude we are observig the effects of shocks and/or winds at such small scales, this could imply that the observed HDO emission is tracing the water abundance jump expected at temperatures ~ 220-250 K, when the activation barrier of the gas phase reactions leading to the formation of water can be overcome. We derive X(HDO) ~ 3 10-6, and a H2O deuteration > 1.5 10-2, suggesting that water deuteration does not decrease as the protostar evolves from the Class 0 to the Class I stage.

preprint2016arXiv

Magnetically-regulated fragmentation of a massive, dense and turbulent clump

Massive stars, multiple stellar systems and clusters are born from the gravitational collapse of massive dense gaseous clumps, and the way these systems form strongly depends on how the parent clump fragments into cores during collapse. Numerical simulations show that magnetic fields may be the key ingredient in regulating fragmentation. Here we present ALMA observations at ~0.25'' resolution of the thermal dust continuum emission at ~278 GHz towards a turbulent, dense, and massive clump, IRAS 16061-5048c1, in a very early evolutionary stage. The ALMA image shows that the clump has fragmented into many cores along a filamentary structure. We find that the number, the total mass and the spatial distribution of the fragments are consistent with fragmentation dominated by a strong magnetic field. Our observations support the theoretical prediction that the magnetic field plays a dominant role in the fragmentation process of massive turbulent clump.

preprint2016arXiv

Molecular gas kinematics within the central 250 pc of the Milky Way

Using spectral-line observations of HNCO, N2H+, and HNC, we investigate the kinematics of dense gas in the central ~250 pc of the Galaxy. We present SCOUSE (Semi-automated multi-COmponent Universal Spectral-line fitting Engine), a line fitting algorithm designed to analyse large volumes of spectral-line data efficiently and systematically. Unlike techniques which do not account for complex line profiles, SCOUSE accurately describes the {l, b, v_LSR} distribution of CMZ gas, which is asymmetric about Sgr A* in both position and velocity. Velocity dispersions range from 2.6 km/s<σ<53.1 km/s. A median dispersion of 9.8 km/s, translates to a Mach number, M_3D>28. The gas is distributed throughout several "streams", with projected lengths ~100-250 pc. We link the streams to individual clouds and sub-regions, including Sgr C, the 20 and 50 km/s clouds, the dust ridge, and Sgr B2. Shell-like emission features can be explained by the projection of independent molecular clouds in Sgr C and the newly identified conical profile of Sgr B2 in {l ,b, v_LSR} space. These features have previously invoked supernova-driven shells and cloud-cloud collisions as explanations. We instead caution against structure identification in velocity-integrated emission maps. Three geometries describing the 3-D structure of the CMZ are investigated: i) two spiral arms; ii) a closed elliptical orbit; iii) an open stream. While two spiral arms and an open stream qualitatively reproduce the gas distribution, the most recent parameterisation of the closed elliptical orbit does not. Finally, we discuss how proper motion measurements of masers can distinguish between these geometries, and suggest that this effort should be focused on the 20 km/s and 50 km/s clouds and Sgr C.

preprint2016arXiv

Physical Properties of Galactic Planck Cold Cores revealed by the Hi-GAL survey

Previous studies of the initial conditions of massive star formation have mainly targeted Infrared-Dark Clouds (IRDCs) toward the inner Galaxy. This is due to the fact that IRDCs were first detected in absorption against the bright mid-IR background, requiring a favourable location to be observed. By selection, IRDCs represent only a fraction of the Galactic clouds capable of forming massive stars and star clusters. Due to their low dust temperatures, IRDCs are bright in the far-IR and millimeter and thus, observations at these wavelengths have the potential to provide a complete sample of star-forming massive clouds across the Galaxy. Our aim is to identify the clouds at the initial conditions of massive star formation across the Galaxy and compare their physical properties as a function of their Galactic location. We have examined the physical properties of a homogeneous galactic cold core sample obtained with the Planck satellite across the Galactic Plane. With the use of Herschel Hi-GAL observations, we have characterized the internal structure of them. By using background-subtracted Herschel images, we have derived the H2 column density and dust temperature maps for 48 Planck clumps. Their basic physical parameters have been calculated and analyzed as a function of location within the Galaxy. These properties have also been compared with the empirical relation for massive star formation derived by Kauffmann & Pillai (2010). Most of the Planck clumps contain signs of star formation. About 25% of them are massive enough to form high mass stars. Planck clumps toward the Galactic center region show higher peak column densities and higher average dust temperatures than those of the clumps in the outer Galaxy. Although we only have seven clumps without associated YSOs, the Hi-GAL data show no apparent differences in the properties of Planck cold clumps with and without star formation.

preprint2016arXiv

The link between turbulence, magnetic fields, filaments, and star formation in the Central Molecular Zone cloud G0.253+0.016

Star formation is primarily controlled by the interplay between gravity, turbulence, and magnetic fields. However, the turbulence and magnetic fields in molecular clouds near the Galactic Center may differ substantially from spiral-arm clouds. Here we determine the physical parameters of the central molecular zone (CMZ) cloud G0.253+0.016, its turbulence, magnetic field and filamentary structure. Using column-density maps based on dust-continuum emission observations with ALMA+Herschel, we identify filaments and show that at least one dense core is located along them. We measure the filament width W_fil=0.17$\pm$0.08pc and the sonic scale λ_sonic=0.15$\pm$0.11pc of the turbulence, and find W_fil~λ_sonic. A strong velocity gradient is seen in the HNCO intensity-weighted velocity maps obtained with ALMA+Mopra, which is likely caused by large-scale shearing of G0.253+0.016, producing a wide double-peaked velocity PDF. After subtracting the gradient to isolate the turbulent motions, we find a nearly Gaussian velocity PDF typical for turbulence. We measure the total and turbulent velocity dispersion, 8.8$\pm$0.2km/s and 3.9$\pm$0.1km/s, respectively. Using magnetohydrodynamical simulations, we find that G0.253+0.016's turbulent magnetic field B_turb=130$\pm$50$μ$G is only ~1/10 of the ordered field component. Combining these measurements, we reconstruct the dominant turbulence driving mode in G0.253+0.016 and find a driving parameter b=0.22$\pm$0.12, indicating solenoidal (divergence-free) driving. We compare this to spiral-arm clouds, which typically have a significant compressive (curl-free) driving component (b>0.4). Motivated by previous reports of strong shearing motions in the CMZ, we speculate that shear causes the solenoidal driving in G0.253+0.016 and show that this reduces the star formation rate (SFR) by a factor of 6.9 compared to typical nearby clouds.

preprint2016arXiv

The sulfur depletion problem: upper limits on the H2S2, HS2, and S2 gas-phase abundances toward the low-mass warm core IRAS16293-2422

A fraction of the missing sulfur in dense clouds and circumstellar regions could be in the form of three species not yet de- tected in the interstellar medium: H2S2, HS2, and S2 according to experimental simulations performed under astrophysically relevant conditions. These S-S bonded molecules can be formed by the energetic processing of H2S-bearing ice mantles on dust grains, and subsequently desorb to the gas phase. The detection of these species could partially solve the sulfur depletion problem, and would help to improve our knowledge of the poorly known chemistry of sulfur in the interstellar medium. To this purpose we performed dedicated ground-based observations toward the low-mass warm core IRAS16293-2422. Observations in the submillimeter regime were obtained with the APEX 12 m telescope during 15 hours of observation, targeting a wide selection of the predicted rotational transitions of the three molecules. The 1σ noise rms values were extracted in the spectral regions where the targeted species should have been detected. These values were a factor of 2-7 lower than those reached by previous observations toward the same source, and allowed us to estimate a 1σ upper limit to their molecular abundances of 8.1E-9 , 1.1E-8 , and 2.9E-7 relative to H2, for H2S2, HS2, and S2, respectively. The upper limit abundances of the three molecules containing the S2 unit are up to two orders of magnitude lower than the H2S abundance in the source, and one order of magnitude lower than the expected abundances from the experimental simulations using ice analogs. Subsequent gas-phase chemistry after desorption could lower the abundances of the three species to undetectable levels in our observations.

preprint2016arXiv

Understanding discs in binary YSOs: detailed modelling of VV CrA

Given that a majority of stars form in multiple systems, in order to fully understand the star- and planet-formation processes we must seek to understand them in multiple stellar systems. With this in mind, we present an analysis of the enigmatic binary T-Tauri system VV Corona Australis, in which both components host discs, but only one is visible at optical wavelengths. We seek to understand the peculiarities of this system by searching for a model for the binary which explains all the available continuum observations of the system. We present new mid-infrared interferometry and near-infrared spectroscopy along with archival millimetre-wave observations, which resolve the binary at 1.3mm for the first time. We compute a grid of pre-main-sequence radiative transfer models and calculate their posterior probabilities given the observed spectral energy distributions and mid-infrared interferometric visibilities of the binary components, beginning with the assumption that the only differences between the two components are their inclination and position angles. Our best-fitting solution corresponds to a relatively low luminosity T-tauri binary, with each component's disc having a large scale height and viewed at moderate inclination ($\sim 50^\circ$), with the infrared companion inclined by $\sim 5^\circ$ degrees more than the primary. Comparing the results of our model to evolutionary models suggests stellar masses $\sim 1.7$M$_\odot$ and an age for the system of 3.5Myr, towards the upper end of previous estimates. Combining these results with accretion indicators from near-IR spectroscopy, we determine an accretion rate of 4.0$\times 10^{-8}$M$_\odot$ yr$^{-1}$ for the primary. We suggest that future observations of VV~CrA and similar systems should prioritise high angular resolution sub-mm and near-IR imaging of the discs and high resolution optical/NIR spectroscopy of the central stars.

preprint2015arXiv

A cluster in the making: ALMA reveals the initial conditions for high-mass cluster formation

G0.253+0.016 is a molecular clump that appears to be on the verge of forming a high mass, Arches-like cluster. Here we present new ALMA observations of its small-scale (~0.07 pc) 3mm dust continuum and molecular line emission. The data reveal a complex network of emission features, the morphology of which ranges from small, compact regions to extended, filamentary structures that are seen in both emission and absorption. The dust column density is well traced by molecules with higher excitation energies and critical densities, consistent with a clump that has a denser interior. A statistical analysis supports the idea that turbulence shapes the observed gas structure within G0.253+0.016. We find a clear break in the turbulent power spectrum derived from the optically thin dust continuum emission at a spatial scale of ~0.1 pc, which may correspond to the spatial scale at which gravity has overcome the thermal pressure. We suggest that G0.253+0.016 is on the verge of forming a cluster from hierarchical, filamentary structures that arise from a highly turbulent medium. Although the stellar distribution within Arches-like clusters is compact, centrally condensed and smooth, the observed gas distribution within G0.253+0.016 is extended, with no high-mass central concentration, and has a complex, hierarchical structure. If this clump gives rise to a high-mass cluster and its stars are formed from this initially hierarchical gas structure, then the resulting cluster must evolve into a centrally condensed structure via a dynamical process.

preprint2015arXiv

A multi-wavelength analysis for interferometric (sub-)mm observations of protoplanetary disks: radial constraints on the dust properties and the disk structure

Theoretical models of grain growth predict dust properties to change as a function of protoplanetary disk radius, mass, age and other physical conditions. We lay down the methodology for a multi-wavelength analysis of (sub-)mm and cm continuum interferometric observations to constrain self-consistently the disk structure and the radial variation of the dust properties. The computational architecture is massively parallel and highly modular. The analysis is based on the simultaneous fit in the uv-plane of observations at several wavelengths with a model for the disk thermal emission and for the dust opacity. The observed flux density at the different wavelengths is fitted by posing constraints on the disk structure and on the radial variation of the grain size distribution. We apply the analysis to observations of three protoplanetary disks (AS 209, FT Tau, DR Tau) for which a combination of spatially resolved observations in the range ~0.88mm to ~10mm is available (from SMA, CARMA, and VLA), finding evidence of a decreasing maximum dust grain size (a_max) with radius. We derive large a_max values up to 1 cm in the inner disk between 15 and 30 AU and smaller grains with a_max~1 mm in the outer disk (R > 80AU). In this paper we develop a multi-wavelength analysis that will allow this missing quantity to be constrained for statistically relevant samples of disks and to investigate possible correlations with disk or stellar parameters.

preprint2015arXiv

A nebular analysis of the central Orion Nebula with MUSE

A nebular analysis of the central Orion Nebula and its main structures is presented. We exploit MUSE integral field observations in the wavelength range 4595-9366 Å to produce the first O, S and N ionic and total abundance maps of a region spanning 6' x 5' with a spatial resolution of 0.2". We use the S$_{23}$ ( = ([SII]$λ$6717,31+[SIII]$λ$9068)/H$β$) parameter, together with [OII]/[OIII] as an indicator of the degree of ionisation, to distinguish between the various small-scale structures. The only Orion Bullet covered by MUSE is HH 201, which shows a double component in the [FeII]$λ$8617 line throughout indicating an expansion, and we discuss a scenario in which this object is undergoing a disruptive event. We separate the proplyds located south of the Bright Bar into four categories depending on their S$_{23}$ values, propose the utility of the S$_{23}$ parameter as an indicator of the shock-contribution to the excitation of line-emitting atoms, and show that the MUSE data is able to identify the proplyds associated with disks and microjets. We compute the second order structure function for the H$α$, [OIII]$λ$5007, [SII]$λ$6731 and [OI]$λ$6300 emission lines to analyse the turbulent velocity field of the region covered with MUSE. We find that the spectral and spatial resolution of MUSE is not able to faithfully reproduce the structure functions of previous works.

preprint2015arXiv

A network of filaments detected by Herschel in the Serpens Core: A laboratory to test simulations of low-mass star formation

Filaments represent a key structure during the early stages of the star formation process. Simulations show filamentary structure commonly formed before and during the formation of cores. Aims. The Serpens Core represents an ideal laboratory to test the state-of-the-art of simulations of turbulent Giant Molecular Clouds. We use Herschel observations of the Serpens Core to compute temperature and column density maps of the region. Among the simulations of Dale et al. (2012), we select the early stages of their Run I, before stellar feedback is initiated, with similar total mass and physical size as the Serpens Core. We derive temperature and column density maps also from the simulations. The observed distribution of column densities of the filaments has been analysed first including and then masking the cores. The same analysis has been performed on the simulations as well. A radial network of filaments has been detected in the Serpens Core. The analysed simulation shows a striking morphological resemblance to the observed structures. The column density distribution of simulated filaments without cores shows only a log-normal distribution, while the observed filaments show a power-law tail. The power-law tail becomes evident in the simulation if one focuses just on the column density distribution of the cores. In contrast, the observed cores show a flat distribution. Even though the simulated and observed filaments are subjectively similar-looking, we find that they behave in very different ways. The simulated filaments are turbulence-dominated regions, the observed filaments are instead self-gravitating structures that will probably fragment into cores.

preprint2015arXiv

An ATCA survey of debris disks at 7 millimeters

We present ATCA continuum observations at a wavelength of 6.8 mm of five debris disks: $β$ Pictoris, q$^1$ Eridani, HD 107146, HD 181327, and HD 95086. These observations provide the detection at the longest wavelengths obtained to date for all these debris disks. By combining our 6.8 mm data with previous detections at shorter sub-millimeter/millimeter wavelengths we measure the long wavelength spectral index of these sources. We then use previous estimates for the temperature of the emitting dust to derive the spectral index of the dust emissivity. Under the assumption that all the detected flux comes from dust only, we constrain the slope of the solid size distribution, assumed to be a power-law. The values that we infer for the slope of the size distribution range between about 3.36 and 3.50. We compare our findings with the case of the Fomalhaut debris disk and use these results to test the predictions of collisional cascades of planetesimal belts.

preprint2015arXiv

An Overview of the 2014 ALMA Long Baseline Campaign

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

preprint2015arXiv

First Results from High Angular Resolution ALMA Observations Toward the HL Tau Region

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations from the 2014 Long Baseline Campaign in dust continuum and spectral line emission from the HL Tau region. The continuum images at wavelengths of 2.9, 1.3, and 0.87 mm have unprecedented angular resolutions of 0.075 arcseconds (10 AU) to 0.025 arcseconds (3.5 AU), revealing an astonishing level of detail in the circumstellar disk surrounding the young solar analogue HL Tau, with a pattern of bright and dark rings observed at all wavelengths. By fitting ellipses to the most distinct rings, we measure precise values for the disk inclination (46.72pm0.05 degrees) and position angle (+138.02pm0.07 degrees). We obtain a high-fidelity image of the 1.0 mm spectral index ($α$), which ranges from $α\sim2.0$ in the optically-thick central peak and two brightest rings, increasing to 2.3-3.0 in the dark rings. The dark rings are not devoid of emission, we estimate a grain emissivity index of 0.8 for the innermost dark ring and lower for subsequent dark rings, consistent with some degree of grain growth and evolution. Additional clues that the rings arise from planet formation include an increase in their central offsets with radius and the presence of numerous orbital resonances. At a resolution of 35 AU, we resolve the molecular component of the disk in HCO+ (1-0) which exhibits a pattern over LSR velocities from 2-12 km/s consistent with Keplerian motion around a ~1.3 solar mass star, although complicated by absorption at low blue-shifted velocities. We also serendipitously detect and resolve the nearby protostars XZ Tau (A/B) and LkHa358 at 2.9 mm.

preprint2015arXiv

Hunting for planets in the HL Tau disk

Recent ALMA images of HL Tau show gaps in the dusty disk that may be caused by planetary bodies. Given the young age of this system, if confirmed, this finding would imply very short timescales for planet formation, probably in a gravitationally unstable disk. To test this scenario, we searched for young planets by means of direct imaging in the L'-band using the Large Binocular Telescope Interferometer mid-infrared camera. At the location of two prominent dips in the dust distribution at ~70AU (~0.5") from the central star we reach a contrast level of ~7.5mag. We did not detect any point source at the location of the rings. Using evolutionary models we derive upper limits of ~10-15MJup at <=0.5-1Ma for the possible planets. With these sensitivity limits we should have been able to detect companions sufficiently massive to open full gaps in the disk. The structures detected at mm-wavelengths could be gaps in the distributions of large grains on the disk midplane, caused by planets not massive enough to fully open gaps. Future ALMA observations of the molecular gas density profile and kinematics as well as higher contrast infrared observations may be able to provide a definitive answer.

preprint2015arXiv

Italian Science Case for ALMA Band 2+3

The Premiale Project "Science and Technology in Italy for the upgraded ALMA Observatory - iALMA" has the goal of strengthening the scientific, technological and industrial Italian contribution to the Atacama Large Millimeter/submillimeter Array (ALMA), the largest ground based international infrastructure for the study of the Universe in the microwave. One of the main objectives of the Science Working Group (SWG) inside iALMA, the Work Package 1, is to develop the Italian contribution to the Science Case for the ALMA Band 2 or Band 2+3 receiver. ALMA Band 2 receiver spans from ~67 GHz (bounded by an opaque line complex of ozone lines) up to 90 GHz which overlaps with the lower frequency end of ALMA Band 3. Receiver technology has advanced since the original definition of the ALMA frequency bands. It is now feasible to produce a single receiver which could cover the whole frequency range from 67 GHz to 116 GHz, encompassing Band 2 and Band 3 in a single receiver cartridge, a so called Band 2+3 system. In addition, upgrades of the ALMA system are now foreseen that should double the bandwidth to 16 GHz. The science drivers discussed below therefore also discuss the advantages of these two enhancements over the originally foreseen Band 2 system.

preprint2015arXiv

Jet multiplicity in the proto-binary system NGC1333-IRAS4A. The detailed CALYPSO IRAM-PdBI view

Owing to the paucity of sub-arcsecond (sub)mm observations required to probe the innermost regions of newly forming protostars, several fundamental questions are still being debated, such as the existence and coevality of close multiple systems. We study the physical and chemical properties of the jets and protostellar sources in the NGC1333-IRAS4A proto-binary system using continuum emission and molecular tracers of shocked gas. We observed NGC1333-IRAS4A in the SiO(6-5), SO(6_5-5_4), and CO(2-1) lines and the continuum emission at 1.3, 1.4, and 3 mm using the IRAM Plateau de Bure Interferometer in the framework of the CALYPSO large program. We clearly disentangle for the first time the outflow emission from the two sources A1 and A2. The two protostellar jets have very different properties: the A1 jet is faster, has a short dynamical timescale (<10^3 yr), and is associated with H2 shocked emission, whereas the A2 jet, which dominates the large-scale emission, is associated with diffuse emission, bends, and emits at slower velocities. The observed bending of the A2 jet is consistent with the change of propagation direction observed at large scale and suggests jet precession on very short timescales (~200-600 yr). In addition, a chemically rich spectrum with emission from several COMs (e.g. HCOOH, CH3OCHO, CH3OCH3) is only detected towards A2. Finally, very high-velocity shocked emission (~50 km s^-1) is observed along the A1 jet. An LTE analysis shows that SiO, SO, and H2CO abundances in the gas phase are enhanced up to (3-4)x10^{-7}, (1.4-1.7)x10^{-6}, and (3-7.9)x10^{-7}, respectively. The intrinsic different properties of the jets and driving sources in NGC1333-IRAS4A suggest different evolutionary stages for the two protostars, with A1 being younger than A2, in a very early stage of star formation previous to the hot-corino phase.

preprint2015arXiv

Molecular depletion times and the CO-to-H2 conversion factor in metal-poor galaxies

Tracing molecular hydrogen content with carbon monoxide in low-metallicity galaxies has been exceedingly difficult. Here we present a new effort, with IRAM 30-m observations of 12CO(1-0) of a sample of 8 dwarf galaxies having oxygen abundances ranging from 12+logO/H=7.7 to 8.4. CO emission is detected in all galaxies, including the most metal-poor galaxy of our sample (0.1 Zsun); to our knowledge this is the largest number of 12CO(1-0) detections ever reported for galaxies with 12+logO/H<=8 (0.2 Zsun) outside the Local Group. We calculate stellar masses (Mstar) and star-formation rates (SFRs), and analyze our results by combining our observations with galaxy samples from the literature. Extending previous results for a correlation of the molecular gas depletion time, tau(dep), with Mstar and specific SFR (sSFR), we find a variation in tau(dep) of a factor of 200 or more (from <50 Myr to 10 Gyr) over a spread of 1000 in sSFR and Mstar. We exploit the variation of tau(dep) to constrain the CO-to-H2 mass conversion factor alpha(CO) at low metallicity, and assuming a power-law variation find alpha(CO) \propto (Z/Zsun)^1.9, similar to results based on dust continuum measurements compared with gas mass. By including HI measurements, we show that the fraction of total gas mass relative to the baryonic mass is higher in galaxies that are metal poor, of low mass, and of high sSFR. Finally, comparisons of the data with star-formation models of the molecular gas phases suggest that, at metallicities Z/Zsun<=0.2, there are some discrepancies with model predictions.

preprint2015arXiv

The EChO science case

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

preprint2015arXiv

The first Galaxy scale hunt for the youngest high-mass protostars

The origin of massive stars is a fundamental open issue in modern astrophysics. Pre-ALMA interferometric studies reveal precursors to early B to late O type stars with collapsing envelopes of 15-20 M$_\odot$ on 1000-3000 AU size-scales. To search for more massive envelopes we selected the most massive nearby young clumps from the ATLASGAL survey to study their protostellar content with ALMA. Our first results using the intermediate scales revealed by the ALMA ACA array providing 3-5" angular resolution, corresponding to $\sim$0.05-0.1 pc size-scales, reveals a sample of compact objects. These massive dense cores are on average two-times more massive than previous studies of similar types of objects. We expect that once the full survey is completed, it will provide a comprehensive view on the origin of the most massive stars.

preprint2015arXiv

The formation of the solar system

The solar system started to form about 4.56 Gyr ago and despite the long intervening time span, there still exist several clues about its formation. The three major sources for this information are meteorites, the present solar system structure and the planet-forming systems around young stars. In this introduction we give an overview of the current understanding of the solar system formation from all these different research fields. This includes the question of the lifetime of the solar protoplanetary disc, the different stages of planet formation, their duration, and their relative importance. We consider whether meteorite evidence and observations of protoplanetary discs point in the same direction. This will tell us whether our solar system had a typical formation history or an exceptional one. There are also many indications that the solar system formed as part of a star cluster. Here we examine the types of cluster the Sun could have formed in, especially whether its stellar density was at any stage high enough to influence the properties of today's solar system. The likelihood of identifying siblings of the Sun is discussed. Finally, the possible dynamical evolution of the solar system since its formation and its future are considered.

preprint2015arXiv

The Pillars of Creation revisited with MUSE: gas kinematics and high-mass stellar feedback traced by optical spectroscopy

Integral field unit (IFU) data of the iconic Pillars of Creation in M16 are presented. The ionisation structure of the pillars was studied in great detail over almost the entire visible wavelength range, and maps of the relevant physical parameters, e.g. extinction, electron density, electron temperature, line-of-sight velocity of the ionised and neutral gas are shown. In agreement with previous authors, we find that the pillar tips are being ionised and photo-evaporated by the massive members of the nearby cluster NGC 6611. They display a stratified ionisation structure where the emission lines peak in a descending order according to their ionisation energies. The IFU data allowed us to analyse the kinematics of the photo-evaporative flow in terms of the stratified ionisation structure, and we find that, in agreement with simulations, the photo-evaporative flow is traced by a blueshift in the position-velocity profile. The gas kinematics and ionisation structure have allowed us to produce a sketch of the 3D geometry of the Pillars, positioning the pillars with respect to the ionising cluster stars. We use a novel method to detect a previously unknown bipolar outflow at the tip of the middle pillar and suggest that it has an embedded protostar as its driving source. Furthermore we identify a candidate outflow in the leftmost pillar. With the derived physical parameters and ionic abundances, we estimate a mass loss rate due to the photo-evaporative flow of 70 M$_{\odot}$ Myr$^{-1}$ which yields an expected lifetime of approximately 3 Myr.

preprint2015arXiv

X-Shooter study of accretion in $ρ$-Ophiucus: very low-mass stars and brown dwarfs

We present new VLT/X-Shooter optical and NIR spectra of a sample of 17 candidate young low-mass stars and BDs in the rho-Ophiucus cluster. We derived SpT and Av for all the targets, and then we determined their physical parameters. All the objects but one have M*<0.6 Msun, and 8 have mass below or close to the hydrogen-burning limit. Using the intensity of various emission lines present in their spectra, we determined the Lacc and Macc for all the objects. When compared with previous works targeting the same sample, we find that, in general, these objects are not as strongly accreting as previously reported, and we suggest that the reason is our more accurate estimate of the photospheric parameters. We also compare our findings with recent works in other slightly older star-forming regions to investigate possible differences in the accretion properties, but we find that the accretion properties for our targets have the same dependence on the stellar and substellar parameters as in the other regions. This leads us to conclude that we do not find evidence for a different dependence of Macc with M* when comparing low-mass stars and BDs. Moreover, we find a similar small (1 dex) scatter in the Macc-M* relation as in some of our recent works in other star-forming regions, and no significant differences in Macc due to different ages or properties of the regions. The latter result suffers, however, from low statistics and sample selection biases in the current studies. The small scatter in the Macc-M* correlation confirms that Macc in the literature based on uncertain photospheric parameters and single accretion indicators, such as the Ha width, can lead to a scatter that is unphysically large. Our studies show that only broadband spectroscopic surveys coupled with a detailed analysis of the photospheric and accretion properties allows us to properly study the evolution of disk accretion rates.

preprint2014arXiv

Absorption Filaments Towards the Massive Clump G0.253+0.016

ALMA HCO+ observations of the infrared dark cloud G0.253+0.016 located in the Central Molecular Zone of the Galaxy are presented. The 89 GHz emission is area-filling, optically thick, and sub-thermally excited. Two types of filaments are seen in absorption against the HCO+ emission. Broad-line absorption filaments (BLAs) have widths of less than a few arcseconds (0.07 - 0.14 pc), lengths of 30 to 50 arcseconds (1.2 - 1.8 pc), and absorption profiles extending over a velocity range larger than 20 km/sec. The BLAs are nearly parallel to the nearby G0.18 non-thermal filaments and may trace HCO+ molecules gyrating about highly ordered magnetic fields located in front of G0.253+0.016 or edge-on sheets formed behind supersonic shocks propagating orthogonal to our line-of-sight in the foreground. Narrow-line absorption filaments (NLAs) have line-widths less than 20 km/sec. Some NLAs are also seen in absorption in other species with high optical depth such as HCN and occasionally in emission where the background is faint. The NLAs, which also trace low-density, sub-thermally excited HCO+ molecules, are mostly seen on the blueshifted side of the emission from G0.253+0.016. If associated with the surface of G0.253+0.016, the kinematics of the NLAs indicate that the cloud surface is expanding. The decompression of entrained, milli-Gauss magnetic fields may be responsible for the re-expansion of the surface layers of G0.253+0.016 as it recedes from the Galactic center following a close encounter with Sgr A.

preprint2014arXiv

ALMA sub-mm maser and dust distribution of VY Canis Majoris

Cool, evolved stars have copious, enriched winds. The structure of these winds and the way they are accelerated is not well known. We need to improve our understanding by studying the dynamics from the pulsating stellar surface to about 10 stellar radii, where radiation pressure on dust is fully effective. Some red supergiants have highly asymmetric nebulae, implicating additional forces. We retrieved ALMA Science Verification data providing images of sub-mm line and continuum emission from VY CMa. This enables us to locate water masers with milli-arcsec precision and resolve the dusty continuum. The 658-, 321- and 325-GHz masers lie in irregular, thick shells at increasing distances from the centre of expansion. For the first time this is confirmed as the stellar position, coinciding with a compact peak offset to the NW of the brightest continuum emission. The maser shells (and dust formation zone) overlap but avoid each other on tens-au scales. Their distribution is broadly consistent with excitation models but the conditions and kinematics appear to be complicated by wind collisions, clumping and asymmetries.

preprint2014arXiv

Brown dwarf disks with ALMA

We present ALMA continuum and spectral line data at 0.89 mm and 3.2 mm for three disks surrounding young brown dwarfs and very low mass stars in the Taurus star forming region. Dust thermal emission is detected and spatially resolved for all the three disks, while CO(J=3-2) emission is seen in two disks. We analyze the continuum visibilities and constrain the disks physical structure in dust. The results of our analysis show that the disks are relatively large, the smallest one with an outer radius of about 70 AU. The inferred disk radii, radial profiles of the dust surface density and disk to central object mass ratios lie within the ranges found for disks around more massive young stars. We derive from our observations the wavelength dependence of the millimeter dust opacity. In all the three disks data are consistent with the presence of grains with at least millimeter sizes, as also found for disks around young stars, and confirm that the early stages of the solid growth toward planetesimals occur also around very low mass objects. We discuss the implications of our findings on models of solids evolution in protoplanetary disks, on the main mechanisms proposed for the formation of brown dwarfs and very low mass stars, as well as on the potential of finding rocky and giant planets around very low mass objects.

preprint2014arXiv

CARMA Large Area Star Formation Survey: Project Overview with Analysis of Dense Gas Structure and Kinematics in Barnard 1

We present details of the CARMA Large Area Star Formation Survey (CLASSy), while focusing on observations of Barnard 1. CLASSy is a CARMA Key Project that spectrally imaged N2H+, HCO+, and HCN (J=1-0 transitions) across over 800 square arcminutes of the Perseus and Serpens Molecular Clouds. The observations have angular resolution near 7" and spectral resolution near 0.16 km/s. We imaged ~150 square arcminutes of Barnard 1, focusing on the main core, and the B1 Ridge and clumps to its southwest. N2H+ shows the strongest emission, with morphology similar to cool dust in the region, while HCO+ and HCN trace several molecular outflows from a collection of protostars in the main core. We identify a range of kinematic complexity, with N2H+ velocity dispersions ranging from ~0.05-0.50 km/s across the field. Simultaneous continuum mapping at 3 mm reveals six compact object detections, three of which are new detections. A new non-binary dendrogram algorithm is used to analyze dense gas structures in the N2H+ position-position-velocity (PPV) cube. The projected sizes of dendrogram-identified structures range from about 0.01-0.34 pc. Size-linewidth relations using those structures show that non-thermal line-of-sight velocity dispersion varies weakly with projected size, while rms variation in the centroid velocity rises steeply with projected size. Comparing these relations, we propose that all dense gas structures in Barnard 1 have comparable depths into the sky, around 0.1-0.2 pc; this suggests that over-dense, parsec-scale regions within molecular clouds are better described as flattened structures rather than spherical collections of gas. Science-ready PPV cubes for Barnard 1 molecular emission are available for download.

preprint2014arXiv

Complex organic molecules in protostellar environments in the SKA era

Molecular complexity builds up at each step of the Sun-like star formation process, starting from simple molecules and ending up in large polyatomic species. Complex organic molecules (COMs; such as methyl formate, HCOOCH$_3$, dymethyl ether, CH$_3$OCH$_3$, formamide, NH$_2$CHO, or glycoaldehyde, HCOCH$_2$OH) are formed in all the components of the star formation recipe (e.g. pre-stellar cores, hot-corinos, circumstellar disks, shocks induced by fast jets), due to ice grain mantle sublimation or sputtering as well as gas-phase reactions. Understanding in great detail the involved processes is likely the only way to predict the ultimate molecular complexity reached in the ISM, as the detection of large molecules is increasingly more difficult with the increase of the number of atoms constituting them. Thanks to the recent spectacular progress of astronomical observations, due to the Herschel (sub-mm and IR), IRAM and SMA (mm and sub-mm), and NRAO (cm) telescopes, an enormous activity is being developed in the field of Astrochemistry, extending from astronomical observatories to chemical laboratories. We are involved in several observational projects providing unbiased spectral surveys (in the 80-300 and 500-2000 GHz ranges) with unprecedented sensitivity of templates of dense cores and protostars. Forests of COM lines have been detected. In this chapter we will focus on the chemistry of both cold prestellar cores and hot shocked regions, (i) reviewing results and open questions provided by mm-FIR observations, and (ii) showing the need of carrying on the observations of COMs at lower frequencies, where SKA will operate. We will also emphasize the importance of analysing the spectra by the light of the experimental studies performed by our team, who is investigating the chemical effects induced by ionising radiation bombarding astrophysically relevant ices.

preprint2014arXiv

First results from the CALYPSO IRAM-PdBI survey. II. Resolving the hot corino in the Class 0 protostar NGC 1333-IRAS2A

We investigate the origin of complex organic molecules (COMs) in the gas phase around the low-mass Class~0 protostar NGC1333-IRAS2A, to determine if the COM emission lines trace an embedded disk, shocks from the protostellar jet, or the warm inner parts of the protostellar envelope. In the framework of the CALYPSO (Continuum And Lines in Young ProtoStellar Objects) IRAM Plateau de Bure survey, we obtained large bandwidth spectra at sub-arcsecond resolution towards NGC 1333-IRAS2A. We identify the emission lines towards the central protostar and perform Gaussian fits to constrain the size of the emitting region for each of these lines, tracing various physical conditions and scales. The emission of numerous COMs such as methanol, ethylene glycol, and methyl formate is spatially resolved by our observations. This allows us to measure, for the first time, the size of the COM emission inside the protostellar envelope, finding that it originates from a region of radius 40-100 AU, centered on the NGC 1333-IRAS2A protostellar object. Our analysis shows no preferential elongation of the COM emission along the jet axis, and therefore does not support the hypothesis that COM emission arises from shocked envelope material at the base of the jet. Down to similar sizes, the dust continuum emission is well reproduced with a single envelope model, and therefore does not favor the hypothesis that COM emission arises from the thermal sublimation of grains embedded in a circumstellar disk. Finally, the typical scale $\sim$60 AU observed for COM emission is consistent with the size of the inner envelope where $T_{\rm{dust}} > 100$ K is expected. Our data therefore strongly suggest that the COM emission traces the hot corino in IRAS2A, i.e., the warm inner envelope material where the icy mantles of dust grains evaporate because they are passively heated by the central protostellar object.

preprint2014arXiv

G0.253+0.016: A centrally condensed, high-mass protocluster

Despite their importance as stellar nurseries and the building blocks of galaxies, very little is known about the formation of the highest mass clusters. The dense clump G0.253+0.016 represents an example of a clump that may form an Arches-like, high-mass cluster. Here we present molecular line maps toward G0.253+0.016 taken as part of the MALT90 molecular line survey, complemented with APEX observations. Combined, these data reveal the global physical properties and kinematics of G0.253+0.016. Recent Herschel data show that while the dust temperature is low (~19 K) toward its centre, the dust temperature on the exterior is higher (~27 K) due to external heating. Our new molecular line data reveal that, overall, the morphology of dense gas detected toward G0.253+0.016 matches very well its IR extinction and dust continuum emission. An anti-correlation between the dust and gas column densities toward its centre indicates that the clump is centrally condensed with a cold, dense interior in which the molecular gas is chemically depleted. The velocity field shows a strong gradient along the clump's major axis, with the blue-shifted side at higher Galactic longitude. The optically thick gas tracers are systematically red-shifted with respect to the optically thin and hot gas tracers, indicating radial motions. The gas kinematics and line ratios support the recently proposed scenario in which G0.253+0.016 results from a tidal compression during a recent pericentre passage near SgrA*. Because G0.253+0.016 represents an excellent example of a clump that may form a high-mass cluster, its detailed study should reveal a wealth of knowledge about the early stages of cluster formation.

preprint2014arXiv

Grain growth in the envelopes and disks of Class I protostars

We present new 3 mm ATCA data of two Class I Young Stellar Objects in the Ophiucus star forming region: Elias29 and WL12. For our analysis we compare them with archival 1.1 mm SMA data. In the (u,v) plane the two sources present a similar behavior: a nearly constant non-zero emission at long baselines, which suggests the presence of an unresolved component and an increase of the fluxes at short baselines, related to the presence of an extended envelope. Our data analysis leads to unusually low values of the spectral index $α_{\rm 1.1-3mm}$, which may indicate that mm-sized dust grains have already formed both in the envelopes and in the disk-like structures at such early stages. To explore the possible scenarios for the interpretation of the sources we perform a radiative transfer modeling using a Monte Carlo code, in order to take into account possible deviations from the Rayleigh-Jeans and optically thin regimes. Comparison between the model outputs and the observations indicates that dust grains may form aggregates up to millimeter size already in the inner regions of the envelopes of Class I YSOs. Moreover, we conclude that the embedded disk-like structures in our two Class Is are probably very compact, in particular in the case of WL12, with outer radii down to tens of AU.

preprint2014arXiv

Millimetre spectral indices of transition disks and their relation to the cavity radius

Transition disks are protoplanetary disks with inner depleted dust cavities and excellent candidates to investigate the dust evolution under the existence of a pressure bump. A pressure bump at the outer edge of the cavity allows dust grains from the outer regions to stop their rapid inward migration towards the star and efficiently grow to millimetre sizes. Dynamical interactions with planet(s) have been one of the most exciting theories to explain the clearing of the inner disk. We look for evidence of the presence of millimetre dust particles in transition disks by measuring their spectral index with new and available photometric data. We investigate the influence of the size of the dust depleted cavity on the disk integrated millimetre spectral index. We present the 3mm photometric observations carried out with PdBI of four transition disks: LkHa330, UXTauA, LRLL31, and LRLL67. We use available values of their fluxes at 345GHz to calculate their spectral index, as well as the spectral index for a sample of twenty transition disks. We compare the observations with two kind of models. In the first set of models, we consider coagulation and fragmentation of dust in a disk in which a cavity is formed by a massive planet located at different positions. The second set of models assumes disks with truncated inner parts at different radius and with power-law dust size distributions, where the maximum size of grains is calculated considering turbulence as the source of destructive collisions. We show that the integrated spectral index is higher for transition disks than for regular protoplanetary disks. For transition disks, the probability that the measured spectral index is positively correlated with the cavity radius is 95%. High angular resolution imaging of transition disks is needed to distinguish between the dust trapping scenario and the truncated disk case.

preprint2014arXiv

Old pre-main-sequence Stars: Disc reformation by Bondi-Hoyle accretion

Young stars show evidence of accretion discs which evolve quickly and disperse with an e-folding time of $\sim$ 3Myr. This is in striking contrast with recent observations that suggest evidence for numerous $>30$ Myr old stars with an accretion disc in large star-forming complexes. We consider whether these observations of apparently old accretors could be explained by invoking Bondi-Hoyle accretion to rebuild a new disc around these stars during passage through a clumpy molecular cloud. We combine a simple Monte Carlo model to explore the capture of mass by such systems with a viscous evolution model to infer the levels of accretion that would be observed. We find that a significant fraction of stars may capture enough material via the Bondi-Hoyle mechanism to rebuild a disc of mass $\gtrsim$ 1 minimum-mass solar nebula, and $\lesssim 10\%$ accrete at observable levels at any given time. A significant fraction of the observed old accretors may be explained with our proposed mechanism. Such accretion may provide a chance for a second epoch of planet formation, and have unpredictable consequences for planetary evolution.

preprint2014arXiv

On the gas content of transitional disks: a VLT/X-Shooter study of accretion and winds

Transitional disks (TDs) are thought to be a late evolutionary stage of protoplanetary disks with dust depleted inner regions. The mechanism responsible for this depletion is still under debate. To constrain the models it is mandatory to have a good understanding of the properties of the gas content of the inner disk. Using X-Shooter broad band -UV to NIR- medium resolution spectroscopy we derive the stellar, accretion, and wind properties of a sample of 22 TDs. The analysis of these properties allows us to put strong constraints on the gas content in a region very close to the star (<0.2 AU) which is not accessible with any other observational technique. We fit the spectra with a self-consistent procedure to derive simultaneously SpT,Av,and mass accretion rates (Macc) of the targets. From forbidden emission lines we derive the wind properties of the targets. Comparing our findings to values for cTTs, we find that Macc and wind properties of 80% of the TDs in our sample, which is strongly biased towards strongly accreting objects, are comparable to those of cTTs. Thus, there are (at least) some TDs with Macc compatible with those of cTTs, irrespective of the size of the dust inner hole.Only in 2 cases Macc are much lower, while the wind properties are similar. We do not see any strong trend of Macc with the size of the dust depleted cavity, nor with the presence of a dusty optically thick disk close to the star. In the TDs in our sample there is a gas rich inner disk with density similar to that of cTTs disks. At least for some TDs, the process responsible of the inner disk clearing should allow for a transfer of gas from the outer disk to the inner region. This should proceed at a rate that does not depend on the physical mechanism producing the gap seen in the dust emission and results in a gas density in the inner disk similar to that of unperturbed disks around stars of similar mass.

preprint2014arXiv

On the structure of the transition disk around TW Hya

For over a decade, the structure of the inner cavity in the transition disk of TW Hydrae has been a subject of debate. Modeling the disk with data obtained at different wavelengths has led to a variety of proposed disk structures. Rather than being inconsistent, the individual models might point to the different faces of physical processes going on in disks, such as dust growth and planet formation. Our aim is to investigate the structure of the transition disk again and to find to what extent we can reconcile apparent model differences. A large set of high-angular-resolution data was collected from near-infrared to centimeter wavelengths. We investigated the existing disk models and established a new self-consistent radiative-transfer model. A genetic fitting algorithm was used to automatize the parameter fitting. Simple disk models with a vertical inner rim and a radially homogeneous dust composition from small to large grains cannot reproduce the combined data set. Two modifications are applied to this simple disk model: (1) the inner rim is smoothed by exponentially decreasing the surface density in the inner ~3 AU, and (2) the largest grains (>100 um) are concentrated towards the inner disk region. Both properties can be linked to fundamental processes that determine the evolution of protoplanetary disks: the shaping by a possible companion and the different regimes of dust-grain growth, respectively. The full interferometric data set from near-infrared to centimeter wavelengths requires a revision of existing models for the TW Hya disk. We present a new model that incorporates the characteristic structures of previous models but deviates in two key aspects: it does not have a sharp edge at 4 AU, and the surface density of large grains differs from that of smaller grains. This is the first successful radiative-transfer-based model for a full set of interferometric data.

preprint2014arXiv

The ATLASGAL survey: a catalog of dust condensations in the Galactic plane

The formation processes and the evolutionary stages of high-mass stars are poorly understood compared to low-mass stars. Large-scale surveys are needed to provide an unbiased census of high column density sites which can potentially host precursors to high-mass stars. Here we use the ATLASGAL survey covering 420 sq. degree of the Galactic plane at 870 $μ$m; and use the MRE-GLC method to identify the population of embedded sources throughout the inner Galaxy. We identify in total 10861 compact sub-millimeter sources with fluxes above 5 sigma. Completeness tests show that our catalogue is 97% complete above 5 sigma and >99% complete above 7$σ$. We correlate this sample with mid-infrared point source catalogues (MSX at 21.3 $μ$m and WISE at 22 $μ$m) and determine a lower limit of ~33% that are associated with embedded protostellar objects. We note that the proportion of clumps associated with mid-infrared sources increases with increasing flux density, achieving a rather constant fraction of ~75% of all clumps with fluxes over 5 Jy/beam being associated with star-formation. Examining the source counts as a function of Galactic longitude we are able to identify the most prominent star forming regions in the Galaxy. From the fraction of the likely massive quiescent clumps (~25%) we estimate a formation time-scale of ~7.5+/-2.5 $\times$ 10$^4$yr for the deeply embedded phase before the emergence of luminous YSOs. Such a short duration for the formation of high-mass stars in massive clumps clearly proves that the earliest phases have to be dynamic with supersonic motions.

preprint2014arXiv

Turbulence sets the initial conditions for star formation in high-pressure environments

Despite the simplicity of theoretical models of supersonically turbulent, isothermal media, their predictions successfully match the observed gas structure and star formation activity within low-pressure (P/k < 10^5 K cm^-3) molecular clouds in the solar neighbourhood. However, it is unknown if these theories extend to clouds in high-pressure (P/k > 10^7 K cm^-3) environments, like those in the Galaxy's inner 200 pc Central Molecular Zone (CMZ) and in the early Universe. Here we present ALMA 3mm dust continuum emission within a cloud, G0.253+0.016, which is immersed in the high-pressure environment of the CMZ. While the log-normal shape and dispersion of its column density PDF is strikingly similar to those of solar neighbourhood clouds, there is one important quantitative difference: its mean column density is 1--2 orders of magnitude higher. Both the similarity and difference in the PDF compared to those derived from solar neighbourhood clouds match predictions of turbulent cloud models given the high-pressure environment of the CMZ. The PDF shows a small deviation from log-normal at high column densities confirming the youth of G0.253+0.016. Its lack of star formation is consistent with the theoretically predicted, environmentally dependent volume density threshold for star formation which is orders of magnitude higher than that derived for solar neighbourhood clouds. Our results provide the first empirical evidence that the current theoretical understanding of molecular cloud structure derived from the solar neighbourhood also holds in high-pressure environments. We therefore suggest that these theories may be applicable to understand star formation in the early Universe.

preprint2014arXiv

X-Shooter spectroscopy of young stellar objects: V - Slow winds in T Tauri stars

Disks around T Tauri stars are known to lose mass, as best shown by the profiles of forbidden emission lines of low ionization species. At least two separate kinematic components have been identified, one characterised by velocity shifts of tens to hundreds km/s (HVC) and one with much lower velocity of few km/s (LVC). The HVC are convincingly associated to the emission of jets, but the origin of the LVC is still unknown. In this paper we analyze the forbidden line spectrum of a sample of 44 mostly low mass young stars in Lupus and $σ$-Ori observed with the X-Shooter ESO spectrometer. We detect forbidden line emission of [OI], [OII], [SII], [NI], and [NII], and characterize the line profiles as LVC, blue-shifted HVC and red-shifted HVC. We focus our study on the LVC. We show that there is a good correlation between line luminosity and both L$_{star}$ and the accretion luminosity (or the mass-accretion rate) over a large interval of values (L$_{star}$ $\sim 10^{-2} - 1$ L$_\odot$; L$_{acc}$ $\sim 10^{-5} - 10^{-1}$ L$_\odot$; $\dot M_{acc}$ $\sim 10^{-11} - 10^{-7}$ M$_\odot$/yr). The lines show the presence of a slow wind ($V_{peak}<20$ km/s), dense ($n_H>10^8$ cm$^{-3}$), warm (T$\sim 5000-10000$ K), mostly neutral. We estimate the mass of the emitting gas and provide a value for the maximum volume it occupies. Both quantities increase steeply with the stellar mass, from $\sim 10^{-12}$ M$_\odot$ and $\sim 0.01$ AU$^3$ for M$_{star}$$\sim 0.1$ M$_\odot$, to $\sim 3 \times 10^{-10}$ M$_\odot$ and $\sim 1$ AU$^3$ for M$_{star}$$\sim 1$ M$_\odot$, respectively. These results provide quite stringent constraints to wind models in low mass young stars, that need to be explored further.

preprint2013arXiv

A Community Science Case for E-ELT HIRES

Building on the experience of the high-resolution community with the suite of VLT high-resolution spectrographs, which has been tremendously successful, we outline here the (science) case for a high-fidelity, high-resolution spectrograph with wide wavelength coverage at the E-ELT. Flagship science drivers include: the study of exo-planetary atmospheres with the prospect of the detection of signatures of life on rocky planets; the chemical composition of planetary debris on the surface of white dwarfs; the spectroscopic study of protoplanetary and proto-stellar disks; the extension of Galactic archaeology to the Local Group and beyond; spectroscopic studies of the evolution of galaxies with samples that, unlike now, are no longer restricted to strongly star forming and/or very massive galaxies; the unraveling of the complex roles of stellar and AGN feedback; the study of the chemical signatures imprinted by population III stars on the IGM during the epoch of reionization; the exciting possibility of paradigm-changing contributions to fundamental physics. The requirements of these science cases can be met by a stable instrument with a spectral resolution of R~100,000 and broad, simultaneous spectral coverage extending from 370nm to 2500nm. Most science cases do not require spatially resolved information, and can be pursued in seeing-limited mode, although some of them would benefit by the E-ELT diffraction limited resolution. Some multiplexing would also be beneficial for some of the science cases. (Abridged)

preprint2013arXiv

A Hi-GAL study of the high-mass star-forming region G29.96-0.02

Context. G29.96-0.02 is a high-mass star-forming cloud observed at 70, 160, 250, 350, and 500 microns as part of the Herschel survey of the Galactic Plane during the Science Demonstration Phase. Aims. We wish to conduct a far-infrared study of the sources associated with this star-forming region by estimating their physical properties and evolutionary stage, and investigating the clump mass function, the star formation efficiency and rate in the cloud. Methods. We have identified the Hi-GAL sources associated with the cloud, searched for possible counterparts at centimeter and infrared wavelengths, fitted their spectral energy distribution and estimated their physical parameters. Results. A total of 198 sources have been detected in all 5 Hi-GAL bands, 117 of which are associated with 24 microns emission and 87 of which are not associated with 24 microns emission. We called the former sources 24 microns-bright and the latter ones 24 microns-dark. The [70-160] color of the 24 microns-dark sources is smaller than that of the 24 microns-bright ones. The 24 microns-dark sources have lower L_bol and L_bol/M_env than the 24 microns-bright ones for similar M_env, which suggests that they are in an earlier evolutionary phase. The G29-SFR cloud is associated with 10 NVSS sources and with extended centimeter continuum emission well correlated with the 70 microns emission. Most of the NVSS sources appear to be early B or late O-type stars. The most massive and luminous Hi-GAL sources in the cloud are located close to the G29-UC region, which suggests that there is a privileged area for massive star formation towards the center of the G29-SFR cloud. Almost all the Hi-GAL sources have masses well above the Jeans mass but only 5% have masses above the virial mass, which indicates that most of the sources are stable against gravitational collapse. The sources with M_env > M_virial and that ...

preprint2013arXiv

Accurate determination of accretion and photospheric parameters in Young Stellar Objects: the case of two candidate old disks in the Orion Nebula Cluster

Current planet formation models are largely based on the observational constraint that protoplanetary disks have lifetime 3Myr. Recent studies, however, report the existence of PMS stars with signatures of accretion (strictly connected with the presence of circumstellar disks)and photometrically determined ages of 30 Myr, or more. Here we present a spectroscopic study of two major age outliers in the ONC. We use broad band, intermediate resolution VLT/X-Shooter spectra combined with an accurate method to determine the stellar parameters and the related age of the targets to confirm their peculiar age estimates and the presence of ongoing accretion.The analysis is based on a multi-component fitting technique, which derives simultaneously SpT, extinction, and accretion properties of the objects. With this method we confirm and quantify the ongoing accretion. From the photospheric parameters of the stars we derive their position on the HRD, and the age given by evolutionary models. Together with other age indicators like the lithium equivalent width we estimate with high accuracy the age of the objects. Our study shows that the two objects analyzed are not older than the typical population of the ONC. Our results show that, while photometric determination of the photospheric parameters are an accurate method to estimate the parameters of the bulk of young stellar populations, those of individual objects with high accretion rates and extinction may be affected by large uncertainties. Broad band spectroscopic determinations should thus be used to confirm the nature of individual objects. Our analysis shows that this method allows us to obtain an accurate determination of the photospheric parameters of accreting YSOs in any nearby star-forming region. We suggest that our detailed, broad- band spectroscopy method should be used to derive accurate properties of candidate old and accreting YSOs.

preprint2013arXiv

ALMA and CARMA observations of Brown Dwarfs disks: testing the models of dust evolution

The first steps toward planet formation involve the coagulation of small microscopic grains into larger and larger pebbles and rocks in gas-rich disks around young stars and brown dwarfs. Observations in the sub-millimeter can trace mm/cm-sized pebbles in the outer disks, and investigate the mechanisms of coagulation/fragmentation and radial migration of these solids. These represent key, yet not fully understood ingredients for our understanding of the formation of planetesimals, the building blocks of planets. Here we present the first results from an observational program using the ALMA and CARMA sub-mm/mm interferometers aimed at characterizing the dust properties and disk structure of young disks around brown dwarfs and very low mass stars. Given the physical conditions expected for these disks, they represent critical test beds for the models of the early stages of planet formation in proto-planetary disks.

preprint2013arXiv

ALMA detection of the rotating molecular disk wind from the young star HD 163296

Disk winds have been postulated as a mechanism for angular momentum release in protostellar systems for decades. HD 163296 is a Herbig Ae star surrounded by a disk and has been shown to host a series of HH knots (HH 409) with bow shocks associated with the farthest knots. Here we present ALMA Science Verification data of CO J=2-1 and J=3-2 emission which are spatially coincident with the blue shifted jet of HH knots, and offset from the disk by -18.6 km/s. The emission has a double corkscrew morphology and extends more than 10" from the disk with embedded emission clumps coincident with jet knots. We interpret this double corkscrew as emission from material in a molecular disk wind, and that the compact emission near the jet knots is being heated by the jet which is moving at much higher velocities. We show that the J=3-2 emission is likely heavily filtered by the interferometer, but the J=2-1 emission suffers less due to the larger beam and measurable angular scales. Excitation analysis suggests temperatures exceeding 900 K in these compact features, with the wind mass, momentum and energy being of order 10^{-5} M_sun, 10^{-4} M_sun km/s and 10^{40} erg respectively. The high mass loss rate suggests that this star is dispersing the disk faster than it is funneling mass onto the star.

preprint2013arXiv

ALMA imaging of the CO snowline of the HD 163296 disk with DCO+

The high spatial and line sensitivity of ALMA opens the possibility of resolving emission from molecules in circumstellar disks. With an understanding of physical conditions under which molecules have high abundance, they can be used as direct tracers of distinct physical regions. In particular, DCO+ is expected to have an enhanced abundance within a few Kelvin of the CO freezeout temperature of 19 K, making it a useful probe of the cold disk midplane. We compare ALMA line observations of HD 163296 to a grid of models. We vary the upper- and lower-limit temperatures of the region in which DCO+ is present as well as the abundance of DCO+ in order to fit channel maps of the DCO+ J=5-4 line. To determine the abundance enhancement compared to the general interstellar medium, we carry out similar fitting to HCO+ J=4-3 and H13CO+ J=4-3 observations. ALMA images show centrally peaked extended emission from HCO+ and H13CO+. DCO+ emission lies in a resolved ring from ~110 to 160 AU. The outer radius approximately corresponds to the size of the CO snowline as measured by previous lower resolution observations of CO lines in this disk. The ALMA DCO+ data now resolve and image the CO snowline directly. In the best fitting models, HCO+ exists in a region extending from the 19 K isotherm to the photodissociation layer with an abundance of 3x10^-10 relative to H2. DCO+ exists within the 19-21 K region of the disk with an abundance ratio [DCO+] / [HCO+] = 0.3. This represents a factor of 10^4 enhancement of the DCO+ abundance within this narrow region of the HD 163296 disk. Such a high enhancement has only previously been seen in prestellar cores. The inferred abundances provide a lower limit to the ionization fraction in the midplane of the cold outer disk (approximately greater than 4x10^-10), and suggest the utility of DCO+ as a tracer of its parent molecule H2D+. Abridged

preprint2013arXiv

ALMA observations of cool dust in a low-metallicity starburst, SBS0335-052

We present Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 0 Band 7 observations of an extremely metal-poor dwarf starburst galaxy in the Local Universe, SBS0335-052 (12+log(O/H)~7.2). With these observations, dust is detected at 870micron (ALMA Band 7), but 87% of the flux in this band is due to free-free emission from the starburst. We have compiled a spectral energy distribution (SED) of SBS0335-052 that spans almost 6 orders of magnitude in wavelength and fit it with a spherical dust shell heated by a single-age stellar population; the best-fit model gives a dust mass of (3.8+/-0.6)x10^4 Msun. We have also constructed a SED including Herschel archival data for IZw18, another low-metallicity dwarf starburst (12+log(O/H)=7.17), and fit it with a similar model to obtain a dust mass of (3.4+/-1.0)x10^2 Msun. Compared with their atomic gas mass, the dust mass of SBS0335-052 far exceeds the prediction of a linear trend of dust-to-gas mass ratio with metallicity, while IZw18 falls far below. We use gas scaling relations to assess a putative missing gas component in both galaxies and find that the missing, possibly molecular, gas in SBS0335-052 is a factor of 6 times higher than the value inferred from the observed HI column density; in IZw18 the missing component is 4 times smaller. Ultimately, despite their similarly low metallicity, the differences in gas and dust column densities in SBS0335-052 and IZw18 suggest that metal abundance does not uniquely define star-formation processes. At some level, self-shielding and the survival of molecules may depend just as much on gas and dust column density as on metallicity. The effects of low metallicity may at least be partially compensated for by large column densities in the interstellar medium.

preprint2013arXiv

Candidate super star cluster progenitor gas clouds possibly triggered by close passage to Sgr A*

Super star clusters are the end product of star formation under the most extreme conditions. As such, studying how their final stellar populations are assembled from their natal progenitor gas clouds can provide strong constraints on star formation theories. An obvious place to look for the initial conditions of such extreme stellar clusters are gas clouds of comparable mass and density, with no star formation activity. We present a method to identify such progenitor gas clouds and demonstrate the technique for the gas in the inner few hundred pc of our Galaxy. The method highlights three clouds in the region with similar global physical properties to the previously identified extreme cloud, G0.253+0.016, as potential young massive cluster (YMC) precursors. The fact that four potential YMC progenitor clouds have been identified in the inner 100 pc of the Galaxy, but no clouds with similar properties have been found in the whole first quadrant despite extensive observational efforts, has implications for cluster formation/destruction rates across the Galaxy. We put forward a scenario to explain how such dense gas clouds can arise in the Galactic centre environment, in which YMC formation is triggered by gas streams passing close to the minimum of the global Galactic gravitational potential at the location of the central supermassive black hole, Sgr A*. If this triggering mechanism can be verified, we can use the known time interval since closest approach to Sgr A* to study the physics of stellar mass assembly in an extreme environment as a function of absolute time.

preprint2013arXiv

CARMA interferometric observations of 2MASS J044427+2512: the first spatially resolved observations of thermal emission of a brown dwarf disk

We present CARMA 1.3 mm continuum data of the disk surrounding the young brown dwarf 2MASS J044427+2512 in the Taurus molecular cloud. The high angular resolution of the CARMA observations (0.16 arcsec) allows us to spatially resolve for the first time the thermal emission from dust around a brown dwarf. We analyze the interferometric visibilities and constrain the disk outer radius adopting disk models with power-law radial profiles of the dust surface density. In the case of a power-law index equal to or lower than 1, we obtain a disk radius in the range of about 15 - 30 AU, while larger disks are inferred for steeper radial profiles. By combining this information on the disk spatial extent with the sub-mm spectral index of this source we find conclusive evidence for mm-sized grains, or larger, in this brown dwarf disk. We discuss the implications of our results on the models of dust evolution in proto-planetary disks and brown dwarf formation.

preprint2013arXiv

Constraints on the radial distribution of the dust properties in the CQ Tau protoplanetary disk

Grain growth in protoplanetary disks is the first step towards the formation of the rocky cores of planets. Models predict that grains grow, migrate, and fragment in the disk and predict varying dust properties as a function of radius, age, and physical properties. High-angular resolution observations at more than one (sub-)mm wavelength are the essential tool for constraining grain growth and migration on the disk midplane. We developed a procedure to analyze self-consistently multi wavelength (sub-)mm continuum interferometric observations of protoplanetary disks to constrain the radial distribution of dust properties. We apply this technique to existing multi frequency continuum mm observations of the disk around CQ Tau, a A8 pre-main sequence star with a well-studied disk. We demonstrate that our models can be used to simultaneously constrain the disk and dust structure. In CQ Tau, the best-fitting model has a radial dependence of the maximum grain size, which decreases from a few cm in the inner disk (<40 AU) to a few mm at 80 AU. Nevertheless, the currently available dataset does not allow us to exclude the possibility of a uniform grain size distribution at a 3sigma level.

preprint2013arXiv

Evidence of a SiO collimated outflow from a massive YSO in IRAS 17233-3606

Studies of molecular outflows in high-mass young stellar objects reveal important information about the formation process of massive stars. We therefore selected the close-by IRAS 17233-3606 massive star-forming region to perform SiO observations with the SMA interferometer in the (5-4) line and with the APEX single-dish telescope in the (5-4) and (8-7) transitions. In this paper, we present a study of one of the outflows in the region, OF1, which shows several properties similar to jets driven by low-mass protostars, such as HH211 and HH212. It is compact and collimated, and associated with extremely high velocity CO emission, and SiO emission at high velocities. We used a state-of-the-art shock model to constrain the pre-shock density and shock velocity of OF1. The model also allowed us to self-consistently estimate the mass of the OF1 outflow. The shock parameters inferred by the SiO modelling are comparable with those found for low-mass protostars, only with higher pre-shock density values, yielding an outflow mass in agreement with those obtained for molecular outflows driven by early B-type young stellar objects. Our study shows that it is possible to model the SiO emission in high-mass star-forming regions in the same way as for shocks from low-mass young stellar objects.

preprint2013arXiv

Explaining millimeter-sized particles in brown dwarf disks

Planets have been detected around a variety of stars, including low-mass objects, such as brown dwarfs. However, such extreme cases are challenging for planet formation models. Recent sub-millimeter observations of disks around brown dwarf measured low spectral indices of the continuum emission that suggest that dust grains grow to mm-sizes even in these very low mass environments. To understand the first steps of planet formation in scaled-down versions of T-Tauri disks, we investigate the physical conditions that can theoretically explain the growth from interstellar dust to millimeter-sized grains in disks around brown dwarf. We modeled the evolution of dust particles under conditions of low-mass disks around brown dwarfs. We used coagulation, fragmentation and disk-structure models to simulate the evolution of dust, with zero and non-zero radial drift. For the non-zero radial drift, we considered strong inhomogeneities in the gas surface density profile that mimic long-lived pressure bumps in the disk. We studied different scenarios that could lead to an agreement between theoretical models and the spectral slope found by millimeter observations. We find that fragmentation is less likely and rapid inward drift is more significant for particles in brown dwarf disks than in T-Tauri disks. We present different scenarios that can nevertheless explain millimeter-sized grains. As an example, a model that combines the following parameters can fit the millimeter fluxes measured for brown dwarf disks: strong pressure inhomogeneities of $\sim$ 40% of amplitude, a small radial extent $\sim$ 15 AU, a moderate turbulence strength $α_{\mathrm{turb}}= 10^{-3}$, and average fragmentation velocities for ices $v_f = 10 m s^{-1}$.

preprint2013arXiv

First results from the Herschel Gould Belt Survey in Taurus

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

preprint2013arXiv

Herschel Reveals Massive Cold Clumps in NGC 7538

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

preprint2013arXiv

On the shape of the mass-function of dense clumps in the Hi-GAL fields. II. Using Bayesian inference to study the clump mass function

Context. Stars form in dense, dusty clumps of molecular clouds, but little is known about their origin, their evolution and their detailed physical properties. In particular, the relationship between the mass distribution of these clumps (also known as the "clump mass function", or CMF) and the stellar initial mass function (IMF), is still poorly understood. Aims. In order to better understand how the CMF evolve toward the IMF, and to discern the "true" shape of the CMF, large samples of bona-fide pre- and proto-stellar clumps are required. Two such datasets obtained from the Herschel infrared GALactic Plane Survey (Hi-GAL) have been described in paper I. Robust statistical methods are needed in order to infer the parameters describing the models used to fit the CMF, and to compare the competing models themselves. Methods. In this paper we apply Bayesian inference to the analysis of the CMF of the two regions discussed in Paper I. First, we determine the Bayesian posterior probability distribution for each of the fitted parameters. Then, we carry out a quantitative comparison of the models used to fit the CMF. Results. We have compared the results from several methods implementing Bayesian inference, and we have also analyzed the impact of the choice of priors and the influence of various constraints on the statistical conclusions for the preferred values of the parameters. We find that both parameter estimation and model comparison depend on the choice of parameter priors. Conclusions. Our results confirm our earlier conclusion that the CMFs of the two Hi-GAL regions studied here have very similar shapes but different mass scales. Furthermore, the lognormal model appears to better describe the CMF measured in the two Hi-GAL regions studied here. However, this preliminary conclusion is dependent on the choice of parameters priors.

preprint2013arXiv

Silicon isotopic abundance toward evolved stars and its application for presolar grains

Galactic chemical evolution (GCE) is important for understanding the composition of the present-day interstellar medium (ISM) and of our solar system. In this paper, we aim to track the GCE by using the 29Si/30Si ratios in evolved stars and tentatively relate this to presolar grain composition. We used the APEX telescope to detect thermal SiO isotopologue emission toward four oxygen-rich M-type stars. Together with the data retrieved from the Herschel science archive and from the literature, we were able to obtain the 29Si/30Si ratios for a total of 15 evolved stars inferred from their optically thin 29SiO and 30SiO emission. These stars cover a range of masses and ages, and because they do not significantly alter 29Si/30Si during their lifetimes, they provide excellent probes of the ISM metallicity (or 29Si/30Si ratio) as a function of time. The 29Si/30Si ratios inferred from the thermal SiO emission tend to be lower toward low-mass oxygen-rich stars (e.g., down to about unity for W Hya), and close to an interstellar or solar value of 1.5 for the higher-mass carbon star IRC+10216 and two red supergiants. There is a tentative correlation between the 29Si/30Si ratios and the mass-loss rates of evolved stars, where we take the mass-loss rate as a proxy for the initial stellar mass or current stellar age. This is consistent with the different abundance ratios found in presolar grains. We found that older objects (up to possibly 10 Gyr old) in our sample trace a previous, lower 29Si/30Si value of about 1. Material with this isotopic ratio is present in two subclasses of presolar grains, providing independent evidence of the lower ratio. Therefore, the 29Si/30Si ratio derived from the SiO emission of evolved stars is a useful diagnostic tool for the study of the GCE and presolar grains.

preprint2013arXiv

Strongly star-forming rotating disks in a complex merging system at z = 4,7 as revealed by ALMA

We performed a kinematical analysis of the [CII] line emission of the BR 1202-0725 system at z~4,7 using ALMA observations. The most prominent sources of this system are a quasar and a submillimeter galaxy, separated by a projected distance of about 24 kpc and characterized by very high SFR, higher than 1000 Msun/yr. However, the ALMA observations reveal that these galaxies apparently have undisturbed rotating disks, which is at variance with the commonly accepted scenario in which strong star formation activity is induced by a major merger. We also detected faint components which, after spectral deblending, were spatially resolved from the main QSO and SMG emissions. The relative velocities and positions of these components are compatible with orbital motions within the gravitational potentials generated by the QSO host galaxy and the SMG, suggesting that they are smaller galaxies in interaction or gas clouds in accretion flows of tidal streams. We did not find any clear spectral evidence for outflows caused by AGN or stellar feedback. This suggests that the high star formation rates might be induced by interactions or minor mergers with these companions, which do not affect the large-scale kinematics of the disks, however. Our kinematical analysis also indicates that the QSO and the SMG have similar Mdyn, mostly in the form of molecular gas, and that the QSO host galaxy and the SMG are seen close to face-on with slightly different disk inclinations: the QSO host galaxy is seen almost face-on (i~15), while the SMG is seen at higher inclinations (i~25). Finally, the ratio between the black hole mass of the QSO, obtained from XShooter spectroscopy, and the Mdyn of the host galaxy is similar to value found in very massive local galaxies, suggesting that the evolution of black hole galaxy relations is probably better studied with dynamical than with stellar host galaxy masses.

preprint2013arXiv

The Herschel and JCMT Gould Belt Surveys: Constraining Dust Properties in the Perseus B1 Clump with PACS, SPIRE, and SCUBA-2

We present Herschel observations from the Herschel Gould Belt Survey and SCUBA-2 science verification observations from the JCMT Gould Belt Survey of the B1 clump in the Perseus molecular cloud. We determined the dust emissivity index using four different techniques to combine the Herschel PACS+SPIRE data at 160 - 500 microns with the SCUBA-2 data at 450 microns and 850 microns. Of our four techniques, we found the most robust method was to filter-out the large-scale emission in the Herschel bands to match the spatial scales recovered by the SCUBA-2 reduction pipeline. Using this method, we find beta ~ 2 towards the filament region and moderately dense material and lower beta values (beta > 1.6) towards the dense protostellar cores, possibly due to dust grain growth. We find that beta and temperature are more robust with the inclusion of the SCUBA-2 data, improving estimates from Herschel data alone by factors of ~ 2 for beta and by ~ 40% for temperature. Furthermore, we find core mass differences of < 30% compared to Herschel-only estimates with an adopted beta = 2, highlighting the necessity of long wavelength submillimeter data for deriving accurate masses of prestellar and protostellar cores.

preprint2013arXiv

The Herschel Gould Belt Survey in Chamaeleon II - Properties of cold dust in disks around young stellar objects

We report on the Herschel Gould Belt survey (HGBS) of ChaII, focusing on the detection of Class I to III young stellar objects (YSOs). We aim at characterizing the circumstellar material around these YSOs and understanding which disk parameters are most likely constrained by the new HGBS data. We recovered 29 out of the 63 known YSOs in ChaII with a detection in at least one of the PACS/SPIRE pass-bands: 3 Class I YSOs (i.e.,100%), 1 Flat source (i.e., 50%), 21 Class II objects (i.e., 55%), 3 Class III objects (i.e, 16%) and the unclassified source IRAS 12522-7640. We explore PACS/SPIRE colors of this sample and present modeling of their SEDs using the RADMC-2D radiative transfer code. We find that YSO colors are confined in specific regions of PACS/SPIRE color-color diagrams. These color ranges are expected to be only marginally contaminated by extragalactic sources and field stars and, hence, provide a useful YSO selection tool when applied altogether. We were able to model the SED of 26 out of the 29 detected YSOs. We discuss the degeneracy/limitations of our SED fitting results and adopt the Bayesian method to estimate the probability of different values for the derived disk parameters. The Cha II YSOs present typical disk inner radii around 0.1 AU. The lower limit to Rc is typically around 50 AU. The lower limits to Mdisk are proportional to the stellar masses with a typical 0.3% ratio, i.e., in the range estimated in the literature for young Class II stars and brown dwarfs across a broad range of stellar masses. The estimated flaring angles, although very uncertain, point towards rather flat disks (1+phi less than 1.2), as found for low-mass M-type YSO samples in other star forming regions. Thus, our results support the idea that disk properties show a dependence on stellar properties.

preprint2013arXiv

Unveiling the gas and dust disk structure in HD 163296 using ALMA observations

Aims: The aim of this work is to study the structure of the protoplanetary disk surrounding the Herbig Ae star HD 163296. Methods: We have used high-resolution and high-sensitivity ALMA observations of the CO(3-2) emission line and the continuum at 850 microns, as well as the 3- dimensional radiative transfer code MCFOST to model the data presented in this work. Results: The CO(3-2) emission unveils for the first time at sub-millimeter frequencies the vertical structure details of a gaseous disk in Keplerian rotation, showing the back- and the front-side of a flared disk. Continuum emission at 850 microns reveals a compact dust disk with a 240 AU outer radius and a surface brightness profile that shows a very steep decline at radius larger than 125 AU. The gaseous disk is more than two times larger than the dust disk, with a similar critical radius but with a shallower radial profile. Radiative transfer models of the continuum data confirms the need for a sharp outer edge to the dust disk. The models for the CO(3-2) channel map require the disk to be slightly more geometrically thick than previous models suggested, and that the temperature at which CO gas becomes depleted (frozen-out) from the outer regions of the disk midplane is T < 20 K, in agreement with previous studies.

preprint2013arXiv

Water vapor in the protoplanetary disk of DG Tau

Water is key in the evolution of protoplanetary disks and the formation of comets and icy/water planets. While high excitation water lines originating in the hot inner disk have been detected in several T Tauri stars (TTSs), water vapor from the outer disk, where most of water ice reservoir is stored, was only reported in the closeby TTS TW Hya. We present spectrally resolved Herschel/HIFI observations of the young TTS DG Tau in the ortho- and para- water ground-state transitions at 557, 1113 GHz. The lines show a narrow double-peaked profile, consistent with an origin in the outer disk, and are ~19-26 times brighter than in TW Hya. In contrast, CO and [C II] lines are dominated by emission from the envelope/outflow, which makes H2O lines a unique tracer of the disk of DG Tau. Disk modeling with the thermo-chemical code ProDiMo indicates that the strong UV field, due to the young age and strong accretion of DG Tau, irradiates a disk upper layer at 10-90 AU from the star, heating it up to temperatures of 600 K and producing the observed bright water lines. The models suggest a disk mass of 0.015-0.1 Msun, consistent with the estimated minimum mass of the solar nebula before planet formation, and a water reservoir of ~1e2-1e3 Earth oceans in vapour, and ~100 times larger in the form of ice. Hence, this detection supports the scenario of ocean delivery on terrestrial planets by impact of icy bodies forming in the outer disk.

preprint2013arXiv

X-Shooter spectroscopy of young stellar objects III. Photospheric and chromospheric properties of Class III objects

We analyzed X-Shooter/VLT spectra of 24 ClassIII sources from three nearby star-forming regions (sigmaOrionis, LupusIII, and TWHya). We determined the effective temperature, surface gravity, rotational velocity, and radial velocity by comparing the observed spectra with synthetic BT-Settl model spectra. We investigated in detail the emission lines emerging from the stellar chromospheres and combined these data with archival X-ray data to allow for a comparison between chromospheric and coronal emissions. Both X-ray and Halpha luminosity as measured in terms of the bolometric luminosity are independent of the effective temperature for early-M stars but decline toward the end of the spectral M sequence. For the saturated early-M stars the average emission level is almost one dex higher for X-rays than for Halpha: log(L_x/L_bol) = -2.85 +- 0.36 vs. log(L_Halpha/L_bol) = -3.72 +- 0.21. When all chromospheric emission lines (including the Balmer series up to H11, CaII HK, the CaII infrared triplet, and several HeI lines) are summed up the coronal flux still dominates that of the chromosphere, typically by a factor 2-5. Flux-flux relations between activity diagnostics that probe different atmospheric layers (from the lower chromosphere to the corona) separate our sample of active pre-main sequence stars from the bulk of field M dwarfs studied in the literature. Flux ratios between individual optical emission lines show a smooth dependence on the effective temperature. The Balmer decrements can roughly be reproduced by an NLTE radiative transfer model devised for another young star of similar age. Future, more complete chromospheric model grids can be tested against this data set.

preprint2013arXiv

X-Shooter spectroscopy of young stellar objects: II. Impact of chromospheric emission on accretion rate estimates

Context. The lack of knowledge of photospheric parameters and the level of chromospheric activity in young low-mass pre-main sequence stars introduces uncertainties when measuring mass accretion rates in accreting (Class II) Young Stellar Objects. A detailed investigation of the effect of chromospheric emission on the estimates of mass accretion rate in young low-mass stars is still missing. This can be undertaken using samples of young diskless (Class III) K and M-type stars. Aims. Our goal is to measure the chromospheric activity of Class III pre main sequence stars to determine its effect on the estimates of accretion luminosity (Lacc) and mass accretion rate (Macc) in young stellar objects with disks. Methods. Using VLT/X-Shooter spectra we have analyzed a sample of 24 non-accreting young stellar objects of spectral type between K5 and M9.5. We identify the main emission lines normally used as tracers of accretion in Class II objects, and we determine their fluxes in order to estimate the contribution of the chromospheric activity to the line luminosity. Results. We have used the relationships between line luminosity and accretion luminosity derived in the literature for Class II objects to evaluate the impact of chromospheric activity on the accretion rate measurements. We find that the typical chromospheric activity would bias the derived accretion luminosity by Lacc,noise< 10-3Lsun, with a strong dependence with the Teff of the objects. The noise on Macc depends on stellar mass and age, and the typical values of log(Macc,noise) range between -9.2 to -11.6Msun/yr. Conclusions. Values of Lacc< 10-3Lsun obtained in accreting low-mass pre main sequence stars through line luminosity should be treated with caution as the line emission may be dominated by the contribution of chromospheric activity.

preprint2013arXiv

X-Shooter spectroscopy of young stellar objects: IV -- Accretion in low-mass stars and sub-stellar objects in Lupus

We present X-Shooter/VLT observations of a sample of 36 accreting low-mass stellar and sub-stellar objects (YSOs) in the Lupus star forming region, spanning a range in mass from ~0.03 to ~1.2Msun, but mostly with 0.1Msun < Mstar < 0.5Msun. Our aim is twofold: firstly, analyse the relationship between excess-continuum and line emission accretion diagnostics, and, secondly, to investigate the accretion properties in terms of the physical properties of the central object. The accretion luminosity (Lacc), and from it the accretion rate (Macc), is derived by modelling the excess emission, from the UV to the near-IR, as the continuum emission of a slab of hydrogen. The flux and luminosity (Ll) of a large number of emission lines of H, He, CaII, etc., observed simultaneously in the range from ~330nm to 2500nm, were computed. The luminosity of all the lines is well correlated with Lacc. We provide empirical relationships between Lacc and the luminosity of 39 emission lines, which have a lower dispersion as compared to previous relationships in the literature. Our measurements extend the Pab and Brg relationships to Lacc values about two orders of magnitude lower than those reported in previous studies. We confirm that different methodologies to measure Lacc and Macc yield significantly different results: Ha line profile modelling may underestimate Macc by 0.6 to 0.8dex with respect to Macc derived from continuum-excess measures. Such differences may explain the likely spurious bi-modal relationships between Macc and other YSOs properties reported in the literature. We derive Macc in the range 2e-12 -- 4e-8 Msun/yr and conclude that Macc is proportional to Mstar^1.8(+/-0.2), with a dispersion lower by a factor of about 2 than in previous studies. A number of properties indicate that the physical conditions of the accreting gas are similar over more than 5 orders of magnitude in Macc.

preprint2012arXiv

ALMA observations of rho-Oph 102: grain growth and molecular gas in the disk around a young Brown Dwarf

We present ALMA continuum and spectral line observations of the young Brown Dwarf rho-Oph 102 at about 0.89 mm and 3.2 mm. We detect dust emission from the disk at these wavelengths and derive an upper limit on the radius of the dusty disk of ~ 40 AU. The derived variation of the dust opacity with frequency in the mm provides evidence for the presence of mm-sized grains in the disk outer regions. This result demonstrates that mm-grains are found even in the low density environments of Brown Dwarf disks and challenges our current understanding of dust evolution in disks. The CO map at 345 GHz clearly reveals molecular gas emission at the location of the Brown Dwarf, indicating a gas-rich disk as typically found for disks surrounding young pre-Main Sequence stars. We derive a disk mass of ~ 0.3-1% of the mass of the central Brown Dwarf, similar to the typical values found for disks around more massive young stars.

preprint2012arXiv

An HST Imaging Survey of Low-Mass Stars in the Chamaeleon I Star Forming region

We present new HST/WFPC2 observations of 20 fields centered around T Tauri stars in the Chamaeleon I star forming region. Images have been obtained in the F631N ([OI]6300A), F656N (Ha) and F673N ([SII]6716A+6731A) narrow-band filters, plus the Johnson V-band equivalent F547M filter. We detect 31 T Tauri stars falling within our fields. We discuss the optical morphology of 10 sources showing evidence of either binarity, circumstellar material, or mass loss. We supplement our photometry with a compilation of optical, infrared and sub-millimeter data from the literature, together with new sub-mm data for three objects, to build the Spectral Energy Distributions (SED) of 19 single sources. Using an SED model fitting tool, we self-consistently estimate a number of stellar and disk parameters, while mass accretion rates are directly derived from our Ha photometry. We find that bolometric luminosities derived from dereddened optical data tend to be underestimated in systems with high alpha(2-24} IR spectral index, suggesting that disks seen nearly edge-on may occasionally be interpreted as low luminosity (and therefore more evolved) sources. On the other hand, the same alpha(2-24) spectral index, a tracer of the amount of dust in the warmer layers of the circumstellar disks, and the mass accretion rate appear to decay with the isocronal stellar age, suggesting that the observed age spread (~0.5-5 Myr) within the cluster is real. Our sample contains a few outliers that may have dissipated their circumstellar disks on shorter time-scale.

preprint2012arXiv

Can habitable planets form in clustered environments?

We present observational evidence of environmental effects on the formation and evolution of planetary systems. We combine catalogues of resolved protoplanetary discs (PPDs) and young stellar objects in the solar neighbourhood to analyse the PPD size distribution as a function of ambient stellar density. By running Kolmogorov-Smirnov tests between the PPD radii at different densities, we find empirical evidence, at the >97% confidence level, for a change in the PPD radius distribution at ambient stellar densities Σ> 10^{3.5} pc^{-2}. This coincides with a simple theoretical estimate for the truncation of PPDs or planetary systems by dynamical encounters. If this agreement is causal, the ongoing disruption of PPDs and planetary systems limits the possible existence of planets in the habitable zone, with shorter lifetimes at higher host stellar masses and ambient densities. Therefore, habitable planets are not likely to be present in long-lived stellar clusters, and may have been ejected altogether to form a population of unbound, free-floating planets. We conclude that, while highly suggestive, our results should be verified through other methods. Our simple model shows that truncations should lead to a measurable depletion of the PPD mass function that can be detected with ALMA observations of the densest nearby and young clusters.

preprint2012arXiv

Cluster-formation in the Rosette molecular cloud at the junctions of filaments

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

preprint2012arXiv

Different evolutionary stages in massive star formation. Centimeter continuum and H2O maser emission with ATCA

We present ATCA observations of the H2O maser line and radio continuum at 18.0GHz and 22.8GHz, toward a sample of 192 massive star forming regions containing several clumps already imaged at 1.2mm. The main aim of this study is to investigate the water maser and centimeter continuum emission (likely tracing thermal free-free emission) in sources at different evolutionary stages, using the evolutionary classifications proposed by Palla et al (1991) and Molinari et al (2008). We used the recently comissioned CABB backend at ATCA obtaining images with 20arcsec resolution in the 1.3cm continuum and H2O maser emission, in all targets. For the evolutionary analysis of the sources we used the millimeter continuum emission from Beltran et al (2006) and the infrared emission from the MSX Point Source Catalogue. We detect centimeter continuum emission in 88% of the observed fields with a typical rms noise level of 0.45mJy/beam. Most of the fields show a single radio continuum source, while in 20% of them we identify multiple components. A total of 214 centimeter continuum sources have been identified, likely tracing optically thin HII regions, with physical parameters typical of both extended and compact HII regions. Water maser emission was detected in 41% of the regions, resulting in a total of 85 distinct components. The low angular (20arcsec) and spectral (14km/s) resolutions do not allow a proper analysis of the water maser emission, but suffice to investigate its association with the continuum sources. We have also studied the detection rate of HII regions in the two types of IRAS sources defined by Palla et (1991) on the basis of the IRAS colours: High and Low. No significant differences are found, with large detection rates (>90%) for both High and Low sources. We classify the millimeter and infrared sources in our fields in three evolutionary stages following the scheme presented by ...

preprint2012arXiv

Fomalhaut debris disk emission at 7 millimeters: constraints on the collisional models of planetesimals

We present new spatially resolved observations of the dust thermal emission at 7 mm from the Fomalhaut debris disk obtained with the Australia Telescope Compact Array. These observations provide the longest wavelength detection of the Fomalhaut debris disk to date. We combined the new data to literature sub-mm data to investigate the spectral index of the dust thermal emission in the sub-millimeter and constrained the $q$-slope of the power-law grain size distribution. We derived a value for $q = 3.48 \pm 0.14$ for grains with sizes around 1 mm. This is consistent with the classical prediction for a collisional cascade at the steady-state. The same value cannot be explained by more recent collisional models of planetesimals in which either the velocity distribution of the large bodies or their tensile strength is a strong function of the body size.

preprint2012arXiv

Herschel Observations of a Potential Core Forming Clump: Perseus B1-E

We present continuum observations of the Perseus B1-E region from the Herschel Gould Belt Survey. These Herschel data reveal a loose grouping of substructures at 160 - 500 micron not seen in previous submillimetre observations. We measure temperature and column density from these data and select the nine densest and coolest substructures for follow-up spectral line observations with the Green Bank Telescope. We find that the B1-E clump has a mass of ~ 100 solar masses and appears to be gravitationally bound. Furthermore, of the nine substructures examined here, one substructure (B1-E2) appears to be itself bound. The substructures are typically less than a Jeans length from their nearest neighbour and thus, may interact on a timescale of ~ 1 Myr. We propose that B1-E may be forming a first generation of dense cores, which could provide important constraints on the initial conditions of prestellar core formation. Our results suggest that B1-E may be influenced by a strong, localized magnetic field, but further observations are still required.

preprint2012arXiv

Herschel observations of B1-bS and B1-bN: two first hydrostatic core candidates in the Perseus star-forming cloud

We report far-IR Herschel observations obtained between 70 $μ$m and 500 $μ$m of two star-forming dusty condensations, B1-bS and B1-bN, in the B1 region of the Perseus star-forming cloud. In the Western part of the Perseus cloud, B1-bS is the only source detected in all of the 6 PACS and SPIRE photometric bands without being visible in the Spitzer map at 24 $μ$m. B1-bN is clearly detected between 100 $μ$m and 250 $μ$m. We have fitted the spectral energy distributions of these sources to derive their physical properties, and find that a simple greybody model fails to reproduce the observed SEDs. At least a two-component model, consisting of a central source surrounded by a dusty envelope is required. The properties derived from the fit, however, suggest that the central source is not a Class 0 object. We then conclude that while B1-bS and B1-bN appear to be more evolved than a pre-stellar core, the best-fit models suggest that their central objects are younger than a Class 0 source. Hence, they may be good candidates to be examples of the first hydrostatic core phase. The projected distance between B1-bS and B1-bN is a few Jeans lengths. If their physical separation is close to this value, this pair would allow the mutual interactions between two forming stars at a very early stage of their evolution to be studied.

preprint2012arXiv

On the shape of the mass-function of dense clumps in the Hi-GAL fields. I. SED determination and global properties of the mass-functions

Stars form in dense, dusty clumps of molecular clouds, but little is known about their origin and evolution. In particular, the relationship between the mass distribution of these clumps (also known as the "clump mass function", or CMF) and the stellar initial mass function (IMF), is still poorly understood. In order to discern the "true" shape of the CMF and to better understand how the CMF may evolve toward the IMF, large samples of bona-fide pre- and proto-stellar clumps are required. The sensitive observations of the Herschel Space Observatory (HSO) are now allowing us to look at large clump populations in various clouds with different physical conditions. We analyse two fields in the Galactic plane mapped by HSO during its science demonstration phase, as part of the more complete and unbiased Herschel infrared GALactic Plane Survey (Hi-GAL). These fields undergo a source-extraction and flux-estimation pipeline, which allows us to obtain a sample with thousands of clumps. Starless and proto-stellar clumps are separated using both color and positional criteria to find those coincident with MIPS 24 micron sources. We describe the probability density functions of the power-law and lognormal models that are used to fit the CMFs, and we then find their best-fit parameters. For the lognormal model we apply several statistical techniques to the data and compare their results. The CMFs of the two SDP fields show very similar shapes, but very different mass scales. This similarity is confirmed by the values of the best-fit parameters of either the power-law or lognormal model. The power-law model leads to almost identical CMF slopes, whereas the lognormal model shows that the CMFs have similar widths. The similar CMF shape but different mass scale represents an evidence that the overall process of star formation in the two regions is very different.

preprint2012arXiv

Recent star formation in the Lupus clouds as seen by Herschel

We present a study of the star formation histories of the Lupus I, III, and IV clouds using the Herschel 70-500 micron maps obtained by the Herschel Gould Belt Survey Key Project. By combining the new Herschel data with the existing Spitzer catalog we obtained an unprecedented census of prestellar sources and young stellar objects in the Lupus clouds, which allowed us to study the overall star formation rate (SFR) and efficiency (SFE). The high SFE of Lupus III, its decreasing SFR, and its large number of pre-main sequence stars with respect to proto- and prestellar sources, suggest that Lupus III is the most evolved cloud, and after having experienced a major star formation event in the past, is now approaching the end of its current star-forming cycle. Lupus I is currently undergoing a large star formation event, apparent by the increasing SFR, the large number of prestellar objects with respect to more evolved objects, and the high percentage of material at high extinction (e.g., above A_V=8 mag). Also Lupus IV has an increasing SFR; however, the relative number of prestellar sources is much lower, suggesting that its star formation has not yet reached its peak.

preprint2012arXiv

Search for starless clumps in the ATLASGAL survey

In this study, we present an unbiased sample of the earliest stages of massive star formation across 20 square-degree of the sky. Within the region 10deg < l < 20deg and |b| < 1deg, we search the ATLASGAL survey at 870 micron for dense gas condensations. These clumps are carefully examined for indications of ongoing star formation using YSOs from the GLIMPSE source catalog as well as sources in the 24 micron MIPSGAL images, to search for starless clumps. We calculate the column densities as well as the kinematic distances and masses for sources where the v_lsr is known from spectroscopic observations. Within the given region, we identify 210 starless clumps with peak column densities > 1 x 10e23 cm^(-2). In particular, we identify potential starless clumps on the other side of the Galaxy. The sizes of the clumps range between 0.1 pc and 3 pc with masses between a few tens of solar masses up to several ten thousands of solar masses. Most of them may form massive stars, but in the 20 square-degree we only find 14 regions massive enough to form stars more massive than 20 solar masses and 3 regions with the potential to form stars more massive than 40 40 solar masses. The slope of the high-mass tail of the clump mass function for clumps on the near side of the Galaxy is 2.2 and, therefore, Salpeter-like. We estimate the lifetime of the most massive starless clumps to be 60000 yr. The sample offers a uniform selection of starless clumps. In the large area surveyed, we only find a few potential precursors of stars in the excess of 40 solar masses. It appears that the lifetime of these clumps is somewhat shorter than their free-fall times, although both values agree within the errors. In addition, these are ideal objects for detailed studies and follow-up observations.

preprint2012arXiv

The effect of local optically thick regions in the long-wave emission of young circumstellar disks

Multi-wavelength observations of protoplanetary disks in the sub-millimeter continuum have measured spectral indices values which are significantly lower than what is found in the diffuse interstellar medium. Under the assumption that mm-wave emission of disks is mostly optically thin, these data have been generally interpreted as evidence for the presence of mm/cm-sized pebbles in the disk outer regions. In this work we investigate the effect of possible local optically thick regions on the mm-wave emission of protoplanetary disks without mm/cm-sized grains. A significant local increase of the optical depth in the disk can be caused by the concentration of solid particles, as predicted to result from a variety of proposed physical mechanisms. We calculate the filling factors and implied overdensities these optically thick regions would need to significantly affect the millimeter fluxes of disks, and we discuss their plausibility. We find that optically thick regions characterized by relatively small filling factors can reproduce the mm-data of young disks without requesting emission from mm/cm-sized pebbles. However, these optically thick regions require dust overdensities much larger than what predicted by any of the physical processes proposed in the literature to drive the concentration of solids. We find that only for the most massive disks it is possible and plausible to imagine that the presence of optically thick regions in the disk is responsible for the low measured values of the mm spectral index. For the majority of the disk population, optically thin emission from a population of large mm-sized grains remains the most plausible explanation. The results of this analysis further strengthen the scenario for which the measured low spectral indices of protoplanetary disks at mm wavelengths are due to the presence of large mm/cm-sized pebbles in the disk outer regions.

preprint2012arXiv

The Herschel view of the on-going star formation in the Vela-C molecular cloud

As part of the Herschel guaranteed time key program 'HOBYS', we present the photometric survey of the star forming region Vela-C, one of the nearest sites of low-to-high-mass star formation in the Galactic plane. Vela-C has been observed with PACS and SPIRE in parallel mode between 70 um and 500 um over an area of about 3 square degrees. A photometric catalogue has been extracted from the detections in each band, using a threshold of 5 sigma over the local background. Out of this catalogue we have selected a robust sub-sample of 268 sources, of which 75% are cloud clumps and 25% are cores. Their Spectral Energy Distributions (SEDs) have been fitted with a modified black body function. We classify 48 sources as protostellar and 218 as starless. For two further sources, we do not provide a secure classification, but suggest they are Class 0 protostars. From SED fitting we have derived key physical parameters. Protostellar sources are in general warmer and more compact than starless sources. Both these evidences can be ascribed to the presence of an internal source(s) of moderate heating, which also causes a temperature gradient and hence a more peaked intensity distribution. Moreover, the reduced dimensions of protostellar sources may indicate that they will not fragment further. A virial analysis of the starless sources gives an upper limit of 90% for the sources gravitationally bound and therefore prestellar. We fit a power law N(logM) prop M^-1.1 to the linear portion of the mass distribution of prestellar sources. This is in between that typical of CO clumps and those of cores in nearby star-forming regions. We interpret this as a result of the inhomogeneity of our sample, which is composed of comparable fractions of clumps and cores.

preprint2012arXiv

The M16 molecular complex under the influence of NGC6611. Herschel's perspective of the heating effect on the Eagle Nebula

We present Herschel images from the HOBYS key program of the Eagle Nebula (M16) in the far-infrared and sub-millimetre, using the PACS and SPIRE cameras at 70μm, 160μm, 250μm, 350μm, 500μm. M16, home to the Pillars of Creation, is largely under the influence of the nearby NGC6611 high-mass star cluster. The Herschel images reveal a clear dust temperature gradient running away from the centre of the cavity carved by the OB cluster. We investigate the heating effect of NGC6611 on the entire M16 star-forming complex seen by Herschel including the diffuse cloud environment and the dense filamentary structures identified in this region. In addition, we interpret the three-dimensional geometry of M16 with respect to the nebula, its surrounding environment, and the NGC6611 cavity. The dust temperature and column density maps reveal a prominent eastern filament running north-south and away from the high-mass star-forming central region and the NGC6611 cluster, as well as a northern filament which extends around and away from the cluster. The dust temperature in each of these filaments decreases with increasing distance from the NGC6611 cluster, indicating a heating penetration depth of \sim 10 pc in each direction in 3 - 6 \times 10^{22} cm-2 column density filaments. We show that in high-mass star-forming regions OB clusters impact the temperature of future star-forming sites, modifying the initial conditions for collapse and effecting the evolutionary criteria of protostars developed from spectral energy distributions. Possible scenarios for the origin of the morphology seen in this region are discussed, including a western equivalent to the eastern filament, which was destroyed by the creation of the OB cluster and its subsequent winds and radiation.

preprint2012arXiv

The Pipe Nebula as seen with Herschel: Formation of filamentary structures by large-scale compression ?

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

preprint2012arXiv

Variations in the Galactic star formation rate and density thresholds for star formation

The conversion of gas into stars is a fundamental process in astrophysics and cosmology. Stars are known to form from the gravitational collapse of dense clumps in interstellar molecular clouds, and it has been proposed that the resulting star formation rate is proportional to either the amount of mass above a threshold gas surface density, or the gas volume density. These star-formation prescriptions appear to hold in nearby molecular clouds in our Milky Way Galaxy's disk as well as in distant galaxies where the star formation rates are often much larger. The inner 500 pc of our Galaxy, the Central Molecular Zone (CMZ), contains the largest concentration of dense, high-surface density molecular gas in the Milky Way, providing an environment where the validity of star-formation prescriptions can be tested. Here we show that by several measures, the current star formation rate in the CMZ is an order-of-magnitude lower than the rates predicted by the currently accepted prescriptions. In particular, the region 1 deg < l < 3.5 deg, |b| < 0.5 deg contains ~10^7 Msun of dense molecular gas -- enough to form 1000 Orion-like clusters -- but the present-day star formation rate within this gas is only equivalent to that in Orion. In addition to density, another property of molecular clouds, such as the amplitude of turbulent motions, must be included in the star-formation prescription to predict the star formation rate in a given mass of molecular gas.

preprint2011arXiv

A 100-parsec elliptical and twisted ring of cold and dense molecular clouds revealed by Herschel around the Galactic Center

Thermal images of cold dust in the Central Molecular Zone of the Milky Way, obtained with the far-infrared cameras on-board the Herschel satellite, reveal a 3x10^7 solar masses ring of dense and cold clouds orbiting the Galactic Center. Using a simple toy-model, an elliptical shape having semi-major axes of 100 and 60 parsecs is deduced. The major axis of this 100-pc ring is inclined by about 40 degrees with respect to the plane-of-the-sky and is oriented perpendicular to the major axes of the Galactic Bar. The 100-pc ring appears to trace the system of stable x_2 orbits predicted for the barred Galactic potential. Sgr A* is displaced with respect to the geometrical center of symmetry of the ring. The ring is twisted and its morphology suggests a flattening-ratio of 2 for the Galactic potential, which is in good agreement with the bulge flattening ratio derived from the 2MASS data.

preprint2011arXiv

A low optical depth region in the inner disk of the HerbigAe star HR5999

Circumstellar disks surrounding young stars are known to be the birthplaces of planets, and the innermost astronomical unit is of particular interest. We present new long-baseline spectro-interferometric observations of the HerbigAe star, HR5999, obtained in the H and K bands with the AMBER instrument at the VLTI, and aim to produce near-infrared images at the sub-AU spatial scale. We spatially resolve the circumstellar material and reconstruct images using the MiRA algorithm. In addition, we interpret the interferometric observations using models that assume that the near-infrared excess is dominated by the emission of a circumstellar disk. We compare the images reconstructed from the VLTI measurements to images obtained using simulated model data. The K-band image reveals three main elements: a ring-like feature located at ~0.65 AU, a low surface brightness region inside, and a central spot. At the maximum angular resolution of our observations (1.3 mas), the ring is resolved while the central spot is only marginally resolved, preventing us from revealing the exact morphology of the circumstellar environment. We suggest that the ring traces silicate condensation, i.e., an opacity change, in a circumstellar disk around HR 5999. We build a model that includes a ring at the silicate sublimation radius and an inner disk of low surface brightness responsible for a large amount of the near-infrared continuum emission. The model successfully fits the SED, visibilities, and closure phases, and provides evidence of a low surface brightness region inside the silicate sublimation radius. This study provides additional evidence that in HerbigAe stars, there is material in a low surface brightness region, probably a low optical depth region, located inside the silicate sublimation radius and of unknown nature.

preprint2011arXiv

Characterizing Precursors to Stellar Clusters with Herschel

Context: Despite their profound effect on the universe, the formation of massive stars and stellar clusters remains elusive. In the past decade, compelling evidence has emerged that suggests Infrared Dark Clouds (IRDCs) may be precursors to stellar clusters. However, the usual method for identifying IRDCs is biased by the requirement that they are seen in absorption against bright mid-IR emission, whereas dust continuum observations allow cold, dense pre-stellar-clusters to be identified anywhere. Aims: We aim to understand what physical properties characterize IRDCs, to explore the population of dust continuum sources that are not IRDCs, and to roughly characterize the star formation activity in dust continuum sources. Results: We present temperature and column density maps in the Hi-GAL l=30 and l=59 SDP fields, as well as a robust algorithm for cirrus subtraction and source identification using Hi-GAL. We report on the fraction of Hi-GAL sources which are mid-IR-dark, -neutral, or -bright in both fields. We find significant trends in column density and temperature between mid-IR-dark and -bright pixels; mid-IR-dark pixels are about 10 K colder and have a factor of 2 higher column density on average than mid-IR-bright pixels. We find that Hi-GAL dust continuum sources span a range of evolutionary states from pre- to star-forming, and that warmer sources are associated with more star formation tracers. There is a trend of increasing temperature with tracer type from mid-IR-dark at the coldest, to outflow/maser sources in the middle, and finally to 8 and 24 micron bright sources at the warmest. Finally, we identify five candidate IRDC-like sources on the far-side of the Galaxy. These are cold (~ 20 K), high column density (N(H2) > 10^22 cm^-2) clouds identified with Hi-GAL which, despite bright surrounding mid-IR emission, show little to no absorption at 8 micron. (abridged)

preprint2011arXiv

Clustering Properties of Far-Infrared Sources in Hi-GAL Science Demonstration Phase Fields

We use a Minimum Spanning Tree algorithm to characterize the spatial distribution of Galactic Far-IR sources and derive their clustering properties. We aim to reveal the spatial imprint of different types of star forming processes, e.g. isolated spontaneous fragmentation of dense molecular clouds, or events of triggered star formation around HII regions, and highlight global properties of star formation in the Galaxy. We plan to exploit the entire Hi-GAL survey of the inner Galactic plane to gather significant statistics on the clustering properties of star forming regions, and to look for possible correlations with source properties such as mass, temperature or evolutionary stage. In this paper we present a pilot study based on the two 2x2 square degree fields centered at longitudes l=30 and l=59 obtained during the Science Demonstration Phase (SDP) of the Herschel mission. We find that over half of the clustered sources are associated with HII regions and infrared dark clouds. Our analysis also reveals a smooth chromatic evolution of the spatial distribution where sources detected at short-wavelengths, likely proto-stars surrounded by warm circumstellar material emitting in the far-infrared, tend to be clustered in dense and compact groups around HII regions while sources detected at long-wavelengths, presumably cold and dusty density enhancements of the ISM emitting in the sub-millimeter, are distributed in larger and looser groups.

preprint2011arXiv

Looking for outflow and infall signatures in high mass star forming regions

(Context) Many physical parameters change with time in star forming regions. Here we attempt to correlate changes in infall and outflow motions in high mass star forming regions with evolutionary stage using JCMT observations. (Aims) From a sample of 45 high mass star forming regions in three phases of evolution, we investigate the presence of established infall and outflow tracers to determine whether there are any trends attributable to the age of the source. (Methods) We obtained JCMT observations of HCO+/H13CO+ J=4-3 to trace large scale infall, and SiO J=8-7 to trace recent outflow activity. We compare the infall and outflow detections to the evolutionary stage of the host source (high mass protostellar objects, hypercompact HII regions and ultracompact HII regions). We also note that the integrated intensity of SiO varies with the full width at half maximum of the H13CO+. (Results) We find a surprising lack of SiO detections in the middle stage (Hypercompact HII regions), which may be due to an observational bias. When SiO is detected, we find that the integrated intensity of the line increases with evolutionary stage. We also note that all of the sources with infall signatures onto Ultracompact HII regions have corresponding outflow signatures as well.

preprint2011arXiv

The Herschel view of massive star formation in G035.39--00.33: Dense and cold filament of W48 undergoing a mini-starburst

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

preprint2011arXiv

The mm-colors of a young binary disk system in the Orion Nebula Cluster

We present new EVLA continuum observations at 7 mm of the 253-1536 binary disk system in the Orion Nebula Cluster. The measured fluxes were combined with data in the sub-mm to derive the millimeter spectral index of each individual disk component. We show how these observations can be used to test the models of dust evolution and early growth of solids in protoplanetary disks. Our analysis indicates that the disk with lower density and higher temperature hosts larger grains than the companion disk. This result is the opposite of what predicted by the dust evolution models. The models and observational results can be reconciled if the viscosity $α$-parameter differs by more than a factor of ten in the two disks, or if the distribution of solids in the disks is strongly affected by radial motions. This analysis can be applied to future high-angular resolution observations of young disks with EVLA and ALMA to provide even stronger observational constraints to the models of dust evolution in protoplanetary disks.

preprint2011arXiv

The molecular distribution of the IRDC G351.77-0.51

Infrared dark clouds are massive, dense clouds seen in extinction against the IR Galactic background. Many of these objects appear to be on the verge of star and star cluster formation. Our aim is to understand the physical properties of IRDCs in very early evolutionary phases. We selected the filamentary IRDC G351.77 - 0.51, which is remarkably IR quiet at 8μm. As a first step, we observed mm dust continuum emission and rotational lines of moderate and dense gas tracers to characterise different condensations along the IRDC and study the velocity field of the filament. Our initial study confirms coherent velocity distribution along the infrared dark cloud ruling out any coincidental projection effects. Excellent correlation between MIR extinction, mm continuum emission and gas distribution is found. Large-scale turbulence and line profiles throughout the filament is indicative of a shock in this cloud. Excellent correlation between line width, and MIR brightness indicates turbulence driven by local star formation.

preprint2011arXiv

Trapping dust particles in the outer regions of protoplanetary disks

In order to explain grain growth to mm sized particles and their retention in outer regions of protoplanetary disks, as it is observed at sub-mm and mm wavelengths, we investigate if strong inhomogeneities in the gas density profiles can slow down excessive radial drift and can help dust particles to grow. We use coagulation/fragmentation and disk-structure models, to simulate the evolution of dust in a bumpy surface density profile which we mimic with a sinusoidal disturbance. For different values of the amplitude and length scale of the bumps, we investigate the ability of this model to produce and retain large particles on million years time scales. In addition, we introduced a comparison between the pressure inhomogeneities considered in this work and the pressure profiles that come from magnetorotational instability. Using the Common Astronomy Software Applications ALMA simulator, we study if there are observational signatures of these pressure inhomogeneities that can be seen with ALMA. We present the favorable conditions to trap dust particles and the corresponding calculations predicting the spectral slope in the mm-wavelength range, to compare with current observations. Finally we present simulated images using different antenna configurations of ALMA at different frequencies, to show that the ring structures will be detectable at the distances of the Taurus Auriga or Ophiucus star forming regions.

preprint2010arXiv

A gallery of bubbles - The nature of the bubbles observed by Spitzer and what ATLASGAL tells us about the surrounding neutral material

We attempt to determine the nature of the bubbles observed by Spitzer in the Galactic plane, mainly to establish if possible their association with massive stars. We take advantage of the very simple morphology of these objects to search for star formation triggered by HII regions, and to estimate the importance of this mode of star formation. We consider a sample of 102 bubbles detected by Spitzer-GLIMPSE, and catalogued by Churchwell et al.(2006). We use mid-infrared and radio-continuum public data to discuss their nature. We use the ATLASGAL survey at 870 micron to search for dense neutral material collected on their borders. Results: We find that 86% of the bubbles contain ionized gas detected by means of its radio-continuum emission at 20-cm. Thus, most of the bubbles observed at 8.0 micron enclose HII regions ionized by O-B2 stars. Ninety-eight percent of the bubbles exhibit 24 micron emission in their central regions. The ionized regions at the center of the 8.0 micron bubbles seem to be devoid of PAHs but contain hot dust. Among the 65 regions for which the angular resolution of the observations is high enough to resolve the spatial distribution of cold dust at 870 micron, we find that 40% are surrounded by cold dust, and that another 28% contain interacting condensations. The former are good candidates for the collect and collapse process, as they display an accumulation of dense material at their borders. The latter are good candidates for the compression of pre-existing condensations by the ionized gas. Eighteen bubbles exhibit associated ultracompact HII regions and/or methanol masers in the direction of dust condensations adjacent to their ionization fronts. Our results suggest that more than a quarter of the bubbles may have triggered the formation of massive objects.

preprint2010arXiv

A Herschel study of the properties of starless cores in the Polaris Flare dark cloud region using PACS and SPIRE

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

preprint2010arXiv

A Herschel study of YSO evolutionary stages and formation timelines in two fields of the Hi-GAL survey

We present a first study of the star-forming compact dust condensations revealed by Herschel in the two 2 \times 2 \degr Galactic Plane fields centered at [l;b] = [30\degr; 0 \degr] and [l;b] = [59\degr; 0 \degr], respectively, and observed during the Science Demonstration Phase for the Herschel infrared Galactic Plane survey (Hi-GAL) Key-Project. Compact source catalogs extracted for the two fields in the five Hi-GAL bands (70, 160, 250, 350 and 500 $μ$m) were merged based on simple criteria of positional association and spectral energy distribution (SED) consistency into a final catalog which contains only coherent SEDs with counterparts in at least three adjacent Herschel bands. These final source lists contain 528 entries for the l = 30\degr field, and 444 entries for the l = 59\degr field. The SED coverage has been augmented with ancillary data at 24 $μ$m and 1.1 mm. SED modeling for the subset of 318 and 101 sources (in the two fields, respectively) for which the distance is known was carried out using both a structured star/disk/envelope radiative transfer model and a simple isothermal grey-body. Global parameters like mass, luminosity, temperature and dust properties have been estimated. The Lbol/Menv ratio spans four orders of magnitudes from values compatible with the pre-protostellar phase to embedded massive zero-age main sequence stars. Sources in the l = 59\degr field have on average lower L/M, possibly outlining an overall earlier evolutionary stage with respect to the sources in the l = 30\degr field. Many of these cores are actively forming high-mass stars, although the estimated core surface densities appear to be an order of magnitude below the 1 g cm$^{-2}$ critical threshold for high-mass star formation.

preprint2010arXiv

Clouds, filaments and protostars: the Herschel Hi-GAL Milky Way

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° x 2° tiles approximately centered at l=30° and l=59°. 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° field; a A_V value between 5 and 10 is found for the l=30° 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.

preprint2010arXiv

Dust grain growth in rho-Ophiuchi protoplanetary disks

We present new ATCA observations at 3.3 mm of 27 young stellar objects in the rho-Oph young cluster. 25 of these sources have been detected. We analyze the sub-millimeter and millimeter SED for a subsample of 17 isolated class II protoplanetary disks and derive constraints on the grain growth and total dust mass in the disk outer regions. All the disks in our sample show a mm slope of the SED which is significantly shallower than the one observed for the ISM at these long wavelengths. This indicates that 1) class II disks in Ophiuchus host grains grown to mm/cm-sizes in their outer regions, 2) formation of mm/cm-sized pebbles is a fast process and 3) a mechanism halting or slowing down the inward radial drift of solid particles is required to explain the data. These findings are consistent with previous results in other star forming regions. We compare the dust properties of this sample with those of a uniformly selected sample in Taurus-Auriga and find no statistical evidence of any difference in terms of grain growth between the two regions. Finally, in our sample the mm slope of the SED is not found to correlate with indicators of grain growth to micron sizes in the surface layers of the inner disk.

preprint2010arXiv

Filamentary structures and compact objects in the Aquila and Polaris clouds observed by Herschel

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

preprint2010arXiv

From filamentary clouds to prestellar cores to the stellar IMF: Initial highlights from the Herschel Gould Belt survey

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

preprint2010arXiv

Herschel observations of embedded protostellar clusters in the Rosette Molecular Cloud

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

preprint2010arXiv

Herschel Observations of the W43 "mini-starburst"

Aims: To explore the infrared and radio properties of one of the closest Galactic starburst regions. Methods: Images obtained with the Herschel Space Observatory at wavelengths of 70, 160, 250, 350, and 500 microns using the PACS and SPIRE arrays are analyzed and compared with radio continuum VLA data and 8 micron images from the Spitzer Space Telescope. The morphology of the far-infrared emission is combined with radial velocity measurements of millimeter and centimeter wavelength transitions to identify features likely to be associated with the W43 complex. Results: The W43 star-forming complex is resolved into a dense cluster of protostars, infrared dark clouds, and ridges of warm dust heated by massive stars. The 4 brightest compact sources with L > 1.5 x 10^4 Lsun embedded within the Z-shaped ridge of bright dust emission in W43 remain single at 4" (0.1 pc) resolution. These objects, likely to be massive protostars or compact clusters in early stages of evolution are embedded in clumps with masses of 10^3 to 10^4 Msun, but contribute only 2% to the 3.6 x 10^6 Lsun far-IR luminosity of W43 measured in a 16 by 16 pc box. The total mass of gas derived from the far-IR dust emission inside this region is ~10^6 Msun. Cometary dust clouds, compact 6 cm radio sources, and warm dust mark the locations of older populations of massive stars. Energy release has created a cavity blowing-out below the Galactic plane. Compression of molecular gas in the plane by the older HII region near G30.684-0.260 and the bipolar structure of the resulting younger W43 HII region may have triggered the current mini-star burst.

preprint2010arXiv

New constraints on dust grain size and distribution in CQ Tau

Grain growth in circumstellar disks is expected to be the first step towards the formation of planetary systems. There is now evidence for grain growth in several disks around young stars. Radially resolved images of grain growth in circumstellar disks are believed to be a powerful tool to constrain the dust evolution models and the initial stage for the formation of planets. In this paper we attempt to provide these constraints for the disk surrounding the young star CQ Tau. This system was already suggested from previous studies to host a population of grains grown to large sizes. We present new high angular resolution (0.3-0.9 arcsec) observations at wavelengths from 850um to 3.6cm obtained at the SMA, IRAM-PdBI and NRAO-VLA interferometers. We perform a combined analysis of the spectral energy distribution and of the high-resolution images at different wavelengths using a model to describe the dust thermal emission from the circumstellar disk. We include a prescription for the gas emission from the inner regions of the system. We detect the presence of evolved dust by constraining the disk averaged dust opacity coefficient beta (computed between 1.3 and 7mm) to be 0.6+/-0.1. This confirms the earlier suggestions that the disk contains dust grains grown to significant sizes and puts this on firmer grounds by tightly constraining the gas contamination to the observed fluxes at mm-cm wavelengths. We report some evidence of radial variations in dust properties, but current resolution and sensitivity are still too low for definitive results.

preprint2010arXiv

Resolved images of self-gravitating circumstellar discs with ALMA

In this paper we present simulated observations of massive self-gravitating circumstellar discs using the Atacama Large Millimetre/sub-millimetre Array (ALMA). Using a smoothed particle hydrodynamics model of a $0.2M_{\odot}$ disc orbiting a $1M_{\odot}$ protostar, with a cooling model appropriate for discs at temperatures below $\sim 160$K and representative dust opacities, we have constructed maps of the expected emission at sub-mm wavelengths. We have then used the CASA ALMA simulator to generate simulated images and visibilities with various array configurations and observation frequencies, taking into account the expected thermal noise and atmospheric opacities. We find that at 345 GHz (870 $μ$m) spiral structures at a resolution of a few AU should be readily detectable in approximately face-on discs out to distances of the Taurus-Auriga star-forming complex.

preprint2010arXiv

Results from DROXO. III. Observation, source list and X-ray properties of sources detected in the "Deep Rho Ophiuchi XMM-Newton Observation"

X-rays from very young stars are powerful probes to investigate the mechanisms at work in the very first stages of the star formation and the origin of X-ray emission in very young stars. We present results from a 500 ks long observation of the Rho Ophiuchi cloud with a XMM-Newton large program named DROXO, aiming at studying the X-ray emission of deeply embedded Young Stellar Objects (YSOs). The data acquired during the DROXO program were reduced with SAS software, and filtered in time and energy to improve the signal to noise of detected sources; light curves and spectra were obtained. We detected 111 sources, 61 of them associated with rho Ophiuchi YSOs as identified from infrared observations with ISOCAM. Specifically, we detected 9 out of 11 Class I, 31 out of 48 Class II and 15 out 16 Class III objects. Six objects out of 21 classified Class III candidates are also detected. At the same time we suggest that 15 Class III candidates that remain undetected at log Lx < 28.3 are not related to the cloud population. The global detection rate is ~64%. We have achieved a flux sensitivity of ~5 x 10^{-15} erg s^{-1} cm^{-2}. The Lx to L_bol ratio shows saturation at a value of ~ -3.5 for stars with T_eff <= 5000 K or 0.7 M_sun as observed in the Orion Nebula. The plasma temperatures and the spectrum absorption show a decline with YSO class, with Class I YSOs being hotter and more absorbed than Class II and III YSOs. In one star (GY 266) with infrared counterpart in 2MASS and Spitzer catalogs we have detected a soft excess in the X-ray spectrum which is best fitted by a cold thermal component less absorbed than the main thermal component of the plasma. Such a soft component hints to the presence of plasma heated by shocks due to jets outside the dense circumstellar material.

preprint2010arXiv

Small-scale structure in the Rosette molecular cloud revealed by Herschel

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

preprint2010arXiv

Star formation in the outer Galaxy: membership and fundamental parameters of the young open cluster NGC 1893

Different environmental conditions can play a crucial role in determining final products of the star formation process and in this context, less favorable activities of star formation are expected in the external regions of our Galaxy. We studied the properties of the young open cluster NGC 1893 located about 12 Kpc from the galactic center, to investigate how different physical conditions can affect the process of star formation. By adopting a multiwavelength approach, we compiled a catalog extending from X-rays to NIR data to derive the cluster membership. In addition, optical and NIR photometric properties are used to evaluate the cluster parameters. We find 415 diskless candidate members plus 1061 young stellar objects with a circumstellar disk or class II candidate members, 125 of which are also Halpha emitters. Considering the diskless candidate members, we find that the cluster distance is 3.6$\pm$0.2 kpc and the mean interstellar reddening is E(B-V)=0.6$\pm$0.1 with evidence of differential reddening in the whole surveyed region. NGC 1893 contains a conspicuous population of pre-main sequence stars together with the well studied main sequence cluster population; we found a disk fraction of about 70% similar to that found in clusters of similar age in the solar neighbour and then, despite expected unfavorable conditions for star formation, we conclude that very rich young clusters can form also in the outer regions of our Galaxy.

preprint2010arXiv

Star formation triggered by HII regions in our Galaxy: First results for N49 from the Herschel infrared survey of the Galactic plane

It has been shown that by means of different physical mechanisms the expansion of HII regions can trigger the formation of new stars of all masses. This process may be important to the formation of massive stars but has never been quantified in the Galaxy. We use Herschel-PACS and -SPIRE images from the Herschel Infrared survey of the Galactic plane, Hi-GAL, to perform this study. We combine the Spitzer-GLIMPSE and -MIPSGAL, radio-continuum and sub-millimeter surveys such as ATLASGAL with Hi-GAL to study Young Stellar Objects (YSOs) observed towards Galactic HII regions. We select a representative HII region, N49, located in the field centered on l=30 degr observed as part of the Hi-GAL Science Demonstration Phase, to demonstrate the importance Hi-GAL will have to this field of research. Hi-GAL PACS and SPIRE images reveal a new population of embedded young stars, coincident with bright ATLASGAL condensations. The Hi-GAL images also allow us, for the first time, to constrain the physical properties of the newly formed stars by means of fits to their spectral energy distribution. Massive young stellar objects are observed at the borders of the N49 region and represent second generation massive stars whose formation has been triggered by the expansion of the ionized region. Hi-GAL enables us to detect a population of young stars at different evolutionary stages, cold condensations only being detected in the SPIRE wavelength range. The far IR coverage of Hi-GAL strongly constrains the physical properties of the YSOs. The large and unbiased spatial coverage of this survey offers us a unique opportunity to lead, for the first time, a global study of star formation triggered by HII regions in our Galaxy.

preprint2010arXiv

Testing the theory of grain growth and fragmentation by millimeter observations of protoplanetary disks

Context. Observations at sub-millimeter and mm wavelengths will in the near future be able to resolve the radial dependence of the mm spectral slope in circumstellar disks with a resolution of around a few AU at the distance of the closest star-forming regions. Aims. We aim to constrain physical models of grain growth and fragmentation by a large sample of (sub-)mm observations of disks around pre-main sequence stars in the Taurus-Auriga and Ophiuchus star-forming regions. Methods. State-of-the-art coagulation/fragmentation and disk-structure codes are coupled to produce steady-state grain size distributions and to predict the spectral slopes at (sub-)mm wavelengths. Results. This work presents the first calculations predicting the mm spectral slope based on a physical model of grain growth. Our models can quite naturally reproduce the observed mm-slopes, but a simultaneous match to the observed range of flux levels can only be reached by a reduction of the dust mass by a factor of a few up to about 30 while keeping the gas mass of the disk the same. This dust reduction can either be due to radial drift at a reduced rate or during an earlier evolutionary time (otherwise the predicted fluxes would become too low) or due to efficient conversion of dust into larger, unseen bodies.

preprint2010arXiv

The (sub-)millimeter SED of protoplanetary disks in the outskirts of the Orion Nebula Cluster

We present the sub-mm/mm SED for a sample of eight young circumstellar disks in the outer regions of the Orion Nebula Cluster. New observations were carried out at 2.9 mm with the CARMA array and for one disk, 216-0939, at 3.3 and 6.8 mm with ATCA. By combining these new millimeter data with literature measurements at sub-millimeter wavelengths we investigate grain growth and measure the dust mass in protoplanetary disks in the Orion Nebula Cluster. These data provide evidence for dust grain growth to at least millimeter-sizes for the first time in a high-mass star forming region. The obtained range in sub-mm/mm spectral index, namely 1.5-3.2, indicates that for disks in the outskirts of the Orion Nebula Cluster (projected distance from the cluster center between about 0.4 pc and 1.5 pc) grain growth to mm sizes occurs in the same manner as disks in regions where only low-mass stars form. Finally, in our sample three disks are more massive than about $0.05\,M_\odot$, confirming that massive disks are present in the outer regions of the Orion Nebula.

preprint2010arXiv

The Aquila prestellar core population revealed by Herschel

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

preprint2010arXiv

The Herschel first look at protostars in the Aquila Rift

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

preprint2010arXiv

The molecular environment of the Galactic star forming region G19.61-0.23

We present single-dish (sub)millimeter observations of gas and dust in the Galactic high-mass star-forming region G19.61-0.23, with the aim of studying the large-scale properties and physical conditions of the molecular gas across the region. The final aim is to compare the large-scale (about 100 pc) properties with the small-scale (about 3 pc) properties and to consider possible implications for extragalactic studies. We have mapped CO isotopologues in the J=1-0 transition using the FCRAO-14m telescope and the J=2-1 transition using the IRAM-30m telescope. We have also used data from the ATLASGAL survey and from the BU-FCRAO Galactic Ring Survey, as well as the Spitzer infrared Galactic plane surveys GLIMPSE and MIPSGAL to characterize the star-formation activity within the molecular clouds. We reveal a population of molecular clumps in the 13CO(1-0) emission. Our analysis of the 13CO suggests that the virial parameter (ratio of kinetic to gravitational energy) varies over an order of magnitude between clumps that are unbound and some that are apparently "unstable". This conclusion is independent of whether they show evidence of ongoing star formation. We find that the majority of ATLASGAL sources have MIPSGAL counterparts with luminosities in the range 10^4 - 5 10^4 Lsun and are presumably forming relatively massive stars. We compare our results with previous extragalactic studies of the nearby spiral galaxies M31 and M33; and the blue compact dwarf galaxy Henize2-10. We find that the main giant molecular cloud surrounding G19.61-0.23 has physical properties typical for Galactic GMCs and which are comparable to the GMCs in M31 and M33. However, the GMC studied here shows smaller surface densities and masses than the clouds identified in Henize2-10 and associated with super star cluster formation.

preprint2010arXiv

U-band study of the accretion properties in the sigma Ori star-forming region

This paper presents the results of an U band survey with FORS1/VLT of a large area in the sigma Orionis star-forming region. We combine the U-band photometry with literature data to compute accretion luminosity and mass accretion rates from the U-band excess emission for all objects (187) detected by Spitzer in the FORS1 field and classified by Hernandez et al. (2007) as likely members of the cluster. The sample stars range in mass from ~0.06 to ~1.2 Msun; 72 of them show evidence of disks and we measure mass accretion rates Macc between <10^{-11} and few 10^{-9} Msun/yr, using the colors of the diskless stars as photospheric templates. Our results confirm the dependence of Macc on the mass of the central object, which is stronger for low-mass stars and flattens out for masses larger than ~0.3 Msun; the spread of Macc for any value of the stellar mass is ~2 orders of magnitude. We discuss the implications of these results in the context of disk evolution models. Finally, we analyze the relation between Macc and the excess emission in the Spitzer bands, and find that at Macc ~10^{-10} Msun/yr the inner disks change from optically thin to optically thick.

preprint2010arXiv

X-shooter observations of the accreting brown dwarf J053825.4-024241

We present the first observations of a probable brown dwarf, obtained with the new spectrograph X-shooter mounted on the UT2@VLT. The target (2MASS J053825.4-024241) is a 0.06 Msun object in the star-formation region sigma Orionis. The X-shooter spectrum covers simultaneously the whole range from UV to NIR (300-2500 nm). The J053825.4-024241 spectrum is rich in emission lines that are typical of accreting young object and clearly shows the Balmer jump. Moreover, many photospheric atomic and molecular absorption lines yield the spectral type and confirm that the object is young. We compute the mass accretion rate from all available observed accretion diagnostics. We find that there is a large spread in the Macc values (up to a factor 40) that is not caused by variability; some of this spread may be intrinsic, i.e., owing to different physical conditions of the emitting region for the same Macc. However, within the large error bars all Macc measurements agree, and the mean value is logMacc ~ -9.86 +- 0.45 Myr. The hydrogen Balmer lines are clearly detected up to n=25. Their ratios suggest that the emitting region is cold (T~2000-3000 K), dense and in thermal equilibrium (LTE), and that the lines are optically thick up to n~21. We briefly discuss the implications of this result for magnetospheric accretion models.

preprint2009arXiv

ATLASGAL - The APEX Telescope Large Area Survey of the Galaxy at 870 microns

(Abridged) Studying continuum emission from interstellar dust is essential to locating and characterizing the highest density regions in the interstellar medium. In particular, the early stages of massive star formation remain poorly understood. Our goal is to produce a large-scale, systematic database of massive pre- and proto-stellar clumps in the Galaxy, to understand how and under what conditions star formation takes place. A well characterized sample of star-forming sites will deliver an evolutionary sequence and a mass function of high-mass, star-forming clumps. This systematic survey at submm wavelengths also represents a preparatory work for Herschel and ALMA. The APEX telescope is ideally located to observe the inner Milky Way. The Large APEX Bolometer Camera (LABOCA) is a 295-element bolometer array observing at 870 microns, with a beam size of 19". Taking advantage of its large field of view (11.4') and excellent sensitivity, we started an unbiased survey of the Galactic Plane, with a noise level of 50-70 mJy/beam: the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL). As a first step, we covered 95 sq. deg. These data reveal 6000 compact sources brighter than 0.25 Jy, as well as extended structures, many of them filamentary. About two thirds of the compact sources have no bright infrared counterpart, and some of them are likely to correspond to the precursors of (high-mass) proto-stars or proto-clusters. Other compact sources harbor hot cores, compact HII regions or young embedded clusters. Assuming a typical distance of 5 kpc, most sources are clumps smaller than 1 pc with masses from a few 10 to a few 100 M_sun. In this introductory paper, we show preliminary results from these ongoing observations, and discuss the perspectives of the survey.

preprint2009arXiv

Dust Properties of Protoplanetary Disks in the Taurus-Auriga Star Forming Region from Millimeter Wavelengths

We present the most sensitive 3 mm-survey to date of protoplanetary disks carried in the Taurus-Auriga star forming region (average rms of about 0.3 mJy), using the IRAM PdBI. With our high detection rate of 17/19, we provide the first detections at wavelengths longer than about 1 mm for 12 sources. This enables us to study statistically the mm SED slopes and dust properties of faint disks and compare them to brighter disks using a uniform analysis method. With these new data and literature measurements at sub-millimeter and millimeter wavelengths, we analyze the dust properties of a sample of 21 isolated disks around T Tauri stars in the Taurus-Auriga star forming region. Together with the information about the disks spatial extension from sub/mm-mm interferometric studies, we derive from the observed sub-mm/mm spectral energy distribution constraints on the dust opacity law at these wavelengths, using two-layer flared disk models and a self-consistent dust model that takes properly into account the variation of the dust opacity with grain growth. We find evidence for the presence in the disk midplane of dust particles that have grown to sizes as large as at least 1 millimeter in all the disks of our sample, confirming what was previously observed on smaller brighter objects. This indicates that the dust coagulation from ISM dust to mm-sized grains is a very fast process in protoplanetary disks, that appears to occur before a young stellar object enters the Class II evolutionary stage. Also, the amount of these large grains in the disk outer regions is stationary throughout all the Class II evolutionary stage, indicating that mechanisms slowing down the dust inward migration are playing an important role in the Taurus-Auriga protoplanetary disks.

preprint2009arXiv

Resolving the molecular environment of Super Star Clusters in Henize 2-10

The rate of star formation both in the Galaxy and in external galaxies should be related to the physical properties of the molecular clouds from which stars form. This is expected for the starbursts found both in irregular galaxies and in some mergers. The dwarf galaxy Henize 2-10 is particularly interesting in this context as it shows a number of newly formed Super Star Clusters (SSCs) associated with a very rich molecular environment. We present a high angular resolution study of the molecular gas associated with the SSCs with the aim of deriving the physical properties of the parent molecular clouds. The final goal is to test the expectation that the formation of SSCs requires exceptionally dense and massive clouds. We have used the Submillimeter Array with an angular resolution of 1.9 X 1.3 to map the J=2-1 transition of CO in Henize 2-10. Supplementary measurements of HCN(J=1-0), 13CO(J=2-1) and millimeter continuum were obtained with the APEX, IRAM-30m and SEST single dish telescopes. Our single dish observations confirm the association of the newly formed SSCs in Henize 2-10 with dense molecular gas. Our interferometric observations resolve the CO(2-1) emission in several giant molecular clouds. Overall the molecular gas accounts for approximately half of the mass in the central regions of Henize 2-10. Although we find indications that the molecular clouds associated with the formation of SSCs in Henize 2-10 are massive and dense, the tracer we used (CO) and the linear resolution of our observations (60 X 80 pc) are still not adequate to test the expectation that exceptionally dense and massive cores are required for SSCs formation.

preprint2009arXiv

Results from DROXO. II. - [Ne II] and X-ray emission from Rho Ophiuchi young stellar objects

The infrared [Ne II] and [Ne III] fine structure lines at 12.81um and 15.55um are predicted to trace the circumstellar disk gas subject to X-ray heating and ionization. We investigate the origin of these lines by comparing observations with models of X-ray irradiated disks and by searching for empirical correlations between the line luminosities and stellar and circumstellar parameters. We measure neon line fluxes and X-ray luminosities for 28 young stellar objects in the Rho Ophiuchi star formation region for which good quality infrared spectra and X-ray data have been obtained, the former with the Spitzer IRS and the latter with the Deep Rho Ophiuchi XMM-Newton Observation. We detect the [Ne II] and the [Ne III] lines in 10 and 1 cases, respectively. Line luminosities show no correlation with X-ray emission. The luminosity of the [Ne II] line for one star, and that of both the [Ne II] and [Ne III] lines for a second star, match the predictions of published models of X-ray irradiated disks; for the remaining 8 objects the [Ne II] emission is 1-3 dex higher than predicted on the basis of their L_X. Class I objects show significantly stronger [Ne II] lines than Class II and Class III ones. A correlation is moreover found between the [Ne II] line emission and the disk mass accretion rates. This might point toward a role of accretion-generated UV emission in the generation of the line or to other mechanisms related to mass inflows from circumstellar disks and envelopes and/or to the associated mass outflows (winds and jets). We conclude that the X-ray luminosity is not the only parameter that determines the [Ne II] emission. Explaining the strong [Ne II] emission of Class I objects likely requires the inclusion in the models of additional physical components such as the envelope, inflows, and outflows.

preprint2009arXiv

Strong near-infrared emission in the sub-AU disk of the Herbig Ae star HD163296: evidence for refractory dust?

We present new long-baseline spectro-interferometric observations of the HerbigAe star HD163296 obtained in the H and K bands with the AMBER instrument at VLTI. The observations cover a range of spatial resolutions between 3 and 12 milli-arcseconds, with a spectral resolution of ~30. With a total of 1481 visibilities and 432 closure phases, they result in the best (u,v) coverage achieved on a young star so far. The circumstellar material is resolved at the sub-AU spatial scale and closure phase measurements indicate a small but significant deviation from point-symmetry. We discuss the results assuming that the near-infrared excess in HD163296 is dominated by the emission of a circumstellar disk. A successful fit to the spectral energy distribution, near-infrared visibilities and closure phases is found with a model where a dominant contribution to the H and K band emissions arises from an optically thin, smooth and point-symmetric region extending from about 0.1 to 0.45 AU. At the latter distance from the star, silicates condense, the disk becomes optically thick and develops a puffed-up rim, whose skewed emission can account for the non-zero closure phases. We discuss the nature of the inner disk emission and tentatively rule out dense molecular gas as well as optically thin atomic or ionized gas as its possible origin. We propose instead that the inner emission traces the presence of very refractory grains in a partially cleared region, extending at least to 0.5 AU. If so, we may be observing the disk of HD163296 just before it reaches the transition disk phase. However, we note that the nature of the refractory grains or even the possibility for any grain to survive at the very high temperatures we require (~2100-2300 K at 0.1 AU from the star) is unclear and should be investigated further.

preprint2005arXiv

First AMBER/VLTI observations of hot massive stars

AMBER is the first near infrared focal instrument of the VLTI. It combines three telescopes and produces spectrally resolved interferometric measures. This paper discusses some preliminary results of the first scientific observations of AMBER with three Unit Telescopes at medium (1500) and high (12000) spectral resolution. We derive a first set of constraints on the structure of the circumstellar material around the Wolf Rayet Gamma2 Velorum and the LBV Eta Carinae.

preprint2005arXiv

First results from a VLBA proper motion survey of H2O masers in low-mass YSOs: the Serpens core and RNO15-FIR

This article reports first results of a long-term observational program aimed to study the earliest evolution of jet/disk systems in low-mass YSOs by means of VLBI observations of the 22.2 GHz water masers. We report here data for the cluster of low-mass YSOs in the Serpens molecular core and for the single object RNO~15-FIR. Towards Serpens SMM1, the most luminous sub-mm source of the Serpens cluster, the water maser emission comes from two small (< 5 AU in size) clusters of features separated by ~25 AU, having line of sight velocities strongly red-shifted (by more than 10 km/s) with respect to the LSR velocity of the molecular cloud. The two maser clusters are oriented on the sky along a direction that is approximately perpendicular to the axis of the radio continuum jet observed with the VLA towards SMM1. The spatial and velocity distribution of the maser features lead us to favor the interpretation that the maser emission is excited by interaction of the receding lobe of the jet with dense gas in the accretion disk surrounding the YSO in SMM1. Towards RNO~15-FIR, the few detected maser features have both positions and (absolute) velocities aligned along a direction that is parallel to the axis of the molecular outflow observed on much larger angular scales. In this case the maser emission likely emerges from dense, shocked molecular clumps displaced along the axis of the jet emerging from the YSO. The protostar in Serpens SMM1 is more massive than the one in RNO~15-FIR. We discuss the case where a high mass ejection rate can generate jets sufficiently powerful to sweep away from their course the densest portions of circumstellar gas. In this case, the excitation conditions for water masers might preferably occur at the interface between the jet and the accretion disk, rather than along the jet axis.

preprint2004arXiv

Accretion in Brown Dwarfs: an Infrared View

This paper presents a study of the accretion properties of 19 very low mass objects (Mstar .01-0.1 Msun) in the regions Chamaeleon I and rho-Oph For 8 objects we obtained high resolution Halpha profiles and determined mass accretion rate Macc and accretion luminosity Lacc. Pabeta is detected in emission in 7 of the 10 rho-Oph objects, but only in one in Cha I. Using objects for which we have both a determination of Lacc from Halpha and a Pabeta detection, we show that the correlation between the Pabeta luminosity and luminosity Lacc, found by Muzerolle et al. (1998) for T Tauri stars in Taurus, extends to objects with mass approx 0.03 Msun; L(Pab) can be used to measure Lacc also in the substellar regime. The results were less conclusive for Brgamma, which was detected only in 2 objects, neither of which had an Halpha estimate of Macc. Using the relation between L(Pab) and Lacc we determined the accretion rate for all the objects in our sample (including those with no Halpha spectrum), more than doubling the number of substellar objects with known Macc. When plotted as a function of the mass of the central object together with data from the literature, our results confirm the trend of lower Macc for lower Mstar, although with a large spread. Some of the spread is probably due to an age effect; our very young objects in rho-Oph have on average an accretion rate at least one order of magnitude higher than objects of similar mass in older regions. As a side product, we found that the width of Halpha measured at 10% peak intensity is not only a qualitative indicator of accretion, but can be used to obtain a quantitative estimate of Macc over a large mass range, from T Tauri stars to brown dwarfs. Finally, we found that some of our objects show evidence of mass-loss.