Search

We present the overview and first results from the North-PHASE Legacy Survey, which follows six young clusters for five years, using the 2 deg$^2$ FoV of the JAST80 telescope from the Javalambre Observatory (Spain). North-PHASE investigates stellar variability on timescales from days to years for thousands of young stars distributed over entire clusters. This allows us to find new YSO, characterise accretion and study inner disk evolution within the cluster context. Each region (Tr37, CepOB3, IC5070, IC348, NGC2264, and NGC1333) is observed in six filters (SDSS griz, u band, and J0660, which covers H$\alpha$), detecting cluster members as well as field variable stars. Tr37 is used to prove feasibility and optimise the variability analysis techniques. In Tr37, variability reveals 50 new YSO, most of them proper motion outliers. North-PHASE independently confirms the youth of astrometric members, efficiently distinguishes accreting and non-accreting stars, reveals the extent of the cluster populations along Tr37/IC1396 bright rims, and detects variability resulting from rotation, dips, and irregular bursts. The proper motion outliers unveil a more complex star formation history than inferred from Gaia alone, and variability highlights previously hidden proper motion deviations in the surrounding clouds. We also find that non-YSO variables identified by North-PHASE cover a different variability parameter space and include long-period variables, eclipsing binaries, RR Lyr, and $\delta$ Scuti stars. These early results also emphasize the power of variability to complete the picture of star formation where it is missed by astrometry., Comment: Accepted by MNRAS

(Abridged) The presence of short-period (< 10 days) planets around main sequence (MS) stars has been associated either with the dust-destruction region or with the magnetospheric gas-truncation radius in the protoplanetary disks that surround them during the pre-MS phase. However, previous analyses have only considered low-mass FGK stars, making it difficult to disentangle the two scenarios. This exploratory study is aimed at testing whether it is the inner dust or gas disk driving the location of short-period, giant planets. By combining TESS and Gaia DR3 data, we identified a sample of 47 intermediate-mass (1.5-3 M$_{\odot}$) MS stars hosting confirmed and firm candidate hot Jupiters. We compared their orbits with the rough position of the inner dust and gas disks, which are well separated around their Herbig stars precursors. We also made a comparison with the orbits of confirmed hot Jupiters around a similarly extracted TESS/Gaia sample of low-mass sources (0.5-1.5 M$_{\odot}$). Our results suggest that the inner gas (and not the dust) disk limits the innermost orbits of hot Jupiters around intermediate-mass stars. These findings also provide tentative support to previous works that have claimed this is indeed the case for low-mass sources. We propose that hot Jupiters could be explained via a combination of the core-accretion paradigm and migration up to the gas-truncation radius, which may be responsible for halting inward migration regardless of the stellar mass regime. Larger samples of intermediate-mass stars with hot Jupiters are necessary to confirm our hypothesis, which implies that massive Herbig stars without magnetospheres (> 3-4 M$_{\odot}$) may be the most efficient in swallowing their newborn planets., Comment: Letter accepted in Astronomy & Astrophysics. 9 figures, 1 table, 3 appendixes

Chance-aligned sources or blended companions can cause false positives in planetary transit detections or simply bias the determination of the candidate properties. In the era of high-precision space-based photometers, the need for high-spatial resolution images has demonstrated to be critical for validating and confirming transit signals. This already applied to the Kepler mission, it is now applicable to the TESS survey and will be critical for PLATO. We present the results of the AstraLux-TESS survey, a catalog of high-spatial resolution images obtained with the AstraLux instrument (Calar Alto) in the context of the TESS Follow-up Observing Program. We use the lucky-imaging technique to obtain high-spatial resolution images from planet candidate hosts included mostly in two relevant regimes: exoplanet candidates belonging to the level-one requirement of the TESS mission (planets with radii $R<4~R_{\oplus}$), and candidates around intermediate-mass stars. Among the 185 planet host candidate stars observed, we found 13 (7%) to be accompanied by additional sources within 2.2 arcsec separation. Among them, six are not associated to sources in the Gaia DR3 catalog, thus contaminating the TESS light curve. We provide upper limits and probabilities to the possible existence of field contaminants through the sensitivity limits of our images. Among the isolated hosts, we can discard hazardous companions (bright enough to mimic a planetary transit signals) for all their planets. The results from this catalog are key for the statistical validation of small planets (prime targets of the TESS mission) and planets around intermediate-mass stars in the main-sequence. These two populations of planets are hard to confirm with the radial velocity technique. Our results also demonstrate the importance of this type of follow-up observations for future transit missions like PLATO, even in the Gaia era., Comment: Accepted for publication in A&A. 10 pages, 5 figures, 4 tables. Language revision included

(Abridged) The stellar metallicity, [M/H], may have important implications for planet formation. In particular, Kama et al. proposed that the deficit of refractory elements in the surfaces of some Herbig Ae/Be stars (HAeBes) may be linked to the presence of disk cavities likely caused by Jovian planets that trap the metal-rich content. This work aims to provide a robust test on the previous proposal by analyzing the largest sample of HAeBes with homogeneously derived [M/H] values, stellar, and circumstellar properties. [M/H] values of 67 HAeBes were derived based on observed spectra and Kurucz synthetic models. Statistical analyses were carried out aiming to test the potential relation between [M/H] and the group I sources from the spectral energy distribution (SED) classification by Meeus et al., associated to the presence of cavities potentially carved by giant planets. Our study robustly confirms that group I sources tend to have a lower [M/H] than that of group II HAeBes. A similar analysis involving SED-based transitional disks does not reveal such a relation, indicating that not all processes capable of creating dust holes have an effect on the stellar abundances. We also show that the observed [M/H] differences are not driven by environmental effects. Finally, group I sources tend to have stronger (sub-) mm continuum emission presumably related to the presence of giant planets. Indeed, literature results indicate that disk substructures probably associated to their presence are up to ten times more frequent in group I HAeBes. We provide indirect evidences suggesting that giant planets are more frequent around group I/low [M/H] stars than around the rest of the HAeBes. However, the direct test requires multiple detections of forming planets in their disks, so far limited to the candidate around the metal depleted ([M/H] = -0.35 +- 0.25) group I HAeBe star AB Aur., Comment: 16 pages, 10 figures, 2 tables. Accepted in A&A

(Pre-)transitional disks show gaps and cavities that can be related with on-going planet formation. According to theory, young embedded planets can accrete material from the circumplanetary and circumstellar disks, so that they could be detected in accretion tracers, like the H$_{\alpha}$ emission line. In this work, we present spectral angular differential imaging AO-assisted observations of five (pre-)transitional disks obtained with SPHERE/ZIMPOL at the Very Large Telescope (VLT). They were obtained in the H$_{\alpha}$ line and the adjacent continuum. We have combined spectral and angular differential imaging techniques to increase the contrast in the innermost regions close to the star, and search for the signature of young accreting protoplanets. As a result, the reduced images do not show any clear H$_{\alpha}$ point source around any of the targets. We report faint H$_{\alpha}$ emissions around TW Hya and HD163296: while the former is most probably an artifact related with a spike, the nature of the latter remains unclear. The spectral and angular differential images yield contrasts of 6--8 magnitudes at separations of $\sim$ 100 mas from the central stars, except in the case of LkCa15, with values of $\sim$3 mag. We have estimated upper limits to the accretion luminosity of potential protoplanets, obtaining that planetary models provide an average value of $L_{\rm acc} \sim 10^{-4}$ $L_{\odot}$ at 200 mas, which is $\sim$2 orders of magnitude higher than the $L_{\rm acc}$ estimated from the extrapolation of the $L_{H_{\alpha}}$ - $L_{acc}$ stellar relationship. We explain the lack of protoplanet detections as a combination of different factors, like e.g. episodic accretion, extinction from the circumstellar and circumplanetray disks, and/or a majority of low-mass, low-accreting planets., Comment: 17 pages, Accepted for publication in A&A

