Search

We use the James Webb Space Telescope (JWST) and its Mid-Infrared Instrument (MIRI) (5-28 um), to study the embedded HH 211 flow. We map a 0.95'x0.22' region, covering the full extent of the blue-shifted lobe, the central protostellar region, and a small portion of the red-shifted lobe. The jet driving source is not detected even at the longest mid-IR wavelengths. The overall morphology of the flow consists of a highly collimated jet, mostly molecular (H2, HD) with an inner atomic ([FeI], [FeII], [SI], [NiII]) structure. The jet shocks the ambient medium, producing several large bow-shocks, rich in forbidden atomic and molecular lines, and is driving an H2 molecular outflow, mostly traced by low-J, v=0 transitions. Moreover, 0-0 S(1) uncollimated emission is also detected down to 2"-3" (~650-1000 au) from the source, tracing a cold (T=200-400 K), less dense and poorly collimated molecular wind. The atomic jet ([FeII] at 26 um) is detected down to ~130 au from source, whereas the lack of H2 emission close to the source is likely due to the large visual extinction. Dust continuum-emission is detected at the terminal bow-shocks, and in the blue- and red-shifted jet, being likely dust lifted from the disk. The jet shows an onion-like structure, with layers of different size, velocity, temperature, and chemical composition. Moreover, moving from the inner jet to the outer bow-shocks, different physical, kinematic and excitation conditions for both molecular and atomic gas are observed. The jet mass-flux rate, momentum, and momentum flux of the warm H2 component are up to one order of magnitude higher than those inferred from the atomic jet component. Our findings indicate that the warm H2 component is the primary mover of the outflow, namely it is the most significant dynamical component of the jet, in contrast to jets from more evolved YSOs, where the atomic component is dominant., Comment: Paper accepted in A&A for publication

The first generation of ELT instruments includes an optical-infrared high-resolution spectrograph, indicated as ELT-HIRES and recently christened ANDES (ArmazoNes high Dispersion Echelle Spectrograph). ANDES consists of three fibre-fed spectrographs ([U]BV, RIZ, YJH) providing a spectral resolution of $\sim$100,000 with a minimum simultaneous wavelength coverage of 0.4-1.8 $\mu$m with the goal of extending it to 0.35-2.4 $\mu$m with the addition of a U arm to the BV spectrograph and a separate K band spectrograph. It operates both in seeing- and diffraction-limited conditions and the fibre feeding allows several, interchangeable observing modes including a single conjugated adaptive optics module and a small diffraction-limited integral field unit in the NIR. Modularity and fibre-feeding allow ANDES to be placed partly on the ELT Nasmyth platform and partly in the Coud\'e room. ANDES has a wide range of groundbreaking science cases spanning nearly all areas of research in astrophysics and even fundamental physics. Among the top science cases, there are the detection of biosignatures from exoplanet atmospheres, finding the fingerprints of the first generation of stars, tests on the stability of Nature's fundamental couplings, and the direct detection of the cosmic acceleration. The ANDES project is carried forward by a large international consortium, composed of 35 Institutes from 13 countries, forming a team of almost 300 scientists and engineers which include the majority of the scientific and technical expertise in the field that can be found in ESO member states., Comment: SPIE astronomical telescope and instrumentation 2024, in press

Studying outflows is important as they may significantly contribute to angular momentum removal from the star/disk system, affecting disk evolution and planet formation. To investigate the different outflow components; the collimated jet, wide-angled molecular outflow, and outflow cavity, of the Class I HH 46/47 outflow system. We focus on their kinematics. We present Near Infrared (NIR) K-band integral field observations of the blue-shifted HH 46/47 outflow base obtained using VLT/SINFONI with an angular resolution of 0".81. Our analysis focuses on [Fe II], H2 1-0 S(1), and, Br-gamma emission. We employ a wavelength recalibration technique based on OH telluric lines to probe the kinematics of the wide-angled flow with an accuracy of 1 km/s - 3 km/s. A velocity gradient of 10 km/s transverse to the outflow direction is confirmed in the wide-angled H2 outflow cavity. The H2 cavity peaks at radial velocities of -15 km/s to -30 km/s, and the atomic jet at v = -210 km/s. The outflow exhibits a layered structure; the high-velocity [Fe II] and Br-gamma jet is surrounded by a wide-angled H2 outflow cavity, which is in turn nested within the continuum emission and CO molecular outflow. The continuum emission and H2 outflow cavity are asymmetric with respect to the jet axis. We propose that the origin of the asymmetries and the velocity gradient detected in the wide-angled H2 cavity, is due to a wide-angled outflow or successive jet bowshocks expanding into an inhomogeneous ambient medium, or the presence of a secondary outflow. We eliminate outflow rotation as an exclusive origin of this velocity gradient due to large specific angular momenta values, J(r)= 3000 - 4000 km/s au calculated from 1" to 2" along the outflow. The observations reveal the complexities inherent in outflow systems, and the risk of attributing transverse velocity gradients solely to rotation., Comment: 20 pages, 12 figures, 2 tables, submitted to A&A

We present the first results of the JWST program PROJECT-J (PROtostellar JEts Cradle Tested with JWST ), designed to study the Class I source HH46 IRS and its outflow through NIRSpec and MIRI spectroscopy (1.66 to 28 micron). The data provide line-images (~ 6.6" in length with NIRSpec, and up to 20" with MIRI) revealing unprecedented details within the jet, the molecular outflow and the cavity. We detect, for the first time, the red-shifted jet within ~ 90 au from the source. Dozens of shock-excited forbidden lines are observed, including highly ionized species such as [Ne III] 15.5 micron, suggesting that the gas is excited by high velocity (> 80 km/s) shocks in a relatively high density medium. Images of H2 lines at different excitations outline a complex molecular flow, where a bright cavity, molecular shells, and a jet-driven bow-shock interact with and are shaped by the ambient conditions. Additional NIRCam 2 micron images resolve the HH46 IRS ~ 110 au binary system and suggest that the large asymmetries observed between the jet and the H2 wide angle emission could be due to two separate outflows being driven by the two sources. The spectra of the unresolved binary show deep ice bands and plenty of gaseous lines in absorption, likely originating in a cold envelope or disk. In conclusion, JWST has unraveled for the first time the origin of the HH46 IRS complex outflow demonstrating its capability to investigate embedded regions around young stars, which remain elusive even at near-IR wavelengths., Comment: 28 pages, 15 figures, Accepted for publication on The Astrophysical Journal (9 April 2024)

EX Lupi is the prototype by which EXor-type outbursts were defined. It has experienced multiple accretion-related bursts and outbursts throughout the last decades, whose study have greatly extended our knowledge about the effects of these types of events. This star experienced a new burst in 2022. We used multi-band photometry to create color-color and color-magnitude diagrams to exclude the possibility that the brightening could be explained by a decrease in extinction. We obtained VLT/X-shooter spectra to determine the Lacc and Macc during the peak of the burst and after its return to quiescence using 2 methods: empirical relationships between line luminosity and Lacc, and a slab model of the whole spectrum. We examined the 130 year light curve of EX Lupi to provide statistics on the number of outbursts experienced during this period of time. Our analysis of the data taken during the 2022 burst confirmed that a change in extinction is not responsible for the brightening. Our two approaches in calculating the Macc were in agreement, and resulted in values that are 2 orders of magnitude above what had previously been estimated, thus suggesting that EX Lupi is a strong accretor even when in quiescence. We determined that in 2022 March the Macc increased by a factor of 7 with respect to the quiescent level. We also found hints that even though the Macc had returned to almost its pre-outburst levels, certain physical properties of the gas had not returned to the quiescent values. We found that the mass accreted during this three month event was 0.8 lunar masses, which is approximately half of what is accreted during a year of quiescence. We calculated that if EX Lupi remains as active as it has been for the past 130 years, during which it has experienced at least 3 outbursts and 10 bursts, then it will deplete the mass of its circumstellar material in less than 160000 yr., Comment: 18 pages, 12 figures, 4 tables. Accepted for publication in A&A

