Search

Giant planets grow and acquire their gas envelope during the disk phase. At the time of the discovery of giant planets in their host disk, it is important to understand the interplay between the host disk and the envelope and circum-planetary disk properties of the planet. Our aim is to investigate the dynamical and physical structure of the gas in the vicinity of a Jupiter-mass planet and study how protoplanetary disk propertiesy, determine the planetary system as well as the accretion rate inside the planet's Hill sphere.We ran global 3D simulations with the grid-based code fargOCA, using a fully radiative equation of state and a dust-to-gas ratio of 0.01. We explored three models. The nominal one features a disk with surface density, Sigma, corresponding to the MMSN at the planet's location, characterised by an alpha viscosity value of 0.004. The second model has a surface density ten times smaller than the nominal one and the same viscosity. In the third model, we also reduced the viscosity value by a factor of 10. During gap formation gas is heated by compression and cools according to opacity, density, and temperature values. In the analysis of our disks, we find that the gas flowing into the Hill sphere is approximately scaled as the product Sigma nu , as expected from viscous transport while The accretion rate is scaled as sqrt(Sigma nu ).Previous studies have shown that pressure or rotationally supported structures are formed around giant planets, depending on the equation of state (EoS) or on the opacity. In the case of a fully radiative EoS and a constant dust to gas ratio of 0.01, we find that low-mass and low-viscosity circum-stellar disks favour the formation of a rotationally supported circum-planetary disk. Gas accretion leading to the doubling time of the planetary system >10^5 years has only been found in the case of a low-viscosity disk., Comment: 17 pages, 11 figures

The dust- and gas-rich protoplanetary disks around young stellar systems play a key role in star and planet formation. While considerable progress has recently been made in probing these disks on large scales of a few tens of astronomical units (au), the central au needs to be more investigated. We aim at unveiling the physical processes at play in the innermost regions of the strongly accreting T Tauri Star S CrA N by means of near-infrared interferometric observations. The K-band continuum emission is well reproduced with an azimuthally-modulated dusty ring. As the star alone cannot explain the size of this sublimation front, we propose that magnetospheric accretion is an important dust-heating mechanism leading to this continuum emission. The differential analysis of the Hydrogen Br$\gamma$ line is in agreement with radiative transfer models combining magnetospheric accretion and disk winds. Our observations support an origin of the Br$\gamma$ line from a combination of (variable) accretion-ejection processes in the inner disk region.

The dynamical interactions between young binaries can perturb the material distribution of their circumstellar disks, and modify the planet formation process. In order to constrain the impact and nature of the binary interaction in the RW Aur system (bound or unbound), we analyzed the circumstellar material at 1.3 mm wavelengths, as observed at multiple epochs by ALMA. We analyzed the disk properties through parametric visibility modeling, and we used this information to constrain the dust morphology and the binary orbital period. We imaged the dust continuum emission of RW Aur with a resolution of 3 au, and we find that the radius enclosing 90% of the flux (R90%) is 19 au and 14 au for RW Aur A and B, respectively. By modeling the relative distance of the disks at each epoch, we find a consistent trend of movement for the disk of RW Aur B moving away from the disk of RW Aur A at an approximate rate of 3 mas/yr (about 0.5 au/yr in sky-projected distance). By combining ALMA astrometry, historical astrometry, and the dynamical masses of each star, we constrain the RW Aur binary stars to be most likely in a high-eccentricity elliptical orbit with a clockwise prograde orientation relative to RW Aur A, although low-eccentricity hyperbolic orbits are not ruled out by the astrometry. Our analysis does not exclude the possibility of a disk collision during the last interaction, which occurred $295_{-74}^{+21}$ yr ago relative to beginning of 2024. Evidence for the close interaction is found in a tentative warp of 6 deg in the inner 3 au of the disk of RW Aur A, in the brightness temperature of both disks, and in the morphology of the gas emission. A narrow ring that peaks at 6 au around RW Aur B is suggestive of captured material from the disk around RW Aur A., Comment: Accepted in A&A. All self-calibrated data and products are publicly available in https://doi.org/10.5281/zenodo.12825068 . Short video summary in https://youtu.be/IOjEW8Vrj9Q

With two directly detected protoplanets, the PDS 70 system is a unique source in which to study the complex interplay between forming planets and their natal environment. The large dust cavity carved by the two giant planets can affect the disk chemistry, and therefore the molecular emission morphology. On the other hand, chemical properties of the gas component of the disk are expected to leave an imprint on the planetary atmospheres. In this work, we reconstruct the emission morphology of a rich inventory of molecular tracers in the PDS 70 disk, and we look for possible chemical signatures of the two actively accreting protoplanets, PDS b and c. We leverage Atacama Large Millimeter/submillimeter Array (ALMA) band 6 high-angular-resolution and deep-sensitivity line emission observations, together with image and $uv$-plane techniques, to boost the detection of faint lines. We robustly detect ring-shaped emission from $^{12}$CO, $^{13}$CO, C$^{18}$O, H$^{13}$CN, HC$^{15}$N, DCN, H$_2$CO, CS, C$_2$H, and H$^{13}$CO$^{+}$ lines in unprecedented detail. Most of the molecular tracers show a peak of the emission inside the millimeter dust peak. We interpret this as the direct impact of the effective irradiation of the cavity wall, as a result of the planet formation process. Moreover, we have found evidence of an O-poor gas reservoir in the outer disk, which is supported by the observations of bright C-rich molecules, the non-detection of SO, and a lower limit on the $\mathrm{CS/SO}$ ratio of $\sim1$. Eventually, we provide the first detection of the c-C$_3$H$_2$ transitions at 218.73 GHz, and the marginal detection of an azimuthal asymmetry in the higher-energy H$_2$CO (3$_{2,1}$-2$_{2,0}$) line, which could be due to accretion heating near PDS 70b., Comment: 25 pages, 10 figures, 3 tables, 8 figures and one table in appendix. Accepted for publication in A&A

Since 2019, GRAVITY has provided direct observations of giant planets and brown dwarfs at separations of down to 95 mas from the host star. Some of these observations have provided the first direct confirmation of companions previously detected by indirect techniques (astrometry and radial velocities). We want to improve the observing strategy and data reduction in order to lower the inner working angle of GRAVITY in dual-field on-axis mode. We also want to determine the current limitations of the instrument when observing faint companions with separations in the 30-150 mas range. To improve the inner working angle, we propose a fiber off-pointing strategy during the observations to maximize the ratio of companion-light-to-star-light coupling in the science fiber. We also tested a lower-order model for speckles to decouple the companion light from the star light. We then evaluated the detection limits of GRAVITY using planet injection and retrieval in representative archival data. We compare our results to theoretical expectations. We validate our observing and data-reduction strategy with on-sky observations; first in the context of brown dwarf follow-up on the auxiliary telescopes with HD 984 B, and second with the first confirmation of a substellar candidate around the star Gaia DR3 2728129004119806464. With synthetic companion injection, we demonstrate that the instrument can detect companions down to a contrast of $8\times 10^{-4}$ ($\Delta \mathrm{K}= 7.7$ mag) at a separation of 35 mas, and a contrast of $3\times 10^{-5}$ ($\Delta \mathrm{K}= 11$ mag) at 100 mas from a bright primary (K<6.5), for 30 min exposure time. With its inner working angle and astrometric precision, GRAVITY has a unique reach in direct observation parameter space. This study demonstrates the promising synergies between GRAVITY and Gaia for the confirmation and characterization of substellar companions., Comment: 16 pages, 14 figures. Submitted to A&A

Detecting protoplanets during their formation stage is an important but elusive goal of modern astronomy. Kinematic detections via the spiral wakes in the gaseous disc are a promising avenue to achieve this goal. We aim to test the applicability to observations in the low and intermediate planet mass regimes of a commonly used semi-analytical model for planet induced spiral waves. In contrast with previous works which proposed to use the semi-analytical model to interpret observations, in this study we analyse for the first time both the structure of the velocity and density perturbations. We run a set of FARGO3D hydrodynamic simulations and compare them with the output of the semi-analytic model in the code wakeflow, which is obtained by solving Burgers' equation using the simulations as an initial condition. We find that the velocity field derived from the analytic theory is discontinuous at the interface between the linear and nonlinear regions. After 0.2 r$_p$ from the planet, the behaviour of the velocity field closely follows that of the density perturbations. In the low mass limit, the analytical model is in qualitative agreement with the simulations, although it underestimates the azimuthal width and the amplitude of the perturbations, predicting a stronger decay but a slower azimuthal advance of the shock fronts. In the intermediate regime, the discrepancy increases, resulting in a different pitch angle between the spirals of the simulations and the analytic model. The implementation of a fitting procedure based on the minimisation of intensity residuals is bound to fail due to the deviation in pitch angle between the analytic model and the simulations. In order to apply this model to observations, it needs to be revisited accounting also for higher planet masses., Comment: 10 pages, 7 figures. Accepted for publication in A&A May 27, 2024

