Search

We have used the multiobject mode of the Near-Infrared Spectrograph (NIRSpec) on board the James Webb Space Telescope (JWST) to obtain low-resolution 1–5 μ m spectra of 22 brown dwarf candidates in the Orion Nebula Cluster, which were selected with archival images from the Hubble Space Telescope. One of the targets was previously classified as a Herbig–Haro (HH) object and exhibits strong emission in H i, H2, and the fundamental band of CO, further demonstrating that HH objects can have bright emission in that CO band. The remaining targets have late spectral types (M6.5 to early L) and are young based on gravity-sensitive features, as expected for low-mass members of the cluster. According to theoretical evolutionary models, these objects should have masses that range from the hydrogen burning limit to 0.003–0.007 M ⊙. Two of the NIRSpec targets were identified as proplyds in earlier analysis of Hubble images. They have spectral types of M6.5 and M7.5, making them two of the coolest and least massive known proplyds. Another brown dwarf shows absorption bands at 3–5 μ m from ices containing H2O, CO2, OCN−, and CO, indicating that it is either an edge-on class II system or a class I protostar. It is the coolest and least massive object that has detections of these ice features. In addition, it appears to be the first candidate for a protostellar brown dwarf that has spectroscopy confirming its late spectral type. [ABSTRACT FROM AUTHOR]

Context. Stars and planets form in regions of enhanced stellar density, subjecting protoplanetary discs to gravitational perturbations from neighbouring stars. Observations in the Taurus star-forming region have uncovered evidence of at least three recent, star-disc encounters that have truncated discs (HV/DO Tau, RW Aurigae, and UX Tau), raising questions about the frequency of such events. Aims. We aim to assess the probability of observing truncating star-disc encounters in Taurus. Methods. We generated a physically motivated dynamical model including binaries and a spatial-kinematic substructure to follow the historical dynamical evolution of the Taurus star-forming region. We used this model to track star-disc encounters and the resulting outer disc truncation over the lifetime of Taurus. Results. A quarter of discs are truncated below 30 au by dynamical encounters, but this truncation mostly occurs in binaries over the course of a few orbital periods, on a timescale ≲0.1 Myr. Nonetheless, some truncating encounters still occur up to the present age of Taurus. Strongly truncating encounters (ejecting ≳10 percent of the disc mass) occur at a rate ∼10 Myr−1, sufficient to explain the encounter between HV and DO Tau ∼0.1 Myr ago. If encounters that eject only ∼1 percent of the disc mass are responsible for RW Aurigae and UX Tau, then they are also expected with encounter rate Γenc ∼ 100–200 Myr−1. However, the observed sample of recent encounters is probably incomplete, since these examples occurred in systems that are not consistent with a random drawing from the mass function. One more observed example would statistically imply additional physics, such as replenishment of the outer disc material. Conclusions. The marginal consistency of the frequency of observed recent star-disc encounters with theoretical expectations underlines the value of future large surveys searching for external structures associated with recent encounters. The outcome of such a survey offers a highly constraining, novel probe of protoplanetary disc physics. [ABSTRACT FROM AUTHOR]

Kinematic studies of protoplanetary discs are a valuable method for uncovering hidden companions. In the first paper of this series, we presented five morphological and kinematic criteria that aid in asserting the binary nature of a protoplanetary disc. In this work, we study the kinematic signatures of circumbinary discs in the residuals of their velocity maps. We show that Doppler-flips, spiral arms, eccentric gas motion, fast flows inside of the cavity, and vortex-like kinematic signatures are commonly observed. Unlike in the planetary mass companion case, Doppler-flips in circumbinary discs are not necessarily centred on a companion, and can extend towards the cavity edge. We then compare the kinematic signatures in our simulations with observations and see similarities to the Doppler-flip signal in HD 100546 and the vortex-like kinematic signatures in HD 142527. Our analysis also reveals kinematic evidence for binarity in several protoplantary discs typically regarded as circumstellar rather than circumbinary, including AB Aurigae and HD 100546. [ABSTRACT FROM AUTHOR]

A consensus prevails with regard to star-disk systems accreting most of their mass and angular momentum during the collapse of a prestellar core. However, recent results have indicated that stars experience post-collapse or late infall, during which the star and its disk are refreshed with material from the protostellar environment through accretion streamers. Apart from adding mass to the star-disk system, infall potentially supplies a substantial amount of angular momentum, as the infalling material is initially not bound to the collapsing prestellar core. We investigate the orientation of infall on star-disk systems by analyzing the properties of accreting tracer particles in three-dimensional magnetohydrodynamical (3D MHD) simulations of a molecular cloud that is (4 pc)3 in volume. In contrast to the traditional picture, where the rotational axis is inherited from the collapse of a coherent pre-stellar core, the orientation of star-disk systems changes substantially throughout the accretion process, thereby extending the possibility of primordial misalignment as the source of large obliquities. In agreement with previous results that show larger contributions of late infall for increasing stellar masses, a misaligned infall is more likely to lead to a prolonged change in orientation for stars of higher final mass. On average, brown dwarfs and very low mass stars are more likely to form and accrete all of their mass as part of a multiple system, while stars with final masses above a few 0.1 M⊙ are more likely to accrete part of their mass as single stars. Finally, we find an overall trend among our sample: the post-collapse accretion phase is more anisotropic than the early collapse phase. This result is consistent with a scenario of Bondi-Hoyle-Littletlon accretion during the post-collapse phase, while the initial collapse is less anisotropic – despite the fact that material is funneled through accretion channels. [ABSTRACT FROM AUTHOR]

PDS 70 is a protoplanetary system that hosts two actively accreting gas giants, namely, PDS 70b and PDS 70 c. The system has a ∼60–100 au dusty ring that has been resolved by the Atacama Large Millimeter/Submillimeter Array (ALMA), along with circumplanetary disks around the two gas giants. Here, we report the first Karl G. Jansky Very Large Array (JVLA) Q - (40–48 GHz), Ka - (29–37 GHz), K - (18–26 GHz), and X - (8–12 GHz) bands' continuum observations, and the complementary ALMA Bands 3 (∼98 GHz) and 4 (∼145 GHz) observations towards PDS 70. The dusty ring appears azimuthally asymmetric in our ALMA images. We obtained firm detections at Ka and K bands without spatially resolving the source; we obtained a marginal detection at Q band, and no detection at X band. The spectral indices (α) are 5 ± 1 at 33–44 GHz and 0.6 ± 0.2 at 22–33 GHz. At 10–22 GHz, the conservative lower limit of α is 1.7. The 33–44 GHz flux density is likely dominated by the optically thin thermal emission of grown dust with ≳1 mm maximum grain sizes, which may be associated with the azimuthally asymmetric substructure induced by planet–disk interaction. Since PDS 70 was not detected at X band, we found it hard to explain the low spectral index at 22–33 GHz only with free–free emission. Hence, we attribute the dominant emission at 22–33 GHz to the emission of spinning nano-meter-sized dust particles, while free–free emission may partly contribute to emission at this frequency range. In some protoplanetary disks, the emission of spinning nano-meter-sized dust particles may resemble the 20–50 GHz excess in the spectra of millimeter-sized dust. The finding of strong continuum emission of spinning nano-meter-sized particles can complicate the procedure of constraining the properties of grown dust. Future high resolution, multifrequency JVLA/Next Generation Very Large Array and Square Kilometer Array observations may shed light on this issue. [ABSTRACT FROM AUTHOR]

Crescent-shaped asymmetries are common in millimetre observations of protoplanetary discs and are usually attributed to vortices or dust overdensities. However, they often appear on a single side of the major axis and roughly symmetric about the minor axis, suggesting a geometric origin. In this work, we interpret such asymmetries as emission from the exposed inner cavity walls of inclined discs and use them to characterize their vertical extent. Here we focus on the discs around CIDA 9 and RY Tau, first modelling their observations in visibility space with a simple geometric prescription for the walls, and then exploring more detailed radiative transfer models. Accounting for the wall emission yields significantly better residuals than purely axisymmetric models, and we estimate the dust scale height of these systems to be 0.4 au at 37 au for CIDA 9 and 0.2 au at 12 au for RY Tau. Finally, we identify crescent-shaped asymmetries in twelve discs, nine of which have constraints on their orientation – in all cases, the asymmetry appears on the far-side of the disc, lending support to the hypothesis that they are due to their inner rims. Modelling this effect in larger samples of discs will help to build a statistical view of their vertical structure. [ABSTRACT FROM AUTHOR]