Observing the spatial distribution and excitation processes of atomic and molecular gas in the inner regions (< 20 au) of young (< 10 Myr) protoplanetary disks helps us to understand the conditions for the formation and evolution of planetary systems. In the framework of the PENELLOPE and ULLYSES projects, we aim to characterize the atomic and molecular component of protoplanetary disks in a sample of 11 Classical T Tauri Stars (CTTs) of the Orion OB1 and $\sigma$-Orionis associations. We analyzed the flux-calibrated optical-forbidden lines and the fluorescent ultraviolet $\rm H_2$ progressions using spectra acquired with ESPRESSO at VLT, UVES at VLT and HST-COS. Line morphologies were characterized through Gaussian decomposition. We then focused on the properties of the narrow low-velocity (FWHM < 40 $km$ $s^{-1}$ and |$v_p$| < 30 $km$ $s^{-1}$) component (NLVC) of the [OI] 630 nm line, compared with the properties of the UV-$\rm H_2$ lines. We found that the [OI]630 NLVC and the UV-$\rm H_2$ lines are strongly correlated in terms of peak velocities, full width at half maximum, and luminosity. The luminosities of the [OI]630 NLVC and UV-$\rm H_2$ correlate with the accretion luminosity with a similar slope, as well as with the luminosity of the CIV 154.8, 155 nm doublet. We discuss such correlations in the framework of the currently suggested excitation processes for the [OI]630 NLVC. Our results can be interpreted in a scenario in which the [OI]630 NLVC and UV-$\rm H_2$ have a common disk origin with a partially overlapped radial extension. We also suggest that the excitation of the [OI] NLVC is mainly induced by stellar FUV continuum photons more than being of thermal origin. This study demonstrates the potential of contemporaneous wide-band high-resolution spectroscopy in linking different tracers of protoplanetary disks., Comment: 24 pages, 15 figures, 6 tables. Accepted for publication in Astronomy & Astrophysics

In the framework of the GIARPS@TNG High-resolution Observations of T Tauri stars (GHOsT) project, we study the accretion properties of 37 Classical T Tauri Stars of the Taurus-Auriga star forming region (SFR) with the aim of characterizing their relation with the properties of the central star, of jets and disk winds, and of the global disk structure, in synergy with complementary ALMA millimiter observations. We derive stellar parameters, optical veiling, accretion luminosity ($\rm L_{acc}$) and mass accretion rate ($\rm \dot M_{acc}$) in a homogeneous and self-consistent way using high-resolution spectra acquired at the Telescopio Nazionale Galileo with the HARPS-N and GIANO spectrographs, and flux-calibrated based on contemporaneous low-resolution spectroscopic and photometric ancillary observations. The $\rm L_{acc}$-$\rm L_{\star}$, $\rm \dot{M}_{acc}$-$\rm M_{\star}$ and $\rm \dot{M}_{acc}$-$\rm M_{disk}$ relationships of the Taurus sample are provided and compared with those of the coeval SFRs of Lupus and Chamaeleon I. Our results demonstrate the potential of contemporaneous optical and near-infrared high-resolution spectroscopy to simultaneously provide precise measurements of stellar and accretion/wind properties of young stars., Comment: 30 pages, 18 figures, 6 tables. Accepted for publication on Astronomy and Astrophysics

The He I 1 micron line is a high excitation line which allows us to probe the innermost regions of protostellar disks, and to trace both accreting and outflowing material. We use X-Shooter observations of a sample of 107 young stars in the Lupus (1-3 Myr) and Upper Scorpius (5-10 Myr) star-forming regions to search for correlations between the line properties, as well as the disk inclination and accretion luminosity. We identified eight distinct profile types in the sample. We fitted Gaussian curves to the line features to measure the maximum velocities traced in absorption, the full-width half-maximum (FWHM) of the line features, and the Gaussian area of the features. We compare the proportion of each profile type in our sample to previous studies in Taurus. We find significant variations between Taurus and Lupus in the proportion of P Cygni and inverse P Cygni profiles, and between Lupus and Upper Scorpius in the number of emission-only and combination profile types. We find that the blue-shifted absorption features appear less blue-shifted at disk inclinations close to edge-on, but no such trend with inclination is observed in sources with only red-shifted features. Higher accretion rates were observed in sources with strong blue-shifted features which, along with the changes in the proportions of each profile type observed in the two regions, indicates that younger sources may drive stronger jets or winds. Overall, we observe variations in the proportion of each profile type and in the line properties which indicates and evolution of accretion and ejection signatures over time, and with source properties. These results confirm past works and models of the He I line, but for a larger sample and for multiple star-forming regions. The work highlights the power of the He I line as a probe of the gas in the innermost regions of the disk., Comment: 24 pages, 21 figures, 4 tables - accepted by A&A on 01/08/2022

The science case on studies of accretion and outflows in low-mass ($<$1.5 $M_{\odot}$) young stellar objects (YSOs) with the new CUBES instrument is presented. We show the need for a high-sensitivity, near-ultraviolet (NUV) spectrograph like CUBES, with a resolving power at least four times that of X-Shooter and combined with UVES via a fibrelink for simultaneous observations. Simulations with the CUBES exposure time calculator and the end-to-end software show that a significant gain in signal-to -noise can be achieved compared to current instruments, for both the spectral continuum and emission lines, including for relatively embedded YSOs. Our simulations also show that the low-resolution mode of CUBES will be able to observe much fainter YSOs (V $\sim$22 mag) in the NUV than we can today, allowing us extend studies to YSOs with background-limited magnitudes. The performance of CUBES in terms of sensitivity in the NUV will provide important new insights into the evolution of circumstellar disks, by studying the accretion, jets/winds and photo-evaporation processes, down to the low-mass brown dwarf regime. CUBES will also open-up new science as it will be able to observe targets that are several magnitudes fainter than those reachable with current instruments, facilitating studies of YSOs at distances of $\sim$ kpc scale. This means a step-change in the field of low-mass star formation, as it will be possible to expand the science case from relatively local star-forming regions to a large swathe of distances within the Milky Way.

We present an optical/near-IR survey of 11 variable young stars (EXors and EXor candidates) aimed at deriving and monitoring their accretion properties. About 30 optical and near-infrared spectra ($\Re$ $\sim$ 1500-2000) were collected between 2014-2019 with the Large Binocular Telescope (LBT). From the spectral analysis we have derived the accretion luminosity and mass accretion rate, the visual extinction, the temperature and density of the permitted line formation region, and the signature of the outflowing matter. Two sources (ASASSN-13db and iPTF15afq) have been observed in outburst and quiescence, three during a high-level of brightness (XZ Tau, PV Cep, and NY Ori), and the others in quiescence. These latter have accretion luminosity and mass accretion rates in line with the values measured in classical T Tauri stars of similar mass. All sources observed more than once present accretion variability. The most extreme case is ASASSN-13db, for which the mass accretion rate decreases by two orders of magnitude from the outburst peak in 2015 to quiescence in 2017. Also, in NY Ori the accretion luminosity decreases by a factor 25 in one year. In 80% of the sample we detect the [OI]630 nm line, a tracer of mass loss. From the variability of the H$\alpha$/[OI]630, ratio, we conclude that mass accretion variations are larger than mass loss variations. From the analysis of the HI recombination lines a correlation is suggested between the density of the line formation region and the level of accretion activity of the source., Comment: To appear in ApJ

The Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) Director's Discretionary Program of low-mass pre-main-sequence stars, coupled with forthcoming data from ALMA and JWST, will provide the foundation to revolutionize our understanding of the relationship between young stars and their protoplanetary disks. A comprehensive evaluation of the physics of disk evolution and planet formation requires understanding the intricate relationships between mass accretion, mass outflow, and disk structure. Here we describe the Outflows and Disks around Young Stars: Synergies for the Exploration of ULLYSES Spectra (ODYSSEUS) Survey and present initial results of the classical T Tauri Star CVSO 109 in Orion OB1b as a demonstration of the science that will result from the survey. ODYSSEUS will analyze the ULLYSES spectral database, ensuring a uniform and systematic approach in order to (1) measure how the accretion flow depends on the accretion rate and magnetic structures, (2) determine where winds and jets are launched and how mass-loss rates compare with accretion, and (3) establish the influence of FUV radiation on the chemistry of the warm inner regions of planet-forming disks. ODYSSEUS will also acquire and provide contemporaneous observations at X-ray, optical, NIR, and millimeter wavelengths to enhance the impact of the ULLYSES data. Our goal is to provide a consistent framework to accurately measure the level and evolution of mass accretion in protoplanetary disks, the properties and magnitudes of inner-disk mass loss, and the influence of UV radiation fields that determine ionization levels and drive disk chemistry., Comment: accepted to ApJ

In this era of spatially resolved observations of planet forming disks with ALMA and large ground-based telescopes such as the VLT, Keck and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPICA is an infrared space mission concept developed jointly by JAXA and ESA to address these questions. The key unique capabilities of SPICA that enable this research are (1) the wide spectral coverage 10-220 micron, (2) the high line detection sensitivity of (1-2) 10-19 W m-2 with R~2000-5000 in the far-IR (SAFARI) and 10-20 W m-2 with R~29000 in the mid-IR (SMI, spectrally resolving line profiles), (3) the high far-IR continuum sensitivity of 0.45 mJy (SAFARI), and (4) the observing efficiency for point source surveys. This paper details how mid- to far-IR infrared spectra will be unique in measuring the gas masses and water/ice content of disks and how these quantities evolve during the planet forming period. These observations will clarify the crucial transition when disks exhaust their primordial gas and further planet formation requires secondary gas produced from planetesimals. The high spectral resolution mid-IR is also unique for determining the location of the snowline dividing the rocky and icy mass reservoirs within the disk and how the divide evolves during the build-up of planetary systems. Infrared spectroscopy (mid- to far-IR) of key solid state bands is crucial for assessing whether extensive radial mixing, which is part of our Solar System history, is a general process occurring in most planetary systems and whether extrasolar planetesimals are similar to our Solar System comets/asteroids. ... (abbreviated), Comment: accepted for publication in PASA

The mass-accretion rate, Macc, is a crucial parameter for the study of the evolution of accretion disks around young low-mass stellar objects (YSOs) and for planet formation studies. The Taurus star forming region (SFR) is rich in pre-main sequence (PMS) stars, most of them of the T Tauri class. A variety of methodologies have been used in the past to measure mass accretion in samples of YSOs in Taurus, but despite being a general benchmark for star formation studies, a comprehensive and systematic analysis of the Taurus T Tauri population, where the stellar and accretion properties are derived homogeneously and simultaneously, is still missing. As part of the GIARPS High-resolution Observations of T Tauri stars (GHOsT) project, here we present a pilot study of the stellar and accretion properties of seven YSOs in Taurus using the spectrograph GIARPS at the Telescopio Nazionale Galileo (TNG). Contemporaneous low-resolution spectroscopic and photometric ancillary observations allow us to perform an accurate flux calibration of the high-resolution spectra. The simultaneity of the high-resolution, wide-band spectroscopic observations, from the optical to the near-infrared (NIR), the veiling measurements in such wide spectral range, and many well-calibrated emission line diagnostics allows us to derive the stellar and accretion properties of the seven YSOs in a homogeneous and self-consistent way. The procedures and methodologies presented here will be adopted in future works for the analysis of the complete GHOsT data set. We discuss the accretion properties of the seven YSOs in comparison with the 90\% complete sample of YSOs in the Lupus SFR and investigate possibilities for the origin of the continuum excess emission in the NIR. }

We present SOFIA/FIFI-LS observations of three Class 0 and one Class I outflows (Cep E, HH 1, HH 212, and L1551 IRS5) in the far-infrared [O I]63mum and [O I]145mum transitions. Spectroscopic [O I]63mum maps enabled us to infer the spatial extent of warm, low-excitation atomic gas within these protostellar outflows. If proper shock conditions prevail, the instantaneous mass-ejection rate is directly connected to the [O I]63mum luminosity. In order to unravel evolutionary trends, we analysed a set of 14 Class 0/I outflow sources that were spatially resolved in the [O I]63mum emission. We compared these data with a sample of 72 Class 0/I/II outflow sources that have been observed with Herschel (WISH, DIGIT, WILL, GASPS surveys) without spatially resolving the [O I]63mum line., Comment: 14 pages, 7 figures

The evolution of young stars and disks is driven by the interplay of several processes, notably accretion and ejection of material. Critical to correctly describe the conditions of planet formation, these processes are best probed spectroscopically. About five-hundred orbits of the Hubble Space Telescope (HST) are being devoted in 2020-2022 to the ULLYSES public survey of about 70 low-mass (M<2Msun) young (age<10 Myr) stars at UV wavelengths. Here we present the PENELLOPE Large Program that is being carried out at the ESO Very Large Telescope (VLT) to acquire, contemporaneous to HST, optical ESPRESSO/UVES high-resolution spectra to investigate the kinematics of the emitting gas, and UV-to-NIR X-Shooter medium-resolution flux-calibrated spectra to provide the fundamental parameters that HST data alone cannot provide, such as extinction and stellar properties. The data obtained by PENELLOPE have no proprietary time, and the fully reduced spectra are made available to the whole community. Here, we describe the data and the first scientific analysis of the accretion properties for the sample of thirteen targets located in the Orion OB1 association and in the sigma-Orionis cluster, observed in Nov-Dec 2020. We find that the accretion rates are in line with those observed previously in similarly young star-forming regions, with a variability on a timescale of days of <3. The comparison of the fits to the continuum excess emission obtained with a slab model on the X-Shooter spectra and the HST/STIS spectra shows a shortcoming in the X-Shooter estimates of <10%, well within the assumed uncertainty. Its origin can be either a wrong UV extinction curve or due to the simplicity of this modelling, and will be investigated in the course of the PENELLOPE program. The combined ULLYSES and PENELLOPE data will be key for a better understanding of the accretion/ejection mechanisms in young stars., Comment: Accepted for publication on Astronomy & Astrophysics. 15 pages + appendix, language edited version

We present SOFIA/FIFI-LS observations of five prototypical, low-mass Class I outflows (HH111, SVS13, HH26, HH34, HH30) in the far-infrared [OI]63mum and [OI]145mum transitions. The obtained spectroscopic [OI]63mum and [OI]145mum maps enable us to study the spatial extent of warm, low-excitation atomic gas within outflows driven by Class I protostars. These [OI] maps may potentially allow us to measure the mass-loss rates ($\dot{M}_\text{jet}$) of this warm component of the atomic jet., Comment: 24 pages, 41 figures, accepted for publication in A&A

