Circumstellar Disks Research Articles

Inverse-problem versus principal component analysis methods for angular differential imaging of circumstellar disks

Context. Circumstellar disk images have highlighted a wide variety of morphological features. Recovering disk images from high-contrast angular differential imaging (ADI) sequences is, however, generally affected by geometrical biases, leading to unreliable inferences of the morphology of extended disk features. Recently, two types of approaches have been proposed to recover more robust disk images from ADI sequences: iterative principal component analysis (I-PCA) and inverse problem (IP) approaches. Aims. We introduce mustard, a new IP-based algorithm specifically designed to address the problem of the flux invariant to rotation in ADI sequences – a limitation inherent to the ADI observing strategy – and discuss the advantages of IP approaches with respect to PCA-based algorithms. Methods. The mustard model relies on the addition of morphological priors on the disk and speckle field to a standard IP approach to tackle rotation-invariant signals in circumstellar disk images. We compared the performance of mustard, I-PCA, and standard PCA on a sample of high-contrast imaging data sets acquired in different observing conditions, after injecting a variety of synthetic disk models at different contrast levels. Results. Mustard significantly improves the recovery of rotation-invariant signals in disk images, especially for data sets obtained in good observing conditions. However, the mustard model inadequately handles unstable ADI data sets and provides shallower detection limits than PCA-based approaches. Conclusions. Mustard has the potential to deliver more robust disk images by introducing a prior to address the inherent ambiguity of ADI observations. However, the effectiveness of the prior is partly hindered by our limited knowledge of the morphological and temporal properties of the stellar speckle halo. In light of this limitation, we suggest that the algorithm could be improved by enforcing a data-driven prior based on a library of reference stars.

The first scattered light images of HD 112810, a faint debris disk in the Sco-Cen association

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 modeled 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 ± 9 au and an inclination angle of 75.7−1.3+1.1. This is in good agreement with the constraints from spectral energy distribution modeling and from a partially resolved ALMA image of the system. No planets are detected, though planets below the detection limit (~2.6 MJ at a projected separation of 118 au) 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.

Constraining the formation history of the TOI-1338/BEBOP-1 circumbinary planetary system

ABSTRACT The recent discovery of multiple planets in the circumbinary system TOI-1338/BEBOP-1 raises questions about how such a system formed. The formation of the system was briefly explored in the discovery paper, but only to answer the question do current pebble accretion models have the potential to explain the origin of the system? We use a global model of circumbinary planet formation that utilizes N-body simulations, including prescriptions for planet migration, gas and pebble accretion, and interactions with a circumbinary disc, to explore the disc parameters that could have led to the formation of the TOI-1338/BEBOP-1 system. With the disc lifetime being the main factor in determining how planets form, we limit our parameter space to those that determine the disc lifetime. These are: the strength of turbulence in the disc, the initial disc mass, and the strength of the external radiation field that launches photoevaporative winds. When comparing the simulated systems to TOI-1338/BEBOP-1, we find that only discs with low levels of turbulence are able to produce similar systems. The radiation environment has a large effect on the types of planetary systems that form, whilst the initial disc mass only has limited impact since the majority of planetary growth occurs early in the disc lifetime. With the most TOI-1338/BEBOP-1 like systems all occupying similar regions of parameter space, our study shows that observed circumbinary planetary systems can potentially constrain the properties of planet forming discs.

New transient ULX candidate in NGC 4254: evidence of circumbinary disc?

ABSTRACT This paper presents the identification of a new transient ultraluminous X-ray source (ULX) candidate (ULX-3) with reaching a peak luminosity of ∼4 × 1039 erg s−1 in NGC 4254 by using archival Chandra, Swift/X-Ray Telescope, Hubble Space Telescope, and JWST observations. From precise astrometric calculations, unique optical, near-infrared, and mid-infrared counterparts were found. The spectral energy distribution (SED) and colour–magnitude diagrams of counterparts of the new ULX candidate were plotted to constrain the nature of the possible donor star. Evidence of a circumbinary disc was found from its SED with two blackbody temperatures of 1000 and 200 K. Moreover, according to the X-ray hardness ratios, ULX-3 exhibits very hard to very soft transitions as seen in some high-mass X-ray binaries with Be-star donors. Moreover, ULX-3 varies by more than two orders of magnitude in the 0.3–10 keV energy band as seen in typical transient ULXs.

