Source author record

L. Podio

L. Podio 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

59works
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

59 published item(s)

preprint2025arXiv

ALMA chemical survey of disk-outflow sources in Taurus (ALMA-DOT) VII: the layered molecular outflow from HL Tau and its relationship with the ringed disk

The ringed disk around HL Tau stands out as the iconic signature of planet formation, but the origin of the substructures is still debated. The HL Tau system also drives a powerful bipolar wind, and we analyze its outermost component traced by CO emission, to determine the relationship of the flow with the disk and its substructures. We use ALMA observations of the ${}^{12}$CO (2-1) line at 1.3 mm, with 0.2 km/s and ~ 0.28" resolution, conducted within the ALMA-DOT project. The channel maps and position-velocity diagrams show a rich structure of concatenated bubble- and arc-shaped features, whose size and distance from the source increase with velocity. The superposition of the features generates the apparent conical shape. The tomographic reconstruction of the morphology and kinematics of the red-shifted lobe suggests the presence of distinct nested shells having higher velocity and steeper velocity gradient for shells closer to the axis, rotating in the same sense of the disk. Such configuration can be justified by different classes of models. In this paper we compare the derived wind parameters with the predictions of magnetohydrodynamic (MHD) disk winds. Under this hypothesis, the launch radii of the three outermost shells are found to be at about the position of three adjacent dust rings in the disk at 58, 72 and 86 au. The wind may be capable of removing angular momentum from the outer disk, and we derive a magnetic lever arm of $λ\sim 4 - 5$, higher than that commonly adopted for MHD winds from these regions. Interpretations are discussed. The arrangement of the wind in nested shells with brighter emission rooted at the location of ring substructures could support the results of non-ideal MHD simulations according to which magnetic instabilities can generate the disk ring-gap system with a connected layered wind, alternatively to the action of yet undetected protoplanets.

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

FAUST VI. VLA 1623--2417 B: a new laboratory for astrochemistry around protostars on 50 au scale

The ALMA interferometer, with its unprecedented combination of high-sensitivity and high-angular resolution, allows for (sub-)mm wavelength mapping of protostellar systems at Solar System scales. Astrochemistry has benefited from imaging interstellar complex organic molecules in these jet-disk systems. Here we report the first detection of methanol (CH3OH) and methyl formate (HCOOCH3) emission towards the triple protostellar system VLA1623-2417 A1+A2+B, obtained in the context of the ALMA Large Program FAUST. Compact methanol emission is detected in lines from Eu = 45 K up to 61 K and 537 K towards components A1 and B, respectively. LVG analysis of the CH3OH lines towards VLA1623-2417 B indicates a size of 0.11-0.34 arcsec (14-45 au), a column density N(CH3OH) = 10^16-10^17 cm-2, kinetic temperature > 170 K, and volume density > 10^8 cm-3. An LTE approach is used for VLA1623-2417 A1, given the limited Eu range, and yields Trot < 135 K. The methanol emission around both VLA1623-2417 A1 and B shows velocity gradients along the main axis of each disk. Although the axial geometry of the two disks is similar, the observed velocity gradients are reversed. The CH3OH spectra from B shows two broad (4-5 km s-1) peaks, which are red- and blue-shifted by about 6-7 km s-1 from the systemic velocity. Assuming a chemically enriched ring within the accretion disk, close to the centrifugal barrier, its radius is calculated to be 33 au. The methanol spectra towards A1 are somewhat narrower (about 4 km s-1), implying a radius of 12-24 au.

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.

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.

preprint2021arXiv

Chemical survey of Class I protostars with the IRAM-30m

Class I protostars are a bridge between Class 0 protostars, and Class II protoplanetary disks. Recent studies show gaps and rings in the dust distribution of disks younger than 1 Myr, suggesting that planet formation may start already at the Class I stage. To understand what chemistry planets will inherit, it is crucial to characterize the chemistry of Class I sources and to investigate how chemical complexity evolves from Class 0 protostars to protoplanetary disks. The goal is twofold: to obtain a census of the molecular complexity in a sample of four Class I protostars, and to compare it with the chemical compositions of earlier and later phases of the Sun-like star formation process. We performed IRAM-30m observations towards Class I objects (L1489-IRS, B5-IRS1, L1455-IRS1, and L1551-IRS5). The column densities of the detected species are derived assuming LTE or LVG. We detected 27 species: C-chains, N-bearing, S-bearing, Si-bearing species, deuterated molecules, and iCOMs. Different spectral profiles are observed: narrow lines towards all the sources, broader lines towards L1551-IRS5, and line wings due to outflows. Narrow c-C3H2 emission originates from the envelope. The iCOMs in L1551-IRS5 reveal the occurrence of hot corino chemistry, with CH3OH and CH3CN lines originating from a compact and warm region. Finally, OCS and H2S seem to probe the circumbinary disks in the L1455-IRS1 and L1551-IRS5 binary systems. The deuteration in terms of elemental D/H in the molecular envelopes and hot corino are derived. In addition, B5 IRS1, L1455-IRS1 and L1551-IRS5 show a low excitation methanol line, suggesting an origin from an extended structure, plausibly UV illuminated. The abundance ratios of iCOMs with respect CH3OH measured towards the L1551-IRS5 hot corino and the deuteration in our sample are comparable to that estimated at earlier stages, as well as to that found in comets.

preprint2021arXiv

Organic chemistry in the protosolar analogue HOPS-108: Environment matters

Hot corinos are compact regions around solar-mass protostellar objects that are very rich in interstellar complex organic molecules (iCOMs). They are believed to represent the very early phases of our Solar System's birth, which was very likely also characterized by rich organic chemistry. While most of the studied hot corinos are either isolated or born in a loose protocluster, our Sun was born in a densely packed star cluster, near massive stars whose ultraviolet radiation must have contributed to shaping the evolution of the surrounding environment. In addition, internal irradiation from energetic particles ($>$10 Mev), whose imprint is seen today in the products of short-lived radionuclides in meteoritic material, is also known to have occurred during the Solar System formation. How did all these conditions affect the chemistry of the proto-Sun and its surroundings is still an open question. To answer this question, we studied HOPS-108, the hot corino located in the protosolar analogue OMC-2 FIR4. The study was carried out with ALMA at 1.3mm with an angular resolution of $\sim$100 AU. We detected 11 iCOMs such as CH$_{3}$OH HCOOCH$_{3}$ and CH$_{3}$OCH$_{3}$. Our results can be summarized as follows: (1) an enhancement of HCOOCH3 with respect to other hot corinos, (2) a [CH$_{3}$OCH$_{3}$]/[HCOOCH$_{3}$] abundance ratio of $\sim$0.2 marginally deviating from the usual trend seen in other sources ([CH$_{3}$OCH$_{3}$]/[HCOOCH$_{3}$] $\sim$1), (3) a [CH$_{2}$DOH]/[CH$_{3}$OH] ratio of 2.5\% which is lower than what is seen in Perseus and Ophiuchus hot corinos ($\sim$7\%-9\%) and similar to that seen in HH212 another source located in Orion. This might result from different physical conditions in the Orion molecular complex compared to other regions.

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

Constraining MHD disk winds with ALMA. Apparent rotation signatures and application to HH212

Large millimeter interferometers are revealing a growing number of rotating outflows, which are suggested to trace magneto-centrifugal disk winds (MHD DWs). However, their impact on disk accretion is not yet well quantified. Here we identify systematic biases in retrieving the true launch zone, magnetic lever arm, and angular momentum flux of an MHD DW from apparent rotation signatures. Synthetic position-velocity cuts are constructed from self-similar MHD DWs over a broad range of parameters, and three different methods are applied for estimating the specific angular momentum. We find that the launch radius inferred using the well-known relation from Anderson et al. (2006) can markedly differ from the true outermost launch radius $r_{out}$ of the DW. The "double-peak separation" and "flow width" methods provide only a strict lower limit to $r_{out}$. This bias is independent of angular resolution and can reach a factor ten. In contrast, the "rotation curve" method gives a good estimate of $r_{out}$ when the flow is well resolved, and an upper limit otherwise. The magnetic lever arm is always underestimated. Only comparison with synthetic predictions can take into account properly all observational effects. As an application, we present a comparison with ALMA observations of HH212 at resolutions from 250 au to 16 au, which represents the most stringent observational test of MHD DW to date. This comparison confirms our predicted biases for the double-peak separation method, and the large $r_{out}\sim40~$au and small magnetic lever arm first suggested by Tabone et al. (2017). We also derive the first accurate analytical expression for the fraction of disk angular momentum extracted by an MHD disk wind of given radial extent, magnetic lever arm, and mass flux. Application to HH212 confirms that MHD DWs are serious candidates for the steady angular momentum extraction process in young disks.

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