The study of resolved stellar populations in the Milky Way and other Local Group galaxies can provide us with a fossil record of their chemo-dynamical and star-formation histories over timescales of many billions of years. In the galactic components and stellar systems of the Milky Way and its satellites, individual stars can be resolved. Therefore, they represent a unique laboratory in which to investigate the details of the processes behind the formation and evolution of the disc and dwarf/irregular galaxies. MOONS at the VLT represents a unique combination of an efficient infrared multi-object spectrograph and a large-aperture 8-m-class telescope which will sample the cool stellar populations of the dense central regions of the Milky Way and its satellites, delivering accurate radial velocities, metallicities, and other chemical abundances for several millions of stars over its lifetime (see Cirasuolo et al., this issue). MOONS will observe up to 1000 targets across a 25-arcminute field of view in the optical and near-infrared (0.6-1.8 micron) simultaneously. A high-resolution (R~19700) setting in the H band has been designed for the accurate determination of stellar abundances such as alpha, light, iron-peak and neutron-capture elements., Comment: 6 pages, 2 figures, Published on the ESO Messenger No. 180

In the framework of the GIARPS High-resolution Observations of T Tauri stars (GHOsT) project, we aim to characterize the atomic and molecular winds in a sample of classical T Tauri stars (CTTs) of the Taurus-Auriga region. We analyzed the flux calibrated [OI] 630 nm and $\rm H_2$ 2.12 $\rm \mu m$ lines in a sample of 36 CTTs observed at the Telescopio Nazionale Galileo with the HARPS and GIANO spectrographs. We decomposed the line profiles into different kinematic Gaussian components and focused on the most frequently detected component, the narrow low-velocity (v$\rm_p < 20$ $\rm km$ $\rm s^{-1}$) component (NLVC). We found that the $\rm H_2$ line is detected in 17 sources ($\sim 50 \%$ detection rate), and [OI] is detected in all sources but one. The NLV components of the $\rm H_2$ and [OI] emission are kinematically linked, with a strong correlation between the peak velocities and the full widths at half maximum of the two lines. Assuming Keplerian broadening, we found that the [OI] NVLC originates from a disk region between 0.05 and 20 au and that of $\rm H_2$ in a region from 2 and 20 au. We did not find any clear correlation between v$\rm_p$ of the $\rm H_2$ and [OI] NVLC and the outer disk inclination. This result is in line with previous studies. Our results suggest that molecular and neutral atomic emission in disk winds originate from regions that might overlap, and that the survival of molecular winds in disks strongly depends on the gas exposure to the radiation from the central star. Our results demonstrate the potential of wide-band high-resolution spectroscopy in linking tracers of different manifestations of the same phenomenon., Comment: 22 pages, 10 figures, 6 tables, accepted for publication in A&A

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

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., Comment: 13 pages, 9 figures

The mechanism for jet formation in the disks of T Tauri stars is poorly understood. Observational benchmarks to launching models can be provided by tracing the physical properties of the kinematic components of the wind and jet in the inner 100 au of the disk surface. In the framework of the GHOsT (GIARPS High-resolution Observations of T Tauri stars) project, we aim to perform a multi-line analysis of the velocity components of the gas in the jet acceleration zone. We analyzed the GIARPS-TNG spectra of six objects in the Taurus-Auriga complex (RY Tau, DG Tau, DL Tau, HN Tau, DO Tau, RW Aur A). Thanks to the combined high-spectral resolution (R=50000-115000) and wide spectral coverage (~400-2400 nm) we observed several O, S+, N, N+, and Fe+ forbidden lines spanning a large range of excitation and ionization conditions. In four objects (DG Tau, HN Tau, DO Tau, RW Aur A), temperature (T_e), electron and total density (n_e, n_H), and fractional ionization (x_e) were derived as a function of velocity through an excitation and ionization model. The abundance of gaseous iron, X(Fe), a probe of the dust content in the jet, was derived in selected velocity channels. The physical parameters vary smoothly with velocity, suggesting a common origin for the different kinematic components. In DG Tau and HN Tau, T_e, x_e, and X(Fe) increase with velocity (roughly from 6000 K, 0.05, 10% X(Fe)_sun to 15000 K, 0.6, 90% X(Fe)_sun). This trend is in agreement with disk-wind models in which the jet is launched from regions of the disk at different radii. In DO Tau and RW Aur A, we infer x_e < 0.1, n_H ~10^6-7 cm^-3, and X(Fe) <~ X(Fe)_sun at all velocities. These findings are tentatively explained by the formation of these jets from dense regions inside the inner, gaseous disk, or as a consequence of their high degree of collimation., Comment: Accepted by A&A

We investigate ongoing accretion activity in young stars in the TW Hydrae association (TWA, ~8-10 Myr), an ideal target to probe the final stages of disk accretion down to brown dwarf masses. Our sample comprises eleven TWA members with infrared excess, amounting to 85% of the total TWA population with disks, with spectral types between M0 and M9, and masses between 0.58 and 0.02 Msol. We employed homogeneous spectroscopic data from 300 to 2500 nm, obtained with X-shooter, to derive individual extinction, stellar parameters, and accretion parameters simultaneously. We then examined Balmer lines and forbidden emission lines to probe the physics of the star-disk interaction environment. We detected signatures of ongoing accretion for 70% of our TWA targets. This implies a fraction of accretors of 13-17% across the entire TWA (accounting for potentially accreting members not included in our survey). The spectral emission associated with these stars reveals a more evolved stage of these accretors compared to younger PMS populations: (i) a large fraction (~50%) exhibit nearly symmetric, narrow H_alpha line profiles; (ii) over 80% exhibit Balmer decrements consistent with moderate accretion and optically thin emission; (iii) less than a third exhibit forbidden line emission in [O I] 6300A, indicative of winds and outflows activity. However, the distribution in accretion rates (Mdot) derived for the TWA sample follows closely that of younger regions for Mstar~0.1-0.3 Msun. An overall correlation between Mdot and Mstar is detected, best reproduced by Mdot~Mstar^(2.1+/-0.5). At least in the lowest Mstar regimes, stars that still retain a disk at ages ~8-10 Myr are found to exhibit statistically similar, albeit moderate, accretion levels as those measured around younger objects. This slow Mdot evolution may be associated with longer evolutionary timescales of disks around low-mass stars., Comment: 27 pages (19 in the main text plus Appendices), 7 tables, 19 figures. Accepted for publication in Astronomy & Astrophysics. Abstract shortened for the arXiv listing

The impact of stellar multiplicity on the evolution of planet-forming disks is still the subject of debate. Here we present and analyze disk structures around ten multiple stellar systems that were included in an unbiased, high spatial resolution survey performed with ALMA of 32 protoplanetary disks in the Taurus star-forming region. At the unprecedented spatial resolution of ~0.12" we detect and spatially resolve the disks around all primary stars, and those around eight secondary and one tertiary star. The dust radii of disks around multiple stellar systems are smaller than those around single stars in the same stellar mass range and in the same region. The disks in multiple stellar systems also show a steeper decay of the millimeter continuum emission at the outer radius than disks around single stars, suggestive of the impact of tidal truncation on the shape of the disks in multiple systems. However, the observed ratio between the dust disk radii and the observed separation of the stars in the multiple systems is consistent with analytic predictions of the effect of tidal truncation only if the eccentricities of the binaries are rather high (typically >0.5), or if the observed dust radii are a factor of two smaller than the gas radii, as is typical for isolated systems. Similar high-resolution studies targeting the gaseous emission from disks in multiple stellar systems are required to resolve this question., Comment: Accepted on Astronomy & Astrophysics; 15 pages plus appendix, revised version after language editing