Planet migration in massive circumbinary discs

ABSTRACT Most stars are in multiple systems, with the majority of those being binaries. A large number of planets have been confirmed in binary stars and, therefore, it is important to understand their formation and dynamical evolution. We perform simulations to investigate the migration of wide-orbit giant planets (semimajor axis 100 au) in massive circumbinary discs (mass $0.1\, {\rm M}_{\odot }$) that are marginally gravitationally unstable, using the three-dimensional Smoothed Particle Hydrodynamic code seren. We vary the binary parameters to explore their effect on planet migration. We find that a planet in a massive circumbinary disc initially undergoes a period of rapid inward migration before switching to a slow outward migration, as it does in a circumstellar disc. However, the presence of the binary enhances planet migration and mass growth. We find that a high binary mass ratio (binary with equal mass stars) results in more enhanced outward planet migration. Additionally, larger binary separation and/or higher binary eccentricity results to a faster outward planet migration and stronger planet growth. We conclude that wide-orbit giant planets attain wider final orbits due to migration around binary stars than around single stars.

Disappearing thermal X-ray emission as a tell-tale signature of merging massive black hole binaries

ABSTRACT The upcoming Laser Interferometer Space Antenna (LISA) is expected to detect gravitational waves (GWs) from massive black hole binaries (MBHB). Finding the electromagnetic (EM) counterparts for these GW events will be crucial for understanding how and where MBHBs merge, measuring their redshifts, constraining the Hubble constant and the graviton mass, and for other novel science applications. However, due to poor GW sky localization, multiwavelength, time-dependent EM models are needed to identify the right host galaxy. We studied merging MBHBs embedded in a circumbinary disc (CBD) using high-resolution two-dimensional simulations, with a Γ-law equation of state, incorporating viscous heating, shock heating, and radiative cooling. We simulate the binary from large separation until after merger, allowing us to model the decoupling of the binary from the CBD. We compute the EM signatures and identify distinct features before, during, and after the merger. Our main result is a multiband EM signature: we find that the MBHB produces strong thermal X-ray emission until 1–2 d prior to the merger. However, as the binary decouples from the CBD, the X-ray-bright minidiscs rapidly shrink in size, become disrupted, and the accretion rate drops precipitously. As a result, the thermal X-ray luminosity drops by orders of magnitude, and the source remains X-ray dark for several days, regardless of any post-merger effects such as GW recoil or mass-loss. Looking for the abrupt spectral change where the thermal X-ray disappears is a tell-tale EM signature of LISA mergers that does not require extensive pre-merger monitoring.

FAUST

Context. Deuterium in H-bearing species is enhanced during the early stages of star formation. However, only a small number of high-spatial-resolution deuteration studies exist towards protostellar objects, leaving the small-scale structures of these objects unrevealed and understudied. Aims. We aim to constrain the deuterium fractionation ratios in a Class 0/I protostellar object in formaldehyde (H2CO), which has abundant deuterated isotopologues in this environment. Methods. We used the Atacama Large Millimeter Array (ALMA) within the context of the Large Program Fifty AU STudy of the chemistry in the disk/envelope system of Solar-like protostars (FAUST) to observe the Class 0/I protobinary system [BHB2007] 11, whose emission components are embedded in circumstellar disks that have radii of 2 to 3 au. The system is surrounded by a complex filamentary structure (the so-called streamers) connected to the larger circumbinary disk. In this work, we present the first study of formaldehyde D-fractionation towards this source with detections of H2CO 3(0,3)–2(0,2), combined with HDCO 4(2,2)–3(2,1), HDCO 4(1,4)–3(1,3) and D2CO 4(0,4)–3(0,3). These observations probe the structures of the protobinary system, enabling us to resolve multiple velocity components associated with the methanol hot spots also uncovered by FAUST data, as well as the colder external envelope. In addition, based on the kinematics seen in our observations of the H2CO emission, we propose the presence of a second large-scale outflow. Results. The results derived from our ALMA observations agree with the current literature in that we only find the deuterated species HDCO and D2CO in the central regions of the core, while undeuterated H2CO is found more ubiquitously. From our radiative transfer modelling, we the column density of H2CO to be in the range of (3-8) × 1014 cm−2 and that of HDCO to be within (0.8−2.9) × 1013 cm−2. The column density for the single detected velocity component of D2CO is within (2.6–1.3) × 1012 cm−2. This yields an average D/H ratio for formaldehyde in [BHB2007] 11 of $0.02_{ - 0.01}^{ + 0.02}$ from HDCO. The results of our kinematic model suggest that the dynamic feature is inconsistent with a streamer-like nature given the flat and outflowing velocity relation; we therefore tentatively conclude that the feature is an asymmetric molecular outflow launched by a wide-angle disk wind.