The nearby oxygen-rich AGB star L |$_2$|  Pup hosts a well-studied nearly edge-on disc. To date, discs around AGB stars have not been chemically studied in detail. By combining a parametrization commonly used for protoplanetary discs and archival ALMA observations, we retrieved an updated density and temperature structure of this disc. This physical model was then used as input to the first chemical model of an AGB disc. The model shows that the physical structure of the disc has a large impact on its chemistry, with certain species showing large changes in column density relative to a radial outflow, indicating that chemistry could be used as a tracer of discs that cannot be directly imaged. Despite its oxygen-rich nature, the daughter species formed within the disc are surprisingly carbon-rich. Two chemical regimes can be distinguished: cosmic-ray induced chemistry in the midplane and photochemistry induced by the interstellar radiation field in the outer regions. Certain complex organic molecules are formed in the midplane. This occurs via gas-phase chemistry only, as the disc is too warm for dust-gas chemistry. The photochemistry in the outer regions leads to the efficient formation of (long) carbon-chains. The predictions of the model allow us to tentatively put the disc's age |$\lesssim 10^5$| yr. Additional observations are necessary to better constrain the physical structure of L |$_2$|  Pup's disc and are essential to test the predictions made by the chemical model. Our exploratory work paves the way for a more general study of the chemistry of AGB discs. [ABSTRACT FROM AUTHOR]

In circumstellar disks around young stars, the gravitational influence of nascent planets produces telltale patterns in density, temperature, and kinematics. To better understand these signatures, we first performed 3D hydrodynamical simulations of a 0.012 M⊙ disk with a Saturn-mass planet orbiting circularly in-plane at 40 au. We tested four different disk thermodynamic prescriptions (in increasing order of complexity: local isothermality, β cooling, two-temperature radiation hydrodynamics, and three-temperature radiation hydrodynamics), finding that β cooling offers a reasonable approximation for the three-temperature approach when the planet is not massive or luminous enough to substantially alter the background temperature and density structure. Thereafter, using the three-temperature scheme, we relaxed this assumption, simulating a range of different planet masses (Neptune-mass, Saturn-mass, and Jupiter-mass) and accretion luminosities (0 and 10−3L⊙) in the same disk. Our investigation revealed that signatures of disk–planet interaction strengthen with increasing planet mass, with circumplanetary flows becoming prominent in the high-planet-mass regime. Accretion luminosity, which adds pressure support around the planet, was found to weaken the midplane Doppler flip, which is potentially visible in optically thin tracers such as C18O, while strengthening the spiral signature, particularly in upper disk layers sensitive to thicker lines, such as those of 12CO. [ABSTRACT FROM AUTHOR]

Context. Most ground-based direct-imaging planet search campaigns use angular differential imaging that distorts the signal from extended sources, such as protoplanetary disks. In the case of the young system PDS 70, for which two planets were detected within the cavity of a protoplanetary disk, obtaining a reliable image of both planets and the 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 try to match the near-infrared imaging and polarimetric data within measurement errors. Sub-millimeter imaging data from ALMA were also modeled. Furthermore, we extracted the spectra of the planets and the disk and compared them Results. With strong constraints, 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 ∼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 sub-micron, micron, and grains estimated by the modeling indicate a shallow grain-size index greater than −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, and it is clearly redder than the disk spectra. Planet c is possibly associated with the sudden flaring of the disk starting at ∼50 au. We found no new planets in the system. If we assume DUSTY models and an age of 5 Myr, this indicates no new planets more massive than 5 outside a 12 au separation. [ABSTRACT FROM AUTHOR]

Context. To understand how planets form in protoplanetary disks, it is necessary to characterize their gas and dust distribution and masses. This requires a combination of high-resolution dust continuum and molecular line interferometric observations, coupled with advanced theoretical models of protoplanetary disk physics, chemical composition, and radiative transfer. Aims. We aim to constrain the gas density and temperature distributions as well as gas masses in several T Tauri protoplanetary disks located in Taurus. We use the 12CO, 13CO, and C18O (2–1) isotopologue emission observed at 0.9″ with the IRAM NOrthern Extended Millimeter Array (NOEMA) as part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution PIs: P. Caselli & Th. Henning). Our sample consists of Class II disks with no evidence of strong radial substructures. We use these data to constrain the thermal and chemical structure of these disks through theoretical models for gas emission. Methods. To fit the combined optically thick and thin CO line data in Fourier space, we developed the DiskCheF code, which includes the parameterized disk physical structure, machine-learning (ML) accelerated chemistry, and the RADMC-3D line radiative transfer module. A key novelty of DiskCheF is the fast and feasible ML-based chemistry trained on the extended grid of the disk physical-chemical models precomputed with the ANDES2 code. This ML approach allows complex chemical kinetics models to be included in a time-consuming disk fitting without the need to run a chemical code. Results. We present a novel approach to incorporate chemistry into disk modeling without the need to explicitly calculate a chemical network every time. Using this new disk modeling tool, we successfully fit the 12CO, 13CO, and C18O (2-1) data from the CI, CY, DL, DM, DN, and IQ Tau disks. The combination of optically thin and optically thick CO lines allows us to simultaneously constrain the disk temperature and mass distribution, and derive the CO-based gas masses. The best-fit disk gas masses range between 0.005 and 0.04 M⊙. These values are in reasonable agreement with the disk dust masses rescaled by a factor of 100 as well as with other indirect gas measurements via, for example, modeling of the wavelength dependence of the dust continuum emission radii, and HD and CO isotopologue emission. [ABSTRACT FROM AUTHOR]

The origin and evolution of CO gas in debris discs has been debated since its initial detection. The gas could have a primordial origin, as a remnant of the protoplanetary disc or a secondary exocometary origin. This paper investigates the origin of gas in two debris discs, HD110058 and HD131488, using Hubble Space Telescope (HST) observations of CI and CO, which play critical roles in the gas evolution. We fitted several electronic transitions of CI and CO rovibronic bands to derive column densities and temperatures for each system, revealing high CO column densities (∼3–4 orders of magnitude higher than β Pictoris), and low CI/CO ratios in both. Using the exogas model, we simulated the radial evolution of the gas in the debris disc assuming a secondary gas origin. We explored a wide range of CO exocometary release rates and α viscosities, which are the key parameters of the model. Additionally, we incorporated photodissociation due to stellar ultraviolet (UV) to the exogas model and found that it is negligible for typical CO-rich discs and host stars, even at a few au due to the high radial optical depths in the Extreme ultraviolet radiation (EUV). We find that the current steady-state secondary release model cannot simultaneously reproduce the CO and CI HST -derived column densities, as it predicts larger CI/CO ratios than observed. Our direct UV measurement of low CI/CO ratios agrees with results derived from recent Atacama Large Millimetre/submillimetre Array findings and may point to vertical layering of CI, additional CI removal, CO shielding processes, or different gas origin scenarios. [ABSTRACT FROM AUTHOR]

Understanding the evolution and dispersal via energetic stellar radiation of protoplanetary discs (PPDs) is a prominent challenge in astrophysics. It has been established that X-ray luminosity from the central protostar can significantly heat the surface of the disc, causing powerful photoevaporative winds that eject a considerable fraction of the disc's mass. Recent work in the field has moreover shown the importance of global PPD simulations that simultaneously take into account non-ideal magnetohydrodynamic (MHD) effects and detailed thermochemistry. In this paper, we combine these two aspects and figure out how they interact. Focus is put on the Hall Effect (HE) and the impact it has on the overall field topology and mass-loss/accretion rates. Utilizing a novel X-ray temperature parametrization, we perform 2D-axisymmetric MHD simulations with the nirvana-iii fluid code, covering all non-ideal effects. We find that, in the aligned orientation, the HE causes prominent inward displacement of the poloidal field lines that increase the accretion rate through a laminar Maxwell stress. We find that outflows are mainly driven by photoevaporation – unless the magnetic field strength is considerable (i.e. β p ≤ 103) or the X-ray luminosity low enough (i.e. log  LX ≤ 29.3). Inferred mass-loss rate are in the range of the expected values 10−8– |$10^{-7}\, M_\odot \, {\rm yr}^{-1}$|⁠. Moreover, we performed pure hydrodynamic (HD) runs and compared them with the equivalent MHD runs. We concluded that the magnetic field does indeed contribute to the mass-loss rate, albeit only discernibly so for low enough LX (i.e. log  LX ≤ 30.8). [ABSTRACT FROM AUTHOR]

