Accretion Variability Research Articles

PENELLOPE

The processes regulating protoplanetary disk evolution are constrained by studying how mass accretion rates scale with stellar and disk properties. The spread in these relations can be used as a constraint to the models of disk evolution, but only if the impact of accretion variability is correctly accounted for. While the effect of variability might be substantial in the embedded phases of star formation, it is often considered limited at later stages. Here we report on the observed large variation in the accretion rate for one target, XX Cha, and we discuss the impact on population studies of classical T Tauri stars. The mass accretion rate determined by fitting the UV-to-near-infrared spectrum in recent X-shooter observations is compared with the one measured with the same instrument 11 years before. XX Cha displays an accretion variability of almost 2 dex between 2010 and 2021. Although the timescales on which this variability happens are uncertain, XX Cha displays an extreme accretion variability for a classical T Tauri star. If such behavior is common among classical T Tauri stars, possibly on longer timescales than previously probed, it could be relevant for discussing the disk evolution models constrained by the observed spread in accretion rates. Finally, we remark that previous studies of accretion variability based on spectral lines may have underestimated the variability of some targets.

Understanding Accretion Variability through TESS Observations of Taurus

Interpreting the short-timescale variability of the accreting, young, low-mass stars known as Classical T Tauri stars remains an open task. Month-long, continuous light curves from the Transiting Exoplanet Survey Satellite (TESS) have become available for hundreds of T Tauri stars. With this vast data set, identifying connections between the variability observed by TESS and short-timescale accretion variability is valuable for characterizing the accretion process. To this end, we obtained short-cadence TESS observations of 14 T Tauri stars in the Taurus star formation region along with simultaneous ground-based, UBVRI-band photometry to be used as accretion diagnostics. In addition, we combine our data set with previously published simultaneous near-UV–near-IR Hubble Space Telescope spectra for one member of the sample. We find evidence that much of the short-timescale variability observed in the TESS light curves can be attributed to changes in the accretion rate, but note significant scatter between separate nights and objects. We identify hints of time lags within our data set that increase at shorter wavelengths, which we suggest may be evidence of longitudinal density stratification of the accretion column. Our results highlight that contemporaneous, multiwavelength observations remain critical for providing context for the observed variability of these stars.

Towards an understanding of YSO variability: a multiwavelength analysis of bursting, dipping, and symmetrically varying light curves of disc-bearing YSOs

ABSTRACT Using simultaneous optical and infrared light curves of disc-bearing young stars in NGC 2264, we perform the first multiwavelength structure function study of YSOs. We find that dippers have larger variability amplitudes than bursters and symmetric variables at all time-scales longer than a few hours. By analysing optical-infrared colour time-series, we also find that the variability in the bursters is systematically less chromatic at all time-scales than the other variability types. We propose a model of YSO variability in which symmetric, bursting, and dipping behaviour is observed in systems viewed at low, intermediate, and high inclinations, respectively. We argue that the relatively short thermal time-scale for the disc can explain the fact that the infrared light curves for bursters are more symmetric than their optical counterparts, as the disc reprocesses the light from all rotational phases. From this model, we find that the accretion variability on to these YSOs roughly follows a random-walk.

Accretion variability from minute to decade timescales in the classical T Tauri star CR Cha