Three-temperature radiation hydrodynamics with PLUTO

In circumstellar disks around T Tauri stars, visible and near-infrared stellar irradiation is intercepted by dust at the disk’s optical surface and reprocessed into thermal infrared. It subsequently undergoes radiative diffusion through the optically thick bulk of the disk. The gas component, overwhelmingly dominated by mass but contributing little to the opacity, is heated primarily by gas-grain collisions. However, in hydrodynamical simulations, typical models for this heating process (local isothermality, β-cooling, and two-temperature radiation hydrodynamics) incorporate simplifying assumptions that limit their ranges of validity. To build on these methods, we developed a “three-temperature” numerical scheme, which self-consistently models energy exchange between gas, dust, and radiation, as a part of the PLUTO radiation-hydrodynamics code. With a range of test problems in 0D, 1D, 2D, and 3D, we demonstrate the efficacy of our method and make the case for its applicability across a wide range of problems in disk physics, including hydrodynamic instabilities and disk-planet interactions.

A resolved rotating disk wind from a young T Tauri star in the Bok globule CB 26

Context. The disk-outflow connection plays a key role in extracting excess angular momentum from a forming protostar. Although indications of jet rotation have been reported for a few objects, observational constraints of outflow rotation are still very scarce. We have previously reported the discovery of a small collimated molecular outflow from the edge-on T Tauri star-disk system in the Bok globule CB 26 that shows a peculiar velocity pattern, reminiscent of an outflow that corotates with the Keplerian disk. However, we could not ultimately exclude possible alternative explanations for the origin of the observed velocity field. Aims. We report new, high angular resolution millimeter-interferometric observations of CB 26 with the aim of revealing the morphology and kinematics of the outflow at the disk-outflow interface to unambiguously discriminate between the possible alternative explanations for the observed peculiar velocity pattern. Methods. The IRAM PdBI array and the 30 m telescope were used to observe HCO+(1–0) and H13CO+(1–0) at 3.3 mm and 12CO(2–1) at 1.3 mm in three configurations plus zerospacing, resulting in spectral line maps with angular resolutions of 3.″5 and 0.″5, respectively. The SMA was used to observe the HCO+(3–2) line at 1.1 mm with an angular resolution of 1.″35. Additional earlier observations of 13CO(1-0) at 2.7mm with an angular resolution of 1.″0, obtained with OVRO, are also used for the analysis. Using a physical model of the disk, which was derived from the dust continuum emission, we employed chemo-dynamical modeling combined with line radiative transfer calculations to constrain kinematic parameters of the system and to construct a model of the CO emission from the disk that allowed us to separate the emission of the disk from that of the outflow. Results. Our observations confirm the disk-wind nature of the rotating molecular outflow from CB 26 - YSO 1. The new high-resolution data reveal an X-shaped morphology of the CO emission close to the disk, and vertical streaks extending from the disk surface with a small half-opening angle of ≈7°, which can be traced out to vertical heights of ≈500 au. We interpret this emission as the combination of the disk atmosphere and a well-collimated disk wind, of which we mainly see the outer walls of the outflow cone. The decomposition of this emission into a contribution from the disk atmosphere and the disk wind allowed us to trace the disk wind down to vertical heights of ≈40 au, where it is launched from the surface of the flared disk at radii of RL ≈ 20–45 au. The disk wind is rotating with the same orientation and speed as the Keplerian disk and the velocity structure of the cone walls along the flow is consistent with angular momentum conservation. The observed CO outflow has a total gas mass of ≈ 10−3 M⊙, a dynamical age of τdyn ≈ 740 yr, and a total momentum flux of ṖCO ≈ 1.0 × 10−5 M⊙ km s−1 yr−1, which is nearly three orders of magnitude larger than the maximum thrust that can be provided by the luminosity of the central star. Conclusions. We conclude that photoevaporation cannot be the main driving mechanism for this outflow, but it must be predominantly a magnetohydrodynamic disk wind. It is thus far the best-resolved rotating disk wind observed to be launched from a circumstellar disk in Keplerian rotation around a low-mass young stellar object (YSO), albeit also the one with the largest launch radius. It confirms the observed trend that disk winds from Class I YSOs with transitional disks have much larger launch radii than jets ejected from Class 0 protostars.