Context. Most ground-based planet search direct imaging campaigns use angular differential imaging, which distorts the signal from extended sources like protoplanetary disks. In the case PDS 70, a young system with two planets found within the cavity of a protoplanetary disk, obtaining a reliable image of both planets and disk is essential to understanding planet-disk interactions. Aims. Our goals are to reveal the true intensity of the planets and disk without self-subtraction effects for the first time, search for new giant planets beyond separations of 0.1" and to study the morphology of the disk shaped by two massive planets. Methods. We present YJHK-band imaging, polarimetry, and spatially resolved spectroscopy of PDS 70 using near-simultaneous reference star differential imaging, also known as star-hopping. We created a radiative transfer model of the system to match the near-infrared imaging and polarimetric data, along with sub-millimeter imaging from ALMA. Furthermore, we extracted the spectra of the planets and the disk and compared them. Results. We find that the disk is quite flared with a scale height of ~15% at the outer edge of the disk at ~90 au, similar to some disks in the literature. The gap inside of ~50 au is estimated to have ~1% of the dust density of the outer disk. The Northeast outer disk arc seen in previous observations is likely the outer lip of the flared disk. Abundance ratios of grains estimated by the modeling indicate a shallow grain-size index > -2.7, instead of the canonical -3.5. There is both vertical and radial segregation of grains. Planet c is well separated from the disk and has a spectrum similar to planet b, clearly redder than the disk spectra. Planet c is possibly associated with the sudden flaring of the disk starting at ~50 au. No new planets > 5 Mj were found., Comment: Accepted to A&A on April 11, 2024. 20 pages, 19 figures

The detection of low-mass planets orbiting the nearest stars is a central stake of exoplanetary science, as they can be directly characterized much more easily than their distant counterparts. Here, we present the results of our long-term astrometric observations of the nearest binary M-dwarf Gliese 65 AB (GJ65), located at a distance of only 2.67 pc. We monitored the relative astrometry of the two components from 2016 to 2023 with the VLTI/GRAVITY interferometric instrument. We derived highly accurate orbital parameters for the stellar system, along with the dynamical masses of the two red dwarfs. The GRAVITY measurements exhibit a mean accuracy per epoch of 50-60 microarcseconds in 1.5h of observing time using the 1.8m Auxiliary Telescopes. The residuals of the two-body orbital fit enable us to search for the presence of companions orbiting one of the two stars (S-type orbit) through the reflex motion they imprint on the differential A-B astrometry. We detected a Neptune-mass candidate companion with an orbital period of p = 156 +/- 1 d and a mass of m = 36 +/- 7 Mearth. The best-fit orbit is within the dynamical stability region of the stellar pair. It has a low eccentricity, e = 0.1 - 0.3, and the planetary orbit plane has a moderate-to-high inclination of i > 30{\deg} with respect to the stellar pair, with further observations required to confirm these values. These observations demonstrate the capability of interferometric astrometry to reach microarcsecond accuracy in the narrow-angle regime for planet detection by reflex motion from the ground. This capability offers new perspectives and potential synergies with Gaia in the pursuit of low-mass exoplanets in the solar neighborhood., Comment: Corresponding authors: G.Bourdarot, P.Kervella, O.Pfuhl. Accepted in A&A Letters

With four companions at separations from 16 to 71 au, HR 8799 is a unique target for direct imaging, presenting an opportunity for the comparative study of exoplanets with a shared formation history. Combining new VLTI/GRAVITY observations obtained within the ExoGRAVITY program with archival data, we perform a systematic atmospheric characterisation of all four planets. We explore different levels of model flexibility to understand the temperature structure, chemistry and clouds of each planet using both petitRADTRANS atmospheric retrievals and fits to self-consistent radiative-convective equilibrium models. Using Bayesian Model Averaging to combine multiple retrievals, we find that the HR 8799 planets are highly enriched in metals, with [M/H] $\gtrsim$1, and have stellar to super-stellar C/O ratios. The C/O ratio increases with increasing separation from $0.55^{+0.12}_{-0.10}$ for d to $0.78^{+0.03}_{-0.04}$ for b, with the exception of the innermost planet which has a C/O ratio of $0.87\pm0.03$. By retrieving a quench pressure and using a disequilibrium chemistry model we derive vertical mixing strengths compatible with predictions for high-metallicity, self-luminous atmospheres. Bayesian evidence comparisons strongly favour the presence of HCN in HR 8799 c and e, as well as CH$_{4}$ in HR 8799 c, with detections at $>5\sigma$ confidence. All of the planets are cloudy, with no evidence for patchiness. The clouds of c, d and e are best fit by silicate clouds lying above a deep iron cloud layer, while the clouds of the cooler HR 8799 b are more likely composed of Na$_{2}$S. With well defined atmospheric properties, future exploration of this system is well positioned to unveil further detail in these planets, extending our understanding of the composition, structure, and formation history of these siblings., Comment: 45 pages, 25 figures, accepted by Astronomy & Astrophysics

Precise mass constraints are vital for the characterisation of brown dwarfs and exoplanets. Here we present how the combination of data obtained by Gaia and GRAVITY can help enlarge the sample of substellar companions with measured dynamical masses. We show how the Non-Single-Star (NSS) two-body orbit catalogue contained in Gaia DR3 can be used to inform high-angular-resolution follow-up observations with GRAVITY. Applying the method presented in this work to eight Gaia candidate systems, we detect all eight predicted companions, seven of which were previously unknown and five are of a substellar nature. Among the sample is Gaia DR3 2728129004119806464 B, which - detected at an angular separation of (34.01 $\pm$ 0.15) mas from the host - is the closest substellar companion ever imaged. This translates to a semi-major axis of (0.938 $\pm$ 0.023) AU. WT 766 B, detected at a greater angular separation, was confirmed to be on an orbit exhibiting an even smaller semi-major axis of (0.676 $\pm$ 0.008) AU. The GRAVITY data were then used to break the host-companion mass degeneracy inherent to the Gaia NSS orbit solutions as well as to constrain the orbital solutions of the respective target systems. Knowledge of the companion masses enabled us to further characterise them in terms of their ages, effective temperatures, and radii via the application of evolutionary models. The inferred ages exhibit a distinct bias towards values younger than what is to be expected based on the literature. The results serve as an independent validation of the orbital solutions published in the NSS two-body orbit catalogue and show that the combination of astrometric survey missions and high-angular-resolution direct imaging holds great promise for efficiently increasing the sample of directly imaged companions in the future, especially in the light of Gaia's upcoming DR4 and the advent of GRAVITY+., Comment: Accepted for publication in Astronomy and Astrophysics

Eruptive stars are a class of young stellar objects that show an abrupt increase in their luminosity. These burst-like episodes are thought to dominate the stellar accretion process during the class 0/class I stage. We present an overview of a survey of seven episodically accreting protostars aimed at studying their potentially complex circumstellar surroundings. The observations were performed with the instrument SPHERE, mounted at the VLT. We observed the eruptive stars in $H$-band with the near-infrared imager IRDIS and used the polarimeter to extract the polarized light scattered from the stars' surroundings. We produced polarized light images for three FUor objects, Z CMa, V960 Mon, and FU Ori, and four EXor objects, XZ Tau, UZ Tau, NY Ori, and EX Lup. We calculated the intrinsic polarization fraction for all the observed stars. In all systems we registered scattered light from around the primary star. FU Ori and V960 Mon are surrounded by complex structures including spiral-like features. In Z CMa, we detected a point source 0.7 arcsec to the northeast of the primary. Based on the astrometric measurements from archival Keck/NIRC2 data, we find this source to be a third member of the system. Further, Z CMa displays an outflow extending thousands of au. Unlike the other EXor objects in our sample, XZ Tau shows bright, extended scattered light structures, also associated with an outflow on a scale of hundreds of au. The other EXors show relatively faint disk-like structures in the immediate vicinity of the coronagraph. Asymmetric arms were only found around FUor objects, while faint disks seem to predominantly occur around EXors. Importantly, for Z CMa the detection of the faint extended structure questions previous interpretations of the system's dynamic state. The streamer which was associated with a fly-by object turned out to be part of a huge outflow extending 6000 au., Comment: Accepted for publication in A&A. 14 pages, 7 figures, 2 tables