No nitrogen fractionation on 600 au scale in the Sun progenitor analogue OMC-2 FIR4

We show the first interferometric maps of the 14N/15N ratio obtained with the Atacama Large Millimeter Array (ALMA) towards the Solar-like forming protocluster OMC-2 FIR4. We observed N2H+, 15NNH+, N15NH+ (1-0), and N2D+ (2-1), from which we derive the isotopic ratios 14N/15N and D/H. The target is one of the closest analogues of the environment in which our Sun may have formed. The ALMA images, having synthesised beam corresponding to ~600 au, show that the emission of the less abundant isotopologues is distributed in several cores of ~10" (i.e. ~0.02~pc or 4000 au) embedded in a more extended N2H+ emission. Overall, our results indicate that: (1) 14N/15N does not change across the region at core scales, and (2) 14N/15N does not depend on temperature variations. Our findings also suggest that the 14N/15N variations found in pristine Solar System objects are likely not inherited from the protocluster stage, and hence their reason has to be found elsewhere.

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

Searching for kinematic evidence of Keplerian disks around Class 0 protostars with CALYPSO

The formation of protoplanetary disks is not well understood. To understand how and when these disks are formed, it is crucial to characterize the kinematics of the youngest protostars at a high angular resolution. Here we study a sample of 16 Class 0 protostars to measure their rotation profile at scales from 50 to 500~au and search for Keplerian rotation. We used high-angular-resolution line observations obtained with the Plateau de Bure Interferometer as part of the CALYPSO large program. From $\mathrm{^{13}CO}$ ($J=2-1$), $\mathrm{C^{18}O}$ ($J=2-1$) and SO ($N_{j}=5_{6}-4_{5}$) moment maps, we find that seven sources show rotation about the jet axis at a few hundred au scales: SerpS-MM18, L1448-C, L1448-NB, L1527, NGC1333- IRAS2A, NGC1333-IRAS4B, and SVS13-B. We analyzed the kinematics of these sources in the $uv$ plane to derive the rotation profiles down to 50~au scales. We find evidence for Keplerian rotation in only two sources, L1527 and L1448-C. Overall, this suggests that Keplerian disks larger than 50 au are uncommon around Class 0 protostars. However, in some of the sources, the line emission could be optically thick and dominated by the envelope emission. Due to the optical thickness of these envelopes, some of the disks could have remained undetected in our observations.

preprint2020arXiv

Seeds of Life in Space (SOLIS) V. Methanol and acetaldehyde in the protostellar jet-driven shocks L1157-B0 and B1

Aim: In the past, observations of protostellar shocks have been able to set constraints on the formation route of formamide (NH2CHO), exploiting its observed spatial distribution and comparison with astrochemical model predictions. In this work, we follow the same strategy to study the case of acetaldehyde (CH3CHO). Method: To this end, we used the data obtained with the IRAM-NOEMA interferometer in the framework of the Large Program SOLIS to image the B0 and B1 shocks along the L1157 blueshifted outflow in methanol (CH3OH) and acetaldehyde line emission. Results: We imaged six CH3OH and eight CH3CHO lines which cover upper level energies up to 30 K. Both species trace the B0 molecular cavity as well as the northern B1 portion, i.e. the regions where the youngest shocks (1000 yr) occurred. The CH$_3$OH and CH$_3$CHO emission peaks towards the B1b clump, where we measured the following column densities and relative abundances: 1.3 x 10^16 cm-2 and 6.5 x 10-6 (methanol), and 7 x 10^13 cm-2 and 3.5 x 10-8 (acetaldehyde). We carried out a non-LTE LVG analysis of the observed CH3OH line: the average kinetic temperature and density of the emitting gas are Tkin = 90 K and nH2 = 4 x 10^5 cm-3, respectively. The CH3OH and CH3CHO abundance ratio towards B1b is 190, varying by less than a factor 3 throughout the whole B0-B1 structure. Conclusions: The comparison of astrochemical model predictions with the observed methanol and acetaldehyde spatial distribution does not allow to distinguish whether acetaldehyde is formed on the grain mantles or rather on the gas-phase, as its gas-phase formation, dominated by the reaction of ethyl radical (CH3CH2) with atomic oxygen, is very fast. Observations of acetaldehyde in younger shocks, e.g. 10^2 yr old, or/and of the ethyl radical, whose frequencies are not presently available, are necessary to settle the issue.

preprint2020arXiv

Seeds of Life in Space (SOLIS) VIII. SiO isotopic fractionation and a new insight in the shocks of L1157-B1

L1157-B1 is one of the outflow shocked regions along the blue-shifted outflow driven by the Class 0 protostar L1157-mm, and is an ideal laboratory to study the material ejected from the grains in very short timescales, i.e. its chemical composition is representative of the composition of the grains. We imaged $^{28}$SiO, $^{29}$SiO and $^{30}$SiO J = 2-1 emission towards L1157-B1 and B0 with the NOrthern Extended Millimeter Array (NOEMA) interferometer as part of the Seeds of Life in Space (SOLIS) large project. We present here a study of the isotopic fractionation of SiO towards L1157-B1. Furthermore, we use the high spectral resolution observations on the main isotopologue, $^{28}$SiO, to study the jet impact on the dense gas. We present here also single-dish observations obtained with the IRAM 30m telescope and Herschel-HIFI. We carried out a non-LTE analysis using a Large Velocity Gradient (LVG) code to model the single-dish observations. From our observations we can show that (i) the (2-1) transition of the main isotopologue is optically thick in L1157-B1 even at high velocities, and (ii) the [$^{29}$SiO/$^{30}$SiO] ratio is constant across the source, and consistent with the solar value of 1.5. We report the first isotopic fractionation maps of SiO in a shocked region and show the absence of a mass dependent fractionation in $^{29}$Si and $^{30}$Si across L1157-B1. A high-velocity bullet in $^{28}$SiO has been identified, showing the signature of a jet impacting on the dense gas. With the dataset presented in this paper, both interferometric and single-dish, we were able to study in great detail the gas shocked at the B1a position and its surrounding gas.

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

Seeds of Life in Space SOLIS. IX. Chemical segregation of $\rm SO_2$ and SO toward the low-mass protostellar shocked region of L1157

We present observations of SO and $\rm SO_2$ lines toward the shocked regions along the L1157 chemically rich outflow, taken in the context of the Seeds Of Life In Space IRAM-NOrthern Extended Millimeter Array Large Program, and supported by data from Submillimeter Array and IRAM-30 m telescope at 1.1--3.6 mm wavelengths. We simultaneously analyze, for the first time, all of the brightest shocks in the blueshifted lobe, namely, B0, B1, and B2. We found the following. (1) SO and $\rm SO_2$ may trace different gas, given that the large(-scale) velocity gradient analysis indicates for $\rm SO_2$ a volume density ($\rm 10^5\text{--}10^6\,cm^{-3}$) denser than that of the gas emitting in SO by a factor up to an order of magnitude. (2) Investigating the 0.1 pc scale field of view, we note a tentative gradient along the path of the precessing jet. More specifically, $\rm χ({SO/SO_2})$ decreases from the B0-B1 shocks to the older B2. (3) At a linear resolution of 500--1400 au, a tentative spatial displacement between the two emitting molecules is detected, with the SO peak closer (with respect to $\rm SO_2$) to the position where the recent jet is impinging on the B1 cavity wall. Our astrochemical modeling shows that the SO and $\rm SO_2$ abundances evolve on timescales less than about 1000 years. Furthermore, the modeling requires high abundances ($2\times10^{-6}$) of both $\rm H_2S/H$ and S/H injected in the gas phase due to the shock occurrence, so pre-frozen OCS only is not enough to reproduce our new observations.

preprint2019arXiv

Investigating the nature of the extended structure around the Herbig star RCrA using integral field and high-resolution spectroscopy