Context. Classical T Tauri stars are pre-main-sequence stars that are surrounded by a circumstellar disk from which they accrete material. The mass accretion process is essential in the formation of Sun-like stars. Although often described with simple and static models, the accretion process is inherently time variable. Aims. We examine the accretion process of the low-mass young stellar object CR Cha on a wide range of timescales from minutes to a decade by analyzing both photometric and spectroscopic observations from 2006, 2018, and 2019. Methods. We carried out a period analysis of the light curves of CR Cha from the TESS mission and the ASAS-SN and the ASAS-3 databases. We studied the color variations of the system using I, J, H, K-band photometry obtained contemporaneously with the TESS observing window. We analyzed the amplitude, timescale, and the morphology of the accretion tracers found in a series of high-resolution spectra obtained in 2006 with the AAT/UCLES, in 2018 with the HARPS, and in 2019 with the ESPRESSO and the FEROS spectrographs. Results. All photometric data reveal periodic variations compatible with a 2.327-day rotational period. In addition, the ASAS-SN and ASAS-3 data indicate a long-term brightening by 0.2 mag between 2001 and 2008, and a slightly lower brightening than 0.1 mag in the 2015–2018 period. The near-infrared photometry indicates a short-term brightening trend during the observations in 2019. The corresponding color variations can be explained either by a changing accretion rate or changes in the inner disk structure. The Hα line profile variability suggests that the amplitude variations of the central peak, likely due to accretion, are most significant on daily or hourly timescales. On yearly timescales, the line morphology also changes significantly. Conclusions. The photometric variability shows that the period of about 2.3 days is stable in the system over decades. Our results show that the amplitude of the variations in the Hα emission increases on timescales from hours to days or weeks, after which it remains similar even at decadal timescales. On the other hand, we found significant morphological variations on yearly or decadal timescales, indicating that the different physical mechanisms responsible for the line profile changes, such as accretion or wind, are present to varying degrees at different times.

Efficiently Imaging Accreting Protoplanets from Space: Reference Star Differential Imaging of the PDS 70 Planetary System Using the HST/WFC3 Archival PSF Library

Abstract Accreting protoplanets provide key insights into how planets assemble from their natal protoplanetary disks. Recently, Zhou et al. (2021) used angular differential imaging (ADI) with Hubble Space Telescope’s Wide Field Camera 3 (HST/WFC3) to recover the young accreting planet PDS 70 b in F656N (Hα) at a signal-to-noise ratio (S/N) of 7.9. In this paper, we demonstrate a promising approach to efficiently imaging accreting planets by applying reference star differential imaging (RDI) to the same data set. We compile a reference library from the database of WFC3 point-spread functions (PSFs) provided by the Space Telescope Science Institute and develop a set of morphology-significance criteria for preselection of reference frames to improve RDI subtraction. RDI with this PSF library results in a detection of PDS 70 b at an S/N of 5.3. Astrometry and photometry of PDS 70 b are calibrated using a forward-modeling method and injection-recovery tests, resulting in a separation of 186 ± 13 mas, a position angle of 142° ± 5°, and a Hα flux of (1.7 ± 0.3) × 10−15 erg s−1 cm−2. The lower detection significance with RDI can be attributed to the ∼100 times lower peak-to-background ratios of the reference PSFs compared to the ADI PSFs. Building a high-quality reference library with WFC3 will provide unique opportunities to study accretion variability on short timescales not limited by roll angle scheduling constraints and efficiently search for actively accreting protoplanets in Hα around targets inaccessible to ground-based adaptive optics systems, such as faint transition disk hosts.

Preventing Anomalous Torques in Circumbinary Accretion Simulations

Abstract Numerical experiments are the primary method of studying the evolution of circumbinary disks due to the strong nonlinearities involved. Many circumbinary simulations also require the use of numerical mass sinks: source terms that prevent gas from unphysically accumulating around the simulated point masses by removing gas at a given rate. However, special care must be taken when drawing physical conclusions from such simulations to ensure that results are not biased by numerical artifacts. We demonstrate how improved sink methods reduce some of these biases, using simulations of aspect ratio 0.1 accretion disks around binaries with mass ratios between 0.1 and 1. We show that sink terms that do not reduce the angular momentum of gas relative to the accreting object (1) reduce the dependence on the sink rate of quantities such as the torque on the binary, distribution of accretion between binary components, and evolution of the binary semimajor axis; (2) reduce the degree to which the sink rate affects the structure of the accretion disks around each binary component; (3) alter the variability of accretion onto the binary, regularizing it in time. We also investigate other potential sources of systematic error, such as gravitational softening and simplifications to the viscous stress tensor. Because of the strong dependence of binary orbital evolution on both the torque and the distribution of mass between binary components, the sink method influences the orbital evolution of the binary at all mass ratios, with mass ratios below ∼0.3 most affected.

Stellar Transits across a Magnetized Accretion Torus as a Mechanism for Plasmoid Ejection