Jets are rarely associated with pre-main-sequence intermediate-mass stars. Optical and near-IR observations of jet-driving sources are often hindered by the presence of a natal envelope. Jets around partly embedded sources are a useful diagnostic to constrain the geometry of the concealed protoplanetary disk. In fact, the jet-driving mechanisms are affected by both spatial anisotropies and episodic variations at the (sub-)au scale from the star. We obtained a rich set of high-contrast VLT/SPHERE observations from 0.6 micron to 2.2 micron of the young intermediate-mass star RY Tau. Given the proximity to the Sun of this source, our images have the highest spatial resolution ever obtained for an atomic jet. Optical observations in polarized light show no sign of the protoplanetary disk detected by ALMA. Instead, we observed a diffuse signal resembling a remnant envelope with an outflow cavity. The jet is detected in four spectral lines. The jet appears to be wiggling and its radial width increasing with the distance is complementary to the shape of the outflow cavity suggesting a strong jet/envelope interaction. Through the estimated tangential velocity, we revealed a possible connection between the launching time of the jet sub-structures and the stellar activity of RY Tau. RY Tau is at an intermediate stage toward the dispersal of the natal envelope. This source shows episodic increases of mass accretion/ejection similarly to other known intermediate-mass stars. The amount of observed jet wiggle is consistent with the presence of a precessing disk warp or misaligned inner disk that would be induced by an unseen planetary/sub-stellar companion at sub-/few-au scales. The high disk mass of RY Tau and of two other jet-driving intermediate-mass stars, HD163296 and MWC480, suggests that massive, full disks are more efficient at launching prominent jets., Comment: 15 pages, 9 figures. Accepted for publication by A\&A

ALMA has observed a plethora of ring-like structures in planet forming discs at distances of 10-100 au from their host star. Although several mechanisms have been invoked to explain the origin of such rings, a common explanation is that they trace new-born planets. Under the planetary hypothesis, a natural question is how to reconcile the apparently high frequency of gap-carving planets at 10-100 au with the paucity of Jupiter mass planets observed around main sequence stars at those separations. Here, we provide an analysis of the new-born planet population emerging from observations of gaps in discs, under the assumption that the observed gaps are due to planets. We use a simple estimate of the planet mass based on the gap morphology, and apply it to a sample of gaps recently obtained by us in a survey of Taurus with ALMA. We also include additional data from recent published surveys, thus analysing the largest gap sample to date, for a total of 48 gaps. The properties of the purported planets occupy a distinctively different region of parameter space with respect to the known exo-planet population, currently not accessible through planet finding methods. Thus, no discrepancy in the mass and radius distribution of the two populations can be claimed at this stage. We show that the mass of the inferred planets conforms to the theoretically expected trend for the minimum planet mass needed to carve a dust gap. Finally, we estimate the separation and mass of the putative planets after accounting for migration and accretion, for a range of evolutionary times, finding a good match with the distribution of cold Jupiters., Comment: 9 pages, 5 figures, re-submitted to MNRAS after referee's comments

We have analysed the [OI]6300 A line in a sample of 131 young stars with discs in the Lupus, Chamaeleon and signa Orionis star forming regions, observed with the X-shooter spectrograph at VLT. The stars have mass accretion rates spanning from 10^{-12} to 10^{-7} Mo/yr. The line profile was deconvolved into a low velocity component (LVC, < 40 km/s) and a high velocity component (HVC, > 40 km/s ), originating from slow winds and high velocity jets, respectively. The LVC is by far the most frequent component, with a detection rate of 77%, while only 30% of sources have a HVC. The [OI]6300 luminosity of both the LVC and HVC, when detected, correlates with stellar and accretion parameters of the central sources (i.e. Lstar , Mstar , Lacc , Macc), with similar slopes for the two components. The line luminosity correlates better with the accretion luminosity than with the stellar luminosity or stellar mass. We suggest that accretion is the main drivers for the line excitation and that MHD disc-winds are at the origin of both components. In the sub-sample of Lupus sources observed with ALMA a relationship is found between the HVC peak velocity and the outer disc inclination angle, as expected if the HVC traces jets ejected perpendicularly to the disc plane. Mass loss rates measured from the HVC span from ~ 10^{-13} to ~10^{-7} Mo/yr. The corresponding Mloss/Macc ratio ranges from ~0.01 to ~0.5, with an average value of 0.07. However, considering the upper limits on the HVC, we infer a ratio < 0.03 in more than 40% of sources. We argue that most of these sources might lack the physical conditions needed for an efficient magneto-centrifugal acceleration in the star-disc interaction region. Systematic observations of populations of younger stars, that is, class 0/I, are needed to explore how the frequency and role of jets evolve during the pre-main sequence phase., Comment: 15 pages, 14 figures, Accepted for publication in A&A

With the purpose of performing a homogeneous determination of elemental abundances for members of the Lupus T association, we analyzed three chemical elements: lithium, iron, and barium. The aims were: to derive the Li abundance for ~90% of known class II stars in the Lupus I, II, III, IV clouds; to perform chemical tagging of a region where few Fe abundance measurements have been obtained in the past, and no determination of the Ba content has been done up to now. We also investigated possible Ba enhancement, as this element has become increasingly interesting in the last years following the evidence of Ba over-abundance in young clusters, the origin of which is still unknown. Using X-shooter@VLT, we analyzed the spectra of 89 cluster members, both class II and III stars. We measured the strength of the Li line and derived the abundance of this element through equivalent width measurements and curves of growth. For six class II stars we also measured the Fe and Ba abundances using the spectral synthesis and the code MOOG. The veiling contribution was taken into account for all three elements. We find a dispersion in the strength of the Li line at low Teff and identify three targets with severe Li depletion. The nuclear age inferred for these highly Li-depleted stars is around 15 Myr, which exceeds the isochronal one. As in other star-forming regions, no metal-rich members are found in Lupus, giving support to a recent hypothesis that the Fe abundance distribution of most of the nearby young regions could be the result of a common and widespread star formation episode involving the Galactic thin disk. We find that Ba is over-abundant by ~0.7 dex with respect to the Sun. Since current theoretical models cannot reproduce this Ba abundance pattern, we investigated whether this unusually large Ba content might be related to effects due to stellar parameters, stellar activity, and accretion., Comment: 15 pages, 14 figures, 3 tables; accepted for publication in A&A; abstract shortened

We used the TNG/GIANO instrument to obtain near-IR high-resolution spectra (R~50000) of XZ Tau and DR Tau, two actively accreting T Tauri stars classified as EXors. The analysis of the observed features provides information on the properties of the inner disk, the accretion columns, and the winds. Both sources display composite HI line profiles, with contributions from both accreting gas and high-velocity winds. These lines become more symmetric and narrower with increasing upper energy, which may be interpreted in terms of two components with different decrements or imputed to self-absorption effects. XZ Tau shows a relatively high state of activity with respect to literature observations. The variation of the HeI 1.08um line blue-shifted absorption suggests that the inner wind has undergone a dramatic change in its velocity structure, connected with a recent accretion event. DR Tau has a more stable wind as its HeI 1.08um absorption is stable in time in spite of strong variability of the emission component. The IR veiling can be interpreted as due to blackbody emission at temperatures of 1600K for XZ Tau and 2300K DR Tau, with emitting areas ~30 times larger than the star. For XZ Tau these conditions are consistent with emission from the inner rim of the dusty disk, while for DR Tau the fairly high temperature may suggest an origin from a thick gaseous disk within the dust sublimation radius. Strong and broad metallic lines, mainly from CI and FeI, are detected in XZ Tau, similar to those observed in other EXor sources during burst phases. DR Tau shows weaker and narrower metallic lines, despite its larger accretion luminosity. This suggests that accretion is not the only driver of metallic line excitation. Our observations show the potential of wide-band, high-resolution near-IR spectroscopy to simultaneously probe the phenomena that occur in the star-disk interaction region., Comment: 15 pages, accepted by A&A