We present a detailed analysis of the extended structure detected around the young and close-by Herbig Ae/Be star RCrA. This is a young triple system with an intermediate mass central binary whose separation is of the order of a few tens of the radii of the individual components, and an M-star companion at about 30 au. Our aim is to understand the nature of the extended structure by means of combining integral-field and high-resolution spectroscopy. We conducted the analysis based on FEROS archival optical spectroscopy data and adaptive optics images and integral-field spectra obtained with SINFONI and SPHERE at the VLT. The observations reveal a complex extended structure that is composed of at least two components: a non-uniform wide cavity whose walls are detected in continuum emission up to 400~au, and a collimated wiggling-jet detected in the emission lines of Helium and Hydrogen. Moreover, the presence of [FeII] emission projected close to the cavity walls suggests the presence of a slower moving wind, most likely a disk wind. The multiple components of the optical forbidden lines also indicate the presence of a high-velocity jet co-existing with a slow wind. We constructed a geometrical model of the collimated jet flowing within the cavity using intensity and velocity maps, finding that its wiggling is consistent with the orbital period of the central binary. The cavity and the jet do not share the same position angle, suggesting that the jet is itself experiencing a precession motion possibly due to the wide M-dwarf companion. We propose a scenario that closely agrees with the general expectation of a magneto-centrifugal-launched jet. These results build upon the extensive studies already conducted on RCrA.

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

Hot methanol from the inner region of the HH 212 protostellar system

The mechanisms leading to the formation of disks around young stellar objects (YSOs) and to the launching of the associated jets are crucial to the understanding of the earliest stages of star and planet formation. HH 212 is a privileged laboratory to study a pristine jet-disk system. Therefore we investigate the innermost region ($<100$ AU) around the HH 212-MM1 protostar through ALMA band\,7 observations of methanol. The 8 GHz bandwidth spectrum towards the peak of the continuum emission of the HH 212 system reveals at least 19 transitions of methanol. Several of these lines (among which several vibrationally excited lines in the v$_{\rm t}=1,2$ states) have upper energies above 500 K. They originate from a compact ($<135$ AU in diameter), hot ($\sim 295$ K) region elongated along the direction of the SiO jet. We performed a fit in the $uv$ plane of various velocity channels of the strongest high-excitation lines. The blue- and red-shifted velocity centroids are shifted roughly symmetrically on either side of the jet axis, indicating that the line-of-sight velocity beyond 0.7 km s$^{-1}$ from systemic is dominated by rotational motions. The velocity increases moving away from the protostar further indicating that the emission of methanol is not associated with a Keplerian disk or rotating-infalling cavity, and it is more likely associated with outflowing gas. We speculate that CH$_3$OH traces a disk wind gas accelerated at the base. The launching region would be at a radius of a few astronomical units from the YSO.

preprint2016arXiv

Phosphorus-bearing molecules in solar-type star forming regions: First PO detection

As part of the Large Program ASAI (Astrochemical Surveys At IRAM), we have used the IRAM 30m telescope to lead a systematic search for the emission of rotational transitions of P-bearing species between 80 and 350 GHz towards L1157-B1, a shock position in the solar-type star forming region L1157. We report the detection of several transitions of PN and, for the first time, of prebiotic molecule PO. None of these species are detected towards the driving protostar of the outflow L1157-mm. Analysis of the line profiles shows that PN arises from the outflow cavity, where SiO, a strong shock tracer, is produced. Radiative transfer analysis yields an abundance of 2.5e-9 and 0.9e-9 for PO and PN, respectively. These results imply a strong depletion (approx 100) of Phosphorus in the quiescent cloud gas. Shock modelling shows that atomic N plays a major role in the chemistry of PO and PN. The relative abundance of PO and PN brings constraints both on the duration of the pre-shock phase, which has to be about 1 Myr, and on the shock parameters. The maximum temperature in the shock has to be larger than 4000K, which implies a shock velocity of 40 km/s.

preprint2016arXiv

Probing the CO and methanol snow lines in young protostars. Results from the CALYPSO IRAM-PdBI survey

Context. "Snow lines", marking regions where abundant volatiles freeze out onto the surface of dust grains, play an important role for planet growth and bulk composition in protoplanetary disks. They can already be observed in the envelopes of the much younger, low-mass Class 0 protostars that are still in their early phase of heavy accretion. Aims. We aim at using the information on the sublimation regions of different kinds of ices to understand the chemistry of the envelope, its temperature and density structure, and the history of the accretion process. Methods. As part of the CALYPSO IRAM Large Program, we have obtained observations of C$^{18}$O, N$_2$H$^+$ and CH$_3$OH towards nearby Class 0 protostars with the IRAM Plateau de Bure interferometer at sub-arcsecond resolution. For four of these sources we have modeled the emission using a chemical code coupled with a radiative transfer module. Results. We observe an anti-correlation of C$^{18}$O and N$_2$H$^+$ in NGC 1333-IRAS4A, NGC 1333-IRAS4B, L1157, and L1448C, with N$_2$H$^+$ forming a ring around the centrally peaked C$^{18}$O emission due to N$_2$H$^+$ being chemically destroyed by CO. The emission regions of models and observations match for a CO binding energy of 1200 K, which is higher than the binding energy of pure CO ices ($\sim$855 K). Furthermore, we find very low CO abundances inside the snow lines in our sources, about an order of magnitude lower than the total CO abundance observed in the gas on large scales in molecular clouds before depletion sets in. Conclusions. The high CO binding energy may hint at CO being frozen out in a polar ice environment like amorphous water ice or in non-polar CO$_2$-rich ice. The low CO abundances are comparable to values found in protoplanetary disks, which may indicate an evolutionary scenario where these low values are already established in the protostellar phase. (Abbr. Version)

preprint2016arXiv

Sub-0.1" optical imaging of the Z CMa jets with SPHERE/ZIMPOL

Crucial information on the mass accretion-ejection in young stars can be obtained from high spatial resolution images of jets in sources with recurrent accretion outbursts. Using the SPHERE/ZIMPOL instrument, we observed the young binary Z CMa that is composed of a Herbig Be star and a FUor object, both driving a jet. We analyse the structure of the two jets in relation with previous accretion events observed in this target. We obtained optical images in the Halpha and OI 6300A lines at the unprecedented angular resolution of ~0.03 arcsec, on which we have performed both continuum subtraction and deconvolution, thereby deriving results that are consistent with each other. Our images reveal extended emission from both sources: a fairly compact and poorly collimated emission SW of the Herbig component and an extended collimated and precessing jet from the FUor component. The compact emission from the Herbig star is compatible with a wide-angle wind and is possibly connected to the recent outburst events shown by this component. The FUor jet is traced down to 70 mas (80 AU) from the source and is highly collimated with a width of 26-48 AU at distances 100-200 AU, which is similar to the width of jets from T Tauri stars. This strongly suggests that the same magneto-centrifugal jet-launching mechanism also operates in FUors. The observed jet wiggle can be modelled as originating from an orbital motion with a period of 4.2 yr around an unseen companion with mass between 0.48 and 1 Msun. The jet mass loss rate Mloss was derived from the OI luminosity and is in the range 1E-8 and 1E-6 Msun/yr. This is the first direct Mloss measurement from a jet in a FUor. If we assume previous Macc estimates obtained through modelling of the accretion disk, the derived range of Mloss would imply a very low mass-ejection efficiency (Mloss/Macc < 0.02), which is lower than that typical of T Tauri stars.

preprint2016arXiv

The small-scale HH34 IRS jet as seen by X-shooter

Very little information has been so far gathered on atomic jets from young embedded low mass sources (class I stars), especially in the inner jet region. We exploit multiwave spectroscopic observations to infer physical conditions of the inner region of HH34 IRS, a prototypical class I jet. We use a deep X-shooter spectrum ( lambda = 350-2300 nm, R between 8000 and 18000) detecting lines with upper energies from 8000 to 31000 cm-1 . Statistical equilibrium and ionization models are adopted to derive the jet main physical parameters. We derive the physical conditions for the extended high velocity jet (HVC, V about -100 km/s) and for the low velocity and compact gas (LVC, V about -20-50 \kms). At the jet base (< 200 AU) the HVC is mostly neutral (x(e)< 0.1) and very dense (n(H) > 5x10^5 cm-3). The LVC gas has the same density and A_V as the HVC, but it is a factor of two colder. Iron abundance in the HVC is close to solar, while it is 3 times subsolar in the LVC, suggesting an origin of the LVC gas from a dusty disk. We infer a jet mass flux rate > 5x10^{-7} Mo/yr. We find that the relationships between accretion luminosity and line luminosity derived for T Tauri stars cannot be directly extended to class I sources surrounded by reflection nebulae since the permitted lines are seen through scattered light. An accretion rate of about 7\,10^{-6} Mo/yr is instead derived from an empirical relationship connecting the mass accretion rate with L([OI])(HVC). In conclusion, the HH34 IRS jet shares many properties with jets from active CTTs stars. It is however denser as a consequence of the larger mass flux rate. These findings suggest that the acceleration and excitation mechanisms in jets are not influenced by evolution and are similar in CTTs and in still embedded and highly accreting sources.