Abstract The close neighborhood of a supermassive black hole contains not only the accreting gas and dust but also stellar-sized objects, such as late-type and early-type stars and compact remnants that belong to the nuclear star cluster. When passing through the accretion flow, these objects perturb it by the direct action of stellar winds, as well as their magnetic and gravitational effects. By performing general-relativistic magnetohydrodynamic simulations, we investigate how the passages of a star can influence the supermassive black hole gaseous environment. We focus on the changes in the accretion rate and the emergence of blobs of plasma in the funnel of an accretion torus. We compare results from 2D and 3D numerical computations that have been started with comparable initial conditions. We find that a quasi-stationary inflow can be temporarily inhibited by a transiting star, and the plasmoids can be ejected along the magnetic field lines near the rotation axis. We observe the characteristic signatures of the perturbing motion in the power spectrum of the accretion variability, which provides an avenue for a multi-messenger detection of these transient events. Finally, we discuss the connection of our results to multiwavelength observations of galactic nuclei, with the emphasis on ten promising sources (Sgr A*, OJ 287, J0849+5108, RE J1034+396, 1ES 1927+65, ESO 253–G003, GSN 069, RX J1301.9+2747, eRO-QPE1, and eRO-QPE2).

Spiral structures in gravito-turbulent gaseous disks

Context.Gravitational instabilities can drive small-scale turbulence and large-scale spiral arms in massive gaseous disks under conditions of slow radiative cooling. These motions affect the observed disk morphology, its mass accretion rate and variability, and could control the process of planet formation via dust grain concentration, processing, and collisional fragmentation.Aims.We study gravito-turbulence and its associated spiral structure in thin gaseous disks subject to a prescribed cooling law. We characterize the morphology, coherence, and propagation of the spirals and examine when the flow deviates from viscous disk models.Methods.We used the finite-volume code PLUTOto integrate the equations of self-gravitating hydrodynamics in three-dimensional spherical geometry. The gas was cooled over longer-than-orbital timescales to trigger the gravitational instability and sustain turbulence. We ran models for various disk masses and cooling rates.Results.In all cases considered, the turbulent gravitational stress transports angular momentum outward at a rate compatible with viscous disk theory. The dissipation of orbital energy happens via shocks in spiral density wakes, heating the disk back to a marginally stable thermal equilibrium. These wakes drive vertical motions and contribute to mix material from the disk with its corona. They are formed and destroyed intermittently, and they nearly corotate with the gas at every radius. As a consequence, large-scale spiral arms exhibit no long-term global coherence, and energy thermalization is an essentially local process.Conclusions.In the absence of radial substructures or tidal forcing, and provided a local cooling law, gravito-turbulence reduces to a local phenomenon in thin gaseous disks.

Probing accretion variability in CR Cha with TESS and high-resolution spectroscopy

Probing accretion variability in CR Cha with TESS and high-resolution spectroscopy

Synthetic Light Curves of Accretion Variability in T Tauri Stars

Abstract Photometric observations of accreting, low-mass, pre-main-sequence stars (i.e., Classical T Tauri stars; CTTS) have revealed different categories of variability. Several of these classifications have been linked to changes in . To test how accretion variability conditions lead to different light-curve morphologies, we used 1D hydrodynamic simulations of accretion along a magnetic field line coupled with radiative transfer models and a simple treatment of rotation to generate synthetic light curves. We adopted previously developed metrics in order to classify observations to facilitate comparisons between observations and our models. We found that stellar mass, magnetic field geometry, corotation radius, inclination, and turbulence all play roles in producing the observed light curves and that no single parameter is entirely dominant in controlling the observed variability. While the periodic behavior of the light curve is most strongly affected by the inclination, it is also a function of the magnetic field geometry and inner disk turbulence. Objects with either pure dipole fields, strong aligned octupole components, or high turbulence in the inner disk all tend to display accretion bursts. Objects with anti-aligned octupole components or aligned, weaker octupole components tend to show light curves with slightly fewer bursts. We did not find clear monotonic trends between the stellar mass and empirical classification. This work establishes the groundwork for more detailed characterization of well-studied targets as more light curves of CTTS become available through missions such as the Transiting Exoplanet Survey Satellite (TESS).