Pre-stellar cores represent the initial conditions in the process of star and planet formation, therefore it is important to study their physical and chemical structure. Because of their volatility, nitrogen-bearing molecules are key to study the dense and cold gas present in pre-stellar cores. The NH_3 rotational transition detected with Herschel-HIFI provides a unique combination of sensitivity and spectral resolution to further investigate physical and chemical processes in pre-stellar cores. Here we present the velocity-resolved Herschel-HIFI observations of the ortho-NH_3(1_0-0_0) line at 572 GHz and study the abundance profile of ammonia across the pre-stellar core L1544 to test current theories of its physical and chemical structure. Recently calculated collisional coefficients have been included in our non-LTE radiative transfer code to reproduce Herschel observations. A gas-grain chemical model, including spin-state chemistry and applied to the (static) physical structure of L1544 is also used to infer the abundance profile of ortho-NH_3 . The hyperfine structure of ortho-NH_3(1_0-0_0) is resolved for the first time in space. All the hyperfine components are strongly self-absorbed. The profile can be reproduced if the core is contracting in quasi-equilibrium, consistent with previous work, and if the NH_3 abundance is slightly rising toward the core centre, as deduced from previous interferometric observations of para-NH_3(1,1). The chemical model overestimates the NH_3 abundance at radii between ~ 4000 and 15000 AU by about two orders of magnitude and underestimates the abundance toward the core centre by more than one order of magnitude. Our observations show that chemical models applied to static clouds have problems in reproducing NH_3 observations., Comment: accepted for publication in A&A Letters

After a quiescence period of about 10 years, the classical EXor source V1118 Ori has undergone an accretion outburst in 2015 September. The maximum brightness (DV > 4 mag) was reached in 2015 December and was maintained for several months. Since 2016 September, the source is in a declining phase. Photometry and low/ high-resolution spectroscopy were obtained with MODS and LUCI2 at the {\it Large Binocular Telescope}, with the facilities at the Asiago 1.22 and 1.82 m telescopes, and with GIANO at the {\it Telescopio Nazionale Galileo}. The spectra are dominated by emission lines of \hi\ and neutral metallic species. From line and continuum analysis we derive the mass accretion rate and its evolution during the outburst. Considering that extinction may vary between 1.5 and 2.9 mag, we obtain m_acc= 0.3$-$2.0 10$^{-8}$ m_sun/yr, in quiescence and m_acc= 0.2$-$1.9 10$^{-6}$ m_sun/yr, at the outburst peak. The Balmer decrement shape has been interpreted by means of line excitation models, finding that from quiescence to outburst peak, the electron density has increased from $\sim$ 2 10$^9$ cm$^{-3}$ to $\sim$ 4 10$^{11}$ cm$^{-3}$. The profiles of the metallic lines are symmetric and narrower than 100 km s$^{-1}$, while \hi\, and \hei\,\,lines show prominent wings extending up to $\pm$ 500 km s$^{-1}$. The metallic lines likely originate at the base of the accretion columns, where neutrals are efficiently shielded against the ionizing photons, while faster ionized gas is closer to the star. Outflowing activity is testified by the detection of a variable P Cyg-like profile of the H$\alpha$ and \hei\, 1.08\,$\mu$m lines., Comment: Accepted by ApJ

[Abridged] We present spectroscopic observations in H$_{2}$O, CO and related species with \textit{Herschel} HIFI and PACS, as well as ground-based follow-up with the JCMT and APEX in CO, HCO$^{+}$ and isotopologues, of a sample of 49 nearby ($d<$500\,pc) candidate protostars. These data are used to study the outflow and envelope properties of these sources. We also compile their continuum SEDs in order to constrain their physical properties. Water emission is dominated by shocks associated with the outflow, rather than the cooler, slower entrained outflowing gas probed by ground-based CO observations. These shocks become less energetic as sources evolve from Class 0 to Class I. The fraction of mass in the outflow relative to the total envelope (i.e. $M_{\mathrm{out}}/M_{\mathrm{env}}$) remains broadly constant between Class 0 and I. The median value ($\sim$1$\%$) is consistent with a core to star formation efficiency on the order of 50$\%$ and an outflow duty cycle on the order of 5$\%$. Entrainment efficiency, as probed by $F_{\mathrm{CO}}/\dot{M}_{\mathrm{acc}}$, is also invariant with source properties and evolutionary stage. The median value (6.3\kms{}) suggests an entrainment efficiency of between 30 and 60$\%$ if the wind is launched at $\sim$1AU. $L$[O\,{\sc i}] is strongly correlated with $L_{\mathrm{bol}}$ but not with $M_{\mathrm{env}}$, while low-$J$ CO is more closely correlated with the latter than the former. This suggests that [O\,{\sc i}] traces the present-day accretion activity while CO traces time-averaged accretion over the dynamical timescale of the outflow. $L$[O\,{\sc i}] does not vary from Class 0 to Class I, unlike CO and H$_{2}$O. This is likely due to the ratio of atomic to molecular gas in the wind increasing as the source evolves, balancing out the decrease in mass accretion rate. Infall signatures are detected in HCO$^{+}$ and H$_{2}$O in a few sources., Comment: Accepted to A&A, version after language editor corrections. 47 pages, 18 figures, 15 tables. The figures have been converted to pdf due to file-size considerations, which may lead to some degradation

In the BHR71 region, two low-mass protostars drive two distinguishable outflows. They constitute an ideal laboratory to investigate the effects of shock chemistry and the mechanisms that led to their formation. We aim to define the morphology of the warm gas component of the BHR 71 outflow and at modelling its shocked component. We present the first far infrared Herschel images of the BHR71 outflow in the CO(14-13), H$_2$O (2$_{21}$-1$_{10}$), H$_2$O (2$_{12}$-1$_{01}$) and [OI] 145 $\mu$m, lines, revealing the presence of several knots of warm, shocked gas associated with fast outflowing gas. In two of these knots we performed a detailed study of the physical conditions by comparing a large set of transitions from several molecules to a grid of shock models. Herschel lines ratios in the outflow knots are quite similar, showing that the excitation conditions of the fast moving gas do not change significantly within the first $\sim$ 0.068 pc of the outflow, apart at the extremity of the southern blue-shifted lobe that is expanding outside the molecular cloud. Rotational diagram, spectral line profile and LVG analysis of the CO lines in knot A show the presence of two gas components: one extended, cold ($T\sim$80 K) and dense ($n$(H$_2$) = 3$\times$10$^5$-4$\times$10$^6$ cm$^{-3}$) and another compact (18 arcsec), warm ($T$ = 1700-2200 K) with slightly lower density ($n$(H$_2$) = (2-6)$\times$10$^4$ cm$^{-3}$). In the two brightest knots (where we performed shock modelling) we found that H$_2$ and CO are well fitted with non-stationary (young) shocks. These models, however, significantly underestimate the observed fluxes of [OI] and OH lines, but are not too far off those of H$_2$O, calling for an additional, possibly dissociative, J-type shock component. Our modelling indirectly suggests that an additional shock component exists, possibly a remnant of the primary jet, Comment: 16 pages, 15 figures