(Abridged) The lifetime of protoplanetary disks around young stars limits the timescale when planets form. A disk dissipation timescale < 10 Myr was inferred from surveys providing the fraction of stars with disks in young stellar clusters with different ages. However, most previous surveys focused on the compact region within ~ 2 pc from the clusters' centers, for which the disk fraction information considering the outer part is practically absent. We aim to test if disk fraction estimates change when inferred from an extended region around the clusters' centers. Gaia EDR3 data and a best-suited, Virtual Observatory (VO)-based tool -Clusterix-, are used to identify member stars for a representative sample of 19 young clusters considering two concentric fields of view (FOV) with radii ~ 20 pc and ~ 2 pc. Our analysis reveals that the inner disk fractions inferred from the compact and the extended regions are equal within ~ 10%, which does not support a previous hypothesis proposing that disk fractions should be significantly larger considering extended regions. A list of member and disk stars in each cluster is provided and stored in a VO-compliant archive. Averaged values and plots characterizing the whole clusters are also provided, including HR diagrams based on Gaia colors and absolute magnitudes. Our results cover the largest fields ever probed when dealing with disk fractions for all clusters analysed, and imply that their complete characterization requires the use of wide FOVs. The resulting database is a benchmark for future detailed studies of young clusters, whose disk fractions must be accurately determined by using multi-wavelength analysis potentially combined with data from coming Gaia releases., Comment: Accepted in Astronomy & Astrophysics. 3 Appendixes. Tables and figures also available online: http://svocats.cab.inta-csic.es/diskfrac

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

Circumstellar discs are essential for high mass star formation, while multiplicity, in particular binarity, appears to be an inevitable outcome since the vast majority of massive stars (> 8 Msun) are found in binaries (up to 100%). We constrain the sizes of the dust and ionised gas (Brgamma) emission of the innermost regions towards a sample of six MYSOs, and provide high-mass binary statistics of young stars at 2-300 au scales using VLTI (GRAVITY, AMBER) observations. We determine the inner radius of the dust emission and place MYSOs with K-band measurements in a size-luminosity diagram for the first time, and compare our findings to T Tauris and Herbig AeBes. We also compare the observed K-band sizes to the sublimation radius predicted by three different disc scenarios. Lastly, we apply binary geometries to trace close binarity among MYSOs. The inner sizes of MYSOs, Herbig AeBe and T Tauri stars appear to follow a universal trend at which the sizes scale with the square-root of the stellar luminosity. The Brgamma emission originates from somewhat smaller and co-planar area compared to the 2.2 {\mu}m continuum emission. We discuss this new finding with respect to disc-wind or jet origin. Finally, we report an MYSO binary fraction of 17-25% at milli-arcsecond separations (2-300 au). The size-luminosity diagram indicates that the inner regions of discs around young stars scale with luminosity independently of the stellar mass. At the targeted scales (2-300 au), the MYSO binary fraction is lower than what was previously reported for the more evolved main sequence massive stars, which, if further confirmed, could implicate the predictions from massive binary formation theories. Lastly, we spatially resolve the crucial star/disc interface in a sample of MYSOs, showing that au-scale discs are prominent in high-mass star formation and similar to their low-mass equivalents., Comment: Accepted for publication in A&A. 26 pages, 16 Figures, 7 Tables, full abstract to appear in A&A

(Abridged:) It has been hypothesized that the location of Herbig Ae/Be stars (HAeBes) within the empirical relation between the inner disk radius (r$_{in}$), inferred from K-band interferometry, and the stellar luminosity (L$_*$), is related to the presence of the innermost gas, the disk-to-star accretion mechanism, the dust disk properties inferred from the spectral energy distributions (SEDs), or a combination of these effects. This work aims to test whether the previously proposed hypotheses do, in fact, serve as a general explanation for the distribution of HAeBes in the size-luminosity diagram. GRAVITY/VLTI spectro-interferometric observations at $\sim $2.2 $ \mu$m have been obtained for five HBes representing two extreme cases concerning the presence of innermost gas and accretion modes. V590 Mon, PDS 281, and HD 94509 show no excess in the near-ultraviolet, Balmer region of the spectra ($\Delta$D$_B$), indicative of a negligible amount of inner gas and disk-to-star accretion, whereas DG Cir and HD 141926 show such strong $\Delta$D$_B$ values that cannot be reproduced from magnetospheric accretion, but probably come from the alternative boundary layer mechanism. Additional data for these and all HAeBes resolved through K-band interferometry have been compiled from the literature and updated using Gaia EDR3 distances, almost doubling previous samples used to analyze the size-luminosity relation. We find no general trend linking the presence of gas inside the dust destruction radius or the accretion mechanism with the location of HAeBes in the size-luminosity diagram. Underlying trends are present and must be taken into account when interpreting the size-luminosity correlation. Still, it is argued that the size-luminosity correlation is most likely to be physically relevant in spite of the previous statistical warning concerning dependencies on distance., Comment: 20 pages, 17 figures. The main data reduction and analysis process ia available at https://github.com/marcosarenal/gravity-data-reduction

Context. Studying gas chemistry in protoplanetary disks is key to understanding the process of planet formation. Sulfur chemistry in particular is poorly understood in interstellar environments, and the location of the main reservoirs remains unknown. Protoplanetary disks in Taurus are ideal targets for studying the evolution of the composition of planet forming systems. Aims. We aim to elucidate the chemical origin of sulfur-bearing molecular emission in protoplanetary disks, with a special focus on H$_2$S emission, and to identify candidate species that could become the main molecular sulfur reservoirs in protoplanetary systems. Methods. We used IRAM 30m observations of nine gas-rich young stellar objects (YSOs) in Taurus to perform a survey of sulfur-bearing and oxygen-bearing molecular species. In this paper we present our results for the CS 3-2 ($\nu_0$ = 146.969 GHz), H$_2$CO 2$_{11}$-1$_{10}$ ($\nu_0$ = 150.498 GHz), and H$_2$S 1$_{10}$-1$_{01}$ ($\nu_0$ = 168,763 GHz) emission lines. Results. We detected H$_2$S emission in four sources out of the nine observed, significantly increasing the number of detections toward YSOs. We also detected H$_2$CO and CS in six out of the nine. We identify a tentative correlation between H$_2$S 1$_{10}$-1$_{01}$ and H$_2$CO 2$_{11}$-1$_{10}$ as well as a tentative correlation between H$_2$S 1$_{10}$-1$_{01}$ and H$_2$O 8$_{18}$-7$_{07}$. By assuming local thermodynamical equilibrium, we computed column densities for the sources in the sample, with N(o-H$_2$S) values ranging between $2.6\times10^{12}$ cm$^{-2}$ and $1.5\times10^{13}$ cm$^{-2}$., Comment: 13 pages, 6 figures, 6 tables