Investigating episodic accretion in a very low-mass young stellar object

Context. Very low-mass Class I protostars have been investigated very little thus far. Variability of these young stellar objects (YSOs) and whether or not they are capable of strong episodic accretion is also left relatively unstudied. Aims. We investigate accretion variability in IRS 54 (YLW52), a Class I very low-mass protostar with a mass of M⋆ ∼ 0.1 − 0.2 M⊙. Methods. We obtained spectroscopic and photometric data with VLT/ISAAC and VLT/SINFONI in the near-infrared (J, H, and K bands) across four epochs (2005, 2010, 2013, and 2014). We used accretion-tracing lines (Paβ and Brγ) and outflow-tracing lines (H2 and [Fe II]) to examine physical properties and kinematics of the object. Results. A large increase in luminosity was found between the 2005 and 2013 epochs of more than 1 magnitude in the K band, followed in 2014 by a steep decrease. Consistently, the mass accretion rate (Ṁacc) rose by an order of magnitude from ∼10−8 M⊙ yr−1 to ∼10−7 M⊙ yr−1 between the two early epochs. The visual extinction (AV) has also increased from ∼15 mag in 2005 to ∼24 mag in 2013. This rise in AV in tandem with the increase in Ṁacc is explained by the lifting up of a large amount of dust from the disc of IRS 54, following the augmented accretion and ejection activity in the YSO, which intersects our line of sight due to the almost edge-on geometry of the disc. Because of the strength and timescales involved in this dramatic increase, this event is believed to have been an accretion burst possibly similar to bursts of EXor-type objects. IRS 54 is the lowest mass Class I source observed to have an accretion burst of this type, and therefore potentially one of the lowest mass EXor-type objects known so far.

Discovery of a Candidate Binary Supermassive Black Hole in a Periodic Quasar from Circumbinary Accretion Variability

Abstract Binary supermassive black holes (BSBHs) are expected to be a generic byproduct from hierarchical galaxy formation. The final coalescence of BSBHs is thought to be the loudest gravitational wave (GW) siren, yet no confirmed BSBH is known in the GW-dominated regime. While periodic quasars have been proposed as BSBH candidates, the physical origin of the periodicity has been largely uncertain. Here we report discovery of a periodicity (P=1607±7 days) at 99.95% significance (with a global p-value of ∼10−3 accounting for the look elsewhere effect) in the optical light curves of a redshift 1.53 quasar, SDSS J025214.67−002813.7. Combining archival Sloan Digital Sky Survey data with new, sensitive imaging from the Dark Energy Survey, the total ∼20-yr time baseline spans ∼4.6 cycles of the observed 4.4-yr (restframe 1.7-yr) periodicity. The light curves are best fit by a bursty model predicted by hydrodynamic simulations of circumbinary accretion disks. The periodicity is likely caused by accretion rate modulation by a milli-parsec BSBH emitting GWs, dynamically coupled to the circumbinary accretion disk. A bursty hydrodynamic variability model is statistically preferred over a smooth, sinusoidal model expected from relativistic Doppler boost, a kinematic effect proposed for PG1302−102. Furthermore, the frequency dependence of the variability amplitudes disfavors Doppler boost, lending independent support to the circumbinary accretion variability hypothesis. Given our detection rate of one BSBH candidate from circumbinary accretion variability out of 625 quasars, it suggests that future large, sensitive synoptic surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time may be able to detect hundreds to thousands of candidate BSBHs from circumbinary accretion with direct implications for Laser Interferometer Space Antenna.

Circumbinary Disks: Accretion and Torque as a Function of Mass Ratio and Disk Viscosity