The sample of planet-forming disks observed by high-contrast imaging campaigns over the last decade is mature enough to enable the demographical analysis of individual star-forming regions. We present the full census of Taurus sources with VLT/SPHERE polarimetric images available. The whole sample sums up to 43 targets (of which 31 have not been previously published) corresponding to one-fifth of the Class II population in Taurus and about half of such objects that are observable. A large fraction of the sample is apparently made up of isolated faint disks (equally divided between small and large self-shadowed disks). Ambient signal is visible in about one-third of the sample. This probes the interaction with the environment and with companions or the outflow activity of the system. The central portion of the Taurus region almost exclusively hosts faint disks, while the periphery also hosts bright disks interacting with their surroundings. The few bright disks are found around apparently older stars. The overall picture is that the Taurus region is in an early evolutionary stage of planet formation. Yet, some objects are discussed individually, as in an intermediate or exceptional stage of the disk evolution. This census provides a first benchmark for the comparison of the disk populations in different star forming regions., Comment: Accepted by A&A

We present SPHERE/IRDIS H-band data for a sample of 23 stars in the Orion Star forming region observed within the DESTINYS (Disk Evolution Study Through Imaging of Nearby Young Stars) program. We use polarization differential imaging in order to detect scattered light from circumstellar dust. From the scattered light observations we characterize the disk orientation, radius and contrast. We analyse the disks in context of the stellar parameters and the environment of the Orion star-forming region. We use ancillary X-shooter spectroscopic observations to characterize the central stars in the systems. We furthermore use a combination of new and archival ALMA mm-continuum observations to characterize the dust masses present in the circumstellar disks. Within our sample we detect extended circumstellar disks in 10 of 23 systems. Of these, three are exceptionally extended (V351 Ori, V599 Ori and V1012 Ori) and show scattered light asymmetries which may indicate perturbations by embedded planets or (in the case of V599 Ori) by an outer stellar companion. Our high resolution imaging observations are also sensitive to close (sub)stellar companions and we detect 9 such objects in our sample of which 5 were previously unknown. We find in particular a possible sub-stellar companion (either a very low mass star or a high mass brown dwarf) 137 au from the star RY Ori. We find a strong anti-correlation between disk detection and multiplicity, with only 2 of our 10 disk detections located in stellar multiple systems. We also find a correlation between scattered light contrast and the millimetre flux suggesting that disks that have a high dust content are typically bright in near-infrared scattered light. Conversely we do not find significant correlations between scattered light contrast of the disks and the stellar mass or age., Comment: Accepted by A&A

We used VLT/SPHERE to observe 20 systems in the Cha I cloud in polarized scattered light in the near-infrared. We combined the scattered light observations with existing literature data on stellar properties and with archival ALMA continuum data to study trends with system age and dust mass. We also connected resolved near-infrared observations with the spectral energy distributions of the systems. In 13 of the 20 systems included in this study we detected resolved scattered light signals from circumstellar dust. For the CR Cha, CT Cha, CV Cha, SY Cha, SZ Cha, and VZ Cha systems we present the first detailed descriptions of the disks in scattered light. The observations found typically smooth or faint disks, often with little substructure, with the notable exceptions of SZ Cha, which shows an extended multiple-ringed disk, and WW Cha, which shows interaction with the cloud environment. New high S/N K- band observations of the HD 97048 system in our survey reveal a significant brightness asymmetry that may point to disk misalignment and subsequent shadowing of outer disk regions, possibly related to the suggested planet candidate in the disk. We resolve for the first time the stellar binary in the CS Cha system. Multiple wavelength observations of the disk around CS Cha have revealed that the system contains small, compact dust grains that may be strongly settled, consistent with numerical studies of circumbinary disks. We find in our sample that there is a strong anti-correlation between the presence of a (close) stellar companion and the detection of circumstellar material with five of our seven nondetections located in binary systems., Comment: Accepted by A&A

Dual-field interferometric observations with VLTI/GRAVITY sometimes require the use of a "binary calibrator", a binary star whose individual components remain unresolved by the interferometer, with a separation between 400 and 2000 mas for observations with the Units Telescopes (UTs), or 1200 to 3000 mas for the Auxiliary Telescopes (ATs). The separation vector also needs to be predictable to within 10 mas for proper pointing of the instrument. Up until now, no list of properly vetted calibrators was available for dual-field observations with VLTI/GRAVITY on the UTs. Our objective is to compile such a list, and make it available to the community. We identify a list of candidates from the Washington Double Star (WDS) catalogue, all with appropriate separations and brightness, scattered over the Southern sky. We observe them as part of a dedicated calibration programme, and determine whether these objects are true binaries (excluding higher multiplicities resolved interferometrically but unseen by imaging), and extract measurements of the separation vectors. We combine these new measurements with those available in the WDS to determine updated orbital parameters for all our vetted calibrators. We compile a list of 13 vetted binary calibrators for observations with VLTI/GRAVITY on the UTs, and provide orbital estimates and astrometric predictions for each of them. We show that our list guarantees that there are always at least two binary calibrators at airmass < 2 in the sky over the Paranal observatory, at any point in time. Any Principal Investigator wishing to use the dual-field mode of VLTI/GRAVITY with the UTs can now refer to this list to select an appropriate calibrator. We encourage the use of "whereistheplanet" to predict the astrometry of these calibrators, which seamlessly integrates with "p2Gravity" for VLTI/GRAVITY dual-field observing material preparation., Comment: 16 pages, 6 figures. Submitted to A&A

Substructures in disc density are ubiquitous in the bright extended discs that are observed with high resolution. These substructures are intimately linked to the physical mechanisms driving planet formation and disc evolution. Surveys of star-forming regions find that most discs are in fact compact, less luminous, and do not exhibit these same substructures. It remains unclear whether compact discs also have similar substructures or if they are featureless. This suggests that different planet formation and disc evolution mechanisms operate in these discs. We investigated evidence of substructure within two compact discs around the stars Sz 65 and Sz 66 using high angular resolution observations with ALMA at 1.3 mm. The two stars form a wide-binary system with 6.36 arcsec separation. The continuum observations achieve a synthesised beam major axis of 0.026 arcsec, equivalent to about 4.0 au, enabling a search for substructure on these spatial scales and a characterisation of the gas and dust disc sizes with high precision. We analysed the data in the image plane through an analysis of reconstructed images, as well as in the uv plane by modelling the visibilities and by an analysis of the 12CO emission line. Comparisons were made with high-resolution observations of compact discs and radially extended discs. We find evidence of substructure in the dust distribution of Sz 65, namely a shallow gap centred at approximately 20 au, with an emission ring exterior to it. Ninety percent of the measured continuum flux is found within 27 au, and the distance for 12CO is 142 au. The observations show that Sz 66 is very compact: 90 per cent of the continuum flux is contained within 16 au, and 48 au for the gas. While the overall prevalence and diversity of substructure in compact discs relative to larger discs is yet to be determined, we find evidence that substructures can exist in compact discs.

Temporal variability in the photometric and spectroscopic properties of protoplanetary disks is common in YSO. However, evidence pointing toward changes in their morphology over short timescales has only been found for a few sources, mainly due to a lack of high cadence observations at mas resolution. We combine GRAVITY multi-epoch observations of HD98922 at mas resolution with PIONIER archival data covering a total time span of 11 years. We interpret the interferometric visibilities and spectral energy distribution with geometrical models and through radiative transfer techniques. We investigated high-spectral-resolution quantities to obtain information on the properties of the HI BrG-line-emitting region. The observations are best fitted by a model of a crescent-like asymmetric dust feature located at 1 au and accounting for 70% of the NIR emission. The feature has an almost constant magnitude and orbits the central star with a possible sub-Keplerian period of 12 months, although a 9 month period is another, albeit less probable, solution. The radiative transfer models show that the emission originates from a small amount of carbon-rich (25%) silicates, or quantum-heated particles located in a low-density region. Among different possible scenarios, we favor hydrodynamical instabilities in the inner disk that can create a large vortex. The high spectral resolution differential phases in the BrG-line show that the hot-gas component is offset from the star and in some cases is located between the star and the crescent feature. The scale of the emission does not favor magnetospheric accretion as a driving mechanism. The scenario of an asymmetric disk wind or a massive accreting substellar or planetary companion is discussed. With this unique observational data set for HD98922, we reveal morphological variability in the innermost 2 au of its disk region., Comment: 45 pages, 20 figures, accepted by and to be published in Astronomy & Astrophysics (A&A)