(Abridged) Herbig Ae/Be stars (HAeBes) have so far been studied based on relatively small samples that are scattered throughout the sky. Their fundamental stellar and circumstellar parameters and statistical properties were derived with heterogeneous approaches before Gaia. Our main goal is to contribute to the study of HAeBes from the largest sample of such sources to date, for which stellar and circumstellar properties have been determined homogeneously from the analysis of the spectral energy distributions (SEDs) and Gaia EDR3 parallaxes and photometry. Multiwavelength photometry was compiled for 209 bona fide HAeBes for which Gaia EDR3 distances were estimated. Using the Virtual Observatory SED Analyser (VOSA), photospheric models were fit to the optical SEDs to derive stellar parameters, and the excesses at infrared (IR) and longer wavelengths were characterized to derive several circumstellar properties. The stellar temperature, luminosity, radius, mass, and age were derived for each star based on optical photometry. In addition, their IR SEDs were classified according to two different schemes, and their mass accretion rates, disk masses, and the sizes of the inner dust holes were also estimated uniformly. The initial mass function fits the stellar mass distribution of the sample within 2 < M_star/M_sun < 12. In this aspect, the sample is therefore representative of the HAeBe regime and can be used for statistical purposes when it is taken into account that the boundaries are not well probed. A statistical analysis was carried out mainly focused on the different circumstellar properties of Herbig Ae (HAe) and Be (HBe) stars. The complete dataset is available online through a Virtual Observatory-compliant archive, representing the most recent reference for statistical studies on the HAeBe regime., Comment: 45 pages, 13 figures, 3 tables. Accepted in A&A

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

Young accreting stars drive outflows that collimate into jets, which can be seen hundreds of au from their driving sources. Accretion and outflow activity cease with system age, and it is believed that magneto-centrifugally launched disk winds are critical agents in regulating accretion through the protoplanetary disk. Protostellar jets are well studied in classical T~Tauri stars ($M_\star\lesssim2\,M_\odot$), while few nearby ($d\lesssim150\,$pc) intermediate-mass stars ($M_\star=2-10\,M_\odot$), known as Herbig Ae/Be stars, have detected jets. We report VLT/MUSE observations of the Herbig~Ae/Be star HD~100546 and the discovery of a protostellar jet. The jet is similar in appearance to jets driven by low-mass stars and compares well with the jet of HD~163296, the only other known optical jet from a nearby Herbig~Ae/Be star. We derive a (one-sided) mass-loss rate in the jet of $\log \dot{M}_{jet} \sim -9.5$ (in $M_\odot$\,yr$^{-1}$) and a ratio of outflow to accretion of roughly $3\times10^{-3}$, which is lower than that of CTTS jets. The discovery of the HD~100546 jet is particularly interesting because the protoplanetary disk around HD~100546 shows a large radial gap, spiral structure, and might host a protoplanetary system. A bar-like structure previously seen in H$\alpha$ with VLT/SPHERE shares the jet position angle, likely represents the base of the jet, and suggests a jet-launching region within about 2\,au. We conclude that the evolution of the disk at radii beyond a few au does not affect the ability of the system to launch jets., Comment: 6 pages, 6 figures, accepted for publication in A&A letters

The intermediate-mass pre-main sequence Herbig Ae/Be stars are key to understanding the differences in formation mechanisms between low- and high-mass stars. The study of the general properties of these objects is hampered by the fact that few and mostly serendipitously discovered sources are known. Our goal is to identify new Herbig Ae/Be candidates to create a homogeneous and well defined catalogue of these objects. We have applied machine learning techniques to 4,150,983 sources with data from Gaia DR2, 2MASS, WISE, and IPHAS or VPHAS+. Several observables were chosen to identify new Herbig Ae/Be candidates based on our current knowledge of this class, which is characterised by infrared excesses, photometric variabilities, and H$\alpha$ emission lines. Classical techniques are not efficient for identifying new Herbig Ae/Be stars mainly because of their similarity with classical Be stars, with which they share many characteristics. By focusing on disentangling these two types of objects, our algorithm has also identified new classical Be stars. We have obtained a large catalogue of 8470 new pre-main sequence candidates and another catalogue of 693 new classical Be candidates with a completeness of $78.8\pm1.4\%$ and $85.5\pm1.2\%$, respectively. Of the catalogue of pre-main sequence candidates, at least 1361 sources are potentially new Herbig Ae/Be candidates according to their position in the Hertzsprung-Russell diagram. In this study we present the methodology used, evaluate the quality of the catalogues, and perform an analysis of their flaws and biases. For this assessment, we make use of observables that have not been accounted for by the algorithm and hence are selection-independent, such as coordinates and parallax based distances. The catalogue of new Herbig Ae/Be stars that we present here increases the number of known objects of the class by an order of magnitude., Comment: Accepted for publication in Astronomy & Astrophysics. Received 14 February 2020 / Accepted 27 April 2020 (18 pages, 9 figures)

Understanding how young stars gain their masses through disk-to-star accretion is of paramount importance in astrophysics. It affects our knowledge about the early stellar evolution, the disk lifetime and dissipation processes, the way the planets form on the smallest scales, or the connection to macroscopic parameters characterizing star-forming regions on the largest ones, among others. In turn, mass accretion rate estimates depend on the accretion paradigm assumed. For low-mass T Tauri stars with strong magnetic fields there is consensus that magnetospheric accretion (MA) is the driving mechanism, but the transfer of mass in massive young stellar objects with weak or negligible magnetic fields probably occurs directly from the disk to the star through a hot boundary layer (BL). The intermediate-mass Herbig Ae/Be (HAeBe) stars bridge the gap between both previous regimes and are still optically visible during the pre-main sequence phase, thus constituting a unique opportunity to test a possible change of accretion mode from MA to BL. This review deals with our estimates of accretion rates in HAeBes, critically discussing the different accretion paradigms. It shows that although mounting evidence supports that MA may extend to late-type HAes but not to early-type HBes, there is not yet a consensus on the validity of this scenario versus the BL one. Based on MA and BL shock modeling, it is argued that the ultraviolet regime could significantly contribute in the future to discriminating between these competing accretion scenarios., Comment: 26 pages, 5 figures, 1 table. Invited Review accepted for publication in "Galaxies" special issue: "Star Formation in the UV", ed. Jorick Vink

This work presents a spectroscopic study of 163 Herbig Ae/Be stars. Amongst these, we present new data for 30 objects. Stellar parameters such as temperature, reddening, mass, luminosity and age are homogeneously determined. Mass accretion rates are determined from $\rm H\alpha$ emission line measurements. Our data is complemented with the X-Shooter sample from previous studies and we update results using Gaia DR2 parallaxes giving a total of 78 objects with homogeneously determined stellar parameters and mass accretion rates. In addition, mass accretion rates of an additional 85 HAeBes are determined. We confirm previous findings that the mass accretion rate increases as a function of stellar mass, and the existence of a different slope for lower and higher mass stars respectively. The mass where the slope changes is determined to be $3.98^{+1.37}_{-0.94}\,\rm M_{\odot}$. We discuss this break in the context of different modes of disk accretion for low- and high mass stars. Because of their similarities with T Tauri stars, we identify the accretion mechanism for the late-type Herbig stars with the Magnetospheric Accretion. The possibilities for the earlier-type stars are still open, we suggest the Boundary Layer accretion model may be a viable alternative. Finally, we investigated the mass accretion - age relationship. Even using the superior Gaia based data, it proved hard to select a large enough sub-sample to remove the mass dependency in this relationship. Yet, it would appear that the mass accretion does decline with age as expected from basic theoretical considerations.

This paper is framed within a large project devoted to studying the presence of circumstellar material around main sequence stars, and looking for exocometary events. The work concentrates on HR 10 (A2 IV/V), known for its conspicuous variability in the circumstellar narrow absorption features of Ca II K and other lines, so far interpreted as $\beta$ Pic-like phenomena, within the falling evaporating body scenario. The main goal of this paper is to carry out a thorough study of HR 10 to find the origin of the observed variability, determine the nature of the star, its absolute parameters, and evolutionary status. Interferometric near-infrared (NIR) observations, multi-epoch high-resolution optical spectra spanning a time baseline of more than 32 years, and optical and NIR photometry, together with theoretical modelling, were used to tackle the above objectives. Our results reveal that HR 10 is a binary. The narrow circumstellar absorption features superimposed on the photospheric Ca II K lines -- and lines of other species -- can be decomposed into two or more components, the two deep ones tracing the radial velocity of the individual stars, which implies that their origin cannot be ascribed to transient exocometary events, their variability being fully explained by the binarity of the object. There does not appear to be transient events associated with potential exocomets. Each individual star holds its own circumstellar shell and there are no traces of a circumbinary envelope. The combined use of the interferometric and radial velocity data leads to a complete spectrometric and orbital solution for the binary, the main parameters being: an orbital period of 747.6 days, eccentricities of the orbits around the centre of mass 0.25 (HR 10-A), 0.21 (HR 10-B) and a mass ratio of q=M$_{\rm B}$/M$_{\rm A}$=0.72-0.84. The stars are slightly off the main sequence, the binary being $\sim530$ Myr old., Comment: 18 pages, 17 figures