Abstract Using numerical hydrodynamics calculations and a novel method for densely sampling parameter space, we measure the accretion and torque on a binary system from a circumbinary disk. In agreement with some earlier studies, we find that the net torque on the binary is positive for mass ratios close to unity, and that accretion always drives the binary toward equal mass. Accretion variability depends sensitively on the numerical sink prescription, but the torque and relative accretion onto each component do not depend on the sink timescale. Positive torque and highly variable accretion occurs only for mass ratios greater than around 0.05. This means that for mass ratios below 0.05, the binary would migrate inward until the secondary accreted sufficient mass, after which it would execute a U-turn and migrate outward. We explore a range of viscosities, from α = 0.03 to α = 0.15, and find that this outward torque is proportional to the viscous torque, so that torque per unit accreted mass is independent of α. Dependence of accretion and torque on mass ratio is explored in detail, densely sampling mass ratios between 0.01 and unity. For mass ratio q > 0.2, accretion variability is found to exhibit a distinct sawtooth pattern, typically with a five-orbit cycle that provides a smoking gun prediction for variable quasars observed over long periods, as a potential means to confirm the presence of a binary.

Pulsed disc accretion driven by Hot jupiters

ABSTRACT We present 2D hydrodynamical simulations of hot Jupiters orbiting near the inner edge of protoplanetary discs. We systemically explore how the accretion rate at the inner disc edge is regulated by a giant planet of different mass, orbital separation, and eccentricity. We find that a massive (with planet-to-star mass ratio ≳0.003) eccentric (ep ≳ 0.1) planet drives a pulsed accretion at the inner edge of the disc, modulated at one or two times the planet’s orbital frequency. The amplitude of accretion variability generally increases with the planet mass and eccentricity, although some non-monotonic dependences are also possible. Applying our simulation results to the T Tauri system CI Tau, where a young hot Jupiter candidate has been detected, we show that the observed luminosity variability in this system can be explained by pulsed accretion driven by an eccentric giant planet.

Analisis Pengaruh Pertumbuhan Ekonomi, IPM, Pengangguran Terbuka dan Angkatan Kerja Terhadap Kemiskinan di Sumatera Utara

<p>Poverty complexity problem is a multidimensional problem, and poverty is related to various aspects of people's lives so that efforts to solve the poverty problem are not easy. This study has purpose to determine and analyze Economic accretion effect, Human development index, Open Unemployment Rate, and Labor Force Participation Rate on Poverty in North Sumatra Province. Secondary data from 2004 - 2017 is used. Data Source is from Board of Central Statistics in North Sumatra Province with testing conducted through the classic judgement exam and statistical exam. With assist of EViews 6.0 statistical data data processing software, data analysis showed that Economic Accretion and idleness Rate variables have a positive and prominent influence on α = 10%, the Human development index variable and the Labor Force Participation Rate are not prominent at α = 10%, on poverty. Regression results are R<sup>2</sup> = 0.945055, that shows independent variable used affects dependent variable by 94.51% and the rest 5.49% is affected by other variables beyond analyzed approach.</p>

Circumbinary Accretion from Finite and Infinite Disks