Workshop Summary: Exoplanet Orbits and Dynamics

Abstract Exoplanetary systems show a wide variety of architectures, which can be explained by different formation and dynamical evolution processes. Precise orbital monitoring is mandatory to accurately constrain their orbital and dynamical parameters. Although major observational and theoretical advances have been made in understanding the architecture and dynamical properties of exoplanetary systems, many outstanding questions remain. This paper aims to give a brief review of a few current challenges in orbital and dynamical studies of exoplanetary systems and a few future prospects for improving our knowledge. Joint data analyses from several techniques are providing precise measurements of orbits and masses for a growing sample of exoplanetary systems, both with close-in orbits and with wide orbits, as well as different evolutionary stages. The sample of young planets detected around stars with circumstellar disks is also growing, allowing for simultaneous studies of planets and their birthplace environments. These analyses will expand with ongoing and future facilities from both ground and space, allowing for detailed tests of formation, evolution, and atmospheric models of exoplanets. Moreover, these detailed analyses may offer the possibility of finding missing components of exoplanetary systems, such as exomoons, or even finding new exotic configurations such as co-orbital planets. In addition to unveiling the architecture of planetary systems, precise measurements of orbital parameters and stellar properties—in combination with more realistic models for tidal interactions and the integration of such models in N-body codes—will improve the inference of the past history of mature exoplanetary systems in close-in orbits. These improvements will allow a better understanding of planetary formation and evolution, placing the solar system in context.

Long-term Evolution of Massive-star Post-common-envelope Circumbinary Disks and the Environments of Fast Luminous Transients

If the envelope of a massive star is not entirely removed during common envelope (CE) interaction with an orbiting compact (e.g., black hole (BH) or neutron star (NS)) companion, the residual bound material eventually cools, forming a centrifugally supported disk around the binary containing the stripped He core. We present a time-dependent height-integrated model for the long-term evolution of post-CE circumbinary disks (CBDs), accounting for mass and angular momentum exchange with the binary, irradiation heating by the He core, and photoevaporation wind mass loss. A large fraction of the CBD’s mass is accreted prior to its outwards viscous spreading and wind dispersal on a timescale of ∼104–105 yr, driving significant orbital migration, even for disks containing ∼10% of the original envelope mass. Insofar that the CBD lifetime is comparable to the thermal (and, potentially, nuclear) timescale of the He core, over which a second mass-transfer episode onto the companion can occur, the presence of the CBD could impact the stability of this key phase. Disruption of the core by the BH/NS would result in a jetted energetic explosion into the dense gaseous CBD (≲1015 cm) and its wind (≳1016 cm), consistent with the environments of luminous fast blue optical transients like AT2018cow. Evolved He cores that undergo core collapse still embedded in their CBD could generate Type Ibn/Icn supernovae. Thousands of dusty wind-shrouded massive-star CBDs may be detectable as extragalactic luminous infrared sources with the Roman Space Telescope; synchrotron radio nebulae powered by the CBD-fed BH/NS may accompany these systems.

Modeling of the High-velocity Jet Powered by the Massive Star MWC 349A