(Abridged) The detection of forming planets in disks around young stars remains elusive, and state-of-the-art observational techniques provide somewhat ambiguous results. It has been reported that the pre-transitional T Tauri star LkCa 15 could host three planets; candidate planet b is in the process of formation, as inferred from its H$\alpha$ emission. However, a more recent work casts doubts on the planetary nature of the previous detections. We have observed LkCa 15 with ISIS/WHT. The spectrograph's slit was oriented towards the last reported position of LkCa 15 b (parallel direction) and 90degr from that (perpendicular). The photocenter and full width half maximum (FWHM) of the Gaussians fitting the spatial distribution at H$\alpha$ and the adjacent continuum were measured. A well-known binary (GU CMa) was used as a calibrator to test the spectro-astrometric performance of ISIS/WHT, recovering consistent photocenter and FWHM signals. However, the photocenter shift predicted for LkCa 15 b is not detected, but the FWHM in H$\alpha$ is broader than in the continuum for both slit positions. Our simulations show that the photocenter and FWHM observations cannot be explained simultaneously by an accreting planet. In turn, both spectro-astrometric observations are naturally reproduced from a roughly symmetric Halpha emitting region centered on the star and extent comparable to the orbit originally attributed to the planet at several au. The extended H$\alpha$ emission around LkCa 15 could be related to a variable disk wind, but additional multi-epoch data and detailed modeling are necessary to understand its physical nature. Spectro-astrometry in H$\alpha$ is able to test the presence of accreting planets and can be used as a complementary technique to survey planet formation in circumstellar disks., Comment: Accepted as a Letter in A&A. 4 appendices, 10 figures, 2 tables

(abridged) The star formation rate (SFR) linearly correlates with the amount of dense gas mass (Mdg) involved in the formation of stars both for distant galaxies and clouds in our Galaxy. Similarly, the mass accretion rate (Macc) and the disk mass (Mdisk) of young, Class II stars are also linearly correlated. We plotted the corresponding observational data together, finding a statistically significant correlation that spans ~ 16 orders of magnitude. This probably represents one of the widest ranges of any empirical correlation known, encompassing galaxies that are several kpc in size, pc-size star-forming clouds within our Galaxy, down to young, pre-main sequence stars with au-size protoplanetary disks. We propose a bottom-up hypothesis suggesting that a relation between Macc and the total circumstellar mass surrounding Class 0/I sources (Mcs; disk+envelope) drives the correlation in protostellar-hosting clouds and cloud-hosting galaxies. This is consistent with the fact that the SFRs derived for clouds over a timescale of 2 Myr can be roughly recovered from the sum of instantaneous Macc values of the protostars embedded within them, implying that galactic SFRs averaged over ~ 10-100 Myr should be constant over this period too. Moreover, the sum of the Mcs values directly participating in the formation of the protostellar population in a cloud likely represents a non-negligible fraction of the Mdg within the cloud. If such fraction is ~ 1-35% of the Mdg associated with star-forming clouds and galaxies, then the global correlation for all scales has a near unity slope and an intercept consistent with the (proto-)stellar accretion timescale, Mcs/Macc. Therefore, an additional critical test of our hypothesis is that the Macc-Mdisk correlation for Class II stars should also be observed between Macc and Mcs for Class 0/I sources with similar slope and intercept., Comment: Accepted in A&A. 5 figures, 4 Appendices

We use Gaia Data Release 2 (DR2) to place 252 Herbig Ae/Be stars in the HR diagram and investigate their characteristics and properties. For all known Herbig Ae/Be stars with parallaxes in Gaia DR2, we collected their atmospheric parameters and photometric and extinction values from the literature. To these data we added near- and mid-infrared photometry, collected H$\alpha$ equivalent widths and line profiles, and their binarity status. In addition, we developed a photometric variability indicator from Gaia's DR2 information. We provide masses, ages, luminosities, distances, variabilities and infrared excesses homogeneously derived for the most complete sample of Herbig Ae/Be stars to date. We find that high mass stars have a much smaller infrared excess and have much lower optical variabilities compared to lower mass stars, with the break at around 7M$_{\odot}$. H$\alpha$ emission is generally correlated with infrared excess, with the correlation being stronger for infrared emission at wavelengths tracing the hot dust closest to the star. The variability indicator as developed by us shows that $\sim$25% of all Herbig Ae/Be stars are strongly variable. We observe that the strongly variable objects display doubly peaked H$\alpha$ line profiles, indicating an edge-on disk. The fraction of strongly variable Herbig Ae stars is close to that found for A-type UX Ori stars. It had been suggested that this variability is in most cases due to asymmetric dusty disk structures seen edge-on. The observation here is in strong support of this hypothesis. Finally, the difference in dust properties occurs at 7M$_{\odot}$, while various properties traced at UV/optical wavelengths differ at 3M$_{\odot}$. The latter has been linked to different accretion mechanisms at work whereas the differing infrared properties and variabilities are related to different (dust-)disk dispersal mechanisms., Comment: 15 pages, 8 figures. Accepted for publication in Astronomy & Astrophysics

MWC758 is a young star surrounded by a transitional disk. Recently, a protoplanet candidate has been detected around MWC758 through high-resolution $L'$-band observations. The candidate is located inside the disk cavity at a separation of $\sim$111 mas from the central star, and at an average position angle of $\sim$165.5 degrees. We have performed simultaneous adaptive optics observations of MWC758 in the H$_{\alpha}$ line and the adjacent continuum using SPHERE/ZIMPOL at the Very Large Telescope (VLT). We aim at detecting accreting protoplanet candidates through spectral angular differential imaging observations. The data analysis does not reveal any H$_{\alpha}$ signal around the target. The derived contrast curve in the B_Ha filter allows us to derive a 5$\sigma$ upper limit of $\sim$7.6 mag at 111 mas, the separation of the previously detected planet candidate. This contrast translates into a H$_{\alpha}$ line luminosity of $L_{\rm H_{\alpha}}\lesssim$ 5$\times$10$^{-5}$ $L_{\odot}$ at 111 mas, and an accretion luminosity of $L_{acc} <$3.7$\times$10$^{-4}\,L_{\odot}$. For the predicted mass range of MWC758b, 0.5-5 $M_{\rm Jup}$, this implies accretion rates of $\dot M \lesssim$ 3.4$\times$(10$^{-8}$-10$^{-9})\,M_{\odot}/yr$, for an average planet radius of 1.1 $R_{\rm Jup}$. Therefore, our estimates are consistent with the predictions of accreting circumplanetary accretion models for $R_{\rm in} = 1 R_{\rm Jup}$. In any case, the non-detection of any H$_{\alpha}$ emitting source in the ZIMPOL images does not allow us to unveil the true nature of the $L'$ detected source., Comment: 5 pages, accepted for publication in A&A letters