We present the results of a study of the stellar and accretion properties of the (almost) complete sample of class II and transitional YSOs in the Lupus I, II, III and IV clouds, based on spectroscopic data acquired with the VLT/X-Shooter spectrograph. Our study combines the dataset from our previous work with new observations of 55 additional objects. We have investigated 92 YSO candidates in total, 11 of which have been definitely identified with giant stars unrelated to Lupus. The stellar and accretion properties of the 81 bona fide YSOs, which represent more than 90% of the whole class~II and transition disc YSO population in the aforementioned Lupus clouds, have been homogeneously and self-consistently derived, allowing for an unbiased study of accretion and its relationship with stellar parameters. The accretion luminosity, Lacc, increases with the stellar luminosity, Lstar, with an overall slope of ~1.6, similar but with a smaller scatter than in previous studies. There is a significant lack of strong accretors below Lstar~0.1Lsun, where Lacc is always lower than 0.01Lstar. We argue that the Lacc-Lstar slope is not due to observational biases, but is a true property of the Lupus YSOs. The logMacc-logMstar correlation shows a statistically significant evidence of a break, with a steeper relation for Mstar<0.2Msun and a flatter slope for higher masses. The bimodality of the Macc-Mstar relation is confirmed with four different evolutionary models used to derive the stellar mass. The bimodal behaviour of the observed relationship supports the importance of modelling self-gravity in the early evolution of the more massive discs, but other processes, such as photo evaporation and planet formation during the YSO's lifetime, may also lead to disc dispersal on different timescales depending on the stellar mass. We also refined the empirical Lacc vs. Lline relationships., Comment: 43 pages, 22 figures

Hydrogen recombination emission lines commonly observed in accreting young stellar objects represent a powerful tracer for the gas conditions in the circumstellar structures. Here we perform a study of the HI decrements and line profiles, from the Balmer and Paschen lines detected in the X-Shooter spectra of a homogeneous sample of 36 T Tauri stars in Lupus, the accretion and stellar properties of which were already derived in a previous work. We aim to obtain information on the gas physical conditions to derive a consistent picture of the HI emission mechanisms in pre-main sequence low-mass stars. We have empirically classified the sources based on their HI line profiles and decrements. We identified four Balmer decrement types (classified as 1, 2, 3, and 4) and three Paschen decrement types (A, B, and C), characterised by different shapes. We first discussed the connection between the decrement types and the source properties and then compared the observed decrements with predictions from recently published local line excitation models. One third of the objects show lines with narrow symmetric profiles, and present similar Balmer and Paschen decrements (straight decrements, types 2 and A). Lines in these sources are consistent with optically thin emission from gas with hydrogen densities of order 10^9 cm^-3 and 5000 11 cm^-3). Type 1 (curved) Balmer decrements are observed only in three sub-luminous sources viewed edge-on, so we speculate that these are actually reddened type 2 decrements. About 20% of the objects present type 3 Balmer decrements (bumpy), which cannot be reproduced with current models., Comment: 29 pages, accepted by A&A

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., Comment: Letter to Astronomy and Astrophysics

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., Comment: 18 pages, 10 figure, Accepted on Astronomy & Astrophysics

In a previous study of the L1157 B1 shocked cavity, a comparison between NH$_3$(1$_0$-$0_0$) and H$_2$O(1$_{\rm 10}$--1$_{\rm 01}$) transitions showed a striking difference in the profiles, with H$_2$O emitting at definitely higher velocities. This behaviour was explained as a result of the high-temperature gas-phase chemistry occurring in the postshock gas in the B1 cavity of this outflow. If the differences in behaviour between ammonia and water are indeed a consequence of the high gas temperatures reached during the passage of a shock, then one should find such differences to be ubiquitous among chemically rich outflows. In order to determine whether the difference in profiles observed between NH$_3$ and H$_2$O is unique to L1157 or a common characteristic of chemically rich outflows, we have performed Herschel-HIFI observations of the NH$_3$(1$_0$-0$_0$) line at 572.5 GHz in a sample of 8 bright low-mass outflow spots already observed in the H$_2$O(1$_{\rm 10}$--1$_{\rm 01}$) line within the WISH KP. We detected the ammonia emission at high-velocities at most of the outflows positions. In all cases, the water emission reaches higher velocities than NH$_3$, proving that this behaviour is not exclusive of the L1157-B1 position. Comparisons with a gas-grain chemical and shock model confirms, for this larger sample, that the behaviour of ammonia is determined principally by the temperature of the gas., Comment: Accepted for publication in the Monthly Notices of the Royal Astronomical Society

Hydrides are simple compounds containing one or a few hydrogen atoms bonded to a heavier atom. They are fundamental precursor molecules in cosmic chemistry and many hydride ions have become observable in high quality for the first time thanks to the Herschel Space Observatory. Ionized hydrides, such as CH+ and OH+, and also HCO+ that affect the chemistry of molecules such as water, provide complementary information on irradiation by far UV (FUV) or X-rays and gas temperature. The targeted lines of CH+, OH+, H2O+, C+ and CH are detected mostly in blue-shifted absorption. H3O+ and SH+ are detected in emission and only toward some high-mass objects. The observed line parameters and correlations suggest two different origins, related to gas entrained by the outflows and to the circumstellar envelope. The column density ratios of CH+/OH+ are estimated from chemical slab models, assuming that the H2 density is given by the specific density model of each object at the beam radius. For the low-mass YSOs the observed ratio can be reproduced for an FUV flux of 2-400 times the ISRF at the location of the molecules. In two high-mass objects, the UV flux is 20-200 times the ISRF derived from absorption lines, and 300-600 ISRF using emission lines. If the FUV flux required for low-mass objects originates at the central protostar, a substantial FUV luminosity, up to 1.5 L_sun, is required. There is no molecular evidence for X-ray induced chemistry in the low-mass objects on the observed scales of a few 1000 AU. For high-mass regions, the FUV flux required to produce the observed molecular ratios is smaller than the unattenuated flux expected from the central object(s) at the Herschel beam radius. This is consistent with an FUV flux reduced by circumstellar extinction or by bloating of the protostar., Comment: 40 pages, 97 figures, accepted by Astronomy and Astrophysics

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%., Comment: Accepted by The Astrophysical Journal

Context: High-resolution studies of class 0 protostars represent the key to constraining protostar formation models. VLA16234-2417 represents the prototype of class 0 protostars, and it has been recently identified as a triple non-coeval system. Aim: We aim at deriving the physical properties of the jets in VLA16234-2417 using tracers of shocked gas. Methods: ALMA Cycle 0 Early Science observations of CO(2-1) in the extended configuration are presented in comparison with previous SMA CO(3-2) and Herschel-PACS [OI}] 63 micron observations. Gas morphology and kinematics were analysed to constrain the physical structure and origin of the protostellar outflows. Results: We reveal a collimated jet component associated with the [OI] 63 micron emission at about 8'' (about 960 AU) from source B. This newly detected jet component is inversely oriented with respect to the large-scale outflow driven by source A, and it is aligned with compact and fast jet emission very close to source B (about 0.3'') rather than with the direction perpendicular to the A disk. We also detect a cavity-like structure at low projected velocities, which surrounds the [OI] 63 micron emission and is possibly associated with the outflow driven by source A. Finally, no compact outflow emission is associated with source W. Conclusions: Our high-resolution ALMA observations seem to suggest there is a fast and collimated jet component associated with source B. This scenario would confirm that source B is younger than A, that it is in a very early stage of evolution, and that it drives a faster, more collimated, and more compact jet with respect to the large-scale slower outflow driven by A. However, a different scenario of a precessing jet driven by A cannot be firmly excluded from the present observations., Comment: Accepted for publication in Astronomy & Astrophysics

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}., Comment: Editors: S. Antoniucci, J.M. Alcala', C. Codella, B. Nisini

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., Comment: 13 pages, 10 figures, accepted for publication by A&A