Context. Two protoplanets have recently been discovered within the PDS 70 protoplanetary disk. JWST/NIRCam offers a unique opportunity to characterize them and their birth environment at wavelengths that are difficult to access from the ground. Aims. We image the circumstellar environment of PDS 70 at 1.87 μm and 4.83 μm, assess the presence of Pa-α emission due to accretion onto the protoplanets, and probe any IR excess indicative of heated circumplanetary material. Methods. We obtained noncoronagraphic JWST/NIRCam images of PDS 70 within the MIRI mid-INfrared Disk Survey (MINDS) program. We leveraged the Vortex Image Processing (VIP) package for data reduction, and we developed dedicated routines for optimal stellar point spread function subtraction, unbiased imaging of the disk, and protoplanet flux measurement in this type of dataset. A radiative transfer model of the disk was used to separate the contributions from the disk and the protoplanets. Results. We redetect both protoplanets and identify extended emission after subtracting a disk model, including a large-scale spiral-like feature. We interpret its signal in the direct vicinity of planet c as tracing the accretion stream that feeds its circumplanetary disk, while the outer part of the feature may rather reflect asymmetric illumination of the outer disk. We also report a bright signal that is consistent with a previously proposed protoplanet candidate enshrouded in dust near the 1:2:4 mean-motion resonance with planets b and c. The 1.87 μm flux of planet b is consistent with atmospheric model predictions, but the flux of planet c is not. We discuss potential origins for this discrepancy, including significant Pa-α line emission. The 4.83 μm fluxes of planets b and c suggest enshrouding dust or heated CO emission from their circumplanetary environment. Conclusions. The use of image-processing methods that are optimized for extended disk signals on high-sensitivity and high-stability from JWST can uniquely identify signatures of planet–disk interactions and enable accurate photometry of protoplanets at wavelengths that are difficult to probe from the ground. Our results indicate that more protoplanets can be identified and characterized in other JWST datasets. [ABSTRACT FROM AUTHOR]

Context. Resolved observations at near-infrared (near-IR) and millimeter wavelengths have revealed a diverse population of planet-forming disks. In particular, near-IR scattered light observations usually target close-by, low-mass star-forming regions. However, disk evolution in high-mass star-forming regions is likely affected by the different environment. Orion is the closest high-mass star-forming region, enabling resolved observations to be undertaken in the near-IR. Aims. We seek to examine planet-forming disks, in scattered light, within the high-mass star-forming region of Orion in order to study the impact of the environment in a higher-mass star-forming region on disk evolution. Methods. 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 used polarization differential imaging in order to detect scattered light from circumstellar dust. From the scattered light observations we characterized the disk orientation, radius, and contrast. We analysed the disks in the context of the stellar parameters and the environment of the Orion star-forming region. We used ancillary X-shooter spectroscopic observations to characterize the central stars in the systems. We furthermore used a combination of new and archival ALMA mm-continuum photometry to characterize the dust masses present in the circumstellar disks. Results. Within our sample, we detect extended circumstellar disks in ten of 23 systems. Of these, three are exceptionally extended (V351 Ori, V599 Ori, and V1012 Ori) and show scattered light asymmetries that 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 nine such objects in our sample, of which six were previously unknown. We find in particular a possible substellar companion (either a very low-mass star or a high-mass brown dwarf) 137 au from the star RY Ori. We find a strong anticorrelation between disk detection and multiplicity, with only two of our ten disk detections located in stellar multiple systems. We also find a correlation between scattered light contrast and the millimeter flux. This trend is not captured by previous studies of a more diversified sample and is due to the absence of extended, self-shadowed disks in our Orion sample. Conversely, we do not find significant correlations between the scattered light contrast of the disks and the stellar mass or age. We investigate the radial extent of the disks and compare this to the estimated far-ultraviolet (FUV) field strength at the system location. While we do not find a direct correlation, we notice that no extended disks are detected above an FUV field strength of ~300 G0. [ABSTRACT FROM AUTHOR]

Context. The past few years have seen a revolution in the study of circumstellar disks. New instrumentation in the near-infrared and (sub)millimeter regimes have allowed us to routinely spatially resolve disks around young stars of nearby star-forming regions. As a result, we have found that substructures with scales of ~10 au in disks are common. We have also revealed a zoo of different morphologies, sizes, and luminosities that is as complex as the diversity of architectures found in evolved exoplanet systems. Aims. We study disk evolutionary trends as they appear in scattered light observations. Scattered light traces the micron-sized particles at the disk surface that are well coupled to the gas. This means that scattered light observations can be used to trace the distribution of the disk gas and its interaction with embedded perturbers. Methods. 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. Results. 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. We also find a correlation between disk mass, as inferred from millimeter observations, and the detection of scattered light signal. Finally, we find a tentative correlation between relative disk-to-star brightness in scattered light and the presence of a dust cavity in the inner (unresolved) disk, as traced by the system spectral energy distribution. At the same time, faint disks in our sample are generally younger than 2 Myr. [ABSTRACT FROM AUTHOR]

We propose a model of deep and prolonged eclipses in young UX Ori stars. Some of these events last for decades and existing models cannot explain them. We show that such eclipses can be caused by the infall of gas and dust clumps from the remnants of the protostellar cloud onto the protoplanetary disk. The perturbation in the disk caused by the infall of a clump leads to a burst of accretional stellar activity and, as a consequence, to a strengthening of the disk wind. If the circumstellar disk is tilted at a small angle to the line of sight, then the dust raised by the wind from the disk surface can cause a dramatic decrease in the stellar brightness that can last for decades. [ABSTRACT FROM AUTHOR]

FU Ori-type stars (FUors) are young stellar objects (YSOs) experiencing luminosity outbursts by a few orders of magnitude, which last for ∼102 yr. A dozen of FUors are known up to date, but many more currently quiescent YSOs could have experienced such outbursts in the last ∼103 yr. To find observational signatures of possible past outbursts, we utilize andes, radmc-3d code as well as casa ALMA simulator to model the impact of the outburst on the physical and chemical structure of typical FU Ori systems and how it translates to the molecular lines' fluxes. We identify several combinations of molecular lines that may trace past FU Ori objects both with and without envelopes. The most promising outburst tracers from an observational perspective are the molecular flux combinations of the N2H+ J  = 3–2, C18O J  = 2–1, H2CO |$(J_{\rm K_a, K_c}) = 4_{04}-3_{03}$|⁠ , and HCN J  = 3–2 lines. We analyse the processes leading to molecular flux changes and show that they are linked with either thermal desorption or enhanced chemical reactions in the molecular layer. Using observed CO, HCN, N2H+, and H2CO line fluxes from the literature, we identify ten nearby disc systems that might have undergone FU Ori outbursts in the past ∼103 yr: [MGM2012] 556, [MGM2012] 371, and [MGM2012] 907 YSOs in L1641, Class II protoplanetary discs around CI Tau, AS 209, and IM Lup and transitional discs DM Tau, GM Aur, LkCa 15, and J1640-2130. [ABSTRACT FROM AUTHOR]

How protoplanetary discs evolve remains an unanswered question. Competing theories of viscosity and magnetohydrodynamic disc winds have been put forward as the drivers of angular momentum transport in protoplanetary discs. These two models predict distinct differences in the disc mass, radius, and accretion rates over time, that could be used to distinguish them. However that expectation is built on models that do not include another important process – photoevaporation, both internally by the host star and externally by neighbouring stars. In this work we produce numerical models of protoplanetary discs including viscosity, magnetohydrodynamic disc winds, and internal and external photoevaporation. We find that even weak levels of external photoevaporation can significantly affect the evolution of protoplanetary discs, influencing the observable features such as disc radii, that might otherwise distinguish between viscous and wind driven discs. Including internal photoevaporation further suppresses differences in evolution between viscous and wind driven discs. This makes it much more difficult than previously anticipated, to use observations of nearby star forming regions to determine whether discs are viscous or wind driven. Interestingly we find that evolved protoplanetary discs in intermediate FUV environments may be the best cases for differentiating whether they evolve through viscosity or magnetohydrodynamic disc winds. Ultimately this work demonstrates the importance of understanding what are the key evolutionary processes and including as many of those as possible when exploring the evolution of protoplanetary discs. [ABSTRACT FROM AUTHOR]

Aims. We aim to investigate the origin of the HI Brγ emission in young stars by using GRAVITY to image the innermost region of circumstellar disks, where important physical processes such as accretion and winds occur. With high spectral and angular resolution, we focus on studying the continuum and the HI Brγ-emitting area of the Herbig star HD 58647. Methods. Using VLTI-GRAVITY, we conducted observations of HD 58647 with both high spectral and high angular resolution. Thanks to the extensive uv coverage, we were able to obtain detailed images of the circumstellar environment at a sub-au scale, specifically capturing the continuum and the Brγ-emitting region. Through the analysis of velocity-dispersed images and photocentre shifts, we were able to investigate the kinematics of the HI Brγ-emitting region. Results. The recovered continuum images show extended emission where the disk major axis is oriented along a position angle of 14°. The size of the continuum emission at 5-σ levels is ~1.5 times more extended than the sizes reported from geometrical fitting (3.69 mas ± 0.02 mas). This result supports the existence of dust particles close to the stellar surface, screened from the stellar radiation by an optically thick gaseous disk. Moreover, for the first time with GRAVITY, the hot gas component of HD 58647 traced by the Brγ has been imaged. This allowed us to constrain the size of the Brγ-emitting region and study the kinematics of the hot gas; we find its velocity field to be roughly consistent with gas that obeys Keplerian motion. The velocity-dispersed images show that the size of the hot gas emission is from a more compact region than the continuum (2.3 mas ± 0.2 mas). Finally, the line phases show that the emission is not entirely consistent with Keplerian rotation, hinting at a more complex structure in the hot gaseous disk. [ABSTRACT FROM AUTHOR]