HD 100546 is one of the few known pre-main-sequence stars that may host a planetary system in its disk. We analyze new VLT/SPHERE/ZIMPOL polarimetric images of HD 100546 with filters in H$\alpha$ and the adjacent continuum. We have probed the disk gap and the surface layers of the outer disk, covering a region < 500 mas (< 55 au at 109 pc) from the star, at an angular resolution of ~ 20 mas. Our data show an asymmetry: the SE and NW regions of the outer disk are more polarized than the SW and NE. This can be explained from a preferential scattering angle close to 90$^o$, consistent with previous polarization images. The outer disk extends from 13 $\pm$ 2 to 45 $\pm$ 9 au, with a position angle and inclination of 137 $\pm$ 5$^o$ and 44 $\pm$ 8$^o$. The comparison with previous estimates suggests that the disk inclination could increase with the stellocentric distance, although the different measurements are still consistent within the error bars. In addition, no direct signature of the innermost candidate companion is detected from polarimetry, confirming recent results based on intensity imagery. We set an upper limit to its mass accretion rate < 10$^{-8}$ M$_{\odot}$/yr for a sub-stellar mass of 15M$_{Jup}$. Finally, we report the first detection (> 3$\sigma$) of a ~ 20 au bar-like structure that crosses the gap through the central region of HD 100546. It is tentatively suggested that the bar could be dust dragged by infalling gas that radially flows from the outer disk. This could represent an exceptional case in which a small-scale radial inflow is observed in a single system. If this scenario is confirmed, it could explain the presence of atomic gas in the inner disk that would otherwise accrete on to the central star on a timescale of a few months/years, as previously indicated from spectro-interferometric data, and could be related with additional (undetected) planets., Comment: 8 pages, 7 figures. Accepted in A&A

We present medium-resolution (R~7000) near-infrared echelle spectroscopic data for 36 MYSOs drawn from the Red MSX Source (RMS) survey. This is the largest sample observed at this resolution at these wavelengths of MYSOs to date. The spectra are characterized mostly by emission from hydrogen recombination lines and accretion diagnostic lines. One MYSO shows photospheric HI absorption, a comparison with spectral standards indicates that the star is an A type star with a low surface gravity, implying that the MYSOs are probably swollen, as also suggested by evolutionary calculations. An investigation of the Br$\gamma$ line profiles shows that most are in pure emission, while 13$\pm$5 % display P Cygni profiles, indicative of outflow, while less than 8$\pm$4 % have inverse P Cygni profiles, indicative of infall. These values are comparable with investigations into the optically bright Herbig Be stars, but not with those of Herbig Ae and T Tauri stars, consistent with the notion that the more massive stars undergo accretion in a different fashion than lower mass objects which are undergoing magnetospheric accretion. Accretion luminosities and rates as derived from the Br$\gamma$ line luminosities agree with results for lower mass sources, providing tentative evidence for massive star formation theories based on scaling of low-mass scenarios. We present Br$\gamma$/Br12 line profile ratios exploiting the fact that optical depth effects can be traced as a function of Doppler shift across the lines. These show that the winds of MYSOs in this sample are nearly equally split between constant, accelerating, and decelerating velocity structures. There are no trends between the types of features we see and bolometric luminosities or near-infrared colours., Comment: 30 pages, 21 figures (with appendices). Accepted for publication in Monthly Notices of the Royal Astronomical Society

The recent progress in instrumentation and telescope development has brought us different ways to observe protoplanetary disks, including interferometers, space missions, adaptive optics, polarimetry, and time- and spectrally-resolved data. While the new facilities have changed the way we can tackle the existing open problems in disk structure and evolution, there is a substantial lack of interconnection between different observing techniques and their user communities. Here, we explore the complementarity of some of the state-of-the-art observing techniques, and how they can be brought together in a collective effort to understand how disks evolve and disperse at the time of planet formation. This paper was born at the "Protoplanetary Discussions" meeting in Edinburgh, 2016. Its goal is to clarify where multi-wavelength observations of disks converge in unveiling disk structure and evolution, and where they diverge and challenge our current understanding. We discuss caveats that should be considered when linking results from different observations, or when drawing conclusions based on limited datasets (in terms of wavelength or sample). We focus on disk properties that are currently being revolutionized by multi-wavelength observations. Specifically: the inner disk radius, holes and gaps and their link to large-scale disk structures, the disk mass, and the accretion rate. We discuss how the links between them, as well as the apparent contradictions, can help us to disentangle the disk physics and to learn about disk evolution., Comment: Accepted for publication in PASA. 37 pages, 9 figures, 1 table. Revised version: corrected problem in Fig 2

This work presents CHARA/VEGA H$\alpha$ spectro-interferometry (R ~ 6000, and $\lambda$/2B ~ 1 mas) of HD 179218 and HD 141569, doubling the sample of Herbig Ae/Be (HAeBe) stars for which this type of observations is available so far. The observed H$\alpha$ emission is spatially unresolved, indicating that the size of the H$\alpha$ emitting region is smaller than ~ 0.21 and 0.12 au for HD 179218 and HD 141529 (~ 15 and 16 R*, respectively). This is smaller than for the two other HAeBes previously observed with the same instrumentation. Two different scenarios have been explored in order to explain the compact line emitting regions. A hot, several thousand K, blackbody disc is consistent with the observations of HD 179218 and HD 141569. Magnetospheric accretion (MA) is able to reproduce the bulk of the H$\alpha$ emission shown by HD 179218, confirming previous estimates from MA shock modelling with a mass accretion rate of 10^-8 Msun/yr, and an inclination to the line of sight between 30 and 50 degr. The H$\alpha$ profile of HD 141569 cannot be fitted from MA due to the high rotational velocity of this object. Putting the CHARA sample together, a variety of scenarios is required to explain the H$\alpha$ emission in HAeBe stars -compact or extended, discs, accretion, and winds-, in agreement with previous Br$\gamma$ spectro-interferometric observations., Comment: 6 pages, 3 figures, 1 table. Accepted in MNRAS

