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.

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.

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.

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.

(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.

(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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.