Abstract We carry out 2D viscous hydrodynamics simulations of circumbinary disk (CBD) accretion using AREPO. We resolve the accretion flow from a large-scale CBD down to the streamers and disks around individual binary components. Extending our recent studies, we consider circular binaries with various mass ratios (0.1 ≤ q b ≤ 1) and study accretion from “infinite,” steady-supply disks and from finite-sized, viscously spreading tori. For “infinite” disks, a global steady state can be reached, and the accretion variability has a dominant frequency ∼0.2Ωb for q b > 0.5 and Ωb for q b < 0.5 (Ωb is the binary angular frequency). We find that the accretion “eigenvalue” l 0—the net angular momentum transfer from the disk to the binary per unit accreted mass—is always positive and falls in the range (0.65– (with a b the binary separation), depending weakly on the mass ratio and viscosity. This leads to binary expansion when q b ≳ 0.3. Accretion from a finite torus can be separated into two phases: an initial transient phase, corresponding to the filling of the binary cavity, followed by a viscous pseudostationary phase, during which the torus viscously spreads and accretes onto the binary. In the viscous phase, the net torque on the binary per unit accreted mass is close to l 0, the value derived for “infinite” disks. We conclude that binaries of similar mass accreting from CBDs gain angular momentum and expand over long timescales. This result significantly impacts the coalescence of supermassive binary black holes and newly formed binary stars. We offer a word of caution against conclusions drawn from simulations of transient accretion onto empty circumbinary cavities.

Molecular and skeletal fingerprints of scleractinian coral biomineralization: From the sea surface to mesophotic depths

Reef-building corals, the major producers of biogenic calcium carbonate, form skeletons in a plethora of morphological forms. Here we studied skeletal modifications of Stylophora pistillata (clade 4) colonies that adapt to increasing depths with decreasing ambient light. The coral show characteristic transitions from spherical morphologies (shallow depths, 5 m deep) to flat and branching geometries (mesophotic depths, 60 m deep). Such changes are typically ascribed to the algal photosymbiont physiological feedback with the coral that host them. We find specific fine-scale skeletal variability in accretion of structure at shallow- and mesophotic depth morphotypes that suggest underlying genomic regulation of biomineralization pathways of the coral host. To explain this, we conducted comparative morphology-based analyses, including optical and electron microscopy, tomography and X-ray diffraction analysis coupled with a comprehensive transcriptomic analysis of S. pistillata. The samples originated from Gulf of Eilat in the Red Sea collected along a depth gradient from shallow to mesophotic depths (5 to 60 m). Additional samples were experimentally transplanted from 5 m to 60 m and from 60 m to 5 m. Interestingly, both morphologically and functionally, transplanted corals partly adapt by exhibiting typical depth-specific properties. In mesophotic depths, we find that the organic matrix fraction is enriched in the coralla, well matching the overrepresentation of transcripts encoding biomineralization “tool-kit" structural extracellularproteins that was observed. These results provide insights into the molecular mechanisms of calcification and skeletal adaptation that repeatedly allowed this coral group to adapt to a range of environments presumably with a rich geological past. Statement of SignificanceUnderstanding the reef coral physiological plasticity under a rapidly changing climate is of crucial importance for the protection of coral reef ecosystems. Most of the reef corals operate near their upper limit of heat tolerance. A possible rescue for some coral species is migration to deeper, cooler mesophotic depths. However, gradually changing environmental parameters (especially light) along the depth gradient pose new adaptative stress on corals with largely unknown influences on the various biological molecular pathways. This work provides a first comprehensive analysis of changes in gene expression, including biomineralization “tool kit” genes, and reports the fine-scale microstructural and crystallographic skeletal details in S. pistillata collected in the Red Sea along a depth gradient spannign 5 to 60 m.

Silicate features in the circumstellar envelopes of the Class I binary driving source of HH 250

Context. The silicate feature near 10 μm is one of the main tools available to study the mineralogy of circumstellar disks and envelopes, providing information on the thermal processing, growth, location, and circulation of dust grains. Aims. We investigate the silicate feature of the two Class I components of HH 250-IRS, a resolved binary system with a separation of 0″53 driving a Herbig-Haro flow. Each component has its own circumstellar envelope, and the system is surrounded by a circumbinary disk. Methods. We carried out low resolution spectroscopy in the 8–13 μm range using VISIR, the thermal infrared imager and spectrograph at ESO’s Very Large Telescope. Results. The silicate features of both sources are clearly different. The northwest (NW) component has a broad, smooth absorption profile lacking structure. We attribute most of it to foreground interstellar dust absorption, but estimate that additional absorption by amorphous silicates takes place in the circumstellar envelope of the young stellar object. The southeast (SE) component shows the silicate feature in emission, with structure longward of 9.5 μm indicating the presence of crystalline dust in the dominant form of forsterite. The apparent lack of an absorption feature caused by foreground dust is probably due to the filling of the band with emission by amorphous silicates in the envelope of the object. Conclusions. Despite their virtually certain coevality, the differences in the components of the HH 250-IRS binary are most likely due to markedly different circumstellar environments. The NW component displays an unevolved envelope, whereas dust growth and crystallization has taken place in the SE component. The weak or absent signatures of enstatite in the latter are fairly unusual among envelopes with crystalline dust, and we tentatively relate it to a possible wide gap or an inner truncation of the disk already hinted at in previous observations by a drop in the L′-band flux, which might indicate that the SE component could actually be a very close binary. We speculate that the clear differences between the silicate feature spectra of both components of HH 250-IRS may be due either to disk evolution sped up by multiplicity, or by accretion variability leading to episodes of crystal formation. Different inclinations with respect to the line of sight may play a role as well, although it is very unlikely that they are the sole element for the differences between both objects.

On the ALMA observability of nascent massive multiple systems formed by gravitational instability

Massive young stellar object (MYSOs) form during the collapse of high-mass pre-stellar cores, where infalling molecular material is accreted through a centrifugally-balanced accretion disc that is subject to efficient gravitational instabilities. In the resulting fragmented accretion disc of the MYSO, gaseous clumps and low-mass stellar companions can form, which will influence the future evolution of massive protostars in the Hertzsprung-Russell diagram. We perform dust continuum radiative transfer calculations and compute synthetic images of disc structures modelled by the gravito-radiation-hydrodynamics simulation of a forming MYSO, in order to investigate the Atacama Large Millimeter/submillimeter Array (ALMA) observability of circumstellar gaseous clumps and forming multiple systems. Both spiral arms and gaseous clumps located at ~a few 100 au from the protostar can be resolved by interferometric ALMA Cycle 7 C43-8 and C43-10 observations at band 6 (1.2 mm), using a maximal 0.015" beam angular resolution and at least 10-30 min exposure time for sources at distances of 1-2 kpc. Our study shows that substructures are observable regardless of their viewing geometry or can be inferred in the case of an edge-viewed disc. The observation probability of the clumps increases with the gradually increasing efficiency of gravitational instability at work as the disc evolves. As a consequence, large discs around MYSOs close to the zero-age-main-sequence line exhibit more substructures than at the end of the gravitational collapse. Our results motivate further observational campaigns devoted to the close surroundings of the massive protostars S255IR-NIRS3 and NGC 6334I-MM1, whose recent outbursts are a probable signature of disc fragmentation and accretion variability.

Exploring the dimming event of RW Aurigae A through multi-epoch VLT/X-shooter spectroscopy

Context. RW Aur A is a classical T Tauri star that has suddenly undergone three major dimming events since 2010. The reason for these dimming events is still not clear. Aims. We aim to understand the dimming properties, examine accretion variability, and derive the physical properties of the inner disc traced by the CO ro-vibrational emission at near-infrared wavelengths (2.3 μm). Methods. We compared two epochs of X-shooter observations, during and after the dimming. We modelled the rarely detected CO bandhead emission in both epochs to examine whether the inner disc properties had changed. The spectral energy distribution was used to derive the extinction properties of the dimmed spectrum and compare the infrared excess between the two epochs. Lines tracing accretion were used to derive the mass accretion rate in both states. Results. The CO originates from a region with physical properties of T = 3000 K, NCO = 1 × 1021 cm−2 and vk sin i = 113 km s−1. The extinction properties of the dimming layer were derived with the effective optical depth ranging from τeff ~2.5−1.5 from the UV to the near-IR. The inferred mass accretion rate Ṁacc is ~1.5 × 10−8 M⊙ yr−1 and ~2 × 10−8 M⊙ yr−1 after and during the dimming respectively. By fitting the spectral energy distribution, additional emission is observed in the infrared during the dimming event from dust grains with temperatures of 500–700 K. Conclusions. The physical conditions traced by the CO are similar for both epochs, indicating that the inner gaseous disc properties do not change during the dimming events. The extinction curve is flatter than that of the interstellar medium, and large grains of a few hundred microns are thus required. When we correct for the observed extinction, the mass accretion rate is constant in the two epochs, suggesting that the accretion is stable and therefore does not cause the dimming. The additional hot emission in the near-IR is located at about 0.5 au from the star and is not consistent with an occulting body located in the outer regions of the disc. The dimming events could be due to a dust-laden wind, a severe puffing-up of the inner rim, or a perturbation caused by the recent star-disc encounter.