The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy community under the lead of the National Radio Astronomy Observatory (NRAO), and is anticipated to be built between 2027 and 2037. Two workshops have been held in 2022 and 2023 with the goal to discuss and consolidate the scientific interests in the ngVLA within the German astronomical community. This community paper constitutes a collection of 48 science ideas which the German community aims to pursue with the ngVLA in the 2030s. This is not a complete list and the ideas are not developed at the level of a "Science Book", such that the present document is mainly meant provide a basis for further discussion within the community. As such, additional contributions are welcome, and will be considered for inclusion in future revisions., Comment: Version 2.0 (status June 18, 2024): 169 pages, comments and future contributions welcome [v2.0: 7 new science cases added, some minor revisions to other chapters]

Context. The study of protoplanetary disks is fundamental to understand their evolution and interaction with the surrounding environment, and to constrain planet formation mechanisms. Aims. We aim at characterising the young binary system HD 34700 A, which shows a wealth of structures. Methods. Taking advantage of the high-contrast imaging instruments SPHERE at the VLT, LMIRCam at the LBT, and of ALMA observations, we analyse this system at multiple wavelengths. We study the rings and spiral arms morphology and the scattering properties of the dust. We discuss the possible causes of all the observed features. Results. We detect for the first time, in the H${\alpha}$ band, a ring extending from $\sim$65 au to ${\sim}$120 au, inside the ring already known from recent studies. These two have different physical and geometrical properties. Based on the scattering properties, the outer ring may consist of grains of typical size $a_{out} > 4 {\mu}m$, while the inner ring of smaller grains ($a_{in} <= 0.4 {\mu m}$). Two extended logarithmic spiral arms stem from opposite sides of the disk. The outer ring appears as a spiral arm itself, with a variable radial distance from the centre and extended substructures. ALMA data confirm the presence of a millimetric dust substructure centred just outside the outer ring, and detect misaligned gas rotation patterns for HD 34700 A and B. Conclusions. The complexity of HD 34700 A, revealed by the variety of observed features, suggests the existence of one or more disk-shaping physical mechanisms. Possible scenarios, compatible with our findings, involve the presence inside the disk of a yet undetected planet of several Jupiter masses and the system interaction with the surroundings by means of gas cloudlet capture or flybys. Further observations with JWST/MIRI or ALMA (gas kinematics) could shed more light on these., Comment: Accepted for publication on A&A. 14 + 5 pages, 9 + 7 figures (text + appendix)

Giant exoplanets have been directly imaged over orders of magnitude of orbital separations, prompting theoretical and observational investigations of their formation pathways. In this paper, we present new VLTI/GRAVITY astrometric data of HIP 65426 b, a cold, giant exoplanet which is a particular challenge for most formation theories at a projected separation of 92 au from its primary. Leveraging GRAVITY's astrometric precision, we present an updated eccentricity posterior that disfavors large eccentricities. The eccentricity posterior is still prior-dependent, and we extensively interpret and discuss the limits of the posterior constraints presented here. We also perform updated spectral comparisons with self-consistent forward-modeled spectra, finding a best fit ExoREM model with solar metallicity and C/O=0.6. An important caveat is that it is difficult to estimate robust errors on these values, which are subject to interpolation errors as well as potentially missing model physics. Taken together, the orbital and atmospheric constraints paint a preliminary picture of formation inconsistent with scattering after disk dispersal. Further work is needed to validate this interpretation. Analysis code used to perform this work is available at https://github.com/sblunt/hip65426., Comment: 28 pages, 14 figures. Revised and resubmitted to AJ

Tension remains between the observed and modeled properties of substellar objects, but objects in binary orbits, with known dynamical masses can provide a way forward. HD 72946 B is a recently imaged brown dwarf companion to the nearby, solar type star. We achieve $\sim100~\mu\mathrm{as}$ relative astrometry of HD 72946 B in the K-band using VLTI/GRAVITY, unprecedented for a benchmark brown dwarf. We fit an ensemble of measurements of the orbit using orbitize! and derive a strong dynamical mass constraint $\mathrm{M_B}=69.5\pm0.5~\mathrm{M_{Jup}}$ assuming a strong prior on the host star mass $\mathrm{M_A}=0.97\pm0.01~\mathrm{M_\odot}$ from an updated stellar analysis. We fit the spectrum of the companion to a grid of self-consistent BT-Settl-CIFIST model atmospheres, and perform atmospheric retrievals using petitRADTRANS. A dynamical mass prior only marginally influences the sampled distribution on effective temperature, but has a large influence on the surface gravity and radius, as expected. The dynamical mass alone does not strongly influence retrieved pressure-temperature or cloud parameters within our current retrieval setup. Independent of cloud prescription and prior assumptions, we find agreement within $\pm2\,\sigma$ between the C/O ratio of the host ($0.52\pm0.05)$ and brown dwarf ($0.43$ to $0.63$), as expected from a molecular cloud collapse formation scenario, but our retrieved metallicities are implausibly high ($0.6-0.8$) in light of an excellent agreement of the data with the solar abundance model grid. Future work on our retrieval framework will seek to resolve this tension. Additional study of low surface-gravity objects is necessary to assess the influence of a dynamical mass prior on atmospheric analysis., Comment: Accepted for publication in the Astrophysical Journal. 32 pages, 14 figures, 7 tables. v2&3 correct errors in co-author's affiliations, figure rendering, and some grant acknowledgements

Embedded planets are potentially the cause of substructures like gaps and cavities observed in several protoplanetary disks. Thus, the substructures observed in the continuum and in line emission encode information about the presence of planets in the system and how they interact with the natal disk. The pre-transitional disk around the star PDS 70 is the first case of two young planets imaged within a dust depleted gap that was likely carved by themselves. We aim to determine the spatial distribution of the gas and dust components in the PDS 70 disk. The axisymmetric substructures observed in the resulting profiles are interpreted in the context of planet-disk interactions. We develop a thermo-chemical forward model for an axisymmetric disk to explain a subset of the Atacama Large Millimeter/Submillimeter Array (ALMA) band 6 observations of three CO isotopologues plus the continuum towards PDS 70. Combining the inferred gas and dust distributions, the model results in a variable gas-to-dust ratio profile throughout the disk that spans two orders of magnitude within the first $130$ au and shows a step gradient towards the outer disk, which is consistent with the presence of a pressure maxima driven by planet-disk interactions. We find a gas density drop factor of ${\sim} 19$ at the location of the planet PDS 70 c with respect to the peak gas density at $75$ au. Combining this value with literature results on the hydrodynamics of planet-disk interactions, we find this gas gap depth to be consistent with independent planet mass estimates from infrared observations. Our findings point towards gas stirring processes taking place in the common gap due to the gravitational perturbation of both planets., Comment: Accepted for publication in A&A

Understanding the chemical link between protoplanetary disks and planetary atmospheres is complicated by the fact that the popular targets in the study of disks and planets are widely separated both in space and time. The 5 Myr PDS 70 systems offers a unique opportunity to directly compare the chemistry of a giant planet's atmosphere to the chemistry of its natal disk. To that end, we derive our current best physical and chemical model for the PDS 70 disk through forward modelling of the $^{12}$CO, C$^{18}$O, and C$_2$H emission radial profiles with the thermochemical code DALI and find a volatile C/O ratio above unity in the outer disk. Using what we know of the PDS 70 disk today, we analytically estimate the properties of the disk as it was 4 Myr in the past when we assume that the giant planets started their formation, and compute a chemical model of the disk at that time. We compute the formation of PDS 70b and PDS 70c using the standard core accretion paradigm and account for the accretion of volatile and refractory sources of carbon and oxygen to estimate the resulting atmospheric carbon-to-oxygen number ratio (C/O) for these planets. Our inferred C/O ratio of the gas in the PDS 70 disk indicates that it is marginally carbon rich relative to the stellar C/O = 0.44 which we derive from an empirical relation between stellar metallicity and C/O. Under the assumption that the disk has been carbon rich for most of its lifetime, we find that the planets acquire a super-stellar C/O in their atmospheres. If the carbon-rich disk is a relatively recent phenomenon (i.e. developed after the formation of the planets at $\sim 1$ Myr) then the planets should have close to the stellar C/O in their atmospheres. This work lays the groundwork to better understand the disk in the PDS 70 system as well as the planet formation scenario that produce its planets., Comment: 18 pages, 14 figures, 5 tables, accepted for publication in A&A