Context. Transition disks have large central cavities that have been resolved by imaging surveys during recent years. Cavities and other substructures in circumstellar disks are often interpreted as signposts to massive companions. Detecting companions at small angular separations is challenging with coronagraphic imaging observations. Aims. We aim to search for stellar and substellar companions in the central regions of transition disks. Such companions could be responsible for the large dust-depleted cavities. We want to determine if these disks might be circumbinary in their nature, similar to the HD 142527 system. Methods. We observed four systems, HD 100453, HD 100546, HD 135344 B, and PDS 70, with the sparse aperture masking mode of VLT/SPHERE, also leveraging the star-hopping method with the adaptive optics system. We extracted the complex visibilities and bispectra from the H2 and H3 imaging data. A binary model was fit to the closure phases to search for companions and estimate detection limits. For validation, we also analyzed four archival datasets of HD 142527 and inferred the orbital elements and atmospheric parameters of its low-mass stellar companion. Results. We have not detected any significant point sources in the four observed systems. With a contrast sensitivity of ≈0.004, we can rule out stellar companions down to ≈2 au and partially explore the substellar regime at separations ≳3–5 au. The analysis of HD 142527 B revealed that its projected orbit is aligned with dust features in the extended inner disk and that the mutual inclination with the outer disk is close to coplanar for one of the two solutions. Atmospheric modeling confirms the low-gravity and slightly reddened spectral appearance (Teff ≈ 3300 K, log 푔 ≈ 3.7, and AV ≈ 0.7). The inferred and derived bulk parameters (log L*/L0 ≈ −0.65, M* ≈ 0.4 M⊙, and R* ≈ 1.46 R⊙) are in agreement with dynamical constraints and evolutionary tracks. Conclusions. In contrast to HD 142527, we find no evidence that a close-in stellar companion is responsible for the resolved disk features of HD 100453, HD 100546, HD 135344 B, and PDS 70. Instead of a dynamical effect by a stellar companion, the formation of giant planets or even low-mass brown dwarfs could be shaping the innermost environment (≲20 au) of these circumstellar disks, as is the case with the planetary system of PDS 70. [ABSTRACT FROM AUTHOR]

Context. Observations of interstellar material infalling onto star- and planet-forming systems have become increasingly common thanks to recent advancements in radio interferometry. These structures have the potential to alter the dynamics of protoplanetary disks significantly by triggering the formation of substructures, inducing shocks, and modifying their physical and chemical properties. Moreover, the protoplanetary disks are replenished with new material, which increases the overall mass budget for planet formation. Aims. In this study, we combine new ALMA band 3 and archival band 6 observations to characterize the dust content and infall rate of a 4000 au arc-like structure that is infalling onto [MGM2012] 512 (hereafter M512), a class I young stellar object located in the Lynds 1641 region of the Orion A molecular cloud. Methods. We detected the extended dust emission from this structure in both ALMA bands. We tested whether the velocity pattern of the streamer is consistent with infalling trajectories by means of analytical streamline models. We measured spectral index maps for the first time and derived a dust opacity-index profile along a streamer. We constrained its grain properties and mass. Results. We find that the arc structure is consistent with infalling motions. We measure a spectral index α ~ 3.2 across the entire structure and a dust opacity index β ~ 1.6. Considering grain properties consistent with the measured β, the structure can host up to 245 M⊕ of dust, which exceeds or is comparable to the mass of the inner unresolved 600 au, which contains the protoplanetary disk of M512. Assuming a typical dust-to-gas ratio of 1% for the streamer, the free-fall timescales (50 kyr) imply total mass-infall rates up to 1.5 × 10−6M⊙ yr−1. M512 has been classified as an outbursting source with multi-epoch WISE photometry. It is thus an interesting case study for exploring the possible connection between infalling streamers and accretion outbursts. Conclusions. M512 is a unique source for which dust continuum emission of an arc-like streamer extending out to 4000 au can be characterized in a dual-band analysis. The dust properties are similar to those in the interstellar medium and imply a high dust mass. A massive streamer like this can affect the evolution of the star- and planet-forming inner system strongly. [ABSTRACT FROM AUTHOR]

Context. Nearby associations are ideal regions to study coeval samples of protoplanetary and debris disks down to late M-type stars. Those aged 5–10 Myr, where most of the disk should have already dissipated forming planets, are of particular interest. Aims. We present the first complete study of both protoplanetary and debris disks in a young region, using the η Chamaeleontis (η Cha) association as a test bench to study the cold disk content. We obtained submillimeter data for the entire core population down to late M-type stars, plus a few halo members. Methods. We performed a continuum submillimeter survey with APEX/LABOCA of all the core populations of the η Cha association. These data were combined with archival multiwavelength photometry to compile a complete spectral energy distribution. The disk properties were derived by modeling protoplanetary and debris disks using RADMC 2D and DMS, respectively. We compute a lower limit of the disk millimeter fraction, which is then compared to the corresponding disk fraction in the infrared for η Cha. We also revisit and refine the age estimate for the region, using the Gaia eDR3 astrometry and photometry for the core sources. Results. We find that protoplanetary disks in η Cha typically have holes with radii on the order of 0.01–0.03 AU, while ring-like emission from the debris disks is located between 20 and 650 au from the central star. The parallaxes and Gaia eDR3 photometry, in combination with the PARSEC and COLIBRI isochrones, enable us to confirm an age of η Cha between 7 and 9 Myr. In general, the disk mass seems insufficient to support accretion over a long time, even for the lowest mass accretors, a clear difference with other regions and also a sign that the mass budget is further underestimated. We do not find a correlation between the stellar masses, accretion rates, and disk masses, although this could be due to sample issues (very few, mostly low-mass objects). We confirm that the presence of inner holes is not enough to stop accretion unless accompanied by dramatic changes to the total disk mass content. Comparing η Cha with other regions at different ages, we find that the physical processes responsible for debris disks (e.g., dust growth, dust trapping) efficiently act in less than 5 Myr. [ABSTRACT FROM AUTHOR]

Context. The vertical shear instability (VSI) is a hydrodynamical instability predicted to produce turbulence in magnetically inactive regions of protoplanetary disks. The regions in which this instability can occur and the physical phenomena leading to its saturation are a current matter of research. Aims. We explore the secondary instabilities triggered by the nonlinear evolution of the VSI and their role in its saturation. We also expand on previous investigations on stability regions by considering temperature stratifications enforced by stellar irradiation and radiative cooling, and including the effects of dust-gas collisions and molecular line emission. Methods. We modeled the gas-dust mixture in a circumstellar disk around a T Tauri star by means of high-resolution axisymmetric radiation-hydrodynamical simulations including stellar irradiation with frequency-dependent opacities, considering different degrees of depletion of small dust grains. Results. The flow pattern produced by the interplay of the axisymmetric VSI modes and the baroclinic torque forms bands of nearly uniform specific angular momentum. In the high-shear regions in between these bands, the Kelvin–Helmholtz instability (KHI) is triggered. A third instability mechanism, consisting of an amplification of eddies by baroclinic torques, forms meridional vortices with Mach numbers up to ∼0.4. Our stability analysis suggests that protoplanetary disks can be VSI-unstable in surface layers up to tens of au for reasonably high gas emissivities. Conclusions. The significant transfer of kinetic energy to small-scale eddies produced by the KHI and possibly even the baroclinic acceleration of eddies limit the maximum energy of the VSI modes, likely leading to the saturation of the VSI. Depending on the gas molecular composition, the VSI can operate at surface layers even in regions where the midplane is stable. This picture is consistent with current observations of disks showing thin midplane millimeter-sized dust layers while appearing vertically extended in optical and near-infrared wavelengths. [ABSTRACT FROM AUTHOR]

Aims. In this work, we study the silicate dust content in the disk of one of the youngest eruptive stars, V900 Mon, at the highest angular resolution, probing down to the inner 10 au of said disk, and study the historical evolution of the system, traced in part by a newly discovered emission clump. Methods. We performed high angular resolution MIR interferometric observations of V900 Mon with MATISSE/VLTI with a spatial coverage ranging from 38 to 130-m baselines, and compared them to archival MIDI/VLTI data. We also mined and re-analyzed archival optical and infrared photometry of the star to study its long-term evolution since its eruption in the 1990s. We complemented our findings with integral field spectroscopy data from MUSE/VLT. Results. The MATISSE/VLTI data suggest a radial variation in the silicate feature in the dusty disk, whereby at large spatial scales (≥10 au) the protostellar disk's emission is dominated by large-sized (≥1μm) silicate grains, while at smaller spatial scales and closer to the star (≤5 au) silicate emission is absent, suggesting self-shielding. We propose that the self-shielding may be the result of small dust grains at the base of the collimated CO outflow previously detected by ALMA. A newly discovered knot in the MUSE/VLT data, located at a projected distance approximately 27 000 au from the star, is co-aligned with the molecular gas outflow at a P.A. of 250°(±5°) consistent with the position angle and inclination of the disk. The knot is seen in emission in Hα [N II], and the [S II] doublet and its kinematic age is about 5150 yr. This ejected material could originate from a previous eruption. [ABSTRACT FROM AUTHOR]