Herbig Ae/Be objects are pre-main sequence stars surrounded by gas- and dust-rich circumstellar discs. These objects are in the throes of star and planet formation, and their characterisation informs us of the processes and outcomes of planet formation processes around intermediate mass stars. Here we analyse the spectral energy distributions of disc host stars observed by the Herschel Open Time Key Programme `Gas in Protoplanetary Systems'. We present Herschel/PACS far-infrared imaging observations of 22 Herbig Ae/Bes and 5 debris discs, combined with ancillary photometry spanning ultraviolet to sub-millimetre wavelengths. From these measurements we determine the diagnostics of disc evolution, along with the total excess, in three regimes spanning near-, mid-, and far-infrared wavelengths. Using appropriate statistical tests, these diagnostics are examined for correlations. We find that the far-infrared flux, where the disc becomes optically thin, is correlated with the millimetre flux, which provides a measure of the total dust mass. The ratio of far-infrared to sub-millimetre flux is found to be greater for targets with discs that are brighter at millimetre wavelengths and that have steeper sub-millimetre slopes. Furthermore, discs with flared geometry have, on average, larger excesses than flat geometry discs. Finally, we estimate the extents of these discs (or provide upper limits) from the observations., Comment: 18 pages, 12 figures main text, 32 figure in Appendix

We present spatially and spectrally resolved Br-gamma emission around the planet-hosting, transitional Herbig Ae/Be star HD 100546. Aiming to gain insight into the physical origin of the line in possible relation to accretion processes, we carried out Br-gamma spectro-interferometry using AMBER/VLTI from three different baselines achieving spatial and spectral resolutions of 2-4 mas and 12000. The Br-gamma visibility is larger than that of the continuum for all baselines. Differential phases reveal a shift between the photocentre of the Br-gamma line -displaced 0.6 mas (0.06 au at 100 pc) NE from the star- and that of the K-band continuum emission -displaced 0.3 mas NE from the star. The photocentres of the redshifted and blueshifted components of the Br-gamma line are located NW and SE from the photocentre of the peak line emission, respectively. Moreover, the photocentre of the fastest velocity bins within the spectral line tends to be closer to that of the peak emission than the photocentre of the slowest velocity bins. Our results are consistent with a Br-gamma emitting region inside the dust inner rim (<0.25 au) and extending very close to the central star, with a Keplerian, disc-like structure rotating counter-clockwise, and most probably flared (25 deg). Even though the main contribution to the Br-gamma line does not come from gas magnetically channelled on to the star, accretion on to HD 100546 could be magnetospheric, implying a mass accretion rate of a few 10^-7 Msun/yr. This value indicates that the observed gas has to be replenished on time-scales of a few months to years, perhaps by planet-induced flows from the outer to the inner disc as has been reported for similar systems., Comment: 5 Figures. Accepted in MNRAS

Herbig Ae/Be stars span a key mass range that links low and high mass stars, and thus provide an ideal window from which to explore their formation. This paper presents VLT/X-Shooter spectra of 91 Herbig Ae/Be stars, HAeBes; the largest spectroscopic study of HAeBe accretion to date. A homogeneous approach to determining stellar parameters is undertaken for the majority of the sample. Measurements of the ultra-violet (UV) are modelled within the context of magnetospheric accretion, allowing a direct determination of mass accretion rates. Multiple correlations are observed across the sample between accretion and stellar properties: the youngest and often most massive stars are the strongest accretors, and there is an almost 1:1 relationship between the accretion luminosity and stellar luminosity. Despite these overall trends of increased accretion rates in HAeBes when compared to classical T Tauri stars, we also find noticeable differences in correlations when considering the Herbig Ae and Herbig Be subsets. This, combined with the difficulty in applying a magnetospheric accretion model to some of the Herbig Be stars, could suggest that another form of accretion may be occurring within the Herbig Be mass range., Comment: 29 pages, 19 figures, Accepted for publication in MNRAS

Correlations between the accretion luminosity and emission line luminosities (L_acc and L_line) of pre-main sequence (PMS) stars have been published for many different spectral lines, which are used to estimate accretion rates. Despite the origin of those correlations is unknown, this could be attributed to direct or indirect physical relations between the emission line formation and the accretion mechanism. This work shows that all (near-UV/optical/near-IR) L_acc-L_line correlations are the result of the fact that the accretion luminosity and the stellar luminosity (L_star) are correlated, and are not necessarily related with the physical origin of the line. Synthetic and observational data are used to illustrate how the L_acc-L_line correlations depend on the L_acc-L_star relationship. We conclude that because PMS stars show the L_acc-L_star correlation immediately implies that L_acc also correlates with the luminosity of all emission lines, for which the L_acc-L_line correlations alone do not prove any physical connection with accretion but can only be used with practical purposes to roughly estimate accretion rates. When looking for correlations with possible physical meaning, we suggest that L_acc/L_star and L_line/L_star should be used instead of L_acc and L_line. Finally, the finding that L_acc has a steeper dependence on L_star for T-Tauri stars than for intermediate-mass Herbig Ae/Be stars is also discussed. That is explained from the magnetospheric accretion scenario and the different photospheric properties in the near-UV., Comment: 9 pages, 5 figures, 2 tables. Accepted in MNRAS

In this paper we investigate the origin of the mid-infrared (IR) hydrogen recombination lines for a sample of 114 disks in different evolutionary stages (full, transitional and debris disks) collected from the {\it Spitzer} archive. We focus on the two brighter {H~{\sc i}} lines observed in the {\it Spitzer} spectra, the {H~{\sc i}}(7-6) at 12.37$\mu$m and the {H~{\sc i}}(9-7) at 11.32$\mu$m. We detect the {H~{\sc i}}(7-6) line in 46 objects, and the {H~{\sc i}}(9-7) in 11. We compare these lines with the other most common gas line detected in {\it Spitzer} spectra, the {[Ne~{\sc iii}]} at 12.81$\mu$m. We argue that it is unlikely that the {H~{\sc i}} emission originates from the photoevaporating upper surface layers of the disk, as has been found for the {[Ne~{\sc iii}]} lines toward low-accreting stars. Using the {H~{\sc i}}(9-7)/{H~{\sc i}}(7-6) line ratios we find these gas lines are likely probing gas with hydrogen column densities of 10$^{10}$-10$^{11}$~cm$^{-3}$. The subsample of objects surrounded by full and transitional disks show a positive correlation between the accretion luminosity and the {H~{\sc i}} line luminosity. These two results suggest that the observed mid-IR {H~{\sc i}} lines trace gas accreting onto the star in the same way as other hydrogen recombination lines at shorter wavelengths. A pure chromospheric origin of these lines can be excluded for the vast majority of full and transitional disks.We report for the first time the detection of the {H~{\sc i}}(7-6) line in eight young (< 20~Myr) debris disks. A pure chromospheric origin cannot be ruled out in these objects. If the {H~{\sc i}}(7-6) line traces accretion in these older systems, as in the case of full and transitional disks, the strength of the emission implies accretion rates lower than 10$^{-10}$M$_{\odot}$/yr. We discuss some advantages of extending accretion indicators to longer wavelengths.

HD 142527 is a young pre-main sequence star with properties indicative of the presence of a giant planet or/and a low-mass stellar companion. We have analyzed an X-Shooter/Very Large Telescope spectrum to provide accurate stellar parameters and accretion rate. The analysis of the spectrum, together with constraints provided by the SED fitting, the distance to the star (140 +- 20 pc) and the use of evolutionary tracks and isochrones, lead to the following set of parameters T_eff = 6550 +- 100 K, log g = 3.75 +- 0.10, L_*/L_sun = 16.3 +- 4.5, M_*/M_sun = 2.0 +- 0.3 and an age of 5.0 +- 1.5 Myr. This stellar age provides further constrains to the mass of the possible companion estimated by Biller et al. (2012), being in-between 0.20 and 0.35 M_sun. Stellar accretion rates obtained from UV Balmer excess modelling, optical photospheric line veiling, and from the correlations with several emission lines spanning from the UV to the near-IR, are consistent to each other. The mean value from all previous tracers is 2 (+- 1) x 10^-7 M_sun yr^-1, which is within the upper limit gas flow rate from the outer to the inner disk recently provided by Cassasus et al. (2013). This suggests that almost all gas transferred between both components of the disk is not trapped by the possible planet(s) in-between but fall onto the central star, although it is discussed how the gap flow rate could be larger than previously suggested. In addition, we provide evidence showing that the stellar accretion rate of HD 142527 has increased by a factor ~ 7 on a timescale of 2-5 years., Comment: Accepted in ApJ. 8 pages, 7 figures

Context. Debris discs are thought to be formed through the collisional grinding of planetesimals, and can be considered as the outcome of planet formation. Understanding the properties of gas and dust in debris discs can help us to comprehend the architecture of extrasolar planetary systems. Herschel Space Observatory far-infrared (IR) photometry and spectroscopy have provided a valuable dataset for the study of debris discs gas and dust composition. This paper is part of a series of papers devoted to the study of Herschel PACS observations of young stellar associations. Aims. This work aims at studying the properties of discs in the Beta Pictoris Moving Group (BPMG) through far-IR PACS observations of dust and gas. Methods. We obtained Herschel-PACS far-IR photometric observations at 70, 100 and 160 microns of 19 BPMG members, together with spectroscopic observations of four of them. Spectroscopic observations were centred at 63.18 microns and 157 microns, aiming to detect [OI] and [CII] emission. We incorporated the new far-IR observations in the SED of BPMG members and fitted modified blackbody models to better characterise the dust content. Results. We have detected far-IR excess emission toward nine BPMG members, including the first detection of an IR excess toward HD 29391.The star HD 172555, shows [OI] emission, while HD 181296, shows [CII] emission, expanding the short list of debris discs with a gas detection. No debris disc in BPMG is detected in both [OI] and [CII]. The discs show dust temperatures in the range 55 to 264 K, with low dust masses (6.6*10^{-5} MEarth to 0.2 MEarth) and radii from blackbody models in the range 3 to 82 AU. All the objects with a gas detection are early spectral type stars with a hot dust component., Comment: 12 pages, 7 figures, 6 tables

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., Comment: 20 pages, 11 figures, accepted for publication in Astronomy & Astrophysics

This work presents X-Shooter/VLT spectra of the prototypical, isolated Herbig Ae stars HD 31648 (MWC 480) and HD 163296 over five epochs separated by timescales ranging from days to months. Each spectrum spans over a wide wavelength range covering from 310 to 2475 nm. We have monitored the continuum excess in the Balmer region of the spectra and the luminosity of twelve ultraviolet, optical and near infrared spectral lines that are commonly used as accretion tracers for T Tauri stars. The observed strengths of the Balmer excesses have been reproduced from a magnetospheric accretion shock model, providing a mean mass accretion rate of 1.11 x 10^-7 and 4.50 x 10^-7 Msun yr^-1 for HD 31648 and HD 163296, respectively. Accretion rate variations are observed, being more pronounced for HD 31648 (up to 0.5 dex). However, from the comparison with previous results it is found that the accretion rate of HD 163296 has increased by more than 1 dex, on a timescale of ~ 15 years. Averaged accretion luminosities derived from the Balmer excess are consistent with the ones inferred from the empirical calibrations with the emission line luminosities, indicating that those can be extrapolated to HAe stars. In spite of that, the accretion rate variations do not generally coincide with those estimated from the line luminosities, suggesting that the empirical calibrations are not useful to accurately quantify accretion rate variability., Comment: 14 pages, 7 Figures, Accepted in ApJ

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., Comment: Accepted for PASP

Circumstellar disc evolution is paramount for the understanding of planet formation. The GASPS program aims at determining the circumstellar gas and solid mass around ~250 pre-main-sequence Herbig Ae and TTauri stars. We aim to understand the origin and nature of the circumstellar matter orbiting 51 Oph, a young (<1 Myr) luminous B9.5 star. We obtained continuum and line observations with the PACS instrument on board the Herschel Space Observatory and continuum data at 1.2 mm with the IRAM 30m telescope. The SED and line fluxes were modelled using the physico-chemo radiative transfer code ProDiMo. We detected a strong emission by OI at 63 microns using the Herschel Space Observatory. The [OI] emission at 145 microns, the [CII] emission at 158 microns, the high-J CO emissions, and the warm water emissions were not detected. Continuum emission was detected at 1.2 mm. The continuum from the near- to the far-infrared and the [OI] emission are well explained by the emission from a compact hydrostatic disc model with a gas mass of 5E-6 MSun, 100 times that of the solid mass. However, this model fails to match the continuum millimeter flux, which hints at a cold outer disc with a mass in solids of 1E-6 MSun or free-free emission from a photoevaporative disc wind. This outer disc can either be devoid of gas and/or is to cold to emit in the [OI] line. A very flat extended disc model (Rout=400 AU) with a fixed vertical structure and dust settling matches all photometric points and most of the [O I] flux. The observations can be explained by an extended flat disc where dust grains have settled. However, a flat gas disc cannot be reproduced by hydrostatic disc models. The low mass of the 51 Oph inner disc in gas and dust may be explained either by the fast dissipation of an initial massive disc or by a very small initial disc mass., Comment: accepted to A&A

This paper summarizes several results concerning the comparison between accretion-related properties of cool (T Tauri; T < 7000 K, M < 1 Msun and hot (Herbig Ae/Be; 7000 < T(K) < 13000; 1 < M(Msun) < 6) pre-main sequence (PMS) stars. This comparison gives insight into the analogies/differences in the physics of the star-disk interaction and in the physical mechanisms driving disk dissipation. Several optical and near-IR line luminosities used for low-mass objects are also valid to estimate typical accretion rates for intermediate-mass stars under similar empirical expressions. In contrast, the Halpha width at 10% of peak intensity is used as an accretion tracer for T Tauris, but is not reliable to estimate accretion rates for Herbig Ae/Bes. This can be explained as a consequence of the different stellar rotation rates that characterize both types of stars. In addition, there are similar trends when the accretion rate is related to the near-IR colours and disk masses, suggesting that viscous accretion disk models are able to explain these trends for both T Tauri and Herbig Ae/Be stars. However, there are two major differences between cool and hot PMS objects. First, the inner gas dissipation timescale, as estimated by relating the accretion rates and the stellar ages, is slightly faster for Herbig Ae/Be stars. This could have implications on the physical mechanism able to form planets around objects more massive than the Sun. Second, the relative position of Herbig Ae/Bes with disks showing signs of inner dust clearing in the accretion rate--disk mass plane contrasts with that of transitional T Tauri stars, when both samples are compared with "classical", non-evolved disks. This latter difference could be pointing to a different physical mechanism driving disk evolution, depending on the stellar regime., Comment: Proceding article of the CoolStars 17 conference (June 2012, Barcelona). Published in Astronomische Nachrichten

Context. HD 172555 is a young A7 star belonging to the Beta Pictoris Moving Group that harbours a debris disc. The Spitzer IRS spectrum of the source showed mid-IR features such as silicates and glassy silica species, indicating the presence of a warm dust component with small grains, which places HD 172555 among the small group of debris discs with such properties. The IRS spectrum also shows a possible emission of SiO gas. Aims. We aim to study the dust distribution in the circumstellar disc of HD 172555 and to asses the presence of gas in the debris disc. Methods. As part of the GASPS Open Time Key Programme, we obtained Herschel-PACS photometric and spectroscopic observations of the source. We analysed PACS observations of HD 172555 and modelled the Spectral Energy Distribution (SED) with a modified blackbody and the gas emission with a two-level population model with no collisional de-excitation. Results. We report for the first time the detection of OI atomic gas emission at 63.18 microns in the HD 172555 circumstellar disc. We detect excesses due to circumstellar dust toward HD 172555 in the three photometric bands of PACS (70, 100, and 160 microns). We derive a large dust particle mass of 4.8e-4 Earth masses and an atomic oxygen mass of 2.5e-2*R^2 Earth masses, where R in AU is the separation between the star and the inner disc. Thus, most of the detected mass of the disc is in the gaseous phase., Comment: 5 pages, 3 figures

We observed a sample of 20 representative Herbig Ae/Be stars and five A-type debris discs with PACS onboard of Herschel. The observations were done in spectroscopic mode, and cover far-IR lines of [OI], [CII], CO, CH+, H2O and OH. We have a [OI]63 micron detection rate of 100% for the Herbig Ae/Be and 0% for the debris discs. [OI]145 micron is only detected in 25%, CO J=18-17 in 45% (and less for higher J transitions) of the Herbig Ae/Be stars and for [CII] 157 micron, we often found spatially variable background contamination. We show the first detection of water in a Herbig Ae disc, HD 163296, which has a settled disc. Hydroxyl is detected as well in this disc. CH+, first seen in HD 100546, is now detected for the second time in a Herbig Ae star, HD 97048. We report fluxes for each line and use the observations as line diagnostics of the gas properties. Furthermore, we look for correlations between the strength of the emission lines and stellar or disc parameters, such as stellar luminosity, UV and X-ray flux, accretion rate, PAH band strength, and flaring. We find that the stellar UV flux is the dominant excitation mechanism of [OI]63 micron, with the highest line fluxes found in those objects with a large amount of flaring and greatest PAH strength. Neither the amount of accretion nor the X-ray luminosity has an influence on the line strength. We find correlations between the line flux of [OI]63 micron and [OI]145 micron, CO J = 18-17 and [OI]6300 \AA, and between the continuum flux at 63 micron and at 1.3 mm, while we find weak correlations between the line flux of [OI]63 micron and the PAH luminosity, the line flux of CO J = 3-2, the continuum flux at 63 micron, the stellar effective temperature and the Brgamma luminosity. (Abbreviated version), Comment: 20 pages, 29 figures, accepted by Astronomy and Astrophysics

We look for trends relating the mass accretion rate (Macc) and the stellar ages (t), spectral energy distributions (SEDs), and disk masses (Mdisk) for a sample of 38 HAeBe stars, comparing them to analogous correlations found for classical T Tauri stars. Our goal is to shed light on the timescale and physical processes that drive evolution of intermediate-mass pre-main sequence objects. Macc shows a dissipation timescale \tau = 1.3^{+1.0}_{-0.5} Myr from an exponential law fit, while a power law yields Macc(t) \propto t^{-\eta}, with \eta = 1.8^{+1.4}_{-0.7}. This result is based on our whole HAeBe sample (1-6 Msun), but the accretion rate decline most probably depends on smaller stellar mass bins. The near-IR excess is higher and starts at shorter wavelengths (J and H bands) for the strongest accretors. Active and passive disks are roughly divided by 2 x 10^{-7} Msun/yr. The mid-IR excess and the SED shape from the Meeus et al. classification are not correlated with Macc. We find Macc \propto Mdisk^{1.1 +- 0.3}. Most stars in our sample with signs of inner dust dissipation typically show accretion rates ten times lower and disk masses three times smaller than the remaining objects. The trends relating Macc with the near-IR excess and Mdisk extend those for T Tauri stars, and are consistent with viscous disk models. The differences in the inner gas dissipation timescale, and the relative position of the stars with signs of inner dust clearing in the Macc-Mdisk plane, could be suggesting a slightly faster evolution, and that a different process - such as photoevaporation - plays a more relevant role in dissipating disks in the HAeBe regime compared to T Tauri stars. Our conclusions must consider the mismatch between the disk mass estimates from mm fluxes and the disk mass estimates from accretion, which we also find in HAeBe stars., Comment: 11 pages, 7 figures, 1 appendix. Accepted in A&A

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., Comment: 4 letter pages, 2 appendix pages, 5 images, to be published in Astronomy and Astrophysics

We present detailed model fits to observations of the disc around the Herbig Ae star HD 163296. This well-studied object has an age of ~ 4 Myr, with evidence of a circumstellar disc extending out to ~ 540AU. We use the radiation thermo-chemical disc code ProDiMo to model the gas and dust in the circumstellar disc of HD 163296, and attempt to determine the disc properties by fitting to observational line and continuum data. These include new Herschel/PACS observations obtained as part of the open-time key program GASPS (Gas in Protoplanetary Systems), consisting of a detection of the [OI]63mic line and upper limits for several other far infrared lines. We complement this with continuum data and ground-based observations of the 12CO 3-2, 2-1 and 13CO J=1-0 line transitions, as well as the H2 S(1) transition. We explore the effects of stellar ultraviolet variability and dust settling on the line emission, and on the derived disc properties. Our fitting efforts lead to derived gas/dust ratios in the range 9-100, depending on the assumptions made. We note that the line fluxes are sensitive in general to the degree of dust settling in the disc, with an increase in line flux for settled models. This is most pronounced in lines which are formed in the warm gas in the inner disc, but the low excitation molecular lines are also affected. This has serious implications for attempts to derive the disc gas mass from line observations. We derive fractional PAH abundances between 0.007 and 0.04 relative to ISM levels. Using a stellar and UV excess input spectrum based on a detailed analysis of observations, we find that the all observations are consistent with the previously assumed disc geometry., Comment: 20 pages, 15 figures. Accepted for publication in A&A

This work aims to derive accretion rates for a sample of 38 HAeBe stars. We apply magnetospheric accretion (MA) shock modelling to reproduce the observed Balmer excesses. We look for possible correlations with the strength of the Halpha, [OI]6300, and Brgamma emission lines. The median mass accretion rate is 2 x 10^-7 Msun yr^-1 in our sample. The model fails to reproduce the large Balmer excesses shown by the four hottest stars (T* > 12000 K). We derive Macc propto M*^5 and Lacc propto L*^1.2 for our sample, with scatter. Empirical calibrations relating the accretion and the Halpha, [OI]6300, and Brgamma luminosities are provided. The slopes in our expressions are slightly shallower than those for lower mass stars, but the difference is within the uncertainties, except for the [OI]6300 line. The Halpha 10% width is uncorrelated with Macc, unlike for the lower mass regime. The mean Halpha width shows higher values as the projected rotational velocities of HAe stars increase, which agrees with MA. The accretion rate variations in the sample are typically lower than 0.5 dex on timescales of days to months, Our data suggest that the changes in the Balmer excess are uncorrelated to the simultaneous changes of the line luminosities. The Balmer excesses and Halpha line widths of HAe stars can be interpreted within the context of MA, which is not the case for several HBes. The steep trend relating Macc and M* can be explained from the mass-age distribution characterizing HAeBe stars. The line luminosities used for low-mass objects are also valid to estimate typical accretion rates for the intermediate-mass regime under similar empirical expressions. However, we suggest that several of these calibrations are driven by the stellar luminosity., Comment: 13 pages, 10 figures. To be published in A&A

We analysed 337 multi-epoch optical spectra of 38 Herbig Ae/Be (HAeBe) stars to gain insights into the variability behaviour of the circumstellar (CS) atomic gas. Equivalent widths (EWs) and line fluxes of the Halpha, [OI]6300, HeI5876 and NaID lines were obtained for each spectrum; the Halpha line width at 10% of peak intensity (W10) and profile shapes were also measured and classified. The mean line strengths and relative variabilities were quantified for each star. Simultaneous optical photometry was used to estimate the line fluxes. We present a homogeneous spectroscopic database of HAeBe stars. The lines are variable in practically all stars and timescales, although 30 % of the objects show a constant EW in [OI]6300, which is also the only line that shows no variability on timescales of hours. The HeI5876 and NaID EW relative variabilities are typically the largest, followed by those in [OI]6300 and Halpha. The EW changes can be larger than one order of magnitude for the HeI5876 line, and up to a factor 4 for Halpha. The [OI]6300 and Halpha EW relative variabilities are correlated for most stars in the sample. The Halpha mean EW and W10 are uncorrelated, as are their relative variabilities. The Halpha profile changes in 70 % of the objects. The massive stars in the sample usually show more stable Halpha profiles with blueshifted self-absorptions and less variable 10% widths. Our data suggest multiple causes for the different line variations, but the [OI]6300 and Halpha variability must share a similar origin in many objects. The physical mechanism responsible for the Halpha line broadening does not depend on the amount of emission; unlike in lower-mass stars, physical properties based on the Halpha luminosity and W10 would significantly differ. Our results provide additional support to previous works that reported different physical mechanisms in Herbig Ae and Herbig Be stars., Comment: 10 pages, 5 figures, 2 appendixes

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., Comment: Comments: 5 pages, 2 Postscript figures, accepted for publication in the A&A Herschel special issue

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., Comment: 5 pages, 4 figures, accepted for publication in the A&A Herschel special issue

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., Comment: 6 pages, 6 figures, accepted for publication in the A&A Herschel special issue

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., Comment: 5 pages, 3 Postscript figures, accepted for publication in the A&A Herschel special issue