Birth environments of young stars have strong imprints on the star itself and their surroundings. We present a detailed analysis of the wealthy circumstellar environment around the young Herbig Ae/Be star TCrA. Our aim is to understand the nature of the stellar system and the extended circumstellar structures as seen in scattered light images. We conduct our analysis combining archival data, and new adaptive optics high-contrast and high-resolution images. The scattered light images reveal the presence of a complex environment composed of a bright forward scattering rim of the disk's surface that is seen at very high inclination, a dark lane of the disk midplane, bipolar outflows, and streamer features likely tracing infalling material from the surrounding birth cloud onto the disk. The analysis of the light curve suggests the star is a binary with a period of 29.6yrs. The comparison of the scattered light images with ALMA continuum and 12CO line emission shows the disk is in keplerian rotation, with the northern side of the outflowing material receding, while the southern side approaching the observer. The disk is itself seen edge-on. The direction of the outflows seen in scattered light is in agreement with the direction of the more distant molecular hydrogen emission-line objects (MHOs) associated to the star. Modeling of the SED using a radiative transfer scheme well agrees with the proposed configuration, as well as the hydrodynamical simulation performed using a Smoothed Particle Hydrodynamics code. We find evidence of streamers of accreting material around TCrA. These streamers connect the filament along which TCrA is forming with the outer parts of the disk, suggesting that the strong misalignment between the inner and outer disk is due to a change in the direction of the angular momentum of the material accreting on the disk during the late phase of star formation., Comment: 17 pages, 14 figures

Aims: We aim to spatially and spectrally resolve the Br-gamma hydrogen emission line with the methods of interferometry in order to examine the kinematics of the hydrogen gas emission region in the inner accretion disk of a sample of solar-like young stellar objects. The goal is to identify trends and categories among the sources of our sample and to discuss whether or not they can be tied to different origin mechanisms associated with Br-gamma emission in T Tauri stars, chiefly and most prominently magnetospheric accretion. Methods: We observed a sample of seven T Tauri stars for the first time with VLTI GRAVITY, recording spectra and spectrally dispersed interferometric quantities across the Br-gamma line in the NIR K-band. We use them to extract the size of the Br-gamma emission region and the photocenter shifts. To assist in the interpretation, we also make use of radiative transfer models of magnetospheric accretion to establish a baseline of expected interferometric signatures if accretion is the primary driver of Br-gamma emission. Results: From among our sample, we find that five of the seven T~Tauri stars show an emission region with a half-flux radius in the range broadly expected for magnetospheric truncation. Two of the five objects also show Br-gamma emission primarily originating from within the corotation radius, while two other objects exhibit extended emission on a scale beyond 10 R$_*$, one of them even beyond the K~band continuum half-flux radius of 11.3 R$_*$. Conclusions: We find strong evidence to suggest that for the two weakest accretors in the sample, magnetospheric accretion is the primary driver of Br-gamma radiation. The results for the remaining sources imply either partial or strong contributions coming from spatially extended emission components in the form of outflows, such as stellar or disk winds., Comment: Version corresponds to the one in A&A production. The author list has been amended with "GRAVITY Collaboration" as first entry in keeping with the convention established for papers published by the GRAVITY consortium in recent years