preprint2016arXiv

Water and acetaldehyde in HH212: The first hot corino in Orion

Aims: Using the unprecedented combination of high resolution and sensitivity offered by ALMA, we aim to investigate whether and how hot corinos, circumstellar disks, and ejected gas are related in young solar-mass protostars. Methods: We observed CH$_3$CHO and deuterated water (HDO) high-excitation ($E_{\rm u}$ up to 335 K) lines towards the Sun-like protostar HH212--MM1. Results: For the first time, we have obtained images of CH$_3$CHO and HDO emission in the inner $\simeq$ 100 AU of HH212. The multifrequency line analysis allows us to contrain the density ($\geq$ 10$^{7}$ cm$^{-3}$), temperature ($\simeq$ 100 K), and CH$_3$CHO abundance ($\simeq$ 0.2--2 $\times$ 10$^{-9}$) of the emitting region. The HDO profile is asymmetric at low velocities ($\leq$ 2 km s$^{-1}$ from $V_{\rm sys}$). If the HDO line is optically thick, this points to an extremely small ($\sim$ 20--40 AU) and dense ($\ge$ 10$^{9}$ cm$^{-3}$) emitting region. Conclusions: We report the first detection of a hot corino in Orion. The HDO asymmetric profile indicates a contribution of outflowing gas from the compact central region, possibly associated with a dense disk wind.

preprint2015arXiv

Disks, Jets and the dawn of planets, Proceedings of the 2nd JEDI meeting