MWC 349A is a massive star with a well-known circumstellar disk rotating following a Keplerian law, and an ionized wind launched from the disk surface. Recent observations with the Atacama Large Millimeter/submillimeter Array (ALMA) carried out toward this system, however, have revealed an additional high-velocity component in the strong, maser emission of hydrogen radio recombination lines (RRLs), suggesting the presence of a high-velocity ionized jet. In this work, we present 3D non-LTE radiative transfer modeling of the emission of the H30α and H26α maser lines, and of their associated radio continuum emission, toward MWC 349A. By using the MORELI code, we reproduce the spatial distribution and kinematics of the high-velocity emission of the H30α and H26α maser lines with a high-velocity ionized jet expanding at a velocity of ∼250 km s−1, surrounded by MWC 349A’s wide-angle ionized wind. The bipolar jet, which is launched from MWC 349A’s disk, is poorly collimated and slightly misaligned with respect to the disk rotation axis. Thanks to the unprecedented sensitivity and spatial accuracy provided by ALMA, we also find that the already known, wide-angle ionized wind decelerates as it expands radially from the ionized disk. We briefly discuss the implications of our findings for understanding the formation and evolution of massive stars. Our results show the huge potential of RRL masers as powerful probes of the innermost ionized regions around massive stars and of their high-velocity jets.

Hydrodynamical simulations of circumbinary accretion: balance between heating and cooling

ABSTRACT Hydrodynamical interaction in circumbinary discs (CBDs) plays a crucial role in various astrophysical systems, ranging from young stellar binaries to supermassive black hole binaries in galactic centres. Most previous simulations of binary-disc systems have adopted locally isothermal equation of state. In this study, we use the grid-based code Athena++ to conduct a suite of two-dimensional viscous hydrodynamical simulations of circumbinary accretion on a Cartesian grid, resolving the central cavity of the binary. The gas thermodynamics is treated by thermal relaxation towards an equilibrium temperature (based on the constant − β cooling ansatz, where β is the cooling time in units of the local Keplerian time). Focusing on equal mass, circular binaries in CBDs with (equilibrium) disc aspect ratio H/R = 0.1, we find that the cooling of the disc gas significantly influences the binary orbital evolution, accretion variability, and CBD morphology, and the effect depends sensitively on the disc viscosity prescriptions. When adopting a constant kinematic viscosity, a finite cooling time (β ≳ 0.1) leads to a binary inspiral as opposed to an outspiral and the CBD cavity becomes more symmetric. When adopting a dynamically varying α-viscosity, binary inspiral only occurs within a narrow range of cooling time (corresponding to β around 0.5).

The dynamical evolution of protoplanetary discs and planets in dense star clusters

ABSTRACT 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 (rp/a ≤ 100, where rp is the periastron distance of the encountering star and a is the semimajor axis of the planet). We also investigate distant perturbers (rp/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.

Bayesian sampling with BeAtlas, a grid of synthetic Be star spectra I. Recovering the fundamental parameters of α Eri and β CMi

ABSTRACT Classical B emission (Be) stars are fast rotating, near-main-sequence B-type stars. The rotation and the presence of circumstellar discs profoundly modify the observables of active Be stars. Our goal is to infer stellar and disc parameters, as well as distance and interstellar extinction, using the currently most favoured physical models for these objects. We present BeAtlas, a grid of $61\, 600$ non-local thermodynamic equilibrium radiative transfer models for Be stars, calculated with the hdust code. The grid was coupled with a Monte Carlo Markov chain (MCMC) code to sample the posterior distribution. We test our method on two well-studied Be stars, α Eri and β CMi, using photometric, polarimetric, and spectroscopic data as input to the code. We recover literature determinations for most of the parameters of the targets, in particular the mass and age of α Eri, the disc parameters of β CMi, and their distances and inclinations. The main discrepancy is that we estimate lower rotational rates than previous works. We confirm previously detected signs of disc truncation in β CMi and note that its inner disc seems to have a flatter density slope than its outer disc. The correlations between the parameters are complex, further indicating that exploring the entire parameter space simultaneously is a more robust approach, statistically. The combination of BeAtlas and Bayesian-MCMC techniques proves successful, and a powerful new tool for the field: The fundamental parameters of any Be star can now be estimated in a matter of hours or days.

Correlation between the optical veiling and accretion properties

Context. Classical T Tauri stars (cTTs) accrete from their circumstellar disk. The material falls onto the stellar surface, producing an accretion shock, which generates veiling in a star’s spectra. In addition, the shock causes a localized accretion spot at the level of the chromosphere. Aims. Our goal is to investigate the accretion, particularly the mass accretion rates (Ṁacc), for the cTTs DK Tau, over two periods of 17 and 29 days, using two different procedures for comparison purposes. Methods. The first method relies on the derivation of the accretion luminosity via accretion-powered emission lines. The second compares the variability of the optical veiling with accretion shock models to determine mass accretion rates. We used observations taken in 2010 and 2012 with the ESPaDOnS (Echelle SpectroPolarimetric Device for the Observation of Stars) spectropolarimeter at the CFHT (Canada-France-Hawaii Telescope). Results. We find peak values of the veiling (at ∼550 nm) ranging from 0.2 to 1.3, with a steeper trend across the wavelength range for higher peak values. When using the accretion-powered emission lines, we find mass accretion rate values ranging from log (Ṁacc[M⊙ yr−1]) = −8.20 to log (Ṁacc[M⊙ yr−1]) = −7.40. This agrees with the values found in the literature, as well as the values calculated using the accretion shock models and the veiling. In addition, we identify a power-law correlation between the values of the accretion luminosity and the optical veiling. For the 2010 observations, using the values of the filling factors (which represent the area of the star covered by an accretion spot) derived from the shock models, we infer that the accretion spot was located between +45° and +75° in latitude. Conclusion. We show that both methods of determining the mass accretion rate yield similar results. We also present a helpful means of confirming the accretion luminosity values by measuring the veiling at a single wavelength in the optical.

Kinematics and stability of high-mass protostellar disk candidates at sub-arcsecond resolution

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 > 104 L⊙) 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.