We present a new fried grid of mass-loss rates for externally far-ultraviolet (FUV) irradiated protoplanetary discs. As a precursor to the new grid, we also explore the microphysics of external photoevaporation, determining the impact of polycyclic aromatic hydrocarbon (PAH) abundance, metallicity, coolant depletion (via freeze out and radial drift), and grain growth (depletion of small dust in the outer disc) on disc mass-loss rates. We find that metallicity variations typically have a small effect on the mass-loss rate, since the impact of changes in heating, cooling and optical depth to the disc approximately cancel out. The new fried grid therefore focuses on (i) expanding the basic physical parameter space (disc mass, radius, UV field, stellar mass), (ii) on enabling variation of the the PAH abundance, and (iii) including an option for grain growth to have occurred or not in the disc. What we suggest is the fiducial model is comparable to the original fried grid. When the PAH-to-dust ratio is lower, or the dust in the wind more abundant, the mass-loss rate can be substantially lower. We demonstrate with a small set of illustrative disc evolutionary calculations that this in turn can have a significant impact on the disc mass/radius/ evolution and lifetime. [ABSTRACT FROM AUTHOR]

Most stars are born in dense stellar environments, where the formation and early evolution of planetary systems may be significantly perturbed by encounters with neighbouring stars. To investigate the fate of circumstellar gas discs and planets around young stars dense stellar environments, we numerically evolve star–disc–planet systems. We use the N -body codes NBODY6++ GPU and SnIPES for the dynamical evolution of the stellar population, and the SPH -based code GaSPH for the dynamical evolution of protoplanetary discs. The secular evolution of a planetary system in a cluster differs from that of a field star. Most stellar encounters are tidal, adiabatic, and nearly-parabolic. The parameters that characterize the impact of an encounter include the orientation of the protoplanetary disc and planet relative to the orbit of the encountering star, and the orbital phase, and the semimajor axis of the planet. We investigate this dependence for close encounters (r p/ a ≤ 100, where r p is the periastron distance of the encountering star and a is the semimajor axis of the planet). We also investigate distant perturbers (r p/ a ≫ 100), which have a moderate effect on the dynamical evolution of the planet and the protoplanetary disc. We find that the evolution of protoplanetary discs in star clusters differs significantly from that of isolated systems. When interpreting the outcome of the planet formation process, it is thus important to consider their birth environments. [ABSTRACT FROM AUTHOR]

Drastic changes in protoplanets' orbits could occur in the early stages of planetary systems through interactions with other planets and their surrounding protoplanetary or debris discs. The resulting planetary system could exhibit orbits with moderate to high eccentricities and/or inclinations, causing planets to perturb one another as well as the disc significantly. The present work studies the evolution of systems composed of an initially inclined planet and a debris disc. We perform N -body simulations of a narrow, self-gravitating debris disc, and a single interior Neptune-like planet. We simulate systems with various initial planetary inclinations, from coplanar to polar configurations considering different separations between the planet and the disc. We find that except when the planet is initially on a polar orbit, the planet–disc system tends to reach a quasi-coplanar configuration with low vertical dispersion in the disc. When present, the Zeipel–Kozai–Lidov oscillations induced by the disc pump the planet's eccentricity and, in turn, affect the disc structure. We also find that the resulting disc morphology in most of the simulations looks very similar in both radial and vertical directions once the simulations are converged. This contrasts strongly with massless disc simulations, where vertical disc dispersion is set by the initial disc-planet inclination and can be high for initially highly inclined planets. The results suggest caution in interpreting an unseen planet's dynamical history based only on the disc's appearance. [ABSTRACT FROM AUTHOR]

Based on observations obtained with the Nordic Optical Telescope (NOT) we investigate the spectral variability of the Herbig Ae star RR Tau. This star belongs to the UX Ori family, characterized by very deep fadings caused by the screening of the star with opaque fragments (clouds) of the protoplanetary discs. At the moments of such minima one observes strong spectral variability due to the fact that the dust cloud occults, for an observer, not only the star but also a part of the region where the emission spectrum originates. We calculated a series of obscuration models to interpret the observed variability of the H α line parameters. We consider two main obscuration scenarios: (1) the dust screen rises vertically above the circumstellar disc, and (2) the screen intersects the line-of-sight moving azimuthally with the disc. In both cases, the model of the emission region consists of a compact magnetosphere and a magnetocentrifugal disc wind. Comparison with observations shows that the first scenario explains well the variability of the radiation flux, the equivalent width, as well as the asymmetry of the H α line during eclipses, while the second scenario explains them only partly. This permits us to suggest that in the case of RR Tau, the main causes of the eclipses are either a structured disc wind, or the charged dust lifted along the field lines of the poloidal component of the magnetic field of the circumstellar disc. [ABSTRACT FROM AUTHOR]

The Orion Nebula Cluster (ONC) hosts protoplanetary disks experiencing external photoevaporation by the cluster's intense UV field. These "proplyds" are comprised of a disk surrounded by an ionization front. We present ALMA Band 3 (3.1 mm) continuum observations of 12 proplyds. Thermal emission from the dust disks and free–free emission from the ionization fronts are both detected, and the high-resolution (0.″057) of the observations allows us to spatially isolate these two components. The morphology is unique compared to images at shorter (sub)millimeter wavelengths, which only detect the disks, and images at longer centimeter wavelengths, which only detect the ionization fronts. The disks are small (r d = 6.4–38 au), likely due to truncation by ongoing photoevaporation. They have low spectral indices (α ≲ 2.1) measured between Bands 7 and 3, suggesting the dust emission is optically thick. They harbor tens of Earth masses of dust as computed from the millimeter flux using the standard method although their true masses may be larger due to the high optical depth. We derive their photoevaporative mass-loss rates in two ways: first, by invoking ionization equilibrium and second, by using the brightness of the free–free emission to compute the density of the outflow. We find decent agreement between these measurements and M ̇ = 0.6–18.4 × 10−7 M ⊙ yr−1. The photoevaporation timescales are generally shorter than the ∼1 Myr age of the ONC, underscoring the known "proplyd lifetime problem." Disk masses that are underestimated due to being optically thick remains one explanation to ease this discrepancy. [ABSTRACT FROM AUTHOR]

Context. A characteristic feature of young stellar objects is their variability, which is caused by a variety of different physical processes. High-resolution interferometric observations in the near- and mid-infrared wavelength ranges spanning multiple epochs allow the detailed study of these processes. Aims. We aim at investigating the expected variations of the interferometric observables connected to changes in the measured photometric fluxes of a typical variable accreting central young stellar object with a circumstellar disk. Methods. We calculated visibilities and closure phases as well as the photometric flux of brightness distributions obtained using 3D Monte Carlo radiative transfer simulations for a model of a circumstellar disk with an accreting central star. Results. Changes in the accretion luminosity of the central object, that is, an accreting pre-main-sequence star, can lead to significant variations in the visibility and closure phase of the star-disk system measured with instruments at the Very Large Telescope Interferometer (VLTI) that can be related to changes in the photometric flux. Taking into account additional effects due to baseline variation, interferometric observations can provide valuable contributions to the understanding of the underlying processes. Additionally, we provide the web application VLTI B-VAR that allows the impact of the hour angle on the visibility and closure phase for customized intensity maps to be estimated. [ABSTRACT FROM AUTHOR]

Context. The abundance and distribution of ice in protoplanetary disks is critical for an understanding of the link between the composition of circumstellar matter and the composition of exoplanets. Edge-on protoplanetary disks are a useful tool for constraining this ice composition and its location in the disk because the spectral signatures of the ice can be observed in absorption against the continuum emission that arises from the warmer regions in the central disk. Aims. The aim of this work is to model ice absorption features in protoplanetary disks and to determine how well the abundance of the main ice species throughout the disk can be determined within the uncertainty of the physical parameter space. The edge-on proto-planetary disk around HH 48 NE, a target of the James Webb Space Telescope Early Release program Ice Age, is used as a reference system. Methods. We used the full anisotropic scattering capabilities of the radiative transfer code RADMC-3D to ray-trace the mid-infrared continuum. Using a constant parameterized ice abundance, we added ice opacities to the dust opacity in regions in which the disk was cold enough for the main carbon, oxygen, and nitrogen carriers to freeze out. Results. The global abundance relative to the dust content of the main ice carriers in HH 48 NE can be determined within a factor of 3 when the uncertainty of the physical parameters is taken into account. Ice features in protoplanetary disks can be saturated at an optical depth of ≲1 due to local saturation. Ices are observed at various heights in the disk model, but in this model, spatial information is lost for features at wavelengths >7 µm when observing with James Webb Space Telescope because the angular resolution decreases towards longer wavelengths. Spatially observed ice optical depths cannot be directly related to column densities, as would be the case for direct absorption against a bright continuum source, because of radiative transfer effects. Vertical snowlines will not be a clear transition because the height of the snow surface increases radially, but their location may be constrained from observations using radiative transfer modeling. Radial snowlines are not really accessible. Not only the ice abundance, but also the inclination, the settling, the grain size distribution, and the disk mass have a strong impact on the observed ice absorption features in disks. Relative changes in the ice abundance can only be inferred from observations if the source structure is well constrained. [ABSTRACT FROM AUTHOR]