This booklet contains a collection of contributions to the meeting of the JEts and Disks at INAF (JEDI) group, which took place at the Capodimonte Observatory during 9-10 April 2015. Scope of the meeting was to bring together the JEDI researchers of the Italian Istituto Nazionale di Astrofisica (INAF) working in the field of circumstellar disks and jets in young stars, to discuss together the different agents affecting the structure and the evolution of disks, namely accretion, jets and winds. More information on the JEDI group and its activities can be found at \texttt{http://www.oa-roma.inaf.it/irgroup/JEDI}.

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

Solving the excitation and chemical abundances in shocks: the case of HH1

We present deep spectroscopic (3600 - 24700 A) X-shooter observations of the bright Herbig-Haro object HH1, one of the best laboratories to study the chemical and physical modifications caused by protostellar shocks on the natal cloud. We observe atomic fine structure lines, HI, and He, recombination lines and H_2, ro-vibrational lines (more than 500 detections in total). Line emission was analyzed by means of Non Local Thermal Equilibiurm codes to derive the electron temperature and density, and, for the first time, we are able to accurately probe different physical regimes behind a dissociative shock. We find a temperature stratification in the range 4000 - 80000 K, and a significant correlation between temperature and ionization energy. Two density regimes are identified for the ionized gas, a more tenuous, spatially broad component (density about 10^3 cm^-3), and a more compact component (density > 10^5 cm^-3) likely associated with the hottest gas. A further neutral component is also evidenced, having temperature lass than 10000 K and density > 10^4 cm^-3. The gas fractional ionization was estimated solving the ionization equilibrium equations of atoms detected in different ionization stages. We find that neutral and fully ionized regions co-exist inside the shock. Also, indications in favor of at least partially dissociative shock as the main mechanism for molecular excitation are derived. Chemical abundances are estimated for the majority of the detected species. On average, abundances of non-refractory/refractory elements are lower than solar of about 0.15/0.5 dex. This testifies the presence of dust inside the medium, with a depletion factor of Iron of about 40%.

preprint2015arXiv

The jet and the disk of the HH 212 low-mass protostar imaged by ALMA: SO and SO2 emission

To investigate the disk formation and jet launch in protostars is crucial to comprehend the earliest stages of star and planet formation. We aim to constrain the properties of the molecular jet and the disk of the HH 212 protostellar system at unprecedented angular scales through ALMA observations of sulfur-bearing molecules, SO 9(8)-8(7), SO 10(11)-10(10), SO2 8(2,6)-7(1,7). SO 9(8)-8(7) and SO2 8(2,6)-7(1,7) show broad velocity profiles. At systemic velocity they probe the circumstellar gas and the cavity walls. Going from low to high blue-/red-shifted velocities the emission traces the wide-angle outflow and the fast (~100-200 km/s) and collimated (~90 AU) molecular jet revealing the inner knots with timescales <50 years. The jet transports a mass loss rate >0.2-2e-6 Msun/yr, implying high ejection efficiency (>0.03-0.3). The SO and SO2 abundances in the jet are ~1e-7-1e-6. SO 10(11)-10(10) emission is compact and shows small-scale velocity gradients indicating that it originates partly from the rotating disk previously seen in HCO+ and C17O, and partly from the base of the jet. The disk mass is >0.002-0.013 Msun, and the SO abundance in the disk is ~1e-8-1e-7. SO and SO2 are effective tracers of the molecular jet in the inner few hundreds AU from the protostar. Their abundances indicate that 1% - 40% of sulfur is in SO and SO2 due to shocks in the jet/outflow and/or to ambipolar diffusion at the wind base. The SO abundance in the disk is 3-4 orders of magnitude larger than in evolved protoplanetary disks. This may be due to an SO enhancement in the accretion shock at the envelope-disk interface or in spiral shocks if the disk is partly gravitationally unstable.

preprint2014arXiv

[OI] disk emission in the Taurus star forming region

The structure of protoplanetary disks is thought to be linked to the temperature and chemistry of their dust and gas. Whether the disk is flat or flaring depends on the amount of radiation that it absorbs at a given radius, and on the efficiency with which this is converted into thermal energy. The understanding of these heating and cooling processes is crucial to provide a reliable disk structure for the interpretation of dust continuum emission and gas line fluxes. Especially in the upper layers of the disk, where gas and dust are thermally decoupled, the infrared line emission is strictly related to the gas heating/cooling processes. We aim to study the thermal properties of the disk in the oxygen line emission region, and to investigate the relative importance of X-ray (1-120 Angstrom) and far-UV radiation (FUV, 912-2070 Angstrom) for the heating balance there. We use [OI] 63 micron line fluxes observed in a sample of protoplanetary disks of the Taurus/Auriga star forming region and compare it to the model predictions presented in our previous work. The data were obtained with the PACS instrument on board the Herschel Space Observatory as part of the Herschel Open Time Key Program GASPS (GAS in Protoplanetary diskS), published in Howard et al. (2013). Our theoretical grid of disk models can reproduce the [OI] absolute fluxes and predict a correlation between [OI] and the sum Lx+Lfuv. The data show no correlation between the [OI] line flux and the X-ray luminosity, the FUV luminosity or their sum. The data show that the FUV or X-ray radiation has no notable impact on the region where the [OI] line is formed. This is in contrast with what is predicted from our models. Possible explanations are that the disks in Taurus are less flaring than the hydrostatic models predict, and/or that other disk structure aspects that were left unchanged in our models are important. ..abridged..

preprint2014arXiv

Astrochemistry at work in the L1157-B1 shock: acetaldehyde formation

The formation of complex organic molecules (COMs) in protostellar environments is a hotly debated topic. In particular, the relative importance of the gas phase processes as compared to a direct formation of COMs on the dust grain surfaces is so far unknown. We report here the first high-resolution images of acetaldehyde (CH$_3$CHO) emission towards the chemically rich protostellar shock L1157-B1, obtained at 2 mm with the IRAM Plateau de Bure interferometer. Six blueshifted CH$_3$CHO lines with $E_{\rm u}$ = 26-35 K have been detected. The acetaldehyde spatial distribution follows the young ($\sim$ 2000 yr) outflow cavity produced by the impact of the jet with the ambient medium, indicating that this COM is closely associated with the region enriched by iced species evaporated from dust mantles and released into the gas phase. A high CH$_3$CHO relative abundance, 2-3 $\times$ 10$^{-8}$, is inferred, similarly to what found in hot-corinos. Astrochemical modelling indicates that gas phase reactions can produce the observed quantity of acetaldehyde only if a large fraction of carbon, of the order of 0.1%, is locked into iced hydrocarbons.

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

Empirical determination of Einstein A-coefficient ratios of bright [Fe II] lines

The Einstein spontaneous rates (A-coefficients) of Fe^+ lines have been computed by several authors, with results that differ from each other up to 40%. Consequently, models for line emissivities suffer from uncertainties which in turn affect the determination of the physical conditions at the base of line excitation. We provide an empirical determination of the A-coefficient ratios of bright [Fe II] lines, which would represent both a valid benchmark for theoretical computations and a reference for the physical interpretation of the observed lines. With the ESO-VLT X-shooter instrument between 3,000 A, and 24,700 A, we obtained a spectrum of the bright Herbig-Haro object HH1. We detect around 100 [Fe II] lines, some of which with a signal-to-noise ratio > 100. Among these latter, we selected those emitted by the same level, whose de-reddened intensity ratio is a direct function of the Einstein A-coefficient ratios. From the same X-shooter spectrum, we got an accurate estimate of the extinction toward HH1 through intensity ratios of atomic species, HI, recombination lines and H_2 ro-vibrational transitions. We provide seven reliable A-ooefficient ratios between bright [Fe II] lines, which are compared with the literature determinations. In particular, the A-coefficient ratios involving the brightest near-infrared lines (12570A/16440A and 13209A/16440A) are better in agreement with the predictions by Quinet et al. (1996) Relativistic Hartree-Fock model. However, none of the theoretical models predicts A-coefficient ratios in agreement with all our determinations. We also show that literature data of near-infrared intensity ratios better agree with our determinations than with theoretical expectations.

preprint2014arXiv

ESO-H$α$~574 and Par-Lup3-4 Jets: Exploring the spectral, kinematical and physical properties

In this paper a comprehensive analysis of VLT / X-Shooter observations of two jet systems, namely ESO-H$α$ 574 a K8 classical T Tauri star and Par-Lup 3-4 a very low mass (0.13~\Msun) M5 star, is presented. Both stars are known to have near-edge on accretion disks. A summary of these first X-shooter observations of jets was given in a 2011 letter. The new results outlined here include flux tables of identified emission lines, information on the morphology, kinematics and physical conditions of both jets and, updated estimates of $\dot{M}_{out}$ / $\dot{M}_{acc}$. Asymmetries in the \eso flow are investigated while the \para jet is much more symmetric. The density, temperature, and therefore origin of the gas traced by the Balmer lines are investigated from the Balmer decrements and results suggest an origin in a jet for \eso while for \para the temperature and density are consistent with an accretion flow. $\dot{M}_{acc}$ is estimated from the luminosity of various accretion tracers. For both targets, new luminosity relationships and a re-evaluation of the effect of reddening and grey extinction (due to the edge-on disks) allows for substantial improvements on previous estimates of $\dot{M}_{acc}$. It is found that log($\dot{M}_{acc}$) = -9.15 $\pm$ 0.45~\Msun yr$^{-1}$ and -9.30 $\pm$ 0.27~\Msun yr$^{-1}$ for \eso and \para respectively. Additionally, the physical conditions in the jets (electron density, electron temperature, and ionisation) are probed using various line ratios and compared with previous determinations from iron lines. The results are combined with the luminosity of the [SII]$λ$6731 line to derive $\dot{M}_{out}$ through a calculation of the gas emissivity based on a 5-level atom model.

preprint2014arXiv

Fast molecular jet from L1157-mm

L1157-mm powers a molecular outflow that is well-known for its shock-induced chemical activity in several hot-spots. We have studied the molecular emission toward L1157-mm searching for a jet component responsible for these spots. We used the IRAM 30m telescope to observe the vicinity of L1157-mm in several lines of SiO. The SiO(5-4) and SiO(6-5) spectra toward L1157-mm present blue and red detached components about 45 km/s away from the ambient cloud. These extremely high-velocity (EHV) components are similar to those found in the L1448 and IRAS 04166+2706 outflows, and probably arise from a molecular jet driven by L1157-mm. Observations of off-center positions indicate that the jet is unresolved in SiO(5-4) (<11"). The EHV jet seen in SiO probably excites L1157-B1 and the other chemically active spots of the L1157 outflow.

preprint2014arXiv

Herschel -PACS observations of [OI] and H2O in Cha II

Gas plays a major role in the dynamical evolution of protoplanetary discs. Its coupling with the dust is the key to our understanding planetary formation. Studying the gas content is therefore a crucial step towards understanding protoplanetary discs evolution. Such a study can be made through spectroscopic observations of emission lines in the far-infrared, where some of the most important gas coolants emit, such as the [OI] 3P1-3 P2 transition at 63.18 microns. We aim at characterising the gas content of protoplanetary discs in the intermediate-aged Chamaeleon II (Cha II) star forming region. We also aim at characterising the gaseous detection fractions within this age range, which is an essential step tracing gas evolution with age in different star forming regions. We obtained Herschel-PACS line scan spectroscopic observations at 63 microns of 19 Cha II Class I and II stars. The observations were used to trace [OI] and o-H2O at 63 microns. The analysis of the spatial distribution of [OI], when extended, can be used to understand the origin of the emission. We have detected [OI] emission toward seven out of the nineteen systems observed, and o-H2O emission at 63.32 microns in just one of them, Sz 61. Cha II members show a correlation between [OI] line fluxes and the continuum at 70 microns, similar to what is observed in Taurus. We analyse the extended [OI] emission towards the star DK Cha and study its dynamical footprints in the PACS Integral Field Unit (IFU). We conclude that there is a high velocity component from a jet combined with a low velocity component with an origin that may be a combination of disc, envelope and wind emission. The stacking of spectra of objects not detected individually in [OI] leads to a marginal 2.6sigma detection that may indicate the presence of gas just below our detection limits for some, if not all, of them.

preprint2014arXiv

Molecular ions in the protostellar shock L1157-B1

We perform a complete census of molecular ions with an abundance larger than 1e-10 in the protostellar shock L1157-B1 by means of an unbiased high-sensitivity survey obtained with the IRAM-30m and Herschel/HIFI. By means of an LVG radiative transfer code the gas physical conditions and fractional abundances of molecular ions are derived. The latter are compared with estimates of steady-state abundances in the cloud and their evolution in the shock calculated with the chemical model Astrochem. We detect emission from HCO+, H13CO+, N2H+, HCS+, and, for the first time in a shock, from HOCO+, and SO+. The bulk of the emission peaks at blueshifted velocity, ~ 0.5-3 km/s with respect to systemic, has a width of ~ 4-8 km/s, and is associated with the outflow cavities (T_kin ~ 20-70 K, n(H2) ~ 1e5 cm-3). Observed HCO+ and N2H+ abundances are in agreement with steady-state abundances in the cloud and with their evolution in the compressed and heated gas in the shock for cosmic rays ionization rate Z = 3e-16 s-1. HOCO+, SO+, and HCS+ observed abundances, instead, are 1-2 orders of magnitude larger than predicted in the cloud; on the other hand they are strongly enhanced on timescales shorter than the shock age (~2000 years) if CO2, S or H2S, and OCS are sputtered off the dust grains in the shock. The performed analysis indicates that HCO+ and N2H+ are a fossil record of pre-shock gas in the outflow cavity, while HOCO+, SO+, and HCS+ are effective shock tracers and can be used to infer the amount of CO2 and sulphur-bearing species released from dust mantles in the shock. The observed HCS+ (and CS) abundance indicates that OCS should be one of the main sulphur carrier on grain mantles. However, the OCS abundance required to fit the observations is 1-2 orders of magnitude larger than observed. Further studies are required to fully understand the chemistry of sulphur-bearing species.

preprint2014arXiv

Physical properties of the jet from DG Tauri on sub-arcsecond scales with HST/STIS

We derive the physical properties at the base of the jet from DG Tau both along and across the flow and as a function of velocity. We analysed seven optical spectra of the DG Tau jet, taken with the Hubble Space Telescope Imaging Spectrograph. The spectra were obtained by placing a long-slit parallel to the jet axis and stepping it across the jet width. The resulting position-velocity diagrams in optical forbidden emission lines allowed access to plasma conditions via calculation of emission line ratios. We find at the base of the jet high electron density, $n_e \sim $ 10$^5$, and very low ionisation, $x_e \sim 0.02-0.05$, which combine to give a total density up to $n_H \sim $ 3 10$^6$. This analysis confirms previous reports of variations in plasma parameters along the jet, (i.e. decrease in density by several orders of magnitude, increase of $x_e$ from 0.05 to a plateau at 0.7 downstream at 2$''$ from the star). Furthermore, a spatial coincidence is revealed between sharp gradients in the total density and supersonic velocity jumps. This strongly suggests that the emission is caused by shock excitation. The position-velocity diagrams indicate the presence of both fast accelerating gas and slower, less collimated material. We derive the mass outflow rate, $\dot{M}_j$, in the blue-shifted lobe in different velocity channels, that contribute to a total of $\dot{M}_j \sim$ 8 $\pm$ 4 10$^{-9}$ M$_\odot$ yr$^{-1}$. We estimate that a symmetric bipolar jet would transport at the low and intermediate velocities probed by rotation measurements, an angular momentum flux of $\dot{L}_j \sim$ 2.9 $\pm$ 1.5 10$^{-6}$ M$_\odot$ yr$^{-1}$ AU km s$^{-1}$. The derived properties of the DG Tau jet are demonstrated to be consistent with magneto-centrifugal theory. However, non-stationary modelling is required in order to explain all of the features revealed at high resolution.

preprint2014arXiv

Probing the gaseous disk of T Tau N with CN 5-4 lines

We present spectrally resolved Herschel/HIFI observations of the young multiple system T Tau in atomic and molecular lines. While CO, H2O, [C II], and SO lines trace the envelope and the outflowing gas up to velocities of 33 km/s with respect to systemic, the CN 5-4 hyperfine structure lines at 566.7, 566.9 GHz show a narrow double-peaked profile centered at systemic velocity, consistent with an origin in the outer region of the compact disk of T Tau N. Disk modeling of the T Tau N disk with the thermo-chemical code ProDiMo produces CN line fluxes and profiles consistent with the observed ones and constrain the size of the gaseous disk (R_out = 110 (+10, -20) AU) and its inclination (i = 25 \pm 5 degree). The model indicates that the CN lines originate in a disk upper layer at 40-110 AU from the star, which is irradiated by the stellar UV field and heated up to temperatures of 50-700 K. With respect to previously observed CN 2-1 millimeter lines, the CN 5-4 lines appear to be less affected by envelope emission, due to their larger critical density and excitation temperature. Hence, high-J CN lines are a unique confusion-free tracer of embedded disks, such as the disk of T Tau N.

preprint2014arXiv

Relating jet structure to photometric variability: the Herbig Ae star HD 163296

Herbig Ae/Be stars are intermediate-mass pre-main sequence stars surrounded by circumstellar dust disks. Some are observed to produce jets, whose appearance as a sequence of shock fronts (knots) suggests a past episodic outflow variability. This "jet fossil record" can be used to reconstruct the outflow history. We present the first optical to near-infrared (NIR) VLT/X-shooter spectra of the jet from the Herbig Ae star HD 163296. We determine physical conditions in the knots, as well as their kinematic "launch epochs". Knots are formed simultaneously on either side of the disk, with a regular interval of ~16 yr. The velocity dispersion versus jet velocity and the energy input are comparable in both lobes. However, the mass loss rate, velocity, and shock conditions are asymmetric. We find Mjet/Macc ~ 0.01-0.1, consistent with magneto-centrifugal jet launching models. No evidence for dust is found in the high-velocity jet, suggesting it is launched within the sublimation radius (<0.5 au). The jet inclination measured from proper motions and radial velocities confirms it is perpendicular to the disk. A tentative relation is found between the structure of the jet and the photometric variability of the source. Episodes of NIR brightening were previously detected and attributed to a dusty disk wind. We report for the first time significant optical fadings lasting from a few days up to a year, coinciding with the NIR brightenings. These are likely caused by dust lifted high above the disk plane; this supports the disk wind scenario. The disk wind is launched at a larger radius than the high-velocity atomic jet, although their outflow variability may have a common origin. No significant relation between outflow and accretion variability could be established. Our findings confirm that this source undergoes periodic ejection events, which may be coupled with dust ejections above the disk plane.

preprint2014arXiv

The ALMA view of the protostellar system HH212 - The wind, the cavity, and the disk

Context: Because it is viewed simply edge-on, the HH212 protostellar system is an ideal laboratory for studying the interplay of infall, outflow, and rotation in the earliest stages of low-mass star formation. Aims: We wish to exploit the unmatched combination of high angular resolution, high sensitivity, high-imaging fidelity, and spectral coverage provided by ALMA to shed light on the complex kinematics of the innermost central regions of HH212. Methods: We mapped the inner 10" (4500 AU) of the HH212 system at about 0.5 arcsec resolution in several molecular tracers and in the 850 $μ$m dust continuum using the ALMA interferometer in band 7 in the extended configuration of the Early Science Cycle 0 operations. Results: Within a single ALMA spectral set-up, we simultaneously identify all the crucial ingredients known to be involved in the star formation recipe namely: (i) the fast, collimated bipolar SiO jet driven by the protostar, (ii) the large-scale swept-up CO outflow, (iii) the flattened rotating and infalling envelope, with bipolar cavities carved by the outflow (in C$^{17}$O(3--2)), and (iv) a rotating wide-angle flow that fills the cavities and surrounding the axial jet (in C$^{34}$S(7--6)). In addition, the compact high-velocity C$^{17}$O emission ($\pm$ 1.9--3.5 km s$^{-1}$ from systemic) shows a velocity gradient along the equatorial plane consistent with a rotating disk of about 0.2 arcsec = 90 AU around a $\simeq 0.3 \pm 0.1 M_{\rm \odot}$ source. The rotating disk is possibly Keplerian. Conclusions: HH212 is the third Class 0 protostar with possible signatures of a Keplerian disk of radius $\geq 30 AU$. The warped geometry in our CS data suggests that this large keplerian disk might result from misaligned magnetic and rotation axes during the collapse phase. The wide-angle CS flow suggests that disk winds may be present in this source.

preprint2014arXiv

The density structure of the L1157 molecular outflow

We present a multiline CS survey towards the brightest bow-shock B1 in the prototypical chemically active protostellar outflow L1157. We made use of (sub-)mm data obtained in the framework of the Chemical HErschel Surveys of Star forming regions (CHESS) and Astrochemical Surveys at IRAM (ASAI) key science programs. We detected $^{12}$C$^{32}$S, $^{12}$C$^{34}$S, $^{13}$C$^{32}$S, and $^{12}$C$^{33}$S emissions, for a total of 18 transitions, with $E_{\rm u}$ up to $\sim$ 180 K. The unprecedented sensitivity of the survey allows us to carefully analyse the line profiles, revealing high-velocity emission, up to 20 km s$^{-1}$ with respect to the systemic. The profiles can be well fitted by a combination of two exponential laws that are remarkably similar to what previously found using CO. These components have been related to the cavity walls produced by the $\sim$ 2000 yr B1 shock and the older ($\sim$ 4000 yr) B2 shock, respectively. The combination of low- and high-excitation CS emission was used to properly sample the different physical components expected in a shocked region. Our CS observations show that this molecule is highlighting the dense, $n_{\rm H_2}$ = 1--5 $\times$ 10$^{5}$ cm$^{-3}$, cavity walls produced by the episodic outflow in L1157. In addition, the highest excitation (E$_u$ $\geq$ 130 K) CS lines provide us with the signature of denser (1--5 $\times$ 10$^{6}$ cm$^{-3}$) gas, associated with a molecular reformation zone of a dissociative J-type shock, which is expected to arise where the precessing jet impacting the molecular cavities. The CS fractional abundance increases up to $\sim$ 10$^{-7}$ in all the kinematical components. This value is consistent with what previously found for prototypical protostars and it is in agreement with the prediction of the abundances obtained via the chemical code Astrochem.

preprint2013arXiv

Broad N2H+ emission towards the protostellar shock L1157-B1

We present the first detection of N2H+ towards a low-mass protostellar outflow, namely the L1157-B1 shock, at about 0.1 pc from the protostellar cocoon. The detection was obtained with the IRAM 30-m antenna. We observed emission at 93 GHz due to the J = 1-0 hyperfine lines. The analysis of the emission coupled with the HIFI CHESS multiline CO observations leads to the conclusion that the observed N2H+(1-0) line originates from the dense (> 10^5 cm-3) gas associated with the large (20-25 arcsec) cavities opened by the protostellar wind. We find a N2H+ column density of few 10^12 cm-2 corresponding to an abundance of (2-8) 10^-9. The N2H+ abundance can be matched by a model of quiescent gas evolved for more than 10^4 yr, i.e. for more than the shock kinematical age (about 2000 yr). Modelling of C-shocks confirms that the abundance of N2H+ is not increased by the passage of the shock. In summary, N2H+ is a fossil record of the pre-shock gas, formed when the density of the gas was around 10^4 cm-3, and then further compressed and accelerated by the shock.

preprint2013arXiv

GASPS - a Herschel survey of gas and dust in Protoplanetary Disks: Summary and Initial Statistics

GASPS is a far-infrared line and continuum survey of protoplanetary and young debris disks using PACS on the Herschel Space Observatory. The survey includes [OI] at 63 microns, as well as 70, 100 and 160um continuum, with the brightest objects also studied in [OI]145um, [CII]157um, H2O and CO. Targets included T Tauri stars and debris disks in 7 nearby young associations, and a sample of isolated Herbig AeBe stars. The aim was to study the global gas and dust content in a wide disk sample, systemically comparing the results with models. In this paper we review the main aims, target selection and observing strategy. We show initial results, including line identifications, sources detected, and a first statistical study. [OI]63um was the brightest line in most objects, by a factor of ~10. Detection rates were 49%, including 100% of HAeBe stars and 43% of T Tauri stars. Comparison with published dust masses show a dust threshold for [OI]63um detection of ~1e-5 M_solar. Normalising to 140pc distance, 32% with mass 1e-6 - 1e-5 M_solar, and a small number with lower mass were also detected. This is consistent with moderate UV excess and disk flaring. In most cases, continuum and line emission is spatially and spectrally unresolved, suggesting disk emission. ~10 objects were resolved, likely from outflows. Detection rates in [OI]145um, [CII]157um and CO J=18-17 were 20-40%, but [CII] was not correlated with disk mass, suggesting it arises instead from a compact envelope. [OI] detection rates in T Tauri associations of ages 0.3-4Myr were ~50%. ~2 stars were detectable in associations of 5-20Myr, with no detections in associations of age >20Myr. Comparing with the total number of young stars, and assuming a ISM-like gas/dust ratio, this indicates that ~18% of stars retain a gas-rich disk of total mass >1M_Jupiter for 1-4Myr, 1-7% keep such disks for 5-10Myr, and none remain beyond 10-20Myr.

preprint2013arXiv

The diagnostic potential of Fe lines applied to protostellar jets

We investigate the diagnostic capabilities of the iron lines for tracing the physical conditions of the shock-excited gas in jets driven by pre-main sequence stars. We have analyzed the 300-2500 nm X-shooter spectra of two jets driven by the pre-main sequence stars ESO-Halpha 574 and Par-Lup 3-4. Both spectra are very rich in [FeII] lines over the whole spectral range; in addition, lines from [FeIII] are detected in the ESO-Hα574 spectrum. NLTE codes along with codes for the ionization equilibrium are used to derive the gas excitation conditions of electron temperature and density, and fractional ionization. The iron gas-phase abundance is provided by comparing the iron lines emissivity with that of [OI] 630 nm. The [FeII] lines indicate ESO-Halpha 574 jet is, on average, colder (T_e = 9000 K), less dense (n_e = 2 10^4 cm^-3) and more ionized (x_e = 0.7) than the Par-Lup 3-4 jet (T_e = 13000 K, n_e = 6 10^4 cm^-3, x_e < 0.4), even if the existence of a higher density component (n_e = 2 10^5 cm^-3) is probed by the [FeIII] and [FeII] ultra-violet lines. Theoretical models suggest that the shock at work in ESO-Halpha 574 is faster and likely more energetic than the Par-Lup 3-4 shock. This latter feature is confirmed by the high percentage of gas-phase iron measured in ESO-Halpha 574 (50-60% of its solar abundance in comparison with less than 30% in Par-Lup 3-4), which testifies that the ESO-Halpha 574 shock is powerful enough to partially destroy the dust present inside the jet. This work demonstrates that a multiline Fe analysis can be effectively used to probe the excitation and ionization conditions of the gas in a jet without any assumption on ionic abundances. The main limitation on the diagnostics resides in the large uncertainties of the atomic data, which, however, can be overcome through a statistical approach involving many lines.

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.

preprint2012arXiv

Herschel/PACS observations of young sources in Taurus: the far-infrared counterpart of optical jets

Observations of the atomic and molecular line emission associated with jets and outflows emitted by young stellar objects can be used to trace the various evolutionary stages they pass through as they evolve to become main sequence stars. To understand the relevance of atomic and molecular cooling in shocks, and how accretion and ejection efficiency evolves with the source evolutionary state, we will study the far-infrared counterparts of bright optical jets associated with Class I and II sources in Taurus (T Tau, DG Tau A, DG Tau B, FS Tau A+B, and RW Aur). We have analysed Herschel/PACS observations of a number of atomic ([OI]63um, 145um, [CII]158um) and molecular (high-J CO, H2O, OH) lines, collected within the OTKP GASPS. To constrain the origin of the detected lines we have compared the FIR emission maps with the emission from optical-jets and millimetre-outflows, and the line fluxes and ratios with predictions from shock and disk models. All of the targets are associated with extended emission in the atomic lines correlated with the direction of the optical jet/mm-outflow. The atomic lines can be excited in fast dissociative J-shocks. The molecular emission, on the contrary, originates from a compact region, that is spatially and spectrally unresolved. Slow C- or J- shocks with high pre-shock densities reproduce the observed H2O and high-J CO lines; however, the disk and/or UV-heated outflow cavities may contribute to the emission. While the cooling is dominated by CO and H2O lines in Class 0 sources, [OI] becomes an important coolant as the source evolves and the environment is cleared. The cooling and mass loss rates estimated for Class II and I sources are one to four orders of magnitude lower than for Class 0 sources. This provides strong evidence to indicate that the outflow activity decreases as the source evolves.

preprint2012arXiv

The outflow history of two Herbig-Haro jets in RCW 36: HH 1042 and HH 1043

Jets around low- and intermediate-mass young stellar objects (YSOs) contain a fossil record of the recent accretion and outflow activity of their parent star-forming systems. We aim to understand whether the accretion/ejection process is similar across the entire stellar mass range of the parent YSOs. To this end we have obtained VLT/X-shooter spectra of HH 1042 and HH 1043, two newly discovered jets in the massive star-forming region RCW 36. HH 1042 is associated with the intermediate-mass YSO 08576nr292. Over 90 emission lines are detected in the spectra. High-velocity (up to 220 km/s) blue- and redshifted emission from a bipolar flow is observed in typical shock tracers. Low-velocity emission from the background cloud is detected in nebular tracers, including lines from high ionization species. We applied combined optical and infrared spectral diagnostic tools in order to derive the physical conditions (density, temperature, and ionization) in the jets. The measured mass outflow rates are Mjet ~ 10^-7 Msun/yr. We measure a high accretion rate for HH 1042 (Macc ~ 10^-6 Msun/yr) and Mjet/Macc ~ 0.1, comparable to low-mass sources and consistent with models for magneto-centrifugal jet launching. The knotted structure and velocity spread in both jets are interpreted as fossil signatures of a variable outflow rate. The mean velocities in both lobes of the jets are comparable, but the variations in Mjet and velocity in the two lobes are not symmetric, suggesting that the launching mechanism on either side of the accretion disk is not synchronized. For HH 1042, we have constructed an interpretative physical model with a stochastic or periodic outflow rate and a description of a ballistic flow as its constituents. The knotted structure and velocity spread can be reproduced qualitatively with the model, indicating that the outflow velocity varies on timescales on the order of 100 yr.

preprint2011arXiv

Detection of warm water vapour in Taurus protoplanetary discs by Herschel

Line spectra of 68 Taurus T Tauri stars were obtained with the Herschel-PACS (Photodetector Array Camera & Spectrometer) instrument as part of the GASPS (Gas Evolution in Protoplanetary Systems) survey of protoplanetary discs. A careful examination of the line scans centred on the [OI] 63.18 microns fine-structure line unveiled a line at 63.32 micron in some of these spectra. We identify this line with a transition of ortho-water. It is detected confidently (i.e., >3 sigma) in eight sources, i.e., 24% of the sub-sample with gas-rich discs. Several statistical tests were used to search for correlations with other disc and stellar parameters such as line fluxes of [OI] 6300 Armstrong and 63.18 microns; X-ray luminosity and continuum levels at 63 microns and 850 microns. Correlations are found between the water line fluxes and the [OI] 63.18 microns line luminosity, the dust continuum, and possibly with the stellar X-ray luminosity. This is the first time that this line of warm water vapour has been detected in protoplanetary discs. We discuss its origins, in particular whether it comes from the inner disc and/or disc surface or from shocks in outflows and jets. Our analysis favours a disc origin, with the observed water vapour line produced within 2-3AU from the central stars, where the gas temperature is of the order of 500-600 K.

preprint2010arXiv

Gas in Protoplanetary Discs (GASPS) 1. First results

Context - Circumstellar discs are ubiquitous around young stars, but rapidly dissipate their gas and dust on timescales of a few Myr. The Herschel space observatory allows for the study of the warm disc atmosphere, using far-infrared spectroscopy to measure gas content and excitation conditions, and far-IR photometry to constrain the dust distribution. Aims - We aim to detect and characterize the gas content of circumstellar discs in four targets as part of the Herschel science demonstration phase. Methods - We carried out sensitive medium resolution spectroscopy and high sensitivity photometry at lambda ~60-190 micron using the Photodetector Array Camera and Spectrometer instrument on the Herschel space observatory. Results - We detect [OI] 63 micron emission from the young stars HD 169142, TW Hydrae, and RECX 15, but not HD 181327. No other lines, including [CII] 158 and [OI] 145, are significantly detected. All four stars are detected in photometry at 70 and 160 micron. Extensive models are presented in associated papers.

preprint2010arXiv

Gas in the protoplanetary disc of HD 169142: Herschel's view

In an effort to simultaneously study the gas and dust components of the disc surrounding the young Herbig Ae star HD 169142, we present far-IR observations obtained with the PACS instrument onboard the Herschel Space Observatory. This work is part of the Open Time Key Project GASPS, which is aimed at studying the evolution of protoplanetary discs. To constrain the gas properties in the outer disc, we observed the star at several key gas-lines, including [OI] 63.2 and 145.5 micron, [CII] 157.7 micron, CO 72.8 and 90.2 micron, and o-H2O 78.7 and 179.5 micron. We only detect the [OI] 63.2 micron line in our spectra, and derive upper limits for the other lines. We complement our data set with PACS photometry and 12/13CO data obtained with the Submillimeter Array. Furthermore, we derive accurate stellar parameters from optical spectra and UV to mm photometry. We model the dust continuum with the 3D radiative transfer code MCFOST and use this model as an input to analyse the gas lines with the thermo-chemical code ProDiMo. Our dataset is consistent with a simple model in which the gas and dust are well-mixed in a disc with a continuous structure between 20 and 200 AU, but this is not a unique solution. Our modelling effort allows us to constrain the gas-to-dust mass ratio as well as the relative abundance of the PAHs in the disc by simultaneously fitting the lines of several species that originate in different regions. Our results are inconsistent with a gas-poor disc with a large UV excess; a gas mass of 5.0 +/- 2.0 times 10^(-3) Msun is still present in this disc, in agreement with earlier CO observations.

preprint2010arXiv

Herschel-PACS observation of the 10 Myr old T Tauri disk TW Hya: Constraining the disk gas mass

Planets are formed in disks around young stars. With an age of ~10 Myr, TW Hya is one of the nearest T Tauri stars that is still surrounded by a relatively massive disk. In addition a large number of molecules has been found in the TW Hya disk, making TW Hya the perfect test case in a large survey of disks with Herschel-PACS to directly study their gaseous component. We aim to constrain the gas and dust mass of the circumstellar disk around TW Hya. We observed the fine-structure lines of [OI] and [CII] as part of the Open-time large program GASPS. We complement this with continuum data and ground-based 12CO 3-2 and 13CO 3-2 observations. We simultaneously model the continuum and the line fluxes with the 3D Monte-Carlo code MCFOST and the thermo-chemical code ProDiMo to derive the gas and dust masses. We detect the [OI] line at 63 micron. The other lines that were observed, [OI] at 145 micron and [CII] at 157 micron, are not detected. No extended emission has been found. Preliminary modeling of the photometric and line data assuming [12CO]/[13CO]=69 suggests a dust mass for grains with radius < 1 mm of ~1.9 times 10^-4 Msun (total solid mass of 3 times 10^-3 Msun) and a gas mass of (0.5--5) times 10^-3 Msun. The gas-to-dust mass may be lower than the standard interstellar value of 100.

preprint2010arXiv

The Herschel view of Gas in Protoplanetary Systems (GASPS). First comparisons with a large grid of models

The Herschel GASPS Key Program is a survey of the gas phase of protoplanetary discs, targeting 240 objects which cover a large range of ages, spectral types, and disc properties. To interpret this large quantity of data and initiate self-consistent analyses of the gas and dust properties of protoplanetary discs, we have combined the capabilities of the radiative transfer code MCFOST with the gas thermal balance and chemistry code ProDiMo to compute a grid of 300 000 disc models (DENT). We present a comparison of the first Herschel/GASPS line and continuum data with the predictions from the DENT grid of models. Our objective is to test some of the main trends already identified in the DENT grid, as well as to define better empirical diagnostics to estimate the total gas mass of protoplanetary discs. Photospheric UV radiation appears to be the dominant gas-heating mechanism for Herbig stars, whereas UV excess and/or X-rays emission dominates for T Tauri stars. The DENT grid reveals the complexity in the analysis of far-IR lines and the difficulty to invert these observations into physical quantities. The combination of Herschel line observations with continuum data and/or with rotational lines in the (sub-)millimetre regime, in particular CO lines, is required for a detailed characterisation of the physical and chemical properties of circumstellar discs.

preprint2009arXiv

Classical T Tauri-like Outflow Activity in the Brown Dwarf Mass Regime

Over the last number of years spectroscopic studies have strongly supported the assertion that protostellar accretion and outflow activity persists to the lowest masses. In this paper we present the results of our latest investigation of brown dwarf (BD) outflow activity and report on the discovery of two new outflows. Here ISO-Oph 32 is shown to drive a blue-shifted outflow with a radial velocity of 10-20 km/s and spectro-astrometric analysis constrains the position angle of this outflow to 240 +/- 7 degrees. The BD candidate ISO-Cha1 217 is found to have a bipolar outflow bright in several key forbidden lines (radial velocity = -20 km/s, +40 km/s) and with a PA of 190-210 degrees. A striking feature of the ISO-Cha1 217 outflow is the strong asymmetry between the red and blue-shifted lobes. This asymmetry is revealed in the relative brightness of the two lobes (red-shifted lobe is brighter), the factor of two difference in radial velocity (the red-shifted lobe is faster) and the difference in the electron density (again higher in the red lobe). Such asymmetries are common in jets from low mass protostars and the observation of a marked asymmetry at such a low mass supports the idea that BD outflow activity is scaled down from low mass protostellar activity. In addition to presenting these new results, a comprehensive comparison is made between BD outflow activity and jets launched by CTTSs. In particular, the application of current methods for investigating the excitation conditions and mass loss rates in CTT jets to BD spectra is explored.

preprint2009arXiv

Physical Structure and Dust Reprocessing in a sample of HH Jets

We investigate the physical structure and dust reprocessing in the shocks along the beam of a number of classical Herbig-Haro jets in the Orion and Lupus molecular cloud. Spectral diagnostic techniques are applied to obtain the jet physical conditions from the ratios between selected forbidden lines. The presence of dust grains is investigated by estimating the gas-phase abundance of calcium with respect to its Solar value. We find the electron density, ne, varies between 0.05-4 10^3 cm^-3, the ionisation fraction, xe, is 0.01-0.7, the temperature, te, ranges between 0.6-3 10^4 K, and the hydrogen density between 0.01-6 10^4 cm^-3. Interestingly, in the HH 111 jet, ne, xe, and te, peak in the High Velocity Interval (HVI) of the strongest working surfaces, confirming the prediction from shocks models. Calcium turns out to be depleted with respect to its Solar value, but its gas-phase abundance is higher than that estimated in the interstellar medium in Orion. The depletion is high (up to 80%) along the low-excited jets, while low or no depletion is measured in the jets which show higher excitation conditions. Moreover, in HH 111 the depletion is lower in the HVI of the faster shock. Our results confirm the shock structure predicted by models and indicate that the shocks occurring along the jets are partially destroying the dust grains and that the efficiency of dust reprocessing strongly depend on shock velocity. However, the high Ca gas-phase abundance estimated in some of the knots is not well justified by existing models of dust reprocessing in shocks, and indicates that the dust must have been partially reprocessed in the region where the flow originates.