(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

Disk winds are an important mechanism for accretion and disk evolution around young stars. The accreting intermediate-mass T-Tauri star RY Tau has an active jet and a previously known disk wind. Archival optical and new near-infrared observations of the RY Tau system show two horn-like components stretching out as a cone from RY Tau. Scattered light from the disk around RY Tau is visible in near-infrared but not seen at optical wavelengths. In the near-infrared, dark wedges that separates the horns from the disk, indicating we may see the scattered light from a disk wind. We use archived ALMA and SPHERE/ZIMPOL I-band observations combined with newly acquired SPEHRE/IRDIS H-band observations and available literature to build a simple geometric model of the RY Tau disk and disk wind. We use Monte Carlo radiative transfer modelling \textit{MCMax3D} to create comparable synthetic observations that test the effect of a dusty wind on the optical effect in the observations. We constrain the grain size and dust mass needed in the disk wind to reproduce the effect from the observations. A model geometrically reminiscent of a dusty disk wind with small micron to sub-micron size grains elevated above the disk can reproduce the optical effect seen in the observations. The mass in the obscuring component of the wind has been constrained to $1\times10^{-9} M_{\odot} \leq M \leq 5\times10^{-8} M_{\odot}$ which corresponds to a lower limit mass loss rate in the wind of about $\sim 1\times10^{-8}M_{\odot}\mathrm{yr}^{-1}$. While an illuminate dust cavity cannot be ruled out without measurements of the gas velocity, we argue that a magnetically launched disk wind is the most likely scenario., Comment: 11 pages, 5 figures

SPHERE+ is a proposed upgrade of the SPHERE instrument at the VLT, which is intended to boost the current performances of detection and characterization for exoplanets and disks. SPHERE+ will also serve as a demonstrator for the future planet finder (PCS) of the European ELT. The main science drivers for SPHERE+ are 1/ to access the bulk of the young giant planet population down to the snow line ($3-10$ au), to bridge the gap with complementary techniques (radial velocity, astrometry); 2/ to observe fainter and redder targets in the youngest ($1-10$\,Myr) associations compared to those observed with SPHERE to directly study the formation of giant planets in their birth environment; 3/ to improve the level of characterization of exoplanetary atmospheres by increasing the spectral resolution in order to break degeneracies in giant planet atmosphere models. Achieving these objectives requires to increase the bandwidth of the xAO system (from $\sim$1 to 3\,kHz) as well as the sensitivity in the infrared (2 to 3\,mag). These features will be brought by a second stage AO system optimized in the infrared with a pyramid wavefront sensor. As a new science instrument, a medium resolution integral field spectrograph will provide a spectral resolution from 1000 to 5000 in the J and H bands. This paper gives an overview of the science drivers, requirements and key instrumental trade-off that were done for SPHERE+ to reach the final selected baseline concept., Comment: To appear in the Proceedings of the SPIE Astronomical Telescopes + Instrumentation (2022), 13 pages, 6 figure

Long term precise radial velocity (RV) monitoring of the nearby star HD206893, as well as anomalies in the system proper motion, have suggested the presence of an additional, inner companion in the system. Here we describe the results of a multi-epoch search for the companion responsible for this RV drift and proper motion anomaly using the VLTI/GRAVITY instrument. Utilizing information from ongoing precision RV measurements with the HARPS spectrograph, as well as Gaia host star astrometry, we report a high significance detection of the companion HD206893c over three epochs, with clear evidence for Keplerian orbital motion. Our astrometry with $\sim$50-100 $\mu$arcsec precision afforded by GRAVITY allows us to derive a dynamical mass of 12.7$^{+1.2}_{-1.0}$ M$_{\rm Jup}$ and an orbital separation of 3.53$^{+0.08}_{-0.06}$ au for HD206893c. Our fits to the orbits of both companions in the system utilize both Gaia astrometry and RVs to also provide a precise dynamical estimate of the previously uncertain mass of the B component, and therefore derive an age of $155\pm15$ Myr. We find that theoretical atmospheric/evolutionary models incorporating deuterium burning for HD206893c, parameterized by cloudy atmospheres provide a good simultaneous fit to the luminosity of both HD206893B and c. In addition to utilizing long-term RV information, this effort is an early example of a direct imaging discovery of a bona fide exoplanet that was guided in part with Gaia astrometry. Utilizing Gaia astrometry is expected to be one of the primary techniques going forward to identify and characterize additional directly imaged planets. Lastly, this discovery is another example of the power of optical interferometry to directly detect and characterize extrasolar planets where they form at ice-line orbital separations of 2-4\,au., Comment: Accepted to A&A

[abridged] We present new scattered light and millimeter observations of the protoplanetary disk around LkH$\alpha\,330$, using SPHERE/VLT and ALMA, respectively. The scattered-light SPHERE observations reveal an asymmetric ring at around 45au from the star in addition to two spiral arms with similar radial launching points at around 90au. The millimeter observations from ALMA (resolution of 0.06''$\times$0.04'') mainly show an asymmetric ring located at 110au from the star. In addition to this asymmetry, there are two faint symmetric rings at 60au and 200au. The $^{12}$CO, $^{13}$CO, and C$^{18}$O lines seem to be less abundant in the inner disk (these observations have a resolution of 0.16''$\times$0.11''). The $^{13}$CO peaks at a location similar to the inner ring observed with SPHERE, suggesting that this line is optically thick and traces variations of disk temperature instead of gas surface-density variations, while the C$^{18}$O peaks slightly further away at around 60au. We compare our observations with hydrodynamical simulations that include gas and dust evolution, and conclude that a 10$M_{\rm{Jup}}$ mass planet at 60au and in an eccentric orbit ($e=0.1$) can qualitatively explain most of the observed structures. A planet in a circular orbit leads to a much narrower concentration in the millimeter emission, while a planet in a more eccentric orbit leads to a very eccentric cavity as well. In addition, the outer spiral arm launched by the planet changes its pitch angle along the spiral due to the eccentricity and when it interacts with the vortex, potentially appearing in observations as two distinct spirals. Our observations and models show that LkH$\alpha\,330$ is an interesting target to search for (eccentric-) planets while they are still embedded in their parental disk, making it an excellent candidate for studies on planet-disk interaction., Comment: Accepted for publication in A&A, minor changes after language edition

Planets have been detected in circumbinary orbits in several different systems, despite the additional challenges faced during their formation in such an environment. We investigate the possibility of planetary formation in the spectroscopic binary CS Cha by analyzing its circumbinary disk. The system was studied with high angular resolution ALMA observations at 0.87mm. Visibilities modeling and Keplerian fitting are used to constrain the physical properties of CS Cha, and the observations were compared to hydrodynamic simulations. Our observations are able to resolve the disk cavity in the dust continuum emission and the 12CO J:3-2 transition. We find the dust continuum disk to be azimuthally axisymmetric (less than 9% of intensity variation along the ring) and of low eccentricity (of 0.039 at the peak brightness of the ring). Under certain conditions, low eccentricities can be achieved in simulated disks without the need of a planet, however, the combination of low eccentricity and axisymmetry is consistent with the presence of a Saturn-like planet orbiting near the edge of the cavity., Comment: 18 pages, 13 figures, accepted in A&A. Watch a 3min summary by the first author in https://youtu.be/6j5ROSj2N-U

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

High dust density in the midplane of protoplanetary disks is favorable for efficient grain growth and can allow fast formation of planetesimals and planets, before disks dissipate. Vertical settling and dust trapping in pressure maxima are two mechanisms allowing dust to concentrate in geometrically thin and high density regions. In this work, we aim to study these mechanisms in the highly inclined protoplanetary disk SSTC2D J163131.2-242627 (Oph163131, i~84deg). We present new high angular resolution continuum and 12CO ALMA observations of Oph163131. The gas emission appears significantly more extended in the vertical and radial direction compared to the dust emission, consistent with vertical settling and possibly radial drift. In addition, the new continuum observations reveal two clear rings. The outer ring, located at ~100 au, is well resolved in the observations, which allows us to put stringent constraints on the vertical extent of millimeter dust particles. We model the disk using radiative transfer and find that the scale height of millimeter sized grains is 0.5au or less at 100au from the central star. This value is about one order of magnitude smaller than the scale height of smaller micron-sized dust grains constrained by previous modeling, which implies that efficient settling of the large grains is occurring in the disk. When adopting a parametric dust settling prescription, we find that the observations are consistent with a turbulent viscosity coefficient of about alpha<=10^-5 at 100au. Finally, we find that the thin dust scale height measured in Oph163131 is favorable for planetary growth by pebble accretion: a 10 M_Earth planet may grow within less than 10 Myr, even in orbits exceeding 50au., Comment: Accepted for publication in ApJ

In this chapter of the Protostars and Planets VII, we review the breakthrough progress that has been made in the field of high-resolution, high-contrast optical and near-infrared imaging of planet-forming disks. These advancements include the direct detection of protoplanets embedded in some disks, and derived limits on planetary masses in others. Morphological substructures, including: rings, spirals, arcs, and shadows, are seen in all imaged infrared-bright disks to date, and are ubiquitous across spectral types. These substructures are believed to be the result of disk evolution processes, and in particular disk-planet interactions. Since small dust grains that scatter light are tightly bound to the disk's gas, these observations closely trace disk structures predicted by hydrodynamical models and serve as observational tests of the predictions of planet formation theories. We argue that the results of current and next-generation high-contrast imaging surveys will, when combined with complementary data from ALMA, lead to a much deeper understanding of the co-evolution of disks and planets, and the mechanisms by which planets form., Comment: Review Chapter for Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Accepted version, before community feedback. 31 pages (37 with references and appendix), 17 figures. This Chapter is dedicated to our co-author Willy Kley. An obituary written by Richard Nelson can be found at the end of the Chapter

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

For several transition disks (TDs), dark regions interpreted as shadows have been observed in scattered light imaging and are hypothesized to originate from misalignments between distinct disk regions. We aim to investigate the presence of misalignments in TDs. We study the inner disk geometries of 20 well-known transition disks with VLTI/GRAVITY observations and use complementary $^{12}$CO and $^{13}$CO molecular line data from ALMA to derive the orientation of the outer disk regions. We fit simple models to the GRAVITY data to derive the inner disks inclination and position angles. The outer disk geometries were derived from Keplerian fits to the ALMA velocity maps and compared to the inner disk constraints. We also predicted the locations of shadows for significantly misaligned systems. Our analysis reveals six disks to exhibit significant misalignments between their inner and outer disks. The predicted shadow positions agree well with the scattered light images of HD100453 and HD142527, and we find supporting evidence for a shadow in the disk around CQ Tau. In the other three targets for which we infer significantly misaligned disks, V1247 Ori, V1366 Ori, and RY Lup, we do not see any evident sign of shadows in the scattered light images. The scattered light shadows observed in DoAr44, HD135344B, and HD139614 are consistent with our observations, yet the underlying morphology is likely too complex to be described by our models and the accuracy achieved by our observations. Whereas we can derive precise constraints on the potential shadow positions for well-resolved inner disks around HAeBe stars, the statistical uncertainties for the marginally resolved inner disks around the TTS of our sample make it difficult to extract conclusive constraints for the presence of shadows in these systems., Comment: Accepted for publication in Astronomy \& Astrophysics

To understand how the multitude of planetary systems that have been discovered come to be, we need to study systems at different evolutionary stages, with different central stars but also in different environments. The most challenging environment for planet formation may be the harsh UV radiation field of nearby massive stars which quickly erodes disks by external photo-evaporation. We have observed the AT Pyx system, located in the head of a cometary globule in the Gum Nebula, to search for signs of ongoing planet formation. We used the extreme adaptive optics imager VLT/SPHERE to observe AT Pyx in polarized light as well as total intensity in the J, H and K-band. Additionally we employed VLT/NACO to observe the system in the L-band. We resolve the disk around AT Pyx in scattered light across multiple wavelengths. We find an extended (>126 au) disk, with an intermediate inclination between 35 deg and 42 deg. The disk shows complex sub-structure and we identify 2 and possibly 3 spiral-like features. Depending on the precise geometry of the disk (which we can not unambiguously infer from our data) the disk may be eccentric with an eccentricity of ~0.16 or partially self-shadowed. The spiral features and possible eccentricity are both consistent with signatures of an embedded gas giant planet equal in mass to Jupiter. Our own observations can rule out brown dwarf companions embedded in the resolved disk, but are not sensitive enough to detect gas giants. AT Pyx is the first disk in a cometray globule in the Gum Nebula which is spatially resolved. By comparison with disks in the Orion Nebula Cluster we note that the extension of the disk may be exceptional for this environment if the external UV radiation field is comparable to other cometary globules in the region. The signposts of ongoing planet formation are intriguing and need to be followed up with higher sensitivity., Comment: 11 pages, 9 figures, accepted for publication in A&A

To date, nearly two hundred planet-forming disks have been imaged with high resolution. Our propensity to study bright and extended objects is however biasing our view of the disk demography. In this work, we contribute to alleviate this bias by analyzing fifteen disks targeted with VLT/SPHERE that look faint in scattered light. Sources were selected based on a low far-IR excess from the spectral energy distribution. The comparison with the ALMA images available for a few sources shows that the scattered light surveyed by these datasets is only detected from a small portion of the disk extent. The mild anti-correlation between the disk brightness and the near-IR excess demonstrates that these disks are self-shadowed: the inner disk rim intercepts much starlight and leaves the outer disk in penumbra. Based on the uniform distribution of the disk brightness in scattered light across all spectral types, self-shadowing would act similarly for inner rims at a different distance from the star. We discuss how the illumination pattern of the outer disk may evolve with time. Some objects in the sample are proposed to be at an intermediate stage toward bright disks from the literature with either no shadow or with sign of azimuthally confined shadows., Comment: Accepted for publication in A\&A

These protoplanetary disks in T Tauri stars play a central role in star and planet formation. We spatially resolve at sub-au scales the innermost regions of a sample of T Tauri's disks to better understand their morphology and composition. We extended our homogeneous data set of 27 Herbig stars and collected near-IR K-band observations of 17 T Tauri stars, spanning effective temperatures and luminosities in the ranges of ~4000-6000 K and ~0.4-10 Lsun. We focus on the continuum emission and develop semi-physical geometrical models to fit the interferometric data and search for trends between the properties of the disk and the central star. The best-fit models of the disk's inner rim correspond to wide rings. We extend the Radius-luminosity relation toward the smallest luminosities (0.4-10 Lsun) and find the R~L^(1/2) trend is no longer valid, since the K-band sizes measured with GRAVITY are larger than the predicted sizes from sublimation radius computation. No clear correlation between the K-band half-flux radius and the mass accretion rate is seen. Having magnetic truncation radii in agreement with the K-band GRAVITY sizes would require magnetic fields as strong as a few kG, which should have been detected, suggesting that accretion is not the main process governing the location of the half-flux radius of the inner dusty disk. Our measurements agree with models that take into account the scattered light. The N-to-K band size ratio may be a proxy for disentangling disks with silicate features in emission from disks with weak and/or in absorption silicate features. When comparing inclinations and PA of the inner disks to those of the outer disks (ALMA) in nine objects of our sample, we detect misalignments for four objects., Comment: Accepted for publication in A&A; 20 pages, 15 figures, 9 tables

We aim to demonstrate that the presence and mass of an exoplanet can now be effectively derived from the astrometry of another exoplanet. We combined previous astrometry of $\beta$ Pictoris b with a new set of observations from the GRAVITY interferometer. The orbital motion of $\beta$ Pictoris b is fit using Markov chain Monte Carlo simulations in Jacobi coordinates. The inner planet, $\beta$ Pictoris c, was also reobserved at a separation of 96\,mas, confirming the previous orbital estimations. From the astrometry of planet b only, we can (i) detect the presence of $\beta$ Pictoris c and (ii) constrain its mass to $10.04^{+4.53}_{-3.10}\,M_{\rm Jup}$. If one adds the astrometry of $\beta$ Pictoris c, the mass is narrowed down to $9.15^{+1.08}_{-1.06}\,M_{\rm Jup}$. The inclusion of radial velocity measurements does not affect the orbital parameters significantly, but it does slightly decrease the mass estimate to $8.89^{+0.75}_{-0.75}\,M_{\rm Jup}$. With a semimajor axis of $2.68\pm0.02$\,au, a period of $1221\pm15$ days, and an eccentricity of $0.32\pm0.02$, the orbital parameters of $\beta$ Pictoris c are now constrained as precisely as those of $\beta$ Pictoris b. The orbital configuration is compatible with a high-order mean-motion resonance (7:1). The impact of the resonance on the planets' dynamics would then be negligible with respect to the secular perturbations, which might have played an important role in the eccentricity excitation of the outer planet., Comment: Accepted in A&A Letter

The formation and evolution of planetary systems impact the primordial accretion disk. HD141569 is the only known pre-main sequence star characterized by a hybrid disk. Observations probed the outer-disk structure showing a complex system of rings and interferometric observations attempted to characterize its inner 5 au region, but derived limited constraints. The goal of this work was to explore with new high-resolution interferometric observations the properties of the dust and gas in the internal regions of HD141569. We observed HD141569 on mas scales with GRAVITY/VLTI in the near-infrared at low and high spectral resolution. We interpreted the visibilities and spectral energy distribution with geometrical models and radiative transfer techniques to constrain the dust emission. We analyzed the high spectral resolution quantities to investigate the properties of the Br-Gamma line emitting region. Thanks to the combination of three different epochs, GRAVITY resolves the inner dusty disk in the K band. Data modeling shows that an IR excess of about 6% is spatially resolved and that the origin of this emission is confined in a ring of material located at a radius of 1 au from the star with a width smaller than 0.3 au. The MCMax modeling suggests that this emission could originate from a small amount of QHPs, while large silicate grain models cannot reproduce at the same time the observational constraints on the properties of near-IR and mid-IR fluxes. The differential phases in the Br-Gamma line clearly show an S-shape that can be best reproduced witha gas disk in Keplerian rotation, confined within 0.09 au. This is also hinted at by the double-peaked Br-Gamma emission line shape. The modeling of the continuum and gas emission shows that the inclination and position angle of these two components are consistent with a system showing relatively coplanar rings on all scales., Comment: Accepted for publication in A&A; 25 pages, 15 figures, 5 tables

[Full abstract in the paper] In recent years, protoplanetary disks with spiral structures have been detected in scattered light, millimeter continuum, and CO gas emission. The mechanisms causing these structures are still under debate. A popular scenario to drive the spiral arms is the one of a planet perturbing the material in the disk. However, if the disk is massive, gravitational instability is usually the favored explanation. Multiwavelength studies could be helpful to distinguish between the two scenarios. So far, only a handful of disks with spiral arms have been observed in both scattered light and millimeter continuum. We aim to perform an in-depth characterization of the protoplanetary disk morphology around WaOph 6 analyzing data obtained at different wavelengths, as well as to investigate the origin of the spiral features in the disk. We present the first near-infrared polarimetric observations of WaOph 6 obtained with SPHERE at the VLT and compare them to archival millimeter continuum ALMA observations. We traced the spiral features in both data sets and estimated the respective pitch angles. We discuss the different scenarios that can give rise to the spiral arms in WaOph 6. We tested the planetary perturber hypothesis by performing hydrodynamical and radiative transfer simulations to compare them with scattered light and millimeter continuum observations., Comment: Accepted for publication in A&A. 17 pages, 9 figures in the main text and 3 figures in the Appendix

Context. Spiral arms in protoplanetary disks could be shown to be the manifestation of density waves launched by protoplanets and propagating in the gaseous component of the disk. At least two point sources have been identified in the L band in the MWC 758 system as planetary mass object candidates. Aims. We used VLT/SPHERE to search for counterparts of these candidates in the H and K bands, and to characterize the morphology of the spiral arms . Methods. The data were processed with now-standard techniques in high-contrast imaging to determine the limits of detection, and to compare them to the luminosity derived from L band observations. Results. In considering the evolutionary, atmospheric, and opacity models we were not able to confirm the two former detections of point sources performed in the L band. In addition, the analysis of the spiral arms from a dynamical point of view does not support the hypothesis that these candidates comprise the origin of the spirals. Conclusions. Deeper observations and longer timescales will be required to identify the actual source of the spiral arms in MWC 758., Comment: Accepted for publication in Astronomy and Astrophysics

Protoplanetary disks drive some of the formation process (e.g., accretion, gas dissipation, formation of structures, etc.) of stars and planets. Understanding such physical processes is one of the main astrophysical questions. HD 163296 is an interesting young stellar object for which infrared and sub-millimeter observations have shown a prominent circumstellar disk with gaps plausibly created by forming planets. This study aims at characterizing the morphology of the inner disk in HD 163296 with multi-epoch near-infrared interferometric observations performed with GRAVITY at the Very Large Telescope Interferometer (VLTI). Our goal is to depict the K-band (lambda_0 ~ 2.2 um) structure of the inner rim with milliarcsecond (sub-au) angular resolution. Our data is complemented with archival PIONIER (H-band; lambda_0 ~ 1.65 um) data of the source. We performed a Gradient Descent parametric model fitting to recover the sub-au morphology of our source. Our analysis shows the existence of an asymmetry in the disk surrounding the central star of HD 163296. We confirm variability of the disk structure in the inner ~2 mas (0.2 au). While variability of the inner disk structure in this source has been suggested by previous interferometric studies, this is the first time that it is confirmed in the H- and K-bands by using a complete analysis of the closure phases and squared visibilities over several epochs. Because of the separation from the star, position changes, and persistence of this asymmetric structure on timescales of several years, we argue that it is a dusty feature (e.g., a vortex or dust clouds), probably, made by a mixing of sillicate and carbon dust and/or refractory grains, inhomogeneously distributed above the mid-plane of the disk., Comment: Accepted to be published in Astronomy and Astrophysics; main-body: 11 pages, 3 figures and 3 tables

Context. Characterizing the evolution of protoplanetary disks is necessary to improve our understanding of planet formation. Constraints on both dust and gas are needed to determine the dominant disk dissipation mechanisms. Aims. We aim to compare the disk dust masses in the Chamaeleon II (Cha II) star-forming region with other regions with ages between 1 and 10Myr. Methods. We use ALMA band 6 observations (1.3 mm) to survey 29 protoplanetary disks in Cha II. Dust mass estimates are derived from the continuum data. Results. Out of our initial sample of 29 disks, we detect 22 sources in the continuum, 10 in 12CO, 3 in 13CO, and none in C18O (J=2-1). Additionally, we detect two companion candidates in the continuum and 12CO emission. Most disk dust masses are lower than 10Mearth, assuming thermal emission from optically thin dust. We compare consistent estimations of the distributions of the disk dust mass and the disk-to-stellar mass ratios in Cha II with six other low mass and isolated star-forming regions in the age range of 1-10Myr: Upper Sco, CrA, IC 348, Cha I, Lupus, and Taurus. When comparing the dust-to-stellar mass ratio, we find that the masses of disks in Cha II are statistically different from those in Upper Sco and Taurus, and we confirm that disks in Upper Sco, the oldest region of the sample, are statistically less massive than in all other regions. Performing a second statistical test of the dust mass distributions from similar mass bins, we find no statistical differences between these regions and Cha II. Conclusions. We interpret these trends, most simply, as a sign of decline in the disk dust masses with time or dust evolution. Different global initial conditions in star-forming regions may also play a role, but their impact on the properties of a disk population is difficult to isolate in star-forming regions lacking nearby massive stars., Comment: Accepted for publication in A&A

We aim to reveal the nature of the reddest known substellar companion HD 206893 B by studying its near-infrared colors and spectral morphology and by investigating its orbital motion. We fit atmospheric models for giant planets and brown dwarfs and perform spectral retrievals with petitRADTRANS and ATMO on the observed GRAVITY, SPHERE, and GPI spectra of HD 206893 B. To recover its unusual spectral features, we include additional extinction by high-altitude dust clouds made of enstatite grains in the atmospheric model fits. We also infer the orbital parameters of HD 206893 B by combining the $\sim 100~\mu\text{as}$ precision astrometry from GRAVITY with data from the literature and constrain the mass and position of HD 206893 C based on the Gaia proper motion anomaly of the system. The extremely red color and the very shallow $1.4~\mu\text{m}$ water absorption feature of HD 206893 B can be fit well with the adapted atmospheric models and spectral retrievals. Altogether, our analysis suggests an age of $\sim 3$-$300~\text{Myr}$ and a mass of $\sim 5$-$30~\text{M}_\text{Jup}$ for HD 206893 B, which is consistent with previous estimates but extends the parameter space to younger and lower-mass objects. The GRAVITY astrometry points to an eccentric orbit ($e = 0.29^{+0.06}_{-0.11}$) with a mutual inclination of $< 34.4~\text{deg}$ with respect to the debris disk of the system. While HD 206893 B could in principle be a planetary-mass companion, this possibility hinges on the unknown influence of the inner companion on the mass estimate of $10^{+5}_{-4}~\text{M}_\text{Jup}$ from radial velocity and Gaia as well as a relatively small but significant Argus moving group membership probability of $\sim 61\%$. However, we find that if the mass of HD 206893 B is $< 30~\text{M}_\text{Jup}$, then the inner companion HD 206893 C should have a mass between $\sim 8$-$15~\text{M}_\text{Jup}$., Comment: 28 pages, 19 figures, accepted for publication in A&A

Stars form by accreting material from their surrounding disks. There is a consensus that matter flowing through the disk is channelled onto the stellar surface by the stellar magnetic field. This is thought to be strong enough to truncate the disk close to the so-called corotation radius where the disk rotates at the same rate as the star. Spectro-interferometric studies in young stellar objects show that Hydrogen is mostly emitted in a region of a few milliarcseconds across, usually located within the dust sublimation radius. Its origin is still a matter of debate and it can be interpreted as coming from the stellar magnetosphere, a rotating wind or a disk. In the case of intermediate-mass Herbig AeBe stars, the fact that the Br gamma emission is spatially resolved rules out that most of the emission comes from the magnetosphere. This is due to the weak magnetic fields (some tenths of G) detected in these sources, resulting in very compact magnetospheres. In the case of T Tauri sources, their larger magnetospheres should make them easier to resolve. However, the small angular size of the magnetosphere (a few tenths of milliarcseconds), along with the presence of winds emitting in Hydrogen make the observations interpretation challenging. Here, we present direct evidence of magnetospheric accretion by spatially resolving the inner disk of the 60 pc T Tauri star TW Hydrae through optical long baseline interferometry. We find that the hydrogen near-infrared emission comes from a region approximately 3.5 stellar radii (R*) across. This region is within the continuum dusty disk emitting region (Rcont = 7 R*) and smaller than the corotation radius which is twice as big. This indicates that the hydrogen emission originates at the accretion columns, as expected in magnetospheric accretion models, rather than in a wind emitted at much larger distance (>1au)., Comment: Publish in Nature

To determine the origin of the spiral structure observed in the dust continuum emission of Elias 2-27 we analyze multi-wavelength continuum ALMA data with a resolution of $\sim$0.2 arcsec ($\sim$23au) at 0.89, 1.3 and 3.3mm. We also study the kinematics of the disk with $^{13}$CO and C$^{18}$O ALMA observations in the $J=$3-2 transition. The spiral arm morphology is recovered at all wavelengths in the dust continuum observations, where we measure contrast and spectral index variations along the spiral arms and detect subtle dust-trapping signatures. We determine that the emission from the midplane is cold and interpret the optical depth results as signatures of a higher disk mass than previous constraints. From the gas data, we search for deviations from Keplerian motion and trace the morphology of the emitting surfaces and the velocity profiles. We find an azimuthally varying emission layer height in the system, large-scale emission surrounding the disk, and strong perturbations in the channel maps, co-located with the spirals. Additionally, we develop multigrain dust and gas SPH simulations of a gravitationally unstable disk and compare them to the observations. Given the large scale emission and highly perturbed gas structure, together with the comparison of continuum observations to theoretical predictions, we propose infall-triggered gravitational instabilities as origin for the observed spiral structure., Comment: Accepted for publication in ApJ, 38 pages, 23 Figures

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

Observations of protoplanetary disks around very low-mass stars and brown dwarfs remain challenging and little is known about their properties. The disk around CIDA1 ($\sim$0.1-0.2$M_\odot$) is one of the very few known disks that host a large cavity (20au radius in size) around a very low-mass star. We present new ALMA observations at Band7 (0.9mm) and Band4 (2.1mm) of CIDA1 with a resolution of $\sim 0.05''\times 0.034''$. These new ALMA observations reveal a very bright and unresolved inner disk, a shallow spectral index of the dust emission ($\sim2$), and a complex morphology of a ring located at 20au. We also present X-Shooter (VLT) observations that confirm the high accretion rate of CIDA1 of $\dot{M}_{\rm acc}$=1.4 $\times~10^{-8}M_\odot$/yr. This high value of $\dot{M}_{\rm acc}$, the observed inner disk, and the large cavity of 20au exclude models of photo-evaporation to explain the observed cavity. When comparing these observations with models that combine planet-disk interaction, dust evolution, and radiative transfer, we exclude planets more massive than 0.5$M_{\rm{Jup}}$ as the potential origin of the large cavity because with these it is difficult to maintain a long-lived and bright inner disk. Even in this planet mass regime, an additional physical process may be needed to stop the particles from migrating inwards and to maintain a bright inner disk on timescales of millions of years. Such mechanisms include a trap formed by a very close-in extra planet or the inner edge of a dead zone. The low spectral index of the disk around CIDA1 is difficult to explain and challenges our current dust evolution models, in particular processes like fragmentation, growth, and diffusion of particles inside pressure bumps., Comment: Accepted for publication in A&A, minor changes after language edition