The potential of VLTI observations for the study of circumstellar disk variability

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.

Early Planet Formation in Embedded Disks (eDisk). VI. Kinematic Structures around the Very-low-mass Protostar IRAS 16253-2429

Precise estimates of protostellar masses are crucial to characterize the formation of stars of low masses down to brown dwarfs (BDs; M * < 0.08 M ☉). The most accurate estimation of protostellar mass uses the Keplerian rotation in the circumstellar disk around the protostar. To apply the Keplerian rotation method to a protostar at the low-mass end, we have observed the Class 0 protostar IRAS 16253-2429 using the Atacama Large Millimeter/submillimeter Array (ALMA) in the 1.3 mm continuum at an angular resolution of 0.″07 (10 au), and in the 12CO, C18O, 13CO (J = 2–1), and SO (J N = 65−54) molecular lines, as part of the ALMA Large Program Early Planet Formation in Embedded Disks project. The continuum emission traces a nonaxisymmetric, disk-like structure perpendicular to the associated 12CO outflow. The position–velocity (PV) diagrams in the C18O and 13CO lines can be interpreted as infalling and rotating motions. In contrast, the PV diagram along the major axis of the disk-like structure in the 12CO line allows us to identify Keplerian rotation. The central stellar mass and the disk radius are estimated to be ∼0.12–0.17 M ☉ and ∼13–19 au, respectively. The SO line suggests the existence of an accretion shock at a ring (r ∼ 28 au) surrounding the disk and a streamer from the eastern side of the envelope. IRAS 16253-2429 is not a proto-BD but has a central stellar mass close to the BD mass regime, and our results provide a typical picture of such very-low-mass protostars.

Post-outburst Evolution of Bona Fide FU Ori-type V2493 Cygnus: A Spectro-photometric Monitoring

We present here the results of 8 yr of our near-simultaneous optical–near-infrared spectro-photometric monitoring of the bona fide FU Ori (FUor)-type candidate V2493 Cygnus (V2493 Cyg) starting from 2013 September to 2021 June. During our optical monitoring period (between 2015 October 16 and 2019 December 30), the V2493 Cyg is slowly dimming with an average dimming rate of ∼26.6 ± 5.6 mmag yr−1 in the V band. Our optical photometric colors show a significant reddening of the source post the second outburst pointing toward a gradual expansion of the emitting region post the second outburst. The mid-infrared colors, on the contrary, exhibit a blueing trend, which can be attributed to the brightening of the disk due to the outburst. Our spectroscopic monitoring shows a dramatic variation of the Hα line as it transitioned from absorption feature to the emission feature and back. Such transition can possibly be explained by the variation in the wind structure in combination with accretion. Combining our time evolution spectra of the Ca ii infrared triplet lines with the previously published spectra of V2493 Cyg, we find that the accretion region has stabilized compared to the early days of the outburst. The evolution of the O i λ7773 Å line also points toward the stabilization of the circumstellar disk post the second outburst.