V1118 Ori is an eruptive variable belonging to the EXor class of Pre-Main Sequence stars whose episodic outbursts are attributed to disk accretion events. Since 2006, V1118 Ori is in the longest quiescence stage ever observed between two subsequent outbursts of its recent history. We present near-infrared photometry of V1118 Ori carried out during the last eight years, along with a complete spectroscopic coverage from 0.35 to 2.5 um. A longterm sampling of V1118 Ori in quiescence has never been done, hence we can benefit from the current circumstance to determine the lowest values (i.e. the zeroes) of the parameters to be used as a reference for evaluating the physical changes typical of more active phases. A quiescence mass accretion rate between 1--3 $\times$ 10$^{-9}$ M$_{\sun}$ yr$^{-1}$ can be derived and the difference with previous determinations is discussed. From line emission and IR colors analysis a visual extinction of 1-2 mag is consistently derived, confirming that V1118 Ori (at least in quiescence) is a low-extinction T Tauri star with a bolometric luminosity of about 2.1 L$_{\sun}$. An anti-correlation exists between the equivalent width of the emission lines and the underlying continuum. We searched the literature for evaluating whether or not such a behaviour is a common feature of the whole class. The anti-correlation is clearly recognizable for all the available EXors in the optical range (H$\beta$ and H$\alpha$ lines), while it is not as much evident in the infrared (Pa$\beta$ and Br$\gamma$ lines). The observed anti-correlation supports the accretion-driven mechanism as the most likely to account for continuum variations., Comment: 6 figures, 5 tables, accepted on ApJ

We present Herschel PACS mapping observations of the [OI]63 micron line towards protostellar outflows in the L1448, NGC1333-IRAS4, HH46, BHR71 and VLA1623 star forming regions. We detect emission spatially resolved along the outflow direction, which can be associated with a low excitation atomic jet. In the L1448-C, HH46 IRS and BHR71 IRS1 outflows this emission is kinematically resolved into blue- and red-shifted jet lobes, having radial velocities up to 200 km/s. In the L1448-C atomic jet the velocity increases with the distance from the protostar, similarly to what observed in the SiO jet associated with this source. This suggests that [OI] and molecular gas are kinematically connected and that this latter could represent the colder cocoon of a jet at higher excitation. Mass flux rates (\.M$_{jet}$(OI)) have been measured from the [OI]63micron luminosity adopting two independent methods. We find values in the range 1-4 10$^{-7}$ Mo/yr for all sources but HH46, for which an order of magnitude higher value is estimated. \.M$_{jet}$(OI) are compared with mass accretion rates (\.M$_{acc}$) onto the protostar and with \.M$_{jet}$ derived from ground-based CO observations. \.M$_{jet}$(OI)/\.M$_{acc}$ ratios are in the range 0.05-0.5, similar to the values for more evolved sources. \.M$_{jet}$(OI) in HH46 IRS and IRAS4A are comparable to \.M$_{jet}$(CO), while those of the remaining sources are significantly lower than the corresponding \.M$_{jet}$(CO). We speculate that for these three sources most of the mass flux is carried out by a molecular jet, while the warm atomic gas does not significantly contribute to the dynamics of the system., Comment: 37 pages and 12 figures, accepted for publication on Astrophysical Journal

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., Comment: 11 pages, 10 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

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.

As part of the POISSON project (Protostellar Optical-Infrared Spectral Survey on NTT), we present the results of the analysis of low-resolution NIR spectra 0.9-2.4 um) of two samples of YSOs in Lupus and Serpens (52 and 17 objects), with masses 0.1-2.0 Msun and ages from 10^5 to a few 10^7 yr. After determining the accretion parameters of the Lup and Ser targets by analysing their HI near-IR emission features, we added the results to those from previous regions (investigated in POISSON with the same methodology). We obtained a final catalogue (143 objects) of mass accretion rates (Macc) derived in a homogeneous fashion and analysed how Macc correlates with M* and how it evolves in time. We derived the accretion luminosity (Lacc) and Macc for Lup and Ser objects from the Br_gamma line by using relevant empirical relationships from the literature that connect HI line luminosity and Lacc. To minimise the biases and also for self-consistency, we re-derived mass and age for each source using the same set of evolutionary tracks. We observe a correlation MaccM*^2.2, similarly to what has previously been observed in several star-forming clouds. The time variation of Macc is roughly consistent with the expected evolution in viscous disks, with an asymptotic decay that behaves as t^-1.6. However, Macc values are characterised by a large scatter at similar ages and are on average higher than the predictions of viscous models. Although part of the scattering may be related to the employed empirical relationship and to uncertainties on the single measurements, the general distribution and decay trend of the Macc points are real. These findings might be indicative of a large variation in the initial mass of the disks, of fairly different viscous laws among disks, of varying accretion regimes, and of other mechanisms that add to the dissipation of the disks, such as photo-evaporation., Comment: 18 pages, 10 figures, accepted by A&A

The HH54 shock is a Herbig-Haro object, located in the nearby Chamaeleon II cloud. Observed CO line profiles are due to a complex distribution in density, temperature, velocity, and geometry. Resolving the HH54 shock wave in the far-infrared cooling lines of CO constrain the kinematics, morphology, and physical conditions of the shocked region. We used the PACS and SPIRE instruments on board the Herschel space observatory to map the full FIR spectrum in a region covering the HH54 shock wave. Complementary Herschel-HIFI, APEX, and Spitzer data are used in the analysis as well. The observed features in the line profiles are reproduced using a 3D radiative transfer model of a bow-shock, constructed with the Line Modeling Engine code (LIME). The FIR emission is confined to the HH54 region and a coherent displacement of the location of the emission maximum of CO with increasing J is observed. The peak positions of the high-J CO lines are shifted upstream from the lower J CO lines and coincide with the position of the spectral feature identified previously in CO(10-9) profiles with HIFI. This indicates a hotter molecular component in the upstream gas with distinct dynamics. The coherent displacement with increasing J for CO is consistent with a scenario where IRAS12500-7658 is the exciting source of the flow, and the 180 K bow-shock is accompanied by a hot (800 K) molecular component located upstream from the apex of the shock and blueshifted by -7 km s$^{-1}$. The spatial proximity of this knot to the peaks of the atomic fine-structure emission lines observed with Spitzer and PACS ([OI]63, 145 $\mu$m) suggests that it may be associated with the dissociative shock as the jet impacts slower moving gas in the HH54 bow-shock., Comment: 6 pages, 5 figures

Context: Outflows are an important part of the star formation process as both the result of ongoing active accretion and one of the main sources of mechanical feedback on small scales. Water is the ideal tracer of these effects because it is present in high abundance in various parts of the protostar. Method: We present \textit{Herschel} HIFI spectra of multiple water-transitions towards 29 nearby Class 0/I protostars as part of the WISH Survey. These are decomposed into different Gaussian components, with each related to one of three parts of the protostellar system; quiescent envelope, cavity shock and spot shocks in the jet and at the base of the outflow. We then constrain the excitation conditions present in the two outflow-related components. Results: Water emission is optically thick but effectively thin, with line ratios that do not vary with velocity, in contrast to CO. The physical conditions of the cavity and spot shocks are similar, with post-shock H$_{2}$ densities of order 10$^{5}-$10$^{8}$\,cm$^{-3}$ and H$_{2}$O column densities of order 10$^{16}-$10$^{18}$\,cm$^{-2}$. H$_{2}$O emission originates in compact emitting regions: for the spot shocks these correspond to point sources with radii of order 10-200\,AU, while for the cavity shocks these come from a thin layer along the outflow cavity wall with thickness of order 1-30\,AU. Conclusions: Water emission at the source position traces two distinct kinematic components in the outflow; J shocks at the base of the outflow or in the jet, and C shocks in a thin layer in the cavity wall. Class I sources have similar excitation conditions to Class 0 sources, but generally smaller line-widths and emitting region sizes. We suggest that it is the velocity of the wind driving the outflow, rather than the decrease in envelope density or mass, that is the cause of the decrease in H$_{2}$O intensity between Class 0 and I., Comment: Abstract is abridged. Accepted to Astronomy & Astrophysics. 50 pages (21 excluding appendices), 15 figures and 7 tables in the main paper, 18 figures and 10 tables in the appendices. Table 4 updated to Corrigendum version