Search

In the early stages of planetary system formation, young exoplanets gravitationally interact with their surrounding environments and leave observable signatures on protoplanetary disks. Among these structures, a pair of nearly symmetric spiral arms can be driven by a giant protoplanet. For the double-spiraled SAO 206462 protoplanetary disk, we obtained three epochs of observations spanning 7 yr using the Very Large Telescope's SPHERE instrument in near-infrared $J$-band polarized light. By jointly measuring the motion of the two spirals at three epochs, we obtained a rotation rate of $-0.85^\circ\pm0.05^\circ~{\rm yr}^{-1}$. This rate corresponds to a protoplanet at $66\pm3$ au on a circular orbit dynamically driving both spirals. The derived location agrees with the gap in ALMA dust-continuum observations, indicating that the spiral driver may also carve the observed gap. What is more, a dust filament at $\sim$63 au observed by ALMA coincides with the predicted orbit of the spiral-arm-driving protoplanet. This double-spiraled system is an ideal target for protoplanet imaging., Comment: 11 pages, 3 figures. Invited paper accepted to special issue (https://www.mdpi.com/journal/universe/special_issues/Y3T2Z3J1HS) of Universe. Data in ancillary folder

We present near- and mid-infrared (IR) broadband imaging observations of the edge-on protoplanetary disk around HH 30 with the James Webb Space Telescope/Near Infrared Camera (NIRCam) and the Mid-Infrared Instrument (MIRI). We combine these observations with archival optical/near-IR scattered light images obtained with the Hubble Space Telescope (HST) and a millimeter-wavelength dust continuum image obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) with the highest spatial resolution ever obtained for this target. Our multiwavelength images clearly reveal the vertical and radial segregation of micron-sized and sub-mm-sized grains in the disk. In the near- and mid-IR, the images capture not only bi-reflection nebulae separated by a dark lane but also diverse dynamical processes occurring in the HH 30 disk, such as spiral- and tail-like structures, a conical outflow, and a collimated jet. In contrast, the ALMA image reveals a flat dust disk in the disk midplane. By performing radiative transfer simulations, we show that grains of about 3 $\mu$m in radius or larger are fully vertically mixed to explain the observed mid-IR scattered light flux and its morphology, whereas millimeter-sized grains are settled into a layer with a scale height of $\gtrsim1$ au at $100$ au from the central star. We also find a tension in the disk inclination angle inferred from optical/near-IR and mm observations with the latter being closer to an exactly edge-on. Finally, we report the first detection of the proper motion of an emission knot associated with the mid-IR collimated jet detected by combining two epochs of our MIRI 12.8-$\mu$m observations., Comment: 29 pages, 23 figures; Accepted for publication in ApJ

PDS 453 is a rare highly inclined disk where the stellar photosphere is seen at grazing incidence on the disk surface. Our goal is take advantage of this geometry to constrain the structure and composition of this disk, in particular the fact that it shows a 3.1 $\mu$m water ice band in absorption that can be related uniquely to the disk. We observed the system in polarized intensity with the VLT/SPHERE instrument, as well as in polarized light and total intensity using the HST/NICMOS camera. Infrared archival photometry and a spectrum showing the water ice band are used to model the spectral energy distribution under Mie scattering theory. Based on these data, we fit a model using the radiative transfer code MCFOST to retrieve the geometry and dust and ice content of the disk. PDS 453 has the typical morphology of a highly inclined system with two reflection nebulae where the disk partially attenuates the stellar light. The upper nebula is brighter than the lower nebula and shows a curved surface brightness profile in polarized intensity, indicating a ring-like structure. With an inclination of 80{\deg} estimated from models, the line-of-sight crosses the disk surface and a combination of absorption and scattering by ice-rich dust grains produces the water ice band. PDS 453 is seen highly inclined and is composed of a mixture of silicate dust and water ice. The radial structure of the disk includes a significant jump in density and scale height at a radius of 70 au in order to produce a ring-like image. The depth of the 3.1 $\mu$m water ice band depends on the amount of water ice, until it saturates when the optical thickness along the line-of-sight becomes too large. Therefore, quantifying the exact amount of water from absorption bands in edge-on disks requires a detailed analysis of the disk structure and tailored radiative transfer modeling., Comment: 11 pages, 11 figures

We present JWST broadband images of the highly inclined protoplanetary disk SSTc2d J163131.2-242627 (Oph163131) from 2.0 to 21$\mu$m. The images show a remarkable evolution in disk structure with wavelength, quite different from previous JWST observations of other edge-on disks. At 2.0 and 4.4$\mu$m, Oph163131 shows two scattering surfaces separated by a dark lane, typical of highly inclined disks. Starting at 7.7$\mu$m however, 1) the two linear nebulosities flanking the dark lane disappear; 2) the brighter nebula tracing the disk upper surface transitions into a compact central source distinctly larger than the JWST PSF and whose intrinsic size increases with wavelength; and 3) patches of extended emission appear at low latitudes, and at surprisingly large radii nearly twice that of the scattered light seen with $HST$ and NIRCam, and of the gas. We interpret the compact central source as thermal emission from the star and the inner disk that is not seen directly, but which instead is able to progressively propagate to greater distances at longer wavelengths. The lack of sharp-edged structures in the extended patchy emission argues against the presence of shocks and suggests photoexcitation or stochastic heating of material smoothly flowing away from the star along the disk surface. Finally, the dark lane thickness decreases significantly between 0.6$\mu$m and 4.4$\mu$m which indicates that the surface layers of Oph163131 lack grains larger than 1$\mu$m., Comment: Accepted for publication in the Astrophysical Journal. 11 pages, 7 figures

Molecular deuteration is a powerful diagnostic tool for probing the physical conditions and chemical processes in astrophysical environments. In this work, we focus on formaldehyde deuteration in the protobinary system NGC\,1333 IRAS\,4A, located in the Perseus molecular cloud. Using high-resolution ($\sim$\,100\,au) ALMA observations, we investigate the [D$_2$CO]/[HDCO] ratio along the cavity walls of the outflows emanating from IRAS\,4A1. Our analysis reveals a consistent decrease in the deuteration ratio (from $\sim$\,60-20\% to $\sim$\,10\%) with increasing distance from the protostar (from $\sim$\,2000\,au to $\sim$\,4000\,au). Given the large measured [D$_2$CO]/[HDCO], both HDCO and D$_2$CO are likely injected by the shocks along the cavity walls into the gas-phase from the dust mantles, formed in the previous prestellar phase. We propose that the observed [D$_2$CO]/[HDCO] decrease is due to the density profile of the prestellar core from which NGC\,1333 IRAS\,4A was born. When considering the chemical processes at the base of formaldehyde deuteration, the IRAS\,4A's prestellar precursor had a predominantly flat density profile within 3000\,au and a decrease of density beyond this radius.

The interplay between T Tauri stars and their circumstellar disks, and how this impacts the onset of planet formation has yet to be established. We studied a seemingly old T Tauri star, PDS 111, and its disk. We analyzed optical, infrared, and sub-millimeter observations obtained with VLT/X-shooter, Mercator/HERMES, TESS, VLT/SPHERE, and ALMA, providing a new view on PDS 111 and its protoplanetary disk. The multi-epoch spectroscopy yields photospheric lines to classify the star, and emission lines to study variability in the hot inner disk and to determine the mass-accretion rate. The SPHERE and ALMA observations are used to characterize the dust distribution of the small and large grains, respectively. PDS 111 is a weak-line T Tauri star with spectral type G2, exhibits strong H$\alpha$ variability and with a low mass-accretion rate of $1-5\times10^{-10}$\,M$_{\odot}$\,yr$^{-1}$. We measured an age of the system of 15.9$^{+1.7}_{-3.7}$ Myr using pre-main sequence tracks. The SPHERE observations show a strongly flaring disk with an asymmetric substructure. The ALMA observations reveal a 30 au cavity in the dust continuum emission with a low contrast asymmetry in the South-West of the disk and a dust disk mass of 45.8\,$M_\oplus$. The $^{12}$CO radial extension is at least three times larger than that of the dust emission. Although the measured age is younger than suggested in literature, PDS 111 still seems relatively old; this provides insight into disk properties at an advanced stage of pre-main sequence evolution. The characteristics of this disk are very similar to its younger counterparts: strongly flaring, an average disk mass, a typical radial extent of the disk gas and dust, and the presence of common substructures. This suggests that disk evolution has not significantly changed the disk properties. These results show similarities with the "Peter Pan disks" around M-dwarfs., Comment: 23 pages, 23 figures, accepted by A&A (abstract shortened)