Context. The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars (M ≳ 8 M⊙) are still not well understood. A growing number of case studies have found massive young stellar objects (MYSOs) to harbour disk-like structures, painting a picture that the formation of high-mass stars may proceed through disk accretion, similar to that of lower-mass stars. However, the properties of such structures have yet to be uniformly and systematically characterised. Aims. The aim of this work is to uniformly study the kinematic properties of a large sample of MYSOs and characterise the stability of possible circumstellar disks against gravitational fragmentation. Methods. We have undertaken a large observational programme (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous (L > 104L⊙) protostellar objects in the 1.37 mm wavelength regime in both continuum and spectral line emission, reaching 0.4″ resolution (800 au at 2 kpc). Results. We present the gas kinematics of the full sample and detect dense gas emission surrounding 15 regions within the CORE sample. Using the dense gas tracer CH3CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. The extent of the CH3CN emission tracing the disk candidates varies from 1800 to 8500 au. Analysing the free-fall to rotational timescales, we find that the sources are rotationally supported. The rotation profiles of some disk candidates are well described by differential rotation while for others the profiles are poorly resolved. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of ~ 10–25 M⊙. Modelling the level population of CH3CN (12K–11K) K = 0–6 lines, we present temperature maps and find median temperature in the range 70–210 K with a diversity in distributions. Radial profiles of the specific angular momentum (j) for the best disk candidates span a range of 1–2 orders of magnitude, on average ~10−3 km s−1 pc, and they follow j ∝ r1.7, which is consistent with a poorly resolved rotating and infalling envelope-disk model. Studying the Toomre stability of the disk candidates, we find almost all (11 out of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations, as a result of their high disk to stellar mass ratio. In particular, disks with masses greater than ~ 10–20% of the mass of their host (proto)stars are Toomre unstable, and more luminous YSOs tend to have disks that are more massive compared to their host star and hence more prone to fragmentation. Conclusions. In this work, we show that most disk structures around high-mass YSOs are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars. [ABSTRACT FROM AUTHOR]

Circumstellar discs likely have a short window when they are self-gravitating and prone to the effects of disc instability, but during this time the seeds of planet formation can be sown. It has long been argued that disc fragmentation can form large gas giant planets at wide orbital separations, but its place in the planet formation paradigm is hindered by a tendency to form especially large gas giants or brown dwarfs. We instead suggest that planet formation can occur early in massive discs, through the gravitational collapse of dust which can form the seeds of giant planets. This is different from the usual picture of self-gravitating discs, in which planet formation is considered through the gravitational collapse of the gas disc into a gas giant precursor. It is familiar in the sense that the core is formed first, and gas is accreted thereafter, as is the case in the core accretion scenario. However, by forming a ∼1 M⊕ seed from the gravitational collapse of dust within a self-gravitating disc there exists the potential to overcome traditional growth barriers and form a planet within a few times 105 yr. The accretion of pebbles is most efficient with centimetre-sized dust, but the accretion of millimetre sizes can also result in formation within a Myr. Thus, if dust can grow to these sizes, planetary seeds formed within very young, massive discs could drastically reduce the time-scale of planet formation and potentially explain the observed ring and gap structures in young discs. [ABSTRACT FROM AUTHOR]

Photoevaporative disc winds play a key role in our understanding of circumstellar disc evolution, especially in the final stages, and they might affect the planet formation process and the final location of planets. The study of transition discs (i.e. discs with a central dust cavity) is central for our understanding of the photoevaporation process and disc dispersal. However, we need to distinguish cavities created by photoevaporation from those created by giant planets. Theoretical models are necessary to identify possible observational signatures of the two different processes, and models to find the differences between the two processes are still lacking. In this paper, we study a sample of transition discs obtained from radiation–hydrodynamic simulations of internally photoevaporated discs, and focus on the dust dynamics relevant for current Atacama Large Millimetre Array observations. We then compared our results with gaps opened by a super-Earth/giant planets, finding that the photoevaporated cavity steepness depends mildly on gap size, and it is similar to that of a |${1}\, {\rm M_J}$| planet. However, the dust density drops less rapidly inside the photoevaporated cavity compared to the planetary case due to the less efficient dust filtering. This effect is visible in the resulting spectral index, which shows a larger spectral index at the cavity edge and a shallower increase inside it with respect to the planetary case. The combination of cavity steepness and spectral index might reveal the true nature of transition discs. [ABSTRACT FROM AUTHOR]

We simulate the formation and collapse of pre-stellar cores at few-au resolution in a set of radiation-magnetohydrodynamic simulations of giant molecular clouds (GMCs) using the grid-based code RAMSES-RT. We adopt, for the first time to our best knowledge, realistic initial/boundary conditions by zooming in on to individual massive pre-stellar cores within the GMC. We identify two distinct modes of fragmentation: 'quasi-spherical' and 'filamentary'. In both modes, the fragments eventually become embedded in a quasi-steady accretion disc or toroid with radii ∼500–5000 au and opening angles H / R ∼ 0.5 − 1. The discs/toroids are Toomre stable but the accreted pre-existing fragments are found orbiting the outer disc, appearing as disc fragmentation. Each core converts nearly 100 per cent of the gas mass into a few massive stars forming near the disc centre. Large and massive discs around high-mass stars are supported by magnetic pressure in the outer disc, at radii >200–1000 au, and turbulent pressure in the inner disc. The most massive core accretes several times more mass than its initial mass, forming a cluster of 8 massive (proto)stars enshrouded by a toroid, suggesting a competitive accretion scenario for the formation of stars above ∼30 M⊙. We also find that the H  ii regions produced by a single massive star remain trapped in the dense circumstellar discs for a few hundred kiloyears, while the dynamic motions of massive stars in wide binaries or multiple systems displace the stars from the densest parts of the disc, allowing UV radiation to escape producing steady or pulsating bipolar H  ii regions. [ABSTRACT FROM AUTHOR]

While dust disks around optically visible, Class II protostars are found to be vertically thin, when and how dust settles to the midplane are unclear. As part of the Atacama Large Millimeter/submillimeter Array large program, Early Planet Formation in Embedded Disks, we analyze the edge-on, embedded, Class I protostar IRAS 04302+2247, also nicknamed the "Butterfly Star." With a resolution of 0.″05 (8 au), the 1.3 mm continuum shows an asymmetry along the minor axis that is evidence of an optically thick and geometrically thick disk viewed nearly edge-on. There is no evidence of rings and gaps, which could be due to the lack of radial substructure or the highly inclined and optically thick view. With 0.″1 (16 au) resolution, we resolve the 2D snow surfaces, i.e., the boundary region between freeze-out and sublimation, for 12CO J = 2–1, 13CO J = 2–1, C18O J = 2–1, H 2CO J = 30,3–20,2, and SO J = 65–54, and constrain the CO midplane snow line to ∼130 au. We find Keplerian rotation around a protostar of 1.6 ± 0.4 M ⊙ using C18O. Through forward ray-tracing using RADMC-3D, we find that the dust scale height is ∼6 au at a radius of 100 au from the central star and is comparable to the gas pressure scale height. The results suggest that the dust of this Class I source has yet to vertically settle significantly. [ABSTRACT FROM AUTHOR]

The process of accretion through circumstellar disks in young stellar objects is an integral part of star formation and the H α emission line is a prominent signature of accretion in low-mass stars. We present the detection and characterization of H α emission line sources in the central region of a distant, low-metallicity young stellar cluster Dolidze 25 (at ∼ 4.5 kpc) using medium-resolution optical spectra (4750–9350 Å) obtained with the multi-unit spectroscopic explorer (MUSE) at the VLT. We have identified 14 potential accreting sources within a rectangular region of ( 2 ′ × 1 ′ ) towards the center of the cluster based on the detection of strong and broad emissions in H α as well as the presence of other emission lines, such as [OI] and H β . Based on their positions in both photometric color–magnitude and color–color diagrams, we have also confirmed that these objects belong to the pre-main sequence phase of star formation. Our results were compared with the disk and diskless members of the cluster previously identified by Guarcello et al. (2021) using near-IR colors, and all sources they had identified as disks were confirmed to be accreting based on the spectroscopic characteristics. [ABSTRACT FROM AUTHOR]

We demonstrate that planet formation via pebble accretion is sensitive to external photoevaporation of the outer disc. In pebble accretion, planets grow by accreting from a flux of solids (pebbles) that radially drift inwards from the pebble production front. If external photoevaporation truncates the outer disc fast enough, it can shorten the time before the pebble production front reaches the disc outer edge, cutting off the supply of pebble flux for accretion, hence limiting the pebble mass reservoir for planet growth. Conversely, cloud shielding can protect the disc from strong external photoevaporation and preserve the pebble reservoir. Because grain growth and drift can occur quickly, shielding even on a short time-scale (<1 Myr) can have a non-linear impact on the properties of planets growing by pebble accretion. For example, a |$10^{-3}\, \mathrm{ M}_{\oplus }$| planetary seed at 25 au stays at 25 au with a lunar mass if the disc is immediately irradiated by a 103 G0 field, but grows and migrates to be approximately Earth-like in both mass and orbital radius if the disc is shielded for just 1 Myr. In NGC 2024, external photoevaporation is thought to happen to discs that are <0.5 Myr old, which coupled with the results here suggests that the exact planetary parameters can be very sensitive to the star-forming environment. Universal shielding for time-scales of at least |${\sim} 1.5\,$|  Myr would be required to completely nullify the environmental impact on planetary architectures. [ABSTRACT FROM AUTHOR]

The 12CO rotational lines in protoplanetary discs are good tracers of the total spatial extension of the gas component, and potentially planet–disc interactions. We present ALMA long baseline observations of the 12CO(2-1) line of 10 protoplanetary discs from the Ophiuchus DIsc Survey Employing ALMA (ODISEA) project, aiming to set constraints on the gas distribution of these sources. The position angle of the gaseous disc can be inferred for five sources using high-velocity channels, which trace the gas in the inner part of the disc. We compare the high-velocity PAs to the orientations inferred from the continuum, representative of the orientation over ∼53 to 256 au in these resolved discs. We find a significant difference in orientation for DoAr 44, which is evidence of a tilted inner disc. Eight discs show evidence of gas inside inner dust cavities or gaps, and the disc of ISO-Oph 196 is not detected in 12CO(2-1), except for the compact signal located inside its dust cavity. Our observations also point out a possible outflow in WLY 2-63. [ABSTRACT FROM AUTHOR]

Planetesimals – asteroids and comets – are the building blocks of planets in protoplanetary discs and the source of dust, ice, and gas in debris discs. Along with planets they comprise the left-over material after star formation that constitutes a planetary system. Planets influence the dynamics of planetesimals, sculpting the orbits of debris belts to produce asymmetries or gaps. We can constrain the architecture of planetary systems, and infer the presence of unseen planetary companions, by high spatial resolution imaging of debris discs. HD 16743 is a relatively young F-type star that hosts a bright edge-on debris disc. Based on far-infrared Herschel observations its disc was thought to be stirred by a planetary companion. Here, we present the first spatially resolved observations at near-infrared and millimetre wavelengths with HST and ALMA, revealing the disc to be highly inclined at |$87{_{.}^{\circ}} 3~^{+1{_{.}^{\circ}} 9}_{-2{_{.}^{\circ}} 5}$| with a radial extent of 157.7 |$^{+2.6}_{-1.5}$| au and an full width at half maximum of 79.4 |$^{+8.1}_{-7.8}$| au (Δ R / R  = 0.5). The vertical scale height of the disc is 0.13 ± 0.02, significantly greater than typically assumed unstirred value of 0.05, and could be indicative of stirring of the dust-producing planetesimals within the disc by bodies at least a few times the mass of Pluto up to 18.3 M⊕ in the single object limit. [ABSTRACT FROM AUTHOR]

We study the formation and early evolution of young stellar objects (YSOs) using three-dimensional non-ideal magnetohydrodynamic (MHD) simulations to investigate the effect of cosmic ray ionization rate and dust fraction (or amount of dust grains) on circumstellar disc formation. Our simulations show that a higher cosmic ray ionization rate and a lower dust fraction lead to (i) a smaller magnetic resistivity of ambipolar diffusion, (ii) a smaller disc size and mass, and (iii) an earlier timing of outflow formation and a greater angular momentum of the outflow. In particular, at a high cosmic ray ionization rate, the discs formed early in the simulation are dispersed by magnetic braking on a time-scale of about 104 yr. Our results suggest that the cosmic ray ionization rate has particularly a large impact on the formation and evolution of discs, while the impact of the dust fraction is not significant. [ABSTRACT FROM AUTHOR]

In the core accretion model, planetesimals grow by mutual collisions and engulfing millimeter-to-centimeter particles, i.e., pebbles. Pebble accretion can significantly increase the accretion efficiency and help explain the presence of planets on wide orbits. However, the pebble supply is typically parameterized as a coherent pebble mass flux, sometimes being constant in space and time. Here we solve the dust advection and diffusion within viciously evolving protoplanetary disks to determine the pebble supply self-consistently. The pebbles are then accreted by planetesimals interacting with the gas disk via gas drag and gravitational torque. The pebble supply is variable with space and decays with time quickly, with a pebble flux below 10 M ⊕ Myr−1 after 1 Myr in our models. As a result, only when massive planetesimals (>0.01 M ⊕) are luckily produced by the streaming instability or the disk has low viscosity (α ∼ 0.0001) can the herd of planetesimals grow over a Mars mass within 2 Myr. By then, planetesimals only capture pebbles about 50 times their mass and as little as 10 times beyond 20 au due to limited pebble supply. Further studies considering multiple dust species in various disk conditions are warranted to fully assess the realistic pebble supply and its influence on planetesimal growth. [ABSTRACT FROM AUTHOR]

The stellar cluster environment is expected to play a central role in the evolution of circumstellar disks. We use thermochemical modeling to constrain the dust and gas masses, disk sizes, UV and X-ray radiation fields, viewing geometries, and central stellar masses of 20 class II disks in the Orion Nebula Cluster (ONC). We fit a large grid of disk models to 350 GHz continuum, CO J = 3 − 2, and HCO+ J = 4 − 3 Atacama Large Millimeter/submillimeter Array observations of each target, and we introduce a procedure for modeling interferometric observations of gas disks detected in absorption against a bright molecular cloud background. We find that the ONC disks are massive and compact, with typical radii <100 au, gas masses ≥10−3 M ⊙, and gas-to-dust ratios ≥100. The interstellar‐medium‐like gas-to-dust ratios derived from our modeling suggest that compact, externally irradiated disks in the ONC are less prone to gas-phase CO depletion than the massive and extended gas disks that are commonly found in nearby low-mass star-forming regions. The presence of massive gas disks indicates that external photoevaporation may have only recently begun operating in the ONC; though it remains unclear whether other cluster members are older and more evaporated than the ones in our sample. Finally, we compare our dynamically derived stellar masses with the stellar masses predicted from evolutionary models and find excellent agreement. Our study has significantly increased the number of dynamical mass measurements in the mass range ≤0.5 M ⊙, demonstrating that the ONC is an ideal region for obtaining large samples of dynamical mass measurements toward low-mass M-dwarfs. [ABSTRACT FROM AUTHOR]

Close encounters between stars in star-forming regions are important as they can perturb or destroy protoplanetary disks, young planetary systems, and multiple-star systems. We simulate simple, virialized, equal-mass N -body star clusters and find that both the rate and the total number of encounters between stars vary by factors of several in statistically identical clusters due to the stochastic/chaotic details of the orbits and stellar dynamics. Encounters tend to "saturate" rapidly in the core of a cluster, with stars there each having many encounters, while more distant stars have none. However, we find that the fraction of stars that has had at least one encounter within a particular distance grows in the same way (scaling with the crossing time and half-mass radius) in all clusters, and we present a new (empirical) way of estimating the fraction of stars that has had at least one encounter at a particular distance. [ABSTRACT FROM AUTHOR]

Most stars are born in stellar clusters, and their protoplanetary disks, which are the birthplaces of planets, can, therefore, be affected by the radiation of nearby massive stars. However, little is known about the chemistry of externally irradiated disks, including whether or not their properties are similar to the so-far better-studied isolated disks. Motivated by this question, we present ALMA Band 6 observations of two irradiated Class II protoplanetary disks in the outskirts of the Orion Nebula Cluster to explore the chemical composition of disks exposed to (external) far-ultraviolet (FUV) radiation fields: the 216-0939 disk and the binary system 253-1536A/B, which are exposed to radiation fields of 10 ^2 –10 ^3 times the average interstellar radiation field. We detect lines from CO isotopologues, HCN, H _2 CO, and C _2 H toward both protoplanetary disks. Based on the observed disk-integrated line fluxes and flux ratios, we do not find significant differences between isolated and irradiated disks. The observed differences seem to be more closely related to the different stellar masses than to the external radiation field. This suggests that these disks are far enough away from the massive Trapezium stars, that their chemistry is no longer affected by external FUV radiation. Additional observations toward lower-mass disks and disks closer to the massive Trapezium stars are required to elucidate the level of external radiation required to make an impact on the chemistry of planet formation in different kinds of disks.

Scattered light imaging of protoplanetary disks provides key insights on the geometry and dust properties in the disk surface. Here, we present James Webb Space Telescope (JWST) 2–21 μ m images of a 1000 au radius edge-on protoplanetary disk surrounding an 0.4 M _⊙ young star in Taurus, Two Micron All Sky Survey (2MASS) J04202144 + 2813491. These observations represent the longest wavelengths at which a protoplanetary disk is spatially resolved in scattered light. We combine these observations with Hubble Space Telescope optical images and Atacama Large Millimeter/submillimeter Array continuum and CO mapping. We find that the changes in the scattered light disk morphology are remarkably small across a factor of 30 in wavelength, indicating that dust in the disk surface layers is characterized by an almost gray opacity law. Using radiative transfer models, we conclude that grains up to ≳10 μ m in size are fully coupled to the gas in this system, whereas grains ≳100 μ m are strongly settled toward the midplane. Further analyses of these observations, and similar ones of other edge-on disks, will provide strong empirical constraints on disk dynamics and evolution and grain growth models. In addition, the 7.7 and 12. μ m JWST images reveal an X-shaped feature located above the warm molecular layer traced by CO line emission. The highest elevations at which this feature is detectable roughly match the maximal extent of the disk in visible wavelength scattered light as well as of an unusual kinematic signature in CO. We propose that these phenomena could be related to a disk wind entraining small dust grains.

Context. Stable circumbinary discs around evolved post-asymptotic giant branch (post-AGB) binary systems composed of gas and dust show many similarities with protoplanetary discs around young stellar objects. These discs can provide constraints on both binary evolution and the formation of macrostructures within circumstellar discs. Here we focus on one post-AGB binary system: IRAS 08544-4431. Aims. We aim to refine the physical model of IRAS 08544-4431 with a radiative transfer treatment and continue the near-infrared and mid-infrared interferometric analysis covering the H, K, L, and N bands. Results from geometric modelling of these data in our previous study constrain the shape of the inner rim of the disc and its radial dust structure. We aim to capture the previously detected amount of over-resolved flux and the radial intensity profile at and beyond the inner dust disc rim to put constraints on the physical processes in the inner disc regions. Methods. We used a three-dimensional Monte Carlo radiative transfer code to investigate the physical structure of the disc by reproducing both the photometry and the multi-wavelength infrared interferometric dataset. We first performed a parametric study to explore the effect of the individual parameters and selected the most important parameters, which were then used in a thorough grid search to fit the structural characteristics. We developed a strategy to identify the models that were best able to reproduce our extensive multi-wavelength dataset. Results. We find a family of models that successfully fit the infrared photometric and interferometric data in all bands. These models show a flaring geometry with efficient settling. Larger grains are present in the inner disc as probed by our infrared interferometric observations. Some over-resolved flux component was recovered in all bands, but the optimised models still fall short in explaining all the over-resolved flux. This suggests that another dusty structure within the system that is not included in our models plays a role. The structure of this over-resolved component is unclear, but it has a colour temperature between 1400 and 3600 K. Conclusions. Multi-wavelength infrared interferometric observations of circumstellar discs allow the inner disc regions to be studied in unprecedented detail. The refined physical models can reproduce most of the investigated features, including the photometric characteristics, the radial extent, and the overall shape of the visibility curves. Our multi-wavelength interferometric observations combined with photometry show that the disc around IRAS 08544-4431 is similar to protoplanetary discs around young stars with similar dust masses and efficient dust growth. The resulting disc geometry is capable of reproducing part of the over-resolved flux, but to fully reproduce the over-resolved flux component, an additional component is needed. Multi-scale high-angular-resolution analysis combining VLTI, VLT/SPHERE, and ALMA data is needed to fully define the structure of the system. [ABSTRACT FROM AUTHOR]

Context. In recent years, a new hot topic has emerged in the star and planet formation field, namely, the interaction between the circumstellar disk and its birth cloud. The birth environments of young stars leave strong imprints on the star itself and their surroundings. In this context, we present a detailed analysis of the rich circumstellar environment around the young Herbig Ae/Be star T CrA. Aims. Our aim is to understand the nature of the stellar system and the extended circumstellar structures, as seen in scattered light images. Methods. We conducted our analysis on the basis of a set of combined archival data and new adaptive optics images at a high contrast and high resolution. Results. The scattered light images reveal the presence of a complex environment around T CrA, composed of a bright, forward-scattering rim of the disk's surface that is seen at very high inclinations, along with a dark lane of the disk midplane, bipolar outflows, and streamer features that are likely tracing infalling material from the surrounding birth cloud onto the disk. The analysis of the light curve suggests that the star is a binary with a period of 29.6 yr, confirming previous assertions based on spectro-astrometry. The comparison of the scattered light images with the ALMA continuum and 12CO (2–1) line emission shows that the disk is in Keplerian rotation and the northern side of the outflowing material is receding, while the southern side is approaching the observer. The overall system lies on different geometrical planes. The orbit of the binary star is perpendicular to the outflows and is seen edge on. The disk is itself seen edge-on, with a position angle of ~7°. 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 with the star. Modeling of the spectral energy distribution using a radiative transfer scheme is in good agreement with the proposed configuration, as well as the hydrodynamical simulation performed using a smoothed particle hydrodynamics code. Conclusions. We find evidence of streamers of accreting material around T CrA. These streamers connect the filament, along which T CrA is forming along 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. This impacts the accretion taking place in the components of the binary, favoring the growth of the primary with respect the secondary, in contrast to the case of aligned disks. [ABSTRACT FROM AUTHOR]

Context. The young massive-star-forming region M17 contains optically visible massive pre-main-sequence stars that are surrounded by circumstellar disks. Such disks are expected to disappear when these stars enter the main sequence. The physical and dynamical structure of these remnant disks are poorly constrained, especially the inner regions where accretion, photo-evaporation, and companion formation and migration may be ongoing. Aims. We aim to constrain the physical properties of the inner parts of the circumstellar disks of massive young stellar objects B243 (6 M⊙) and B331 (12 M⊙), two systems for which the central star has been detected and characterized previously despite strong dust extinction. Methods. Two-dimensional radiation thermo-chemical modelling with PRODIMO of double-peaked hydrogen lines of the Paschen and Brackett series observed with X-shooter was used to probe the properties of the inner disk of the target sources. The model was modified to treat these lines. Additionally, the dust structure was studied by fitting the optical and near-infrared spectral energy distribution. Results. B243 features a hot gaseous inner disk with dust at the sublimation radius at ~3 AU. The disk appears truncated at roughly 6.5 AU; a cool outer disk of gas and dust may be present, but it cannot be detected with our data. B331 also has a hot gaseous inner disk. A gap separates the inner disk from a colder dusty outer disk starting at up to ~100 AU. In both sources the inner disk extends to almost the stellar surface. Chemistry is essential for the ionization of hydrogen in these disks. Conclusions. The lack of a gap between the central objects and these disks suggests that they accrete through boundary-layer accretion. This would exclude the stars having a strong magnetic field. Their structures suggest that both disks are transitional in nature, that is to say they are in the process of being cleared, either through boundary-layer accretion, photo-evaporation, or through companion activity. [ABSTRACT FROM AUTHOR]

The size of dust grains, a , is key to the physical and chemical processes in circumstellar discs, but observational constraints of grain size remain challenging. (Sub)millimetre continuum observations often show a per cent-level polarization parallel to the disc minor axis, which is generally attributed to scattering by |${\sim}100\, \mu{\rm m}$| -sized spherical grains (with a size parameter x ≡ 2 |$\pi$| a /λ < 1, where λ is the wavelength). Larger spherical grains (with x greater than unity) would produce opposite polarization direction. However, the inferred size is in tension with the opacity index β that points to larger mm/cm-sized grains. We investigate the scattering-produced polarization by large irregular grains with a range of x greater than unity with optical properties obtained from laboratory experiments. Using the radiation transfer code, RADMC-3D, we find that large irregular grains still produce polarization parallel to the disc minor axis. If the original forsterite refractive index in the optical is adopted, then all samples can produce the typically observed level of polarization. Accounting for the more commonly adopted refractive index using the DSHARP dust model, only grains with x of several (corresponding to ∼mm-sized grains) can reach the same polarization level. Our results suggest that grains in discs can have sizes in the millimetre regime, which may alleviate the tension between the grain sizes inferred from scattering and other means. Additionally, if large irregular grains are not settled to the mid-plane, their strong forward scattering can produce asymmetries between the near and far side of an inclined disc, which can be used to infer their presence. [ABSTRACT FROM AUTHOR]