The exploration of outflows in protobinary systems presents a challenging yet crucial endeavour, offering valuable insights into the dynamic interplay between protostars and their evolution. In this study, we examine the morphology and dynamics of jets and outflows within the IRAS\,4A protobinary system. This analysis is based on ALMA observations of SiO(5--4), H$_2$CO(3$_{0,3}$--2$_{0,3}$), and HDCO(4$_{1,4}$--3$_{1,3}$) with a spatial resolution of $\sim$150\,au. Leveraging an astrochemical approach involving the use of diverse tracers beyond traditional ones has enabled the identification of novel features and a comprehensive understanding of the broader outflow dynamics. Our analysis reveals the presence of two jets in the redshifted emission, emanating from IRAS\,4A1 and IRAS\,4A2, respectively. Furthermore, we identify four distinct outflows in the region for the first time, with each protostar, 4A1 and 4A2, contributing to two of them. We characterise the morphology and orientation of each outflow, challenging previous suggestions of bends in their trajectories. The outflow cavities of IRAS\,4A1 exhibit extensions of 10$''$ and 13$''$ with position angles (PA) of 0$^{\circ}$ and -12$^{\circ}$, respectively, while those of IRAS\,4A2 are more extended, spanning 18$''$ and 25$''$ with PAs of 29$^{\circ}$ and 26$^{\circ}$. We propose that the misalignment of the cavities is due to a jet precession in each protostar, a notion supported by the observation that the more extended cavities of the same source exhibit lower velocities, indicating they may stem from older ejection events.

V892 Tau is a young binary star surrounded by a circumbinary disc which show hints of interaction with the low-mass nearby star V892 Tau NE. The goal of this paper is to constrain the orbit of V892 Tau NE and to determine the resulting circumbinary disc dynamics. We present new ALMA observations of the V892 Tau circumbinary disc at a twice higher angular and spectral resolution. We model the data with V892 Tau as a triple system and perform a grid of hydrodynamical simulations testing several orbits of the companion. The simulation outputs are then post-processed to build synthetic maps that we compare to the observations. The 12CO emission of the disc shows clear non-Keplerian features such as spiral arms. When comparing the data with our synthetic observations, we interpret these features as ongoing interactions with the companion. Our simulations indicate that an eccentricity of 0.5 of the companion is needed to reproduce the observed disc extent and that a mutual inclination of approximately 60{\deg} with the inner binary reproduces the measured disc tilt. In order to explain most of the features of the circumbinary disc, we propose that V892 Tau NE follows a misaligned eccentric orbit, with an eccentricity between 0.2 and 0.5 and a mutual inclination between 30{\deg} and 60{\deg}. Such a misaligned companion suggests the disc is oscillating and precessing with time, stabilising in an intermediate plane with a non-zero mutual inclination with the inner binary. Given that orbital configuration, we show that the stability of future planets is compromised in the second half of the disc once the gas has dissipated., Comment: Published in A&A, 29 pages, 15 figures

GG Tau is one of the most studied multiple young stellar systems: GG Tau A is a hierarchical triple surrounded by a massive disc and its companion, GG Tau B, is also a binary. Despite numerous observational attempts, an understanding of the geometry of the GG Tau A system is still elusive. We provide new astrometric measures of the system and we run a set of hydrodynamical simulations with two representative orbits to test how they impact a disc composed of dust and gas. We test the dynamical evolution of the two scenarios on short and long timescales. We obtain synthetic flux emission from our simulations and we compare them with 1300 $\mu$m ALMA dust continuum emission and 1.67 $\mu$m SPHERE dust scattering images to infer the most likely orbital arrangement. We extend the analysis of the binary orbital parameters using six new epochs from archival data, showing that the current measurements alone are not capable of breaking the degeneracy between families of coplanar and misaligned orbits. We found that the time-scale for the onset of the disc eccentricity growth, $\tau_{ecc}$, is a fundamental time-scale for the morphology of the system. Results from numerical simulations show that the best match between is obtained with the misaligned configuration ($\Delta\theta= 30^\circ$) on timescales shorter than $\tau_{ecc}$. The results exhibit an almost circular cavity and dust ring. However, for both scenarios, the cavity size and its eccentricity quickly grow for timescales longer than $\tau_{ecc}$ and the models do not reproduce the observed morphology anymore. This implies that either the age of the system is shorter than $\tau_{ecc}$ or that the disc eccentricity growth is not triggered or dissipated. This finding raises questions on the future evolution of the GG Tau A system and, more in general, on the time evolution of eccentric binaries and their circumbinary discs., Comment: 14 pages, 12 figures

High-order multiple (triple and beyond) systems are relatively common. Their interaction with circumstellar and circumbinary material can have a large impact on the formation and evolution of planetary systems and depends on their orbital properties. GG\,Tau and UX\,Tau are two pre-main sequence high-order multiple systems in which the tightest pair has a projected separation of $\approx5$--20\,au. Characterizing precisely their orbits is crucial to establish their long-term stability, to predict the dynamics and evolution of circumstellar matter, and to evaluate the potential for planet formation in such systems. We combine existing astrometric measurements with previously unpublished high-resolution observations of the GG\,Tau\,Ab and UX\,Tau\,B pairs and perform Keplerian orbital fits. For GG\,Tau\,Ab the data presented here represent the first detection of orbital motion. For both systems they yield dramatic increases in orbital coverage ($\gtrsim60\%$ and $\approx100\%$ for UX\,Tau\,B and GG\,Tau\,Ab, for orbital periods of $\approx32$ and $\approx8$\,yr, respectively) and allow us to obtain well-constrained orbital fits, including dynamical masses with $\lesssim10\%$ and $\lesssim7\%$ random and systematic uncertainties. We find that both GG\,Tau\,A and UX\,Tau\,A--B likely form stable hierarchical systems, although one possible deprojection solution for GG\,Tau is strongly misaligned and could experience von Zeipel-Lidov-Kozai oscillations. We further find that the UX\,Tau\,B orbit is much more eccentric than the GG\,Tau\,Ab one, possibly explaining the lack of circumstellar material in the former. The newly-determined orbits revive the question of the dynamical fate of gas and dust in these two hierarchical systems and should spur new dedicated simulations to assess the long-term evolution of the systems and the dynamical perturbations imposed by the close binaries they host., Comment: Accepted for publication in Astronomy & Astrophysics

We present MIRI MRS observations of the large, multi-gapped protoplanetary disk around the T-Tauri star AS 209. The observations reveal hundreds of water vapor lines from 4.9 to 25.5 $\mu$m towards the inner $\sim1$ au in the disk, including the first detection of ro-vibrational water emission in this disk. The spectrum is dominated by hot ($\sim800$ K) water vapor and OH gas, with only marginal detections of CO$_2$, HCN, and a possible colder water vapor component. Using slab models with a detailed treatment of opacities and line overlap, we retrieve the column density, emitting area, and excitation temperature of water vapor and OH, and provide upper limits for the observable mass of other molecules. Compared to MIRI spectra of other T-Tauri disks, the inner disk of AS 209 does not appear to be atypically depleted in CO$_2$ nor HCN. Based on \textit{Spitzer IRS} observations, we further find evidence for molecular emission variability over a 10-year baseline. Water, OH, and CO$_2$ line luminosities have decreased by factors 2-4 in the new MIRI epoch, yet there are minimal continuum emission variations. The origin of this variability is yet to be understood., Comment: Accepted for publication in ApJ

Extreme adaptive optics instruments have revealed exquisite details on debris discs, allowing to extract the optical properties of the dust particles such as the phase function, the degree of polarisation and the spectral reflectance. These are three powerful diagnostic tools to understand the physical properties of the dust : the size, shape and composition of the dust particles. This can inform us on the population of parent bodies, also called planetesimals, which generate those particles through collisions. It is however very rare to be able to combine all those three observables for the same system, as this requires different high-contrast imaging techniques to suppress the starlight and reveal the faint scattered light emission from the dust. Due to its brightness, the ring detected around the A-type star HR 4796 is a notable exception, with both unpolarised and polarised images covering near-infrared wavelengths. Here, we show how measurements of dust particles in the laboratory can reproduce the observed near-infrared photo-polarimetric properties of the HR 4796 disc. Experimental characterisation of dust allows to bypass the current limitations of dust models to reproduce simultaneously the phase function, the degree of polarisation and the spectral reflectance., Comment: 4 pages, 4 figures, Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics

Polarisation is a powerful remote-sensing tool to study the nature of particles scattering the starlight. It is widely used to characterise interplanetary dust particles in the Solar System and increasingly employed to investigate extrasolar dust in debris discs' systems. We aim to measure the scattering properties of the dust from the debris ring around HD 181327 at near-infrared wavelengths. We obtained high-contrast polarimetric images of HD 181327 in the H band with the SPHERE / IRDIS instrument on the Very Large Telescope (ESO). We complemented them with archival data from HST / NICMOS in the F110W filter reprocessed in the context of the Archival Legacy Investigations of Circumstellar Environments (ALICE) project. We developed a combined forward-modelling framework to simultaneously retrieve the scattering phase function in polarisation and intensity. We detected the debris disc around HD 181327 in polarised light and total intensity. We measured the scattering phase function and the degree of linear polarisation of the dust at 1.6 micron in the birth ring. The maximum polarisation is 23.6% +/- 2.6% and occurs between a scattering angle of 70 deg and 82 deg. We show that compact spherical particles made of a highly refractive and relatively absorbing material in a differential power-law size distribution of exponent $-3.5$ can simultaneously reproduce the polarimetric and total intensity scattering properties of the dust. This type of material cannot be obtained with a mixture of silicates, amorphous carbon, water ice, and porosity, and requires a more refracting component such as iron-bearing minerals. We reveal a striking analogy between the near-infrared polarisation of comets and that of HD 181327. The methodology developed here combining VLT/SPHERE and HST/NICMOS may be applicable in the future to combine the polarimetric capabilities of SPHERE with the sensitivity of JWST., Comment: 21 pages, 22 figures, accepted in Astronomy & Astrophysics

While the most exciting explanation of the observed dust asymmetries in protoplanetary disks is the presence of protoplanets, other mechanisms can also form the dust features. This paper presents dual-wavelength Atacama Large Millimeter/submillimeter Array (ALMA) observations of a large asymmetric dusty ring around the M-type star CIDA 9A. We detect a dust asymmetry in both 1.3 mm and 3.1 mm data. To characterize the asymmetric structure, a parametric model is used to fit the observed visibilities. We report a tentative azimuthal shift of the dust emission peaks between the observations at the two wavelengths. This shift is consistent with a dust trap caused by a vortex, which may be formed by an embedded protoplanet or other hydrodynamical instabilities, such as a dead zone. Deep high-spatial observations of dust and molecular gas are needed to constrain the mechanisms that formed the observed millimeter cavity and dust asymmetry in the protoplanetary disk around CIDA 9A., Comment: Accepted for publication in ApJ, 12 pages, 7 figures with appendix

Diverse protoplanetary disk morphology can result from planet-disk interaction, suggesting planetary presence. To date, most scattered light imaging campaigns have probed polarized light, which is only a fraction of the total light and not very sensitive to planets. To observe and characterize protoplanetary disk systems in the near-infrared in both polarized and total intensity light, we carried out an unprecedented study of scattering properties of disks, as well as of any planetary companions. Using SPHERE with star-hopping at the Very Large Telescope, we observed 29 disk hosts and their reference stars in $K_s$-band polarized light. We extracted disks in total intensity by adopting the data imputation concept with sequential non-negative matrix factorization (DI-sNMF). We obtained high-quality disk images in total intensity for 15 systems and in polarized light for 23. For well-recovered disks in polarized light and total intensity, we parameterized the polarization fraction phase functions using scaled beta distribution: the peak of polarization fraction tentatively correlates with the peak scattering angle, which could be reproduced using certain compact dust, yet more detailed modeling studies are needed. We investigated the empirical DI-sNMF detectability of disks using logistic regression: total intensity detectability of disks primarily depends on host star brightness. For disks with SPHERE data in $Y$-/$J$-/$H$-band, we summarized their polarized color at ~90 deg scattering angle: most of disks are blue in polarized $J-K_s$ color, and they are relatively redder as stellar luminosity increases, indicating larger scatterers. High-quality disk imagery in both total intensity and polarized light thus allows for disk characterization in polarization fraction, and reduces the confusion between disk and planetary signals., Comment: 25 pages, 16 figures, 3 tables, A&A accepted. Data files in FITS format will be publicly available

Context: Circumstellar debris disks provide insight into the formation and early evolution of planetary systems. Resolved belts in particular help to locate planetesimals in exosystems, and can hint at the presence of disk-sculpting exoplanets. Aims: We study the circumstellar environment of HD 112810 (HIP 63439), a mid-F type star in the Sco-Cen association with a significant infrared excess indicating the presence of a circumstellar debris disk. Methods: We collected five high-contrast observations of HD 112810 with VLT/SPHERE. We identified a debris disk in scattered light, and found that the debris signature is robust over a number of epochs and a variety of reduction techniques. We modelled the disk, accounting for self-subtraction and assuming that it is optically thin. Results: We find a single-belt debris disk, with a radius of 118$\pm$9au and an inclination angle of ${75.7}^{+1.1}_{-1.3}$$\deg$. This is in good agreement with the constraints from SED modelling and from a partially-resolved ALMA image of the system. No planets are detected, though planets below the detection limit ($\sim$2.6M$_\textrm{J}$ at a projected separation of 118au) could be present and could have contributed to sculpting the ring of debris. Conclusions: HD 112810 adds to the growing inventory of debris disks imaged in scattered light. The disk is faint, but the radius and the inclination of the disk are promising for follow-up studies of the dust properties., Comment: A&A accepted. 13 pages, 6 figures + appendix

Understanding the driving forces behind spiral arms in protoplanetary disks remains a challenge due to the faintness of young giant planets. MWC 758 hosts such a protoplanetary disk with a two-armed spiral pattern that is suggested to be driven by an external giant planet. We present new thermal infrared observations that are uniquely sensitive to redder (i.e., colder or more attenuated) planets than past observations at shorter wavelengths. We detect a giant protoplanet, MWC 758c, at a projected separation of ~100 au from the star. The spectrum of MWC 758c is distinct from the rest of the disk and consistent with emission from a planetary atmosphere with Teff = 500 +/- 100 K for a low level of extinction (AV<30), or a hotter object with a higher level of extinction. Both scenarios are commensurate with the predicted properties of the companion responsible for driving the spiral arms. MWC 758c provides evidence that spiral arms in protoplanetary disks can be caused by cold giant planets or by those whose optical emission is highly attenuated. MWC 758c stands out both as one of the youngest giant planets known, and also as one of the coldest and/or most attenuated. Furthermore, MWC 758c is among the first planets to be observed within a system hosting a protoplanetary disk., Comment: Published in Nature Astronomy

We study the kinematics of the AS 209 disk using the J=2-1 transitions of $^{12}$CO, $^{13}$CO, and C$^{18}$O. We derive the radial, azimuthal, and vertical velocity of the gas, taking into account the lowered emission surface near the annular gap at ~1.7 (200 au) within which a candidate circumplanetary disk-hosting planet has been reported previously. In $^{12}$CO and $^{13}$CO, we find a coherent upward flow arising from the gap. The upward gas flow is as fast as $150~{\rm m~s}^{-1}$ in the regions traced by $^{12}$CO emission, which corresponds to about 50% of the local sound speed or $6\%$ of the local Keplerian speed. Such an upward gas flow is difficult to reconcile with an embedded planet alone. Instead, we propose that magnetically driven winds via ambipolar diffusion are triggered by the low gas density within the planet-carved gap, dominating the kinematics of the gap region. We estimate the ambipolar Elsasser number, Am, using the HCO$^+$ column density as a proxy for ion density and find that Am is ~0.1 at the radial location of the upward flow. This value is broadly consistent with the value at which numerical simulations find ambipolar diffusion drives strong winds. We hypothesize the activation of magnetically-driven winds in a planet-carved gap can control the growth of the embedded planet. We provide a scaling relationship which describes the wind-regulated terminal mass: adopting parameters relevant to 100 au from a solar-mass star, we find the wind-regulated terminal mass is about one Jupiter mass, which may help explain the dearth of directly imaged super-Jovian-mass planets., Comment: This paper has been accepted for publication in the Astrophysical Journal (ApJ)

We report results of a spectropolarimetric and photometric monitoring of the weak-line T Tauri star LkCa 4 within the SPIRou Legacy Survey large programme, based on data collected with SPIRou at the Canada-France-Hawaii Telescope and the TESS space probe between October 2021 and January 2022. We applied Zeeman-Doppler Imaging to our spectropolarimetric and photometric data to recover a surface brightness distribution compatible with TESS photometry, as well as the large-scale magnetic topology of the star. As expected from the difference in wavelength between near-infrared and optical data, the recovered surface brightness distribution is less contrasted than the previously published one based on ESPaDOnS data, but still features mid-latitude dark and bright spots. The large-scale magnetic field is consistent in shape and strength with the one derived previously, with a poloidal component resembling a 2.2 kG dipole and a toroidal component reaching 1.4 kG and encircling the star at the equator. Our new data confirm that the surface differential rotation of LkCa 4 is about 10 times weaker than that of the Sun, and significantly different from zero. Using our brightness reconstruction and Gaussian Process Regression, we were able to filter the radial velocity activity jitter down to a precision of 0.45 and 0.38 km $\rm s^{-1}$ (from an amplitude of 6.10 km $\rm s^{-1}$), respectively, yielding again no evidence for a close-in massive planet orbiting the star., Comment: 16 pages, 18 figures, accepted for publication in MNRAS

We investigate the nature of an unusually faint member of the $\epsilon$ Cha Association ($D\sim100$ pc, age $\sim5$ Myr), the nearest region of star formation of age $<$8 Myr. This object, 2MASS J11550336-7919147 (2M1155$-$79B), is a wide ($\sim$580 AU) separation, comoving companion to low-mass (M3) $\epsilon$ Cha Association member 2MASS J11550485-7919108 (2M1155$-$79A). We present near-infrared spectra of both components, along with analysis of photometry from Gaia EDR3, 2MASS, VHS, and WISE. The near-IR spectrum of 2M1155$-$79B displays strong He I 1.083 emission, a sign of active accretion and/or accretion-driven winds from a circumstellar disk. Analysis of WISE archival data reveals that the mid-infrared excess previously associated with 2M1155$-$79A instead originates from the disk surrounding 2M1155$-$79B. Based on these results, as well as radiative transfer modeling of its optical/IR spectral energy distribution, we conclude that 2M1155$-$79B is most likely a young, late-M, star that is partially obscured by, and actively accreting from, a nearly edge-on circumstellar disk. This would place 2M1155$-$79B among the rare group of nearby ($D\lesssim100$ pc), young (age $<$10 Myr) mid-M stars that are orbited by and accreting from highly inclined protoplanetary disks. Like these systems, the 2M1155$-$79B system is a particularly promising subject for studies of star and planet formation around low-mass stars., Comment: 11 pages, 6 figures, accepted for publication in AJ

Reference-star differential imaging (RDI) is a promising technique in high-contrast imaging that is thought to be more sensitive to exoplanets and disks than angular differential imaging (ADI) at short angular separations (i.e., <0.3"). However, it is unknown whether the performance of RDI on ground-based instruments can be improved by using all the archival data to optimize the subtraction of stellar contributions. We characterize the performance of RDI on SPHERE/IRDIS data in direct imaging of exoplanets and disks. We made use of all the archival data in H23 obtained by SPHERE/IRDIS in the past five years to build a master reference library and perform RDI. In the point-source detection, RDI can outperform ADI at small angular separations (<0.4") if the observing conditions are around the median conditions of our master reference library. On average, RDI has a gain of ~0.8 mag over ADI at 0.15" separation for observations under median conditions. We demonstrate that including more reference targets in the master reference library can indeed help to improve the performance of RDI. In disk imaging, RDI can reveal more disk features and provide a more robust recovery of the disk morphology. We resolve 33 disks in total intensity (19 planet-forming disks and 14 debris disks), and 4 of them can only be detected with RDI. Two disks are resolved in scattered light for the first time. Three disks are detected in total intensity for the first time. The master reference library we built in this work can be easily implemented into legacy or future SPHERE surveys to perform RDI, achieving better performance than that of ADI. To obtain optimal RDI gains over ADI, we recommend future observations be carried out under seeing conditions of 0.6"-0.8"., Comment: Accepted for publication in Astronomy and Astrophysics; 19 pages, 14 figures, 3 tables

We review the role of stellar flybys and encounters in shaping planet-forming discs around young stars, based on the published literature on this topic in the last 30 years. Since most stars $\leq~2$ Myr old harbour protoplanetary discs, tidal perturbations affect planet formation. First, we examine the probability of experiencing flybys or encounters: More than 50\% of stars with planet-forming discs in a typical star forming environment should experience a close stellar encounter or flyby within 1000 au. Second, we detail the dynamical effects of flybys on planet-forming discs. Prograde, parabolic, disc-penetrating flybys are the most destructive. Grazing and penetrating flybys in particular lead to the capture of disc material by the secondary to form a highly misaligned circumsecondary disc with respect to the disc around the primary. One or both discs may undergo extreme accretion and outburst events, similar to the ones observed in FU Orionis-type stars. Warps and broken discs are distinct signatures of retrograde flybys. Third, we review some recently observed stellar systems with discs where a stellar flyby or an encounter is suspected -- including UX Tau, RW Aur, AS 205, Z CMa, and FU Ori. Finally, we discuss the implications of stellar flybys for planet formation and exoplanet demographics, including possible imprints of a flyby in the Solar System in the orbits of trans-Neptunian objects and the Sun's obliquity., Comment: 16 pages, 7 figures, 1 table, 255 references. Invited review article accepted in EPJ+

We report the discovery of a circumplanetary disk (CPD) candidate embedded in the circumstellar disk of the T Tauri star AS 209 at a radial distance of about 200 au (on-sky separation of 1."4 from the star at a position angle of $161^\circ$), isolated via $^{13}$CO $J=2-1$ emission. This is the first instance of CPD detection via gaseous emission capable of tracing the overall CPD mass. The CPD is spatially unresolved with a $117\times82$ mas beam and manifests as a point source in $^{13}$CO, indicating that its diameter is $\lesssim14$ au. The CPD is embedded within an annular gap in the circumstellar disk previously identified using $^{12}$CO and near-infrared scattered light observations, and is associated with localized velocity perturbations in $^{12}$CO. The coincidence of these features suggests that they have a common origin: an embedded giant planet. We use the $^{13}$CO intensity to constrain the CPD gas temperature and mass. We find that the CPD temperature is $\gtrsim35$ K, higher than the circumstellar disk temperature at the radial location of the CPD, 22 K, suggesting that heating sources localized to the CPD must be present. The CPD gas mass is $\gtrsim 0.095 M_{\rm Jup} \simeq 30 M_{\rm Earth}$ adopting a standard $^{13}$CO abundance. From the non-detection of millimeter continuum emission at the location of the CPD ($3\sigma$ flux density $\lesssim26.4~\mu$Jy), we infer that the CPD dust mass is $\lesssim 0.027 M_{\rm Earth} \simeq 2.2$ lunar masses, indicating a low dust-to-gas mass ratio of $\lesssim9\times10^{-4}$. We discuss the formation mechanism of the CPD-hosting giant planet on a wide orbit in the framework of gravitational instability and pebble accretion., Comment: Accepted for publication in the ApJ Letters (July 7, 2022), 19 pages, 13 figures, interactive figures (Figure 7, 8, 9) are available at http://jaehanbae.com/as209/

While protoplanetary disks are often treated as isolated systems in planet formation models, observations increasingly suggest that vigorous interactions between Class II disks and their environments are not rare. DO Tau is a T Tauri star that has previously been hypothesized to have undergone a close encounter with the HV Tau system. As part of the DESTINYS ESO Large Programme, we present new VLT/SPHERE polarimetric observations of DO Tau and combine them with archival HST scattered light images and ALMA observations of CO isotopologues and CS to map a network of complex structures. The SPHERE and ALMA observations show that the circumstellar disk is connected to arms extending out to several hundred au. HST and ALMA also reveal stream-like structures northeast of DO Tau, some of which are at least several thousand au long. These streams appear not to be gravitationally bound to DO Tau, and comparisons with previous Herschel far-IR observations suggest that the streams are part of a bridge-like structure connecting DO Tau and HV Tau. We also detect a fainter redshifted counterpart to a previously known blueshifted CO outflow. While some of DO Tau's complex structures could be attributed to a recent disk-disk encounter, they might be explained alternatively by interactions with remnant material from the star formation process. These panchromatic observations of DO Tau highlight the need to contextualize the evolution of Class II disks by examining processes occurring over a wide range of size scales., Comment: 42 pages, 20 figures, accepted by ApJ, reduced data available at https://zenodo.org/record/6408903 (typos from v1 fixed)

We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circum-binary VLA 1623A disk as well as the VLA 1623B disk. We found that the minor axis of the circum-binary disk of VLA 1623A is misaligned by about 12 degrees with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circum-binary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be $5-16$ au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circum-binary disk. The origin of its opposite gradient is also discussed., Comment: 27 pages, 21 figures, 2 Tables, Accepted for publication in ApJ

We report results of a spectropolarimetric and photometric monitoring of the weak-line T Tauri star V410 Tau based on data collected mostly with SPIRou, the near-infrared (NIR) spectropolarimeter recently installed at the Canada-France-Hawaii Telescope, as part of the SPIRou Legacy Survey large programme, and with TESS between October and December 2019. Using Zeeman-Doppler Imaging (ZDI), we obtained the first maps of photospheric brightness and large-scale magnetic field at the surface of this young star derived from NIR spectropolarimetric data. For the first time, ZDI is also simultaneously applied to high-resolution spectropolarimetric data and very-high-precision photometry. V410 Tau hosts both dark and bright surface features and magnetic regions similar to those previously imaged with ZDI from optical data, except for the absence of a prominent dark polar spot. The brightness distribution is significantly less contrasted than its optical equivalent, as expected from the difference in wavelength. The large-scale magnetic field (~410 G), found to be mainly poloidal, features a dipole of ~390 G, again compatible with previous studies at optical wavelengths. NIR data yield a surface differential rotation slightly weaker than that estimated in the optical at previous epochs. Finally, we measured the radial velocity of the star and filtered out the stellar activity jitter using both ZDI and Gaussian Process Regression down to a precision of ~0.15 and 0.08 $\mathrm{km\,s^{-1}}$ RMS, respectively, confirming the previously published upper limit on the mass of a potential close-in massive planet around V410 Tau., Comment: 19 pages, 27 figures, accepted for publication in MNRAS

UV photochemistry in the surface layers of protoplanetary disks dramatically alters their composition relative to previous stages of star formation. The abundance ratio CN/HCN has long been proposed to trace the UV field in various astrophysical objects, however to date the relationship between CN, HCN, and the UV field in disks remains ambiguous. As part of the ALMA Large Program MAPS (Molecules with ALMA at Planet-forming Scales), we present observations of CN N=1-0 transitions at 0.3'' resolution towards five disk systems. All disks show bright CN emission within $\sim$50-150 au, along with a diffuse emission shelf extending up to 600 au. In all sources we find that the CN/HCN column density ratio increases with disk radius from about unity to 100, likely tracing increased UV penetration that enhances selective HCN photodissociation in the outer disk. Additionally, multiple millimeter dust gaps and rings coincide with peaks and troughs, respectively, in the CN/HCN ratio, implying that some millimeter substructures are accompanied by changes to the UV penetration in more elevated disk layers. That the CN/HCN ratio is generally high (>1) points to a robust photochemistry shaping disk chemical compositions, and also means that CN is the dominant carrier of the prebiotically interesting nitrile group at most disk radii. We also find that the local column densities of CN and HCN are positively correlated despite emitting from vertically stratified disk regions, indicating that different disk layers are chemically linked. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.

During the main phase of evolution of a protoplanetary disk, accretion regulates the inner-disk properties, such as the temperature and mass distribution, and in turn, the physical conditions associated with planet formation. The driving mechanism behind accretion remains uncertain; however, one promising mechanism is the removal of a fraction of angular momentum via a magnetohydrodynamic (MHD) disk wind launched from the inner tens of astronomical units of the disk. This paper utilizes CO isotopologue emission to study the unique molecular outflow originating from the HD 163296 protoplanetary disk obtained with the Atacama Large Millimeter/submillimeter Array. HD~163296 is one of the most well-studied Class II disks and is proposed to host multiple gas-giant planets. We robustly detect the large-scale rotating outflow in the 12CO J=2-1 and the 13CO J=2-1 and J=1-0 transitions. We constrain the kinematics, the excitation temperature of the molecular gas, and the mass-loss rate. The high ratio of the rates of ejection to accretion (5 - 50), together with the rotation signatures of the flow, provides solid evidence for an MHD disk wind. We find that the angular momentum removal by the wind is sufficient to drive accretion through the inner region of the disk; therefore, accretion driven by turbulent viscosity is not required to explain HD~163296's accretion. The low temperature of the molecular wind and its overall kinematics suggest that the MHD disk wind could be perturbed and shocked by the previously observed high-velocity atomic jet. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: Accepted ApJ July 30th 2021 This paper is part of the MAPS special issue of the Astrophysical Journal Supplement

Deuterium fractionation is dependent on various physical and chemical parameters. Thus, the formation location and thermal history of material in the solar system is often studied by measuring its D/H ratio. This requires knowledge about the deuteration processes operating during the planet formation era. We aim to study these processes by radially resolving the DCN/HCN (at 0.3" resolution) and N$_2$D$^+$/N$_2$H$^+$ (0.3 to 0.9") column density ratios toward the five protoplanetary disks observed by the Molecules with ALMA at Planet-forming scales (MAPS) Large Program. DCN is detected in all five sources, with one newly reported detection. N$_2$D$^+$ is detected in four sources, two of which are newly reported detections. We derive column density profiles that allow us to study the spatial variation of the DCN/HCN and N$_2$D$^+$/N$_2$H$^+$ ratios at high resolution. DCN/HCN varies considerably for different parts of the disks, ranging from $10^{-3}$ to $10^{-1}$. In particular, the inner disk regions generally show significantly lower HCN deuteration compared with the outer disk. In addition, our analysis confirms that two deuterium fractionation channels are active, which can alter the D/H ratio within the pool of organic molecules. N$_2$D$^+$ is found in the cold outer regions beyond $\sim$50 au, with N$_2$D$^+$/N$_2$H$^+$ ranging between $10^{-2}$ and 1 across the disk sample. This is consistent with the theoretical expectation that N$_2$H$^+$ deuteration proceeds via the low-temperature channel only. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: Accepted for publication in the Astrophysical Journal Supplement. 55 pages, 30 figures. Replacement of earlier version with updated references

Constraining dust properties of planet-forming disks via high angular resolution observations is fundamental to understanding how solids are trapped in substructures and how dust growth may be favored or accelerated therein. We use ALMA dust continuum observations of the Molecules with ALMA at Planet-forming Scales (MAPS) disks and explore a large parameter space to constrain the radial distribution of solid mass and maximum grain size in each disk, including or excluding dust scattering. In the nonscattering model, the dust surface density and maximum grain size profiles decrease from the inner disks to the outer disks, with local maxima at the bright ring locations, as expected from dust trapping models. The inferred maximum grain sizes from the inner to outer disks decrease from ~1 cm to 1 mm. For IM Lup, HD 163296, and MWC 480 in the scattering model, two solutions are compatible with their observed inner disk emission: one solution corresponding to a maximum grain size of a few millimeters (similar to the nonscattering model), and the other corresponding to a few hundred micrometer sizes. Based on the estimated Toomre parameter, only IM Lup -- which shows a prominent spiral morphology in millimeter dust -- is found to be gravitationally unstable. The estimated maximum Stokes number in all the disks lies between 0.01 and 0.3, and the estimated turbulence parameters in the rings of AS 209 and HD 163296 are close to the threshold where dust growth is limited by turbulent fragmentation. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.

We observed HCO$^+$ $J=1-0$ and H$^{13}$CO$^+$ $J=1-0$ emission towards the five protoplanetary disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480 as part of the MAPS project. HCO$^+$ is detected and mapped at 0.3\arcsec\,resolution in all five disks, while H$^{13}$CO$^+$ is detected (SNR$>6 \sigma$) towards GM Aur and HD 163296 and tentatively detected (SNR$>3 \sigma$) towards the other disks by a matched filter analysis. Inside a radius of $R\sim 100$ au, the HCO$^+$ column density is flat or shows a central dip. At outer radii ($\gtrsim 100$ au), the HCO$^+$ column density decreases outwards, while the column density ratio of HCO$^+$/CO is mostly in the range of $\sim 10^{-5}-10^{-4}$. We derived the HCO$^+$ abundance in the warm CO-rich layer, where HCO$^+$ is expected to be the dominant molecular ion. At $R\gtrsim 100$ au, the HCO$^+$ abundance is $\sim 3 \times 10^{-11} - 3\times 10^{-10}$, which is consistent with a template disk model with X-ray ionization. At the smaller radii, the abundance decreases inwards, which indicates that the ionization degree is lower in denser gas, especially inside the CO snow line, where the CO-rich layer is in the midplane. Comparison of template disk models with the column densities of HCO$^+$, N$_2$H$^+$, and N$_2$D$^+$ indicates that the midplane ionization rate is $\gtrsim 10^{-18}$ s$^{-1}$ for the disks around IM Lup, AS 209, and HD 163296. We also find hints of an increased HCO$^+$ abundance around the location of dust continuum gaps in AS 209, HD 163296, and MWC 480. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: accepted to ApJS, 33 pages, 20 figures

Small organic molecules, such as C2H, HCN, and H2CO, are tracers of the C, N, and O budget in protoplanetary disks. We present high angular resolution (10-50 au) observations of C2H, HCN, and H2CO lines in five protoplanetary disks from the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program. We derive column density and excitation temperature profiles for HCN and C2H, and find that the HCN emission arises in a temperate (20-30 K) layer in the disk, while C2H is present in relatively warmer (20-60 K) layers. In the case of HD 163296, we find a decrease in column density for HCN and C2H inside one of the dust gaps near 83 au, where a planet has been proposed to be located. We derive H2CO column density profiles assuming temperatures between 20 and 50 K, and find slightly higher column densities in the colder disks around T Tauri stars than around Herbig Ae stars. The H2CO column densities rise near the location of the CO snowline and/or millimeter dust edge, suggesting an efficient release of H2CO ices in the outer disk. Finally, we find that the inner 50 au of these disks are rich in organic species, with abundances relative to water that are similar to cometary values. Comets could therefore deliver water and key organics to future planets in these disks, similar to what might have happened here on Earth. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: 23 pages, 10 figures, 4 tables. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement

The precursors to larger, biologically-relevant molecules are detected throughout interstellar space, but determining the presence and properties of these molecules during planet formation requires observations of protoplanetary disks at high angular resolution and sensitivity. Here we present 0.3" observations of HC$_3$N, CH$_3$CN, and $c$-C$_3$H$_2$ in five protoplanetary disks observed as part of the Molecules with ALMA at Planet-forming Scales (MAPS) Large Program. We robustly detect all molecules in four of the disks (GM Aur, AS 209, HD 163296 and MWC 480) with tentative detections of $c$-C$_3$H$_2$ and CH$_3$CN in IM Lup. We observe a range of morphologies -- central peaks, single or double rings -- with no clear correlation in morphology between molecule nor disk. Emission is generally compact and on scales comparable with the millimetre dust continuum. We perform both disk-integrated and radially-resolved rotational diagram analysis to derive column densities and rotational temperatures. The latter reveals 5-10 times more column density in the inner 50-100 au of the disks when compared with the disk-integrated analysis. We demonstrate that CH$_3$CN originates from lower relative heights in the disks when compared with HC$_3$N, in some cases directly tracing the disk midplane. Finally, we find good agreement between the ratio of small to large nitriles in the outer disks and comets. Our results indicate that the protoplanetary disks studied here are host to significant reservoirs of large organic molecules, and that this planet- and comet-building material can be chemically similar to that in our own Solar System. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement Series., Comment: 24 pages, 11 figures, 5 tables. Accepted for publication in ApJSS. Updated to cross-reference other MAPS publications

Sulfur-bearing molecules play an important role in prebiotic chemistry and planet habitability. They are also proposed probes of chemical ages, elemental C/O ratio, and grain chemistry processing. Commonly detected in diverse astrophysical objects, including the Solar System, their distribution and chemistry remain, however, largely unknown in planet-forming disks. We present CS ($2-1$) observations at $\sim0."3$ resolution performed within the ALMA-MAPS Large Program toward the five disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. CS is detected in all five disks, displaying a variety of radial intensity profiles and spatial distributions across the sample, including intriguing apparent azimuthal asymmetries. Transitions of C$_2$S and SO were also serendipitously covered but only upper limits are found. For MWC 480, we present complementary ALMA observations at $\sim0."5$, of CS, $^{13}$CS, C$^{34}$S, H$_2$CS, OCS, and SO$_2$. We find a column density ratio N(H$_{2}$CS)/N(CS)$\sim2/3$, suggesting that a substantial part of the sulfur reservoir in disks is in organic form (i.e., C$_x$H$_y$S$_z$). Using astrochemical disk modeling tuned to MWC 480, we demonstrate that $N$(CS)/$N$(SO) is a promising probe for the elemental C/O ratio. The comparison with the observations provides a super-solar C/O. We also find a depleted gas-phase S/H ratio, suggesting either that part of the sulfur reservoir is locked in solid phase or that it remains in an unidentified gas-phase reservoir. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: Accepted for publication in The Astrophysical Journal Supplement (27 pages, 13 figures, 5 tables)

Planets form and obtain their compositions in dust and gas-rich disks around young stars, and the outcome of this process is intimately linked to the disk chemical properties. The distributions of molecules across disks regulate the elemental compositions of planets, including C/N/O/S ratios and metallicity (O/H and C/H), as well as access to water and prebiotically relevant organics. Emission from molecules also encodes information on disk ionization levels, temperature structures, kinematics, and gas surface densities, which are all key ingredients of disk evolution and planet formation models. The Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program was designed to expand our understanding of the chemistry of planet formation by exploring disk chemical structures down to 10 au scales. The MAPS program focuses on five disks - around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480 - in which dust substructures are detected and planet formation appears to be ongoing. We observed these disks in 4 spectral setups, which together cover ~50 lines from over 20 different species. This paper introduces the ApJS MAPS Special Issue by presenting an overview of the program motivation, disk sample, observational details, and calibration strategy. We also highlight key results, including discoveries of links between dust, gas, and chemical sub-structures, large reservoirs of nitriles and other organics in the inner disk regions, and elevated C/O ratios across most disks. We discuss how this collection of results is reshaping our view of the chemistry of planet formation., Comment: Accepted for publication in the ApJS MAPS Special Issue. v2 has updated MAPS references and a correction to Fig. 3

Emission substructures in gas and dust are common in protoplanetary disks. Such substructures can be linked to planet formation or planets themselves. We explore the observed gas substructures in AS 209 using thermochemical modeling with RAC2D and high-spatial resolution data from the Molecules with ALMA at Planet-forming Scales(MAPS) program. The observations of C$^{18}$O J=2-1 emission exhibit a strong depression at 88 au overlapping with the positions of multiple gaps in millimeter dust continuum emission. We find that the observed CO column density is consistent with either gas surface-density perturbations or chemical processing, while C$_2$H column density traces changes in the C/O ratio rather than the H$_2$ gas surface density. However, the presence of a massive planet (> 0.2 M$_{Jup}$) would be required to account for this level of gas depression, which conflicts with constraints set by the dust emission and the pressure profile measured by gas kinematics. Based on our models, we infer that a local decrease of CO abundance is required to explain the observed structure in CO, dominating over a possible gap-carving planet present and its effect on the H$_2$ surface density. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement Series., Comment: 21 pages, 18 figures

Here we present high resolution (15-24 au) observations of CO isotopologue lines from the Molecules with ALMA on Planet-forming Scales (MAPS) ALMA Large Program. Our analysis employs $^{13}$CO and C$^{18}$O ($J$=2-1), (1-0), and C$^{17}$O (1-0) line observations of five protoplanetary disks. We retrieve CO gas density distributions, using three independent methods: (1) a thermo-chemical modeling framework based on the CO data, the broadband spectral energy distribution, and the mm-continuum emission; (2) an empirical temperature distribution based on optically thick CO lines; and (3) a direct fit to the C$^{17}$O hyperfine lines. Results from these methods generally show excellent agreement. The CO gas column density profiles of the five disks show significant variations in the absolute value and the radial shape. Assuming a gas-to-dust mass ratio of 100, all five disks have a global CO-to-H$_2$ abundance of 10-100 times lower than the ISM ratio. The CO gas distributions between 150-400 au match well with models of viscous disks, supporting the long-standing theory. CO gas gaps appear to be correlated with continuum gap locations, but some deep continuum gaps do not have corresponding CO gaps. The relative depths of CO and dust gaps are generally consistent with predictions of planet-disk interactions, but some CO gaps are 5-10 times shallower than predictions based on dust gaps. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. 36 pages, 21 figures, accepted for publication in ApJS

Gas mass remains one of the most difficult protoplanetary disk properties to constrain. With much of the protoplanetary disk too cold for the main gas constituent, H2, to emit, alternative tracers such as dust, CO, or the H2 isotopolog HD are used. However, relying on disk mass measurements from any single tracer requires assumptions about the tracer's abundance relative to \hh\ and the disk temperature structure. Using new Atacama Large Millimeter/submillimeter Array (ALMA) observations from the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program as well as archival ALMA observations, we construct a disk physical/chemical model of the protoplanetary disk GM Aur. Our model is in good agreement with the spatially resolved CO isotopolog emission from eleven rotational transitions with spatial resolution ranging from 0.15'' to 0.46'' (24-73 au at 159 pc) and the spatially unresolved HD J=1-0 detection from Herschel. Our best-fit model favors a cold protoplanetary disk with a total gas mass of approximately 0.2 solar masses, a factor of 10 reduction in CO gas inside roughly 100 au and a factor of 100 reduction outside of 100 au. Despite its large mass, the disk appears to be on the whole gravitationally stable based on the derived Toomre Q parameter. However, the region between 70 and 100 au, corresponding to one of the millimeter dust rings, is close to being unstable based on the calculated Toomre Q of <1.7. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: 18 pages, 12 figures

The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a unique opportunity to study the vertical distribution of gas, chemistry, and temperature in the protoplanetary disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. By using the asymmetry of molecular line emission relative to the disk major axis, we infer the emission height ($z$) above the midplane as a function of radius ($r$). Using this method, we measure emitting surfaces for a suite of CO isotopologues, HCN, and C$_2$H. We find that $^{12}$CO emission traces the most elevated regions with $z/r > 0.3$, while emission from the less abundant $^{13}$CO and C$^{18}$O probes deeper into the disk at altitudes of $z/r \lesssim 0.2$. C$_2$H and HCN have lower opacities and SNRs, making surface fitting more difficult, and could only be reliably constrained in AS 209, HD 163296, and MWC 480, with $z/r \lesssim 0.1$, i.e., relatively close to the planet-forming midplanes. We determine peak brightness temperatures of the optically thick CO isotopologues and use these to trace 2D disk temperature structures. Several CO temperature profiles and emission surfaces show dips in temperature or vertical height, some of which are associated with gaps and rings in line and/or continuum emission. These substructures may be due to local changes in CO column density, gas surface density, or gas temperatures, and detailed thermo-chemical models are necessary to better constrain their origins and relate the chemical compositions of elevated disk layers with those of planet-forming material in disk midplanes. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: 31 pages, 20 figures, accepted for publication in ApJS, MAPS cross-references updated

The concentric gaps and rings commonly observed in protoplanetary disks in millimeter continuum emission have lent the impression that planet formation generally proceeds within orderly, isolated systems. While deep observations of spatially resolved molecular emission have been comparatively limited, they are increasingly suggesting that some disks interact with their surroundings while planet formation is underway. We present an analysis of complex features identified around GM Aur in $^{12}$CO $J=2-1$ images at a spatial resolution of $\sim40$ au. In addition to a Keplerian disk extending to a radius of $\sim550$ au, the CO emission traces flocculent spiral arms out to radii of $\sim$1200 au, a tail extending $\sim1800$ au southwest of GM Aur, and diffuse structures extending from the north side of the disk up to radii of $\sim1900$ au. The diffuse structures coincide with a "dust ribbon" previously identified in scattered light. The large-scale asymmetric gas features present a striking contrast with the mostly axisymmetric, multi-ringed millimeter continuum tracing the pebble disk. We hypothesize that GM Aur's complex gas structures result from late infall of remnant envelope or cloud material onto the disk. The morphological similarities to the SU Aur and AB Aur systems, which are also located in the L1517 cloud, provide additional support to a scenario in which interactions with the environment are playing a role in regulating the distribution and transport of material in all three of these Class II disk systems. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: 34 pages, 21 figures, in press at ApJS, cross-references updated

Constraining the distribution of gas and dust in the inner 20 au of protoplanetary disks is difficult. At the same time, this region is thought to be responsible for most planet formation, especially around the water ice line at 3-10 au. Under the assumption that the gas is in a Keplerian disk, we use the exquisite sensitivity of the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA large program to construct radial surface brightness profiles with a ~3 au effective resolution for the CO isotopologue J=2-1 lines using the line velocity profile. IM Lup reveals a central depression in 13CO and C18O that is ascribed to a pileup of ~500 $M_\oplus$ of dust in the inner 20 au, leading to a gas-to-dust ratio of around <10. This pileup is consistent with efficient drift of grains ($\gtrsim$ 100 $M_\oplus$ Myr$^{-1}$) and a local gas-to-dust ratio that suggests that the streaming instability could be active. The CO isotopologue emission in the GM Aur disk is consistent with a small (~15 au), strongly depleted gas cavity within the ~40 au dust cavity. The radial surface brightness profiles for both the AS 209 and HD 163296 disks show a local minimum and maximum in the C18O emission at the location of a known dust ring (~14 au) and gap (~10 au), respectively. This indicates that the dust ring has a low gas-to-dust ratio ($>$ 10) and that the dust gap is gas-rich enough to have optically thick C18O., Comment: 22 pages, 14 figures, accepted by ApJS. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. Update with correct references to other MAPS papers

The elemental composition of the gas and dust in a protoplanetary disk influences the compositions of the planets that form in it. We use the Molecules with ALMA at Planet-forming Scales (MAPS) data to constrain the elemental composition of the gas at the locations of potentially forming planets. The elemental abundances are inferred by comparing source-specific gas-grain thermochemical models, with variable C/O ratios and small-grain abundances, from the DALI code with CO and C2H column densities derived from the high-resolution observations of the disks of AS 209, HD 163296, and MWC 480. Elevated C/O ratios (~2.0), even within the CO ice line, are necessary to match the inferred C2H column densities, over most of the pebble disk. Combined with constraints on the CO abundances in these systems, this implies that both the O/H and C/H ratios in the gas are substellar by a factor of 4-10, with the O/H depleted by a factor of 20-50, resulting in the high C/O ratios. This necessitates that even within the CO ice line, most of the volatile carbon and oxygen is still trapped on grains in the midplane. Planets accreting gas in the gaps of the AS 209, HD 163296, and MWC 480 disks will thus acquire very little carbon and oxygen after reaching the pebble isolation mass. In the absence of atmosphere-enriching events, these planets would thus have a strongly substellar O/H and C/H and superstellar C/O atmospheric composition., Comment: 19 pages, 8 figures This paper is part of the MAPS special issue of the Astrophysical Journal Supplement. Updates references for other MAPS papers

We explore the dynamical structure of the protoplanetary disks surrounding HD 163296 and MWC 480 as part of the Molecules with ALMA at Planet Forming Scales (MAPS) large program. Using the $J = 2-1$ transitions of $^{12}$CO, $^{13}$CO and C$^{18}$O imaged at spatial resolutions of $\sim 0.^{\prime \prime}15$ and with a channel spacing of $200$ ${\rm m\,s^{-1}}$, we find perturbations from Keplerian rotation in the projected velocity fields of both disks ($\lesssim\!5\%$ of the local Keplerian velocity), suggestive of large-scale (10s of au in size), coherent flows. By accounting for the azimuthal dependence on the projection of the velocity field, the velocity fields were decomposed into azimuthally averaged orthogonal components, $v_{\phi}$, $v_r$ and $v_z$. Using the optically thick $^{12}$CO emission as a probe of the gas temperature, local variations of $\approx\! 3$ K ($\approx\! 5 \%$ relative changes) were observed and found to be associated with the kinematic substructures. The MWC 480 disk hosts a suite of tightly wound spiral arms. The spirals arms, in conjunction with the highly localized perturbations in the gas velocity structure (kinematic planetary signatures), indicate a giant planet, $\sim\! 1$ $M_{\rm Jup}$, at a radius of $\approx 245$ au. In the disk of HD 163296, the kinematic substructures were consistent with previous studies of Pinte et al. (2018a) and Teague et al. (2018a) advocating for multiple $\sim\! 1$ $M_{\rm Jup}$ planets embedded in the disk. These results demonstrate that molecular line observations that characterize the dynamical structure of disks can be used to search for the signatures of embedded planets. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: 32 pages, 21 figures, accepted for publication in ApJS, MAPS cross-references updated

The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a detailed, high resolution (${\sim}$10-20 au) view of molecular line emission in five protoplanetary disks at spatial scales relevant for planet formation. Here, we present a systematic analysis of chemical substructures in 18 molecular lines toward the MAPS sources: IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. We identify more than 200 chemical substructures, which are found at nearly all radii where line emission is detected. A wide diversity of radial morphologies - including rings, gaps, and plateaus - is observed both within each disk and across the MAPS sample. This diversity in line emission profiles is also present in the innermost 50 au. Overall, this suggests that planets form in varied chemical environments both across disks and at different radii within the same disk. Interior to 150 au, the majority of chemical substructures across the MAPS disks are spatially coincident with substructures in the millimeter continuum, indicative of physical and chemical links between the disk midplane and warm, elevated molecular emission layers. Some chemical substructures in the inner disk and most chemical substructures exterior to 150 au cannot be directly linked to dust substructure, however, which indicates that there are also other causes of chemical substructures, such as snowlines, gradients in UV photon fluxes, ionization, and radially-varying elemental ratios. This implies that chemical substructures could be developed into powerful probes of different disk characteristics, in addition to influencing the environments within which planets assemble. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: 62 pages, 31 figures, accepted for publication in ApJS, MAPS cross-references updated, corrected Figure 21, updated gas disk sizes (Table 2, Figures 15-16) from associated Erratum

Understanding the temperature structure of protoplanetary disks is key to interpreting observations, predicting the physical and chemical evolution of the disk, and modeling planet formation processes. In this study, we constrain the two-dimensional thermal structure of the disk around Herbig Ae star HD 163296. Using the thermo-chemical code RAC2D, we derive a thermal structure that reproduces spatially resolved ALMA observations (~0.12 arcsec (13 au) - 0.25 arcsec (26 au)) of CO J = 2-1, 13CO J = 1-0, 2-1, C18O J = 1-0, 2-1, and C17O J = 1-0, the HD J = 1-0 flux upper limit, the spectral energy distribution (SED), and continuum morphology. The final model incorporates both a radial depletion of CO motivated by a time scale shorter than typical CO gas-phase chemistry (0.01 Myr) and an enhanced temperature near the surface layer of the the inner disk (z/r <= 0.21). This model agrees with the majority of the empirically derived temperatures and observed emitting surfaces derived from the J = 2-1 CO observations. We find an upper limit for the disk mass of 0.35 Msun, using the upper limit of the HD J = 1-0 and J = 2-1 flux. With our final thermal structure, we explore the impact that gaps have on the temperature structure constrained by observations of the resolved gaps. Adding a large gap in the gas and small dust additionally increases gas temperature in the gap by only 5-10%. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: 15 pages + 11 pages of appendix, accepted to ApJS, part of MAPS collaboration

The Molecules with ALMA at Planet-forming Scales large program (MAPS LP) surveyed the chemical structures of five protoplanetary disks across more than 40 different spectral lines at high angular resolution (0.15" and 0.30" beams for Bands 6 and 3, respectively) and sensitivity (spanning 0.3 - 1.3 mJy/beam and 0.4 - 1.9 mJy/beam for Bands 6 and 3, respectively). In this article, we describe our multi-stage workflow -- built around the CASA tclean image deconvolution procedure -- that we used to generate the core data product of the MAPS LP: the position-position-velocity image cubes for each spectral line. Owing to the expansive nature of the survey, we encountered a range of imaging challenges; some are familiar to the sub-mm protoplanetary disk community, like the benefits of using an accurate CLEAN mask, and others less well-known, like the incorrect default flux scaling of the CLEAN residual map first described in Jorsater & van Moorsel 1995 (the "JvM effect"). We distill lessons learned into recommended workflows for synthesizing image cubes of molecular emission. In particular, we describe how to produce image cubes with accurate fluxes via the "JvM correction," a procedure that is generally applicable to any image synthesized via CLEAN deconvolution but is especially critical for low S/N emission. We further explain how we used visibility tapering to promote a common, fiducial beam size and contextualize the interpretation of signal to noise ratio when detecting molecular emission from protoplanetary disks. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement., Comment: ApJS accepted. Part of the MAPS ALMA large program series: http://www.alma-maps.info/. Updated bibliography with MAPS LP arXiv references

PDS70 is a unique system in which two protoplanets, PDS70b and c, have been discovered within the dust-depleted cavity of their disk, at $\sim$22 and 34au respectively, by direct imaging at infrared wavelengths. Subsequent detection of the planets in the H$\alpha$ line indicates that they are still accreting material through circumplanetary disks. In this Letter, we present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the dust continuum emission at 855$\mu$m at high angular resolution ($\sim$20mas, 2.3au) that aim to resolve the circumplanetary disks and constrain their dust masses. Our observations confirm the presence of a compact source of emission co-located with PDS70c, spatially separated from the circumstellar disk and less extended than $\sim$1.2au in radius, a value close to the expected truncation radius of the cicumplanetary disk at a third of the Hill radius. The emission around PDS70c has a peak intensity of $\sim$86$\pm$16 $\mu \mathrm{Jy}~\mathrm{beam}^{-1}$ which corresponds to a dust mass of $\sim$0.031M$_{\oplus}$ or $\sim$0.007M$_{\oplus}$, assuming that it is only constituted of 1 $\mu$m or 1 mm sized grains, respectively. We also detect extended, low surface brightness continuum emission within the cavity near PDS70b. We observe an optically thin inner disk within 18au of the star with an emission that could result from small micron-sized grains transported from the outer disk through the orbits of b and c. In addition, we find that the outer disk resolves into a narrow and bright ring with a faint inner shoulder., Comment: ApJ Letters, in press; 19 pages, 9 figures

IRAS~04158+2805 has long been thought to be a very low mass T-Tauri star (VLMS) surrounded by a nearly edge-on, extremely large disc. Recent observations revealed that this source hosts a binary surrounded by an extended circumbinary disc with a central dust cavity. In this paper, we combine ALMA multi-wavelength observations of continuum and $^{12}$CO line emission, with H$\alpha$ imaging and Keck astrometric measures of the binary to develop a coherent dynamical model of this system. The system features an azimuthal asymmetry detected at the western edge of the cavity in Band~7 observations and a wiggling outflow. Dust emission in ALMA Band 4 from the proximity of the individual stars suggests the presence of marginally resolved circumstellar discs. We estimate the binary orbital parameters from the measured arc of the orbit from Keck and ALMA astrometry. We further constrain these estimates using considerations from binary-disc interaction theory. We finally perform three SPH gas + dust simulations based on the theoretical constraints; we post-process the hydrodynamic output using radiative transfer Monte Carlo methods and directly compare the models with observations. Our results suggest that a highly eccentric $e\sim 0.5\textrm{--}0.7$ equal mass binary, with a semi-major axis of $\sim 55$ au, and small/moderate orbital plane vs. circumbinary disc inclination $\theta\lesssim 30^\circ$ provides a good match with observations. A dust mass of $\sim 1.5\times 10^{-4} {\rm M_\odot}$ best reproduces the flux in Band 7 continuum observations. Synthetic CO line emission maps qualitatively capture both the emission from the central region and the non-Keplerian nature of the gas motion in the binary proximity., Comment: 19 pages, 14 figures, 1 table, accepted for publication in MNRAS

We obtained polarimetric differential imaging of a gas-rich debris disk around HD 141569A with SPHERE in the H-band to compare the scattering properties of the innermost ring at 44 au with former observations in total intensity with the same instrument. In polarimetric imaging, we observed that the intensity of the ring peaks in the south-east, mostly in the forward direction, whereas in total intensity imaging, the ring is detected only at the south. This noticeable characteristic suggests a non-uniform dust density in the ring. We implemented a density function varying azimuthally along the ring and generated synthetic images both in polarimetry and in total intensity, which are then compared to the actual data. We find that the dust density peaks in the south-west at an azimuthal angle of $220^{\circ} \sim 238^{\circ}$ with a rather broad width of $61^{\circ} \sim 127^{\circ}$. Although there are still uncertainties that remain in the determination of the anisotropic scattering factor, the implementation of an azimuthal density variation to fit the data proved to be robust. Upon elaborating on the origin of this dust density distribution, we conclude that it could be the result of a massive collision when we account for the effect of the high gas mass that is present in the system on the dynamics of grains. Using the outcome of this modelization, we further measured the polarized scattering phase function for the observed scattering angle between 33$^{\circ}$ and 147$^{\circ}$ as well as the spectral reflectance of the southern part of the ring between 0.98 $\mu$m and 2.1 $\mu$m. We tentatively derived the grain properties by comparing these quantities with MCFOST models and assuming Mie scattering. Our preliminary interpretation indicates a mixture of porous sub-micron sized astro-silicate and carbonaceous grains., Comment: 21 pages, 18 figures, Accepted for publication in Astronomy & Astrophysics journal

We present high-resolution $^{12}$CO and $^{13}$CO 2-1 ALMA observations, as well as optical and near-infrared spectroscopy, of the highly-inclined protoplanetary disk around SSTC2D J163131.2-242627. The spectral type we derive for the source is consistent with a $\rm 1.2 \, M_{\odot}$ star inferred from the ALMA observations. Despite its massive circumstellar disk, we find little to no evidence for ongoing accretion on the star. The CO maps reveal a disk that is unusually compact along the vertical direction, consistent with its appearance in scattered light images. The gas disk extends about twice as far away as both the submillimeter continuum and the optical scattered light. CO is detected from two surface layers separated by a midplane region in which CO emission is suppressed, as expected from freeze-out in the cold midplane. We apply a modified version of the Topographically Reconstructed Distribution method presented by Dutrey et al. 2017 to derive the temperature structure of the disk. We find a temperature in the CO-emitting layers and the midplane of $\sim$33 K and $\sim$20 K at $\rm R<200$ au, respectively. Outside of $\rm R>200$ au, the disk's midplane temperature increases to $\sim$30 K, with a nearly vertically isothermal profile. The transition in CO temperature coincides with a dramatic reduction in the sub-micron and sub-millimeter emission from the disk. We interpret this as interstellar UV radiation providing an additional source of heating to the outer part of the disk., Comment: 27 pages, 18 figures, 1 table