Showing 1–50 of 122 results for author: Vorobyov, E I

Primordial dust rings, hidden dust mass, and the first generation of planetesimals in gravitationally unstable protoplanetary disks

Authors: Eduard I. Vorobyov, Aleksandr M. Skliarevskii, Manuel Guedel, Tamara Molyarova

Abstract: Aims. A new mechanism of dust accumulation and planetesimal formation in a gravitationally unstable disk with suppressed magnetorotational instability is studied and compared with the classical dead zone in a layered disk model. Methods. We use numerical hydrodynamics simulations in the thin-disk limit FEOSAD code to model the formation and long-term evolution of gravitationally unstable disks, in… ▽ More Aims. A new mechanism of dust accumulation and planetesimal formation in a gravitationally unstable disk with suppressed magnetorotational instability is studied and compared with the classical dead zone in a layered disk model. Methods. We use numerical hydrodynamics simulations in the thin-disk limit FEOSAD code to model the formation and long-term evolution of gravitationally unstable disks, including dust dynamics and growth. Results. We found that in gravitationally unstable disks with a radially varying strength of gravitational instability a region of low mass and angular momentum transport forms in the inner several astronomical units. This region is characterized by low effective α_GI and is similar in characteristics to the dead zone in the layered disk model. As the disk forms and evolves, the GI-induced dead zone accumulates a massive dust ring, which is susceptible to the development of the streaming instability. The model and observationally inferred dust masses and radii may differ significantly in gravitationally unstable disks with massive inner dust rings. Conclusions. The early occurrence of the GI-induced dust ring followed by the presumed development of the streaming instability suggest that this mechanism may form the first generation of planetesimals in the inner terrestrial zone of the disk. The proposed mechanism, however, crucially depends on the susceptibility of the disk to gravitational instability and requires that the magnetorotational instability be suppressed. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

Metallicity dependence of dust growth in a protoplanetary disk

Authors: Ryoki Matsukoba, Eduard I. Vorobyov, Takashi Hosokawa

Abstract: In the context of planet formation, growth from micron-sized grains to kilometer-sized planetesimals is a crucial question. Since the dust growth rate depends on the amount of dust, realizing planet formation scenarios based on dust growth is challenging in environments with low metallicity, i.e. less dust. We investigate dust growth during disk evolution, particularly focusing on the relationship… ▽ More In the context of planet formation, growth from micron-sized grains to kilometer-sized planetesimals is a crucial question. Since the dust growth rate depends on the amount of dust, realizing planet formation scenarios based on dust growth is challenging in environments with low metallicity, i.e. less dust. We investigate dust growth during disk evolution, particularly focusing on the relationship with metallicity. We perform two-dimensional thin-disk hydrodynamic simulations to track the disk evolution over 300 kyr from its formation. The dust motion is solved separately from the gas motion, with its distribution changing due to drag forces from the gas. Dust size growth is also accounted for, with the magnitude of the drag force varying according to the dust size. We employ three models with metallicities of 1.0, 0.1, and 0.01 ${\rm Z}_{\odot}$, i.e. dust-to-gas mass ratios of 10$^{-2}$, 10$^{-3}$, and 10$^{-4}$, respectively. In the disks with the metallicities $\ge0.1$ ${\rm Z}_{\odot}$, the dust radii reach cm sizes, consistent with estimations from the dust growth timescale. Conversely, for the metallicity of 0.01 ${\rm Z}_{\odot}$, the maximum dust size is only 10$^{-2}$ cm, with almost no growth observed across the entire disk scale ($\sim$100 au). At the metallicities $\ge0.1$ ${\rm Z}_{\odot}$, the decoupling between grown dust and gas leads to non-uniform dust-to-gas mass ratios. However, deviations from the canonical value of this ratio have no impact on the gravitational instability of the disk. The formation of dust rings is confirmed in the innermost part of the disk ($\sim$10--30 au). The dust rings where the dust-to-gas mass ratio is enhanced, and the Stokes number reaches $\sim$0.1, are suitable environments for the streaming instability. We conjecture that planetesimal formation occurs through the streaming instability in these dust rings. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: 13 pages, 12 figures, submitted to A&A

Luminosity Outbursts in Interacting Protoplanetary Systems

Authors: Aleksandr M. Skliarevskii, Eduard I. Vorobyov

Abstract: FU Orionis type objects (fuors) are characterized by rapid (tens to hundreds years) episodic outbursts, during which the luminosity increases by orders of magnitude. One of the possible causes of such events is a close encounter between stars and protoplanetary disks. Numerical simulations show that the fuor-like outburst ignition requires a very close encounter ranging from a few to a few tens of… ▽ More FU Orionis type objects (fuors) are characterized by rapid (tens to hundreds years) episodic outbursts, during which the luminosity increases by orders of magnitude. One of the possible causes of such events is a close encounter between stars and protoplanetary disks. Numerical simulations show that the fuor-like outburst ignition requires a very close encounter ranging from a few to a few tens of au. In contrast, the observed stellar objects in fuor binaries are usually hundreds of au apart. Simple mathematical estimates show that if such a close approach took place, the binary stellar components would have an unrealistic relative velocity, at least an order of magnitude greater than the observed velocity dispersion in young stellar clusters. Thus, the bursts are either triggered with a certain delay after passage of the periastron or their ignition does not necessary require a close encounter and hence the outburst is not caused by the primordial gravitational perturbation of the protoplanetary disk. In this work, an encounter of a star surrounded by a protoplanetary disk with a diskless external stellar object was modeled using numerical hydrodynamics simulations. We showed that even fly-bys with a relatively large periastron (at least 500 au) can result in fuor-like outbursts. Moreover, the delay between the periastron passage and the burst ignition can reach several kyr. It was shown for the first time by means of numerical modeling that the perturbation of the disk caused by the external object can trigger a cascade process, which includes the development of the thermal instability in the innermost disk followed by the magneto-rotational instability ignition. Because of the sequential development of these instabilities, the rapid increase in the accretion rate occurs, resulting in the luminosity increase by more than two orders of magnitude. △ Less

Submitted 21 February, 2024; originally announced February 2024.

Comments: Accepted for publication by Astronomy Reports

Effect of Dust Evaporation and Thermal Instability on Temperature Distribution in a Protoplanetary Disk

Authors: Ya. N. Pavlyuchenkov, V. V. Akimkin, A. P. Topchieva, E. I. Vorobyov

Abstract: The thermal instability of accretion disks is widely used to explain the activity of cataclysmic variables, but its development in protoplanetary disks has been studied in less detail. We present a semi-analytical stationary model for calculating the midplane temperature of a gas and dust disk around a young star. The model takes into account gas and dust opacities, as well as the evaporation of d… ▽ More The thermal instability of accretion disks is widely used to explain the activity of cataclysmic variables, but its development in protoplanetary disks has been studied in less detail. We present a semi-analytical stationary model for calculating the midplane temperature of a gas and dust disk around a young star. The model takes into account gas and dust opacities, as well as the evaporation of dust at temperatures above 1000 K. Using this model, we calculate the midplane temperature distributions of the disk under various assumptions about the source of opacity and the presence of dust. We show that when all considered processes are taken into account, the heat balance equation in the region r<1 au has multiple temperature solutions. Thus, the conditions for thermal instability are met in this region. To illustrate the possible influence of instability on the accretion state in a protoplanetary disk, we consider a viscous disk model with alpha parameterization of turbulent viscosity. We show that in such a model the disk evolution is non-stationary, with alternating phases of accumulation of matter in the inner disk and its rapid accretion onto the star, leading to an episodic accretion pattern. These results indicate that this instability needs to be taken into account in evolutionary models of protoplanetary disks. △ Less

Submitted 28 July, 2023; originally announced July 2023.

Comments: Published in Astronomy Reports Vol. 67, No. 5, pp. 470-482 (2023)

Formation of a wide-orbit giant planet in a gravitationally unstable subsolar-metallicity protoplanetary disc

Authors: Ryoki Matsukoba, Eduard I. Vorobyov, Takashi Hosokawa, Manuel Guedel

Abstract: Direct imaging observations of planets revealed that wide-orbit ($>10$ au) giant planets exist even around subsolar-metallicity host stars and do not require metal-rich environments for their formation. A possible formation mechanism of wide-orbit giant planets in subsolar-metallicity environments is the gravitational fragmentation of massive protoplanetary discs. Here, we follow the long-term evo… ▽ More Direct imaging observations of planets revealed that wide-orbit ($>10$ au) giant planets exist even around subsolar-metallicity host stars and do not require metal-rich environments for their formation. A possible formation mechanism of wide-orbit giant planets in subsolar-metallicity environments is the gravitational fragmentation of massive protoplanetary discs. Here, we follow the long-term evolution of the disc for 1 Myr after its formation, which is comparable to disc lifetime, by way of a two-dimensional thin-disc hydrodynamic simulation with the metallicity of 0.1 ${\rm Z}_{\odot}$. We find a giant protoplanet that survives until the end of the simulation. The protoplanet is formed by the merger of two gaseous clumps at $\sim$0.5 Myr after disc formation, and then it orbits $\sim$200 au from the host star for $\sim$0.5 Myr. The protoplanet's mass is $\sim$10 ${\rm M}_{\rm J}$ at birth and gradually decreases to 1 ${\rm M}_{\rm J}$ due to the tidal effect from the host star. The result provides the minimum mass of 1 ${\rm M}_{\rm J}$ for protoplanets formed by gravitational instability in a subsolar-metallicity disc. We anticipate that the mass of a protoplanet experiencing reduced mass loss thanks to the protoplanetary contraction in higher resolution simulations can increase to $\sim$10 ${\rm M}_{\rm J}$. We argue that the disc gravitational fragmentation would be a promising pathway to form wide-orbit giant planets with masses of $\ge1$ ${\rm M}_{\rm J}$ in subsolar-metallicity environments. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 11 pages, 8 figures, submitted to MNRAS

What governs the spin distribution of very young < 1 Myr low mass stars

Authors: L. Gehrig, E. I. Vorobyov

Abstract: We compute the evolution and rotational periods of young stars, using the MESA code, starting from a stellar seed, and take protostellar accretion, stellar winds, and the magnetic star-disk interaction into account. Furthermore, we add a certain fraction of the energy of accreted material into the stellar interior as additional heat and combine the resulting effects on stellar evolution with the s… ▽ More We compute the evolution and rotational periods of young stars, using the MESA code, starting from a stellar seed, and take protostellar accretion, stellar winds, and the magnetic star-disk interaction into account. Furthermore, we add a certain fraction of the energy of accreted material into the stellar interior as additional heat and combine the resulting effects on stellar evolution with the stellar spin model. For different combinations of parameters, stellar periods at an age of 1~Myr range between 0.6~days and 12.9~days. Thus, during the relatively short time period of 1~Myr, a significant amount of stellar angular momentum can already be removed by the interaction between the star and its accretion disk. The amount of additional heat added into the stellar interior, the accretion history, and the presence of disk and stellar winds have the strongest impact on the stellar spin evolution during the first million years. The slowest stellar rotations result from a combination of strong magnetic fields, a large amount of additional heat, and effective winds. The fastest rotators combine weak magnetic fields and ineffective winds or result from a small amount of additional heat added to the star. Scenarios that could lead to such configurations are discussed. Different initial rotation periods of the stellar seed, on the other hand, quickly converges and do not affect the stellar period at all. Our model matches up to 90\% of the observed rotation periods in six young ($\lesssim 3$~Myr) clusters. Based on these intriguing results, we motivate to combine our model with a hydrodynamic disk evolution code to self-consistently include several important aspects such as episodic accretion events, magnetic disk winds, internal, and external photo-evaporation. (Shortened) △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: 17 pages, 15 figures, accepted for publication in A&A

Journal ref: A&A 673, A54 (2023)

Predicting Stellar Mass Accretion: An Optimized Echo State Network Approach in Time Series Modeling

Authors: Gianfranco Bino, Shantanu Basu, Ramit Dey, Sayantan Auddy, Lyle Muller, Eduard I. Vorobyov

Abstract: Modeling the dynamics of the formation and evolution of protostellar disks as well as the history of stellar mass accretion typically involve the numerical solution of complex systems of coupled differential equations. The resulting mass accretion history of protostars is known to be highly episodic due to recurrent instabilities and also exhibits short timescale flickering. By leveraging the stro… ▽ More Modeling the dynamics of the formation and evolution of protostellar disks as well as the history of stellar mass accretion typically involve the numerical solution of complex systems of coupled differential equations. The resulting mass accretion history of protostars is known to be highly episodic due to recurrent instabilities and also exhibits short timescale flickering. By leveraging the strong predictive abilities of neural networks, we extract some of the critical temporal dynamics experienced during the mass accretion including periods of instability. Particularly, we utilize a novel form of the Echo-State Neural Network (ESN), which has been shown to efficiently deal with data having inherent nonlinearity. We introduce the use of Optimized-ESN (Opt-ESN) to make model-independent time series forecasting of mass accretion rate in the evolution of protostellar disks. We apply the network to multiple hydrodynamic simulations with different initial conditions and exhibiting a variety of temporal dynamics to demonstrate the predictability of the Opt-ESN model. The model is trained on simulation data of $\sim 1-2$ Myr, and achieves predictions with a low normalized mean square error ($\sim 10^{-5}$ to $10^{-3}$) for forecasts ranging between 100 and 3800 yr. This result shows the promise of the application of machine learning based models to time-domain astronomy. △ Less

Submitted 13 April, 2023; v1 submitted 7 February, 2023; originally announced February 2023.

Comments: 16 pages, 13 figures

Journal ref: Open Journal of Astrophysics, Vol. 6, 2023

Formation of pebbles in (gravito-)viscous protoplanetary disks with various turbulent strengths

Authors: Eduard I. Vorobyov, Vardan G. Elbakyan, Anders Johansen, Michiel Lambrechts, Aleksandr M. Skliarevskii, O. P. Stoyanovskaya

Abstract: Aims. Dust plays a crucial role in the evolution of protoplanetary disks. We study the dynamics and growth of initially sub-$μm$ dust particles in self-gravitating young protoplanetary disks with various strengths of turbulent viscosity. We aim to understand the physical conditions that determine the formation and spatial distribution of pebbles when both disk self-gravity and turbulent viscosity… ▽ More Aims. Dust plays a crucial role in the evolution of protoplanetary disks. We study the dynamics and growth of initially sub-$μm$ dust particles in self-gravitating young protoplanetary disks with various strengths of turbulent viscosity. We aim to understand the physical conditions that determine the formation and spatial distribution of pebbles when both disk self-gravity and turbulent viscosity can be concurrently at work. Methods. We perform the thin-disk hydrodynamics simulations of self-gravitating protoplanetary disks over an initial time period of 0.5 Myr using the FEOSAD code. Turbulent viscosity is parameterized in terms of the spatially and temporally constant $α$-parameter, while the effects of gravitational instability on dust growth is accounted for by calculating the effective parameter $α_{\rm GI}$. We consider the evolution of dust component including momentum exchange with gas, dust self-gravity, and also a simplified model of dust growth. Results. We find that the level of turbulent viscosity strongly affects the spatial distribution and total mass of pebbles in the disk. The $α=10^{-2}$ model is viscosity-dominated, pebbles are completely absent, and dust-to-gas mass ratio deviates from the reference 1:100 value no more than by 30\% throughout the disk extent. On the contrary, the $α=10^{-3}$ model and, especially, the $α=10^{-4}$ model are dominated by gravitational instability. The effective parameter $α+α_{\rm GI}$ is now a strongly varying function of radial distance. As a consequence, a bottle neck effect develops in the innermost disk regions, which makes gas and dust accumulate in a ring-like structure. Abridged. △ Less

Submitted 2 December, 2022; originally announced December 2022.

Comments: Accepted for publication in Astronomy and Astrophysics

The influence of metallicity on a combined stellar and disk evolution

Authors: L. Gehrig, T. Steindl, E. I. Vorobyov, R. Guadarrama, K. Zwintz

Abstract: The effects of an accretion disk are crucial to understanding the evolution of young stars. During the combined evolution, stellar and disk parameters influence each other, motivating a combined stellar and disk model. This makes a combined numerical model, evolving the disk alongside the star, the next logical step in the progress of studying early stellar evolution. We aim to understand the effe… ▽ More The effects of an accretion disk are crucial to understanding the evolution of young stars. During the combined evolution, stellar and disk parameters influence each other, motivating a combined stellar and disk model. This makes a combined numerical model, evolving the disk alongside the star, the next logical step in the progress of studying early stellar evolution. We aim to understand the effects of metallicity on the accretion disk and the stellar spin evolution during the T~Tauri phase. We combine the numerical treatment of a hydrodynamic disk with stellar evolution, including a stellar spin model, allowing a self-consistent calculation of the back-reactions between the individual components. We present the self-consistent theoretical evolution of T-Tauri stars coupled to a stellar disk. We find that disks in low metallicity environments are heated differently and have shorter lifetimes, compared to their solar metallicity counterparts. Differences in stellar radii, the contraction rate of the stellar radius, and the shorter disk lifetimes result in faster rotation of low metallicity stars. We present an additional explanation for the observed short disk lifetimes in low metallicity clusters. A combination of our model with previous studies (e.g., a metallicity-based photo-evaporation) could help to understand disk evolution and dispersal at different metallicities. Furthermore, the stellar spin evolution model includes several important effects, previously ignored (e.g., the stellar magnetic field strength and a realistic calculation of the disk lifetime) and we motivate to include our results as initial or input parameters for further spin evolution models, covering the stellar evolution towards and during the main sequence. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: accepted for publication in A&A

Journal ref: A&A 669, A84 (2023)

Episodic accretion and mergers during growth of massive protostars

Authors: Vardan G. Elbakyan, Sergei Nayakshin, Dominique M. -A. Meyer, Eduard I. Vorobyov

Abstract: 3D simulations of high mass young stellar object (HMYSO) growth show that their circumstellar discs fragment onto multiple self-gravitating objects. Accretion of these by HMYSO may explain episodic accretion bursts discovered recently. We post-process results of a previous 3D simulation of a HMYSO disc with a 1D code that resolves the disc and object dynamics down to the stellar surface. We find t… ▽ More 3D simulations of high mass young stellar object (HMYSO) growth show that their circumstellar discs fragment onto multiple self-gravitating objects. Accretion of these by HMYSO may explain episodic accretion bursts discovered recently. We post-process results of a previous 3D simulation of a HMYSO disc with a 1D code that resolves the disc and object dynamics down to the stellar surface. We find that burst-like deposition of material into the inner disc seen in 3D simulations by itself does not always signify powerful accretion bursts. Only high density post-collapse clumps crossing the inner computational boundary may result in observable bursts. The rich physics of the inner disc has a significant impact on the expected accretion bursts: (1) In the standard turbulent viscosity discs, migrating objects can stall at a migration trap at the distance of a few au from the star. However, in discs powered by magnetised winds, the objects are able to cross the trap and produce bursts akin to those observed so far. (2) Migrating objects may interact with and modify the thermal (hydrogen ionisation) instability of the inner disc, which can be responsible for longer duration and lower luminosity bursts in HMYSOs. (3) If the central star is bloated to a fraction of an au by a previous episode of high accretion rate, or if the migrating object is particularly dense, a merger rather than a disc-mediated accretion burst results; (4) Object disruption bursts may be super-Eddington, leading to episodic feedback on HMYSO surroundings via powerful outflows. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: 19 pages, 16 figures, accepted to MNRAS

The burst mode of accretion in massive star formation with stellar inertia

Authors: D. M. -A. Meyer, E. I. Vorobyov, V. G. Elbakyan, S. Kraus, S. -Y. Liu, S. Nayakshin, A. M. Sobolev

Abstract: The burst mode of accretion in massive star formation is a scenario linking the initial gravitational collapse of parent pre-stellar cores to the properties of their gravitationally unstable discs and of their accretion-driven bursts. In this study, we present a series of high-resolution 3D radiation-hydrodynamics numerical simulations for young massive stars formed out of collapsing 100 Mo molecu… ▽ More The burst mode of accretion in massive star formation is a scenario linking the initial gravitational collapse of parent pre-stellar cores to the properties of their gravitationally unstable discs and of their accretion-driven bursts. In this study, we present a series of high-resolution 3D radiation-hydrodynamics numerical simulations for young massive stars formed out of collapsing 100 Mo molecular cores spinning with several values of the ratio of rotational-to-gravitational energies beta=5%-9%. The models include the indirect gravitational potential caused by disc asymmetries. We find that this modifies the barycenter of the disc, causing significant excursions of the central star position, which we term stellar wobbling. The stellar wobbling slows down and protracts the development of gravitational instability in the disc, reducing the number and magnitude of the accretion-driven bursts undergone by the young massive stars, whose properties are in good agreement with that of the burst monitored from the massive protostar M17 MIR. Including stellar wobbling is therefore important for accurate modeling disc structures. Synthetic ALMA interferometric images in the millimeter waveband show that the outcomes of efficient gravitational instability such as spiral arms and gaseous clumps can be detected for as long as the disc is old enough and has already entered the burst mode of accretion. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: Accepted to MNRAS, see https://academic.oup.com/mnras/advance-article-abstract/doi/10.1093/mnras/stac2956/6761717?utm_source=advanceaccess&utm_campaign=mnras&utm_medium=email

The effect of metallicity on the abundances of molecules in protoplanetary disks

Authors: R. Guadarrama, Eduard I. Vorobyov, Christian Rab, Manuel Güdel

Abstract: We study the influence of different metallicities on the physical, thermal, and chemical properties of protoplanetary disks, and in particular on the formation and destruction of carbon-based molecules. With the thermo-chemical code ProoDIMO we investigate the impact of lower metallicities on the radiation field, disk temperature, and the abundance of different molecules (H$_2$O, CH$_4$, CO, CO… ▽ More We study the influence of different metallicities on the physical, thermal, and chemical properties of protoplanetary disks, and in particular on the formation and destruction of carbon-based molecules. With the thermo-chemical code ProoDIMO we investigate the impact of lower metallicities on the radiation field, disk temperature, and the abundance of different molecules (H$_2$O, CH$_4$, CO, CO$_2$, HCN, CN, HCO$^+$ and N$_2$H$^+$). We use a fiducial disk model as a reference model and produce two models with lower metallicity. The resulting influence on different chemical species is studied by analyzing their abundance distribution throughout the disk and their vertical column density. Furthermore, the formation and destruction reactions of the chemical species are studied. The results show a relation between the metallicity of the disk and the strength of the stellar radiation field inside the disk. As the metallicity decreases the radiation field is able to penetrate deeper regions of the disk. As a result, there is a stronger radiation field overall in the disk with lower metallicity which also heats up the disk. This triggers a series of changes in the chemical formation and destruction efficiencies for different chemical species. In most cases, the available species abundances change and have greater values compared to scaled-down abundances by constant factors. Metallicity has a clear impact on the snowline of the molecules studied here as well. As metallicity decreases the snowlines are pushed further out and existing snow rings shrink in size. △ Less

Submitted 19 August, 2022; originally announced August 2022.

Comments: 19 pages, 18 figures

Journal ref: A&A 667, A28 (2022)

Protostellar-disc fragmentation across all metallicities

Authors: Ryoki Matsukoba, Kei E. I. Tanaka, Kazuyuki Omukai, Eduard I. Vorobyov, Takashi Hosokawa

Abstract: Cosmic metallicity evolution possibly creates the diversity of star formation modes at different epochs. Gravitational fragmentation of circumstellar discs provides an important formation channel of multiple star systems, including close binaries. We here study the nature of disc fragmentation, systematically performing a suite of two-dimensional radiation-hydrodynamic simulations, in a broad rang… ▽ More Cosmic metallicity evolution possibly creates the diversity of star formation modes at different epochs. Gravitational fragmentation of circumstellar discs provides an important formation channel of multiple star systems, including close binaries. We here study the nature of disc fragmentation, systematically performing a suite of two-dimensional radiation-hydrodynamic simulations, in a broad range of metallicities, from the primordial to the solar values. In particular, we follow relatively long-term disc evolution over 15 kyr after the disc formation, incorporating the effect of heating by the protostellar irradiation. Our results show that the disc fragmentation occurs at all metallicities $1$--$0$ $Z_{\odot}$, yielding self-gravitating clumps. Physical properties of the clumps, such as their number and mass distributions, change with the metallicity due to different gas thermal evolution. For instance, the number of clumps is the largest for the intermediate metallicity range of $10^{-2}$--$10^{-5}$ $Z_{\odot}$, where the dust cooling is effective exclusively in a dense part of the disc and causes the fragmentation of spiral arms. The disc fragmentation is more modest for $1$--$0.1$ $Z_{\odot}$ thanks to the disc stabilization by the stellar irradiation. Such metallicity dependence agrees with the observed trend that the close binary fraction increases with decreasing metallicity in the range of $1$--$10^{-3}$ $Z_{\odot}$. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 17 pages, 16 figures, submitted to MNRAS

Millimeter-sized Dust Grains Appear Surviving the Water-sublimating Temperature in the Inner 10 au of the FU Ori Disk

Authors: Hauyu Baobab Liu, An-Li Tsai, Wen Ping Chen, Jin Zhong Liu, Xuan Zhang, Shuo Ma, Vardan Elbakyan, Joel D. Green, Antonio S. Hales, Michihiro Takami, Sebastián Pérez, Eduard I. Vorobyov, Yao-Lun Yang

Abstract: Previous observations have shown that the $\lesssim$10 au, $\gtrsim$400 K hot inner disk of the archetypal accretion outburst young stellar object, FU Ori, is dominated by viscous heating. To constrain dust properties in this region, we have performed radio observations toward this disk using the Karl G. Jansky Very Large Array (JVLA) in 2020 June-July, September, and November. We also performed c… ▽ More Previous observations have shown that the $\lesssim$10 au, $\gtrsim$400 K hot inner disk of the archetypal accretion outburst young stellar object, FU Ori, is dominated by viscous heating. To constrain dust properties in this region, we have performed radio observations toward this disk using the Karl G. Jansky Very Large Array (JVLA) in 2020 June-July, September, and November. We also performed complementary optical photometric monitoring observations. We found that the dust thermal emission from the hot inner disk mid-plane of FU Ori has been approximately stationary and the maximum dust grain size is $\gtrsim$1.6 mm in this region. If the hot inner disk of FU Ori which is inward of the 150-170 K water snowline is turbulent (e.g., corresponding to a Sunyaev & Shakura viscous $α_{t}\gtrsim$0.1), or if the actual maximum grain size is still larger than the lower limit we presently constrain, then as suggested by the recent analytical calculations and the laboratory measurements, water-ice free dust grains may be stickier than water-ice coated dust grains in protoplanetary disks. Additionally, we find that the free-free emission and the Johnson B and V bands magnitudes of these binary stars are brightening in 2016-2020. The optical and radio variability might be related to the dynamically evolving protostellar or disk accretion activities. Our results highlight that hot inner disks of outbursting objects are important laboratories for testing models of dust grain growth. Given the active nature of such systems, to robustly diagnose the maximum dust grain sizes, it is important to carry out coordinated multi-wavelength radio observations. △ Less

Submitted 20 October, 2021; v1 submitted 18 October, 2021; originally announced October 2021.

Comments: 27 pages, 11 figures, accepted to ApJ

Time-dependent, long-term hydrodynamic simulations of the inner protoplanetary disk I: The importance of stellar magnetic torques

Authors: D. Steiner, L. Gehrig, B. W. Ratschiner, F. Ragossnig, E. I. Vorobyov, M. Güdel, E. A. Dorfi

Abstract: We conduct simulations of the inner regions of protoplanetary disks (PPDs) to investigate the effects of protostellar magnetic fields on their long-term evolution. We use an inner boundary model that incorporates the influence of a stellar magnetic field. The position of the inner disk is dependent on the mass accretion rate as well as the magnetic field strength. We use this model to study the re… ▽ More We conduct simulations of the inner regions of protoplanetary disks (PPDs) to investigate the effects of protostellar magnetic fields on their long-term evolution. We use an inner boundary model that incorporates the influence of a stellar magnetic field. The position of the inner disk is dependent on the mass accretion rate as well as the magnetic field strength. We use this model to study the response of a magnetically truncated inner disk to an episodic accretion event. Additionally, we vary the protostellar magnetic field strength and investigate the consequences of the magnetic field on the long-term behavior of PPDs. We use the fully implicit 1+1D TAPIR code which solves the axisymmetric hydrodynamic equations self-consistently. Our model allows us to investigate disk dynamics close to the star and to conduct long-term evolution simulations simultaneously and includes the radial radiation transport in the stationary diffusion limit. We include stellar magnetic torques, the influence of a pressure gradient, and a variable inner disk radius in the TAPIR code to describe the innermost disk region in a more self-consistent manner and can show that this approach alters the disk dynamics considerably compared to a simplified diffusive evolution equation, especially during outbursts. The influences of a prescribed stellar magnetic field, local pressure gradients, and a variable inner disk radius result in a more consistent description of the gas dynamics in the innermost regions of PPDs. Combining magnetic torques acting on the innermost disk regions with the long-term evolution of PPDs yields previously unseen results, whereby the whole disk structure is affected over its entire lifetime. Additionally, we want to emphasize that a combination of our 1+1D model with more sophisticated multi-dimensional codes could improve the understanding of PPDs even further. △ Less

Submitted 18 September, 2021; originally announced September 2021.

Journal ref: A&A 655, A110 (2021)

Pulsational instability of pre-main sequence models from accreting protostars. I. Constraining the input physics for accretion with spectroscopic parameters and stellar pulsations

Authors: T. Steindl, K. Zwintz, T. G. Barnes, M. Muellner, E. I. Vorobyov

Abstract: Context: The pre-main sequence evolution is often simplified by choosing classical initial models. These have large initial radii and sufficient uniform contraction to make them fully convective. Contrary to that, real stars are born as small protostellar seeds in collapsing molecular clouds and obtain their final mass by means of accretion. Aims: We aim to constrain the input physics of accretion… ▽ More Context: The pre-main sequence evolution is often simplified by choosing classical initial models. These have large initial radii and sufficient uniform contraction to make them fully convective. Contrary to that, real stars are born as small protostellar seeds in collapsing molecular clouds and obtain their final mass by means of accretion. Aims: We aim to constrain the input physics of accretion on protostellar seeds with observed spectroscopic parameters and stellar pulsations of young stellar objects and pre-main sequence stars. Methods: We conducted a literature search for spectroscopic samples of young stellar objects and pre-main sequence stars including all previously known pulsators. The sample size of pulsating pre-main sequence stars is increased by analysing TESS observations and presenting discoveries in CoRoT data. We employ MESA and GYRE to calculate evolutionary tracks of accreting protostellar seeds in a constant accretion scenario, the subsequent pre-main sequence evolution, and their pulsation properties. The results are then compared with observations to constrain the input physics. Results: We discuss 16 formerly unknown pulsating pre-main sequence stars and candidates that are either of SPB, $δ$ Scuti,$γ$ Doradus or $δ$ Scuti - $γ$ Doradus hybrid type. We find that evolutionary tracks with a mass accretion rate of $5\times10^{-6} M_\odot/{\rm yr}$ and fraction of injected accretion energy of $β=0.1$ provide the best results in enveloping the spectroscopic parameters of pre-main sequence stars in a constant accretion scenario. The calculated instability regions constrain the atmospheric boundary conditions to Eddington Gray atmospheres; we discuss the future potential for additional constraints by instability regions that are dependent on radial order. We present a possible candidate for pulsations in M-type young stellar objects. △ Less

Submitted 15 July, 2021; originally announced July 2021.

Comments: 21 pages plus 12 for the appendix, 17 figures plus 27 for the appendix, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 654, A36 (2021)

Massive compact disks around FU Orionis-type young eruptive stars revealed by ALMA

Authors: Á. Kóspál, F. Cruz-Sáenz de Miera, J. A. White, P. Ábrahám, L. Chen, T. Csengeri, R. Dong, M. M. Dunham, O. Fehér, J. D. Green, J. Hashimoto, Th. Henning, M. Hogerheijde, T. Kudo, H. B. Liu, M. Takami, E. I. Vorobyov

Abstract: FU Orionis-type objects (FUors) are low-mass pre-main sequence stars undergoing a temporary, but significant increase of mass accretion rate from the circumstellar disk onto the protostar. It is not yet clear what triggers the accretion bursts and whether the disks of FUors are in any way different from disks of non-bursting young stellar objects. Motivated by this, we conducted a 1.3 mm continuum… ▽ More FU Orionis-type objects (FUors) are low-mass pre-main sequence stars undergoing a temporary, but significant increase of mass accretion rate from the circumstellar disk onto the protostar. It is not yet clear what triggers the accretion bursts and whether the disks of FUors are in any way different from disks of non-bursting young stellar objects. Motivated by this, we conducted a 1.3 mm continuum survey of ten FUors and FUor-like objects with ALMA, using both the 7 m array and the 12 m array in two different configurations to recover emission at the widest possible range of spatial scales. We detected all targeted sources and several nearby objects as well. To constrain the disk structure, we fit the data with models of increasing complexity from 2D Gaussian to radiative transfer, enabling comparison with other samples modeled in a similar way. The radiative transfer modeling gives disk masses that are significantly larger than what is obtained from the measured millimeter fluxes assuming optically thin emission, suggesting that the FUor disks are optically thick at this wavelength. In comparison with samples of regular Class II and Class I objects, the disks of FUors are typically a factor of 2.9-4.4 more massive and a factor of 1.5-4.7 smaller in size. A significant fraction of them (65-70%) may be gravitationally unstable. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Comments: 34 pages, 8 tables, 29 figures, accepted for publication in the Astrophysical Journal Supplement Series

Accretion bursts in high-mass protostars: a new testbed for models of episodic accretion

Authors: Vardan G. Elbakyan, Sergei Nayakshin, Eduard I. Vorobyov, Alessio Caratti o Garatti, Jochen Eislöffel

Abstract: It is well known that low mass young stellar objects (LMYSOs) gain a significant portion of their final mass through episodes of very rapid accretion, with mass accretion rates up to $\dot M_* \sim 10^{-4} M_{\odot}$~yr$^{-1}$. Recent observations of high mass young stellar objects (HMYSO) with masses $M_* \gtrsim 10 M_{\odot}$ uncovered outbursts with accretion rates exceeding… ▽ More It is well known that low mass young stellar objects (LMYSOs) gain a significant portion of their final mass through episodes of very rapid accretion, with mass accretion rates up to $\dot M_* \sim 10^{-4} M_{\odot}$~yr$^{-1}$. Recent observations of high mass young stellar objects (HMYSO) with masses $M_* \gtrsim 10 M_{\odot}$ uncovered outbursts with accretion rates exceeding $\dot M_*\sim 10^{-3}M_{\odot}$~yr$^{-1}$. Here we examine which scenarios proposed in the literature so far to explain accretion bursts of LMYSOs can apply to the episodic accretion in HMYSOs. We utilise a 1D time dependent models of protoplanetary discs around HMYSOs to study burst properties. We find that discs around HMYSOs are much hotter than those around their low mass cousins. As a result, much more extended regions of the disc are prone to the thermal hydrogen ionisation and MRI activation instabilities. The former in particular is found to be ubiquitous in a very wide range of accretion rates and disc viscosity parameters. The outbursts triggered by these instabilities, however, always have too low $\dot M_*$, and are one to several orders of magnitude too long compared to those observed from HMYSOs so far. On the other hand, bursts generated by tidal disruptions of gaseous giant planets formed by the gravitational instability of the protoplanetary discs yield properties commensurate with observations, provided that the clumps are in the post-collapse configuration with planet radius $R_{\rm p} \gtrsim 10 $ Jupiter radii. Furthermore, if observed bursts are caused by disc ionisation instabilities then they should be periodic phenomena with the duration of the quiescent phase comparable to that of the bursts. This may yield potentially observable burst periodicity signatures in the jets, the outer disc, or the surrounding diffuse material of massive HMYSOs. (abridged) △ Less

Submitted 18 June, 2021; v1 submitted 16 June, 2021; originally announced June 2021.

Comments: 8 pages, 6 figures, Accepted to A&A Letters

Journal ref: A&A 651, L3 (2021)

Restoration of the Parameters of a Gas-Dust Disk Based on Its Synthetic Images

Authors: A. M. Skliarevskii, Ya. N. Pavlyuchenkov, E. I. Vorobyov

Abstract: The topic of the present study is combining a dynamic model of a protoplanetary disk with the computations of radiation transfer for obtaining synthetic spectra and disk images suitable for immediate comparison of the model with observations. Evolution of the disk was computed using the FEOSAD hydrodynamic model, which includes a self-consistent calculation of the dynamics of dust and gas in the 2… ▽ More The topic of the present study is combining a dynamic model of a protoplanetary disk with the computations of radiation transfer for obtaining synthetic spectra and disk images suitable for immediate comparison of the model with observations. Evolution of the disk was computed using the FEOSAD hydrodynamic model, which includes a self-consistent calculation of the dynamics of dust and gas in the 2D thin disk approximation. Radiation transfer was simulated by the open code RADMC-3D. Three phases of disk evolution were considered: a young gravitationally unstable disk, a disk during an accretion luminosity burst, and an evolved disk. For these stages, the influence of various processes upon the disk's thermal structure was analyzed, as well as the differences between the temperatures obtained in the initial dynamic model and in the model with a detailed calculation of the radiation transfer. It is shown that viscous heating in the inner regions and adiabatic heating in the disk spirals can be important sources of heating. On the basis of the calculated spectral energy distributions, using SED-fitter software package used for the observations, physical parameters of the model disks were reconstructed. A significant spread between reconstructed parameters and initial characteristics of the disk indicates verification necessity of the models within the framework of spatially resolved observations of disks in the different spectral ranges △ Less

Submitted 21 April, 2021; originally announced April 2021.

Journal ref: Astron. Rep. 65, 170-183 (2021)

Gravitoviscous protoplanetary disks with a dust component. V. The dynamic model for freeze-out and sublimation of volatiles

Authors: Tamara Molyarova, Eduard I. Vorobyov, Vitaly Akimkin, Aleksandr Skliarevskii, Dmitri Wiebe, Manuel Güdel

Abstract: The snowlines of various volatile species in protoplanetary disks are associated with abrupt changes in gas composition and dust physical properties. Volatiles may affect dust growth, as they cover grains with icy mantles that can change the fragmentation velocity of the grains. In turn, dust coagulation, fragmentation, and drift through the gas disk can contribute to the redistribution of volatil… ▽ More The snowlines of various volatile species in protoplanetary disks are associated with abrupt changes in gas composition and dust physical properties. Volatiles may affect dust growth, as they cover grains with icy mantles that can change the fragmentation velocity of the grains. In turn, dust coagulation, fragmentation, and drift through the gas disk can contribute to the redistribution of volatiles between the ice and gas phases. Here we present the hydrodynamic model FEOSAD for protoplanetary disks with two dust populations and volatile dynamics. We compute the spatial distributions of major volatile molecules (H$_2$O, CO$_2$, CH$_4$, and CO) in the gas, on small and grown dust, and analyze the composition of icy mantles over the initial 0.5 Myr of disk evolution. We show that most of ice arrives to the grown dust through coagulation with small grains. Spiral structures and dust rings forming in the disk, as well as photodissociation in the outer regions, lead to the formation of complex snowline shapes and multiple snowlines for each volatile species. During the considered disk evolution, the snowlines shift closer to the star, with their final position being a factor $4-5$ smaller than that at the disk formation epoch. We demonstrate that volatiles tend to collect in the vicinity of their snowlines, both in the ice and gas phases, leading to the formation of thick icy mantles potentially important for dust dynamics. The dust size is affected by a lower fragmentation velocity of bare grains in the model with a higher turbulent viscosity. △ Less

Submitted 10 March, 2021; originally announced March 2021.

Comments: 28 pages, 14 figures, accepted to ApJ

Simulations of Dynamical Gas-Dust Circumstellar Disks: Going Beyond the Epstein Regime

Authors: Olga P. Stoyanovskaya, Fedor A. Okladnikov, Eduard I. Vorobyov, Yaroslav N. Pavlyuchenkov, Vitaliy V. Akimkin

Abstract: In circumstellar disks, the size of dust particles varies from submicron to several centimeters, while planetesimals have sizes of hundreds of kilometers. Therefore, various regimes for the aerodynamic drag between solid bodies and gas can be realized in these disks, depending on the grain sizes and velocities: Epstein, Stokes, and Newton, as well as transitional regimes between them. For small bo… ▽ More In circumstellar disks, the size of dust particles varies from submicron to several centimeters, while planetesimals have sizes of hundreds of kilometers. Therefore, various regimes for the aerodynamic drag between solid bodies and gas can be realized in these disks, depending on the grain sizes and velocities: Epstein, Stokes, and Newton, as well as transitional regimes between them. For small bodies moving in the Epstein regime, the time required to establish the constant relative velocity between the gas and bodies can be much less than the dynamical time scale for the problem - the time for the rotation of the disk about the central body. In addition, the dust may be concentrated in individual regions of the disk, making it necessary to take into account the transfer of momentum between the dust and gas. It is shown that, for a system of equations for gas and monodisperse dust, a semi-implicit first-order approximation scheme in time in which the interphase interaction is calculated implicitly, while other forces, such as the pressure gradient and gravity are calculated explicitly, is suitable for stiff problems with intense interphase interactions and for computations of the drag in non-linear regimes. The piece-wise drag coefficient widely used in astrophysical simulations has a discontinuity at some values of the Mach and Knudsen numbers that are realized in a circumstellar disk. A continuous drag coefficient is presented, which corresponds to experimental dependences obtained for various drag regimes. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: 24 pages, 13 figures

Journal ref: Astronomy Reports 64 (2020) 107-125

Distinguishing between different mechanisms of FU-Orionis-type luminosity outbursts

Authors: Eduard I. Vorobyov, Vardan G. Elbakyan, Hauyu Baobab Liu, Michihiro Takami

Abstract: Aims. Accretion and luminosity bursts triggered by three distinct mechanisms: the magnetorotational instability in the inner disk regions, clump infall in gravitationally fragmented disks and close encounters with an intruder star, were studied to determine the disk kinematic characteristics that can help to distinguish between these burst mechanisms. Methods. Numerical hydrodynamics simulations i… ▽ More Aims. Accretion and luminosity bursts triggered by three distinct mechanisms: the magnetorotational instability in the inner disk regions, clump infall in gravitationally fragmented disks and close encounters with an intruder star, were studied to determine the disk kinematic characteristics that can help to distinguish between these burst mechanisms. Methods. Numerical hydrodynamics simulations in the thin-disk limit were employed to model the bursts in disk environments that are expected for each burst mechanism. Results. We found that the circumstellar disks featuring accretion bursts can bear kinematic features that are distinct for different burst mechanisms, which can be useful when identifying the burst origin. The disks in the stellar encounter and clump-infall models are characterized by tens of percent deviations from the Keplerian rotation, whie the disks in the MRI models are characterized only a few percent deviation, which is mostly caused by the gravitational instability that fuels the MRI bursts. Velocity channel maps also show distinct kinks and wiggles, which are caused by gas disk flows that are peculiar to each considered burst mechanism. The deviations of velocity channels in the burst-hosting disks from a symmetric pattern typical of Keplerian disks are strongest for the clump-infall and collision models, and carry individual features that may be useful for the identification of the corresponding burst mechanism. The considered burst mechanisms produce a variety of light curves with the burst amplitudes varying in the Δm=2.5-3.7 limits, except for the clump-infall model where Δm can reach 5.4, although the derived numbers may be affected by a small sample and boundary conditions. Conclusions. Burst triggering mechanisms are associated with distinct kinematic features in the burst-hosting disks that may be used for their identification. Abridged. △ Less

Submitted 5 January, 2021; originally announced January 2021.

Comments: Accepted for publication by Astronomy & Astrophysics

Journal ref: A&A 647, A44 (2021)

Parameter study for the burst mode of accretion in massive star formation

Authors: D. M. -A. Meyer, E. I. Vorobyov, V. G. Elbakyan, J. Eisloeffel, A. M. Sobolev, M. Stoehr

Abstract: It is now a widely held view that, in their formation and early evolution, stars build up mass in bursts. The burst mode of star formation scenario proposes that the stars grow in mass via episodic accretion of fragments migrating from their gravitationally-unstable circumstellar discs and it naturally explains the existence of observed pre-main-sequence bursts from high mass protostars. We presen… ▽ More It is now a widely held view that, in their formation and early evolution, stars build up mass in bursts. The burst mode of star formation scenario proposes that the stars grow in mass via episodic accretion of fragments migrating from their gravitationally-unstable circumstellar discs and it naturally explains the existence of observed pre-main-sequence bursts from high mass protostars. We present a parameter study of hydrodynamical models of massive young stellar objects (MYSOs) that explores the initial masses of the collapsing clouds (Mc = 60-200Mo) and ratio of rotational-to-gravitational energies (beta = 0:005-0:33). An increase in Mc and/or beta produces protostellar accretion discs that are more prone to develop gravitational instability and to experience bursts. We find that all MYSOs have bursts even if their pre-stellar core is such that beta <= 0.01. Within our assumptions, the lack of stable discs is therefore a major difference between low- and high-mass star formation mechanisms. All our disc masses and disk-to-star mass ratios Md=M* > 1 scale as a power-law with the stellar mass. Our results confirm that massive protostars accrete about 40-60% of their mass in the burst mode. The distribution of time periods between two consecutive bursts is bimodal: there is a short duration (~ 1-10 yr) peak corresponding to the short, faintest bursts and a long duration peak (at ~ 10^3-10^4 yr) corresponding to the long, FU-Orionis-type bursts appearing in later disc evolution, i.e., around 30 kyr after disc formation. We discuss this bimodality in the context of the structure of massive protostellar jets as potential signatures of accretion burst history. △ Less

Submitted 10 November, 2020; originally announced November 2020.

Comments: 22 pages, 11 figures. Accepted at MNRAS

Disk fragmentation and intermittent accretion onto supermassive stars

Authors: Ryoki Matsukoba, Eduard I. Vorobyov, Kazuyuki Sugimura, Sunmyon Chon, Takashi Hosokawa, Kazuyuki Omukai

Abstract: Supermassive stars (SMSs) with $\sim10^{4-5}~\mathrm{M}_{\odot}$ are candidate objects for the origin of supermassive black holes observed at redshift $z$>6. They are supposed to form in primordial-gas clouds that provide the central stars with gas at a high accretion rate, but their growth may be terminated in the middle due to the stellar ionizing radiation if the accretion is intermittent and i… ▽ More Supermassive stars (SMSs) with $\sim10^{4-5}~\mathrm{M}_{\odot}$ are candidate objects for the origin of supermassive black holes observed at redshift $z$>6. They are supposed to form in primordial-gas clouds that provide the central stars with gas at a high accretion rate, but their growth may be terminated in the middle due to the stellar ionizing radiation if the accretion is intermittent and its quiescent periods are longer than the Kelvin-Helmholtz (KH) timescales at the stellar surfaces. In this paper, we examine the role of the ionizing radiation feedback based on the accretion history in two possible SMS-forming clouds extracted from cosmological simulations, following their evolution with vertically-integrated two-dimensional hydrodynamic simulations with detailed thermal and chemical models. The consistent treatment of the gas thermal evolution is crucial for obtaining the realistic accretion history, as we demonstrate by performing an additional run with a barotropic equation of state, in which the fluctuation of the accretion rate is artificially suppressed. We find that although the accretion becomes intermittent due to the formation of spiral arms and clumps in gravitationally unstable disks, the quiescent periods are always shorter than the KH timescales, implying that SMSs can form without affected by the ionizing radiation. △ Less

Submitted 4 November, 2020; originally announced November 2020.

Comments: 14 pages, 11 figures, accepted for publication in MNRAS

ALMA and VLA Observations of EX Lupi in its Quiescent State

Authors: Jacob Aaron White, Á. Kóspál, A. G. Hughes, P. Ábrahám, V. Akimkin, A. Banzatti, L. Chen, F. Cruz-Sáenz de Miera, A. Dutrey, M. Flock, S. Guilloteau, A. S. Hales, T. Henning, K. Kadam, D. Semenov, A. Sicilia-Aguilar, R. Teague, E. I. Vorobyov

Abstract: Extreme outbursts in young stars may be a common stage of pre-main-sequence stellar evolution. These outbursts, caused by enhanced accretion and accompanied by increased luminosity, can also strongly impact the evolution of the circumstellar environment. We present ALMA and VLA observations of EX Lupi, a prototypical outburst system, at 100 GHz, 45 GHz, and 15 GHz. We use these data, along with ar… ▽ More Extreme outbursts in young stars may be a common stage of pre-main-sequence stellar evolution. These outbursts, caused by enhanced accretion and accompanied by increased luminosity, can also strongly impact the evolution of the circumstellar environment. We present ALMA and VLA observations of EX Lupi, a prototypical outburst system, at 100 GHz, 45 GHz, and 15 GHz. We use these data, along with archival ALMA 232 GHz data, to fit radiative transfer models to EX Lupi's circumstellar disk in its quiescent state following the extreme outburst in 2008. The best fit models show a compact disk with a characteristic dust radius of 45 au and a total mass of 0.01 M$_{\odot}$. Our modeling suggests grain growth to sizes of at least 3 mm in the disk, possibly spurred by the recent outburst, and an ice line that has migrated inward to $0.2-0.3$ au post-outburst. At 15 GHz, we detected significant emission over the expected thermal disk emission which we attribute primarily to stellar (gyro)synchrotron and free-free disk emission. Altogether, these results highlight what may be a common impact of outbursts on the circumstellar dust. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Comments: Accepted to ApJ, 15 pages, 8 figures

The effect of luminosity outbursts on the protoplanetary disk dynamics

Authors: E. I. Vorobyov, V. Elbakyan, Michihiro Takami, Hauyu Baobab Liu

Abstract: Aims: Response of a protoplanetary disk to luminosity bursts of various duration is studied with the purpose to determine the effect of the bursts on the strength and sustainability of gravitational instability in the disk. A special emphasis is paid to the spatial distribution of gas and grown dust (from 1 mm to a few cm) during and after the burst. Methods: Numerical hydrodynamics simulations us… ▽ More Aims: Response of a protoplanetary disk to luminosity bursts of various duration is studied with the purpose to determine the effect of the bursts on the strength and sustainability of gravitational instability in the disk. A special emphasis is paid to the spatial distribution of gas and grown dust (from 1 mm to a few cm) during and after the burst. Methods: Numerical hydrodynamics simulations using the FEOSAD code were employed to study the dynamics of gas and dust in the thin-disk limit. Dust-to-gas friction including back reaction and dust growth were also considered. Bursts of various duration (from 100 to 500 yr) were initiated in accordance with a thermally ignited magnetorotational instability. Luminosity curves for constant- and declining-magnitude bursts were adopted to represent two typical limiting cases for FU-Orionis-type eruptions. Results: The short-term effect of the burst is to reduce the strength of gravitational instability by heating and expanding the disk. The longest bursts with duration comparable to the revolution period of the spiral can completely dissolve the original two-armed spiral pattern in the gas disk by the end of the burst, while the shortest bursts only weaken the spiral pattern. The reaction of grown dust to the burst is somewhat different. The spiral-like initial distribution with deep cavities in the inter-armed regions transforms into a ring-like distribution with deep gaps. This transformation is most expressed for the longest-duration bursts. The long-term effect of the burst depends on the initial disk conditions at the onset of the burst. In some cases, vigorous disk fragmentation sets in several thousand years after the burst, which was absent in the model without the bursts (abridged). △ Less

Submitted 3 September, 2020; originally announced September 2020.

Comments: Accepted for publication by Astronomy & Astrophysics

Journal ref: A&A 643, A13 (2020)

Accretion bursts in low-metallicity protostellar disks

Authors: E. I. Vorobyov, V. G. Elbakyan, K. Omukai, T. Hosokawa, R. Matsukoba, M. Guedel

Abstract: The early evolution of protostellar disks with metallicities in the $Z=1.0-0.01~Z_\odot$ range was studied with a particular emphasis on the strength of gravitational instability and the nature of protostellar accretion in low-metallicity systems. Numerical hydrodynamics simulations in the thin-disk limit were employed that feature separate gas and dust temperatures, and disk mass-loading from the… ▽ More The early evolution of protostellar disks with metallicities in the $Z=1.0-0.01~Z_\odot$ range was studied with a particular emphasis on the strength of gravitational instability and the nature of protostellar accretion in low-metallicity systems. Numerical hydrodynamics simulations in the thin-disk limit were employed that feature separate gas and dust temperatures, and disk mass-loading from the infalling parental cloud cores. Models with cloud cores of similar initial mass and rotation pattern, but distinct metallicity were considered to distinguish the effect of metallicity from that of initial conditions. The early stages of disk evolution in low-metallicity models are characterized by vigorous gravitational instability and fragmentation. Disk instability is sustained by continual mass-loading from the collapsing core. The time period that is covered by this unstable stage is much shorter in the $Z=0.01~Z_\odot$ models as compared to their higher metallicity counterparts thanks to the higher mass infall rates caused by higher gas temperatures (that decouple from lower dust temperatures) in the inner parts of collapsing cores. Protostellar accretion rates are highly variable in the low-metallicity models reflecting a highly dynamical nature of the corresponding protostellar disks. The low-metallicity systems feature short, but energetic episodes of mass accretion caused by infall of inward-migrating gaseous clumps that form via gravitational fragmentation of protostellar disks. These bursts seem to be more numerous and last longer in the $Z=0.1~Z_\odot$ models in comparison to the $Z=0.01~Z_\odot$ case. Variable protostellar accretion with episodic bursts is not a particular feature of solar metallicity disks. It is also inherent to gravitationally unstable disks with metallicities up to 100 times lower than solar. △ Less

Submitted 30 June, 2020; originally announced June 2020.

Comments: Accepted for publication by Astronomy & Astrophysics

Thermal evolution of protoplanetary disks: from $β$-cooling to decoupled gas and dust temperatures

Authors: Eduard I. Vorobyov, Ryoki Matsukoba, Kazuyuki Omukai, Manuel Guedel

Abstract: Aims: We explore the long-term evolution of young protoplanetary disks with different approaches to computing the thermal structure determined by various cooling and heating processes in the disk and its surroundings. Methods: Numerical hydrodynamics simulations in the thin-disk limit were complemented with three thermal evolution schemes: a simplified $β$-cooling approach with and without irradia… ▽ More Aims: We explore the long-term evolution of young protoplanetary disks with different approaches to computing the thermal structure determined by various cooling and heating processes in the disk and its surroundings. Methods: Numerical hydrodynamics simulations in the thin-disk limit were complemented with three thermal evolution schemes: a simplified $β$-cooling approach with and without irradiation, in which the rate of disk cooling is proportional to the local dynamical time, a fiducial model with equal dust and gas temperatures calculated taking viscous heating, irradiation, and radiative cooling into account, and also a more sophisticated approach allowing decoupled dust and gas temperatures. Results: We found that the gas temperature may significantly exceed that of dust in the outer regions of young disks thanks to additional compressional heating caused by the infalling envelope material in the early stages of disk evolution and slow collisional exchange of energy between gas and dust in low-density disk regions. The outer envelope however shows an inverse trend with the gas temperatures dropping below that of dust. The global disk evolution is only weakly sensitive to temperature decoupling. Nevertheless, separate dust and gas temperatures may affect the chemical composition, dust evolution, and disk mass estimates. Constant-$β$ models without stellar and background irradiation fail to reproduce the disk evolution with more sophisticated thermal schemes because of intrinsically variable nature of the $β$-parameter. Constant-$β$ models with irradiation can better match the dynamical and thermal evolution, but the agreement is still incomplete. Conclusions: Models allowing separate dust and gas temperatures are needed when emphasis is placed on the chemical or dust evolution in protoplanetary disks, particularly in sub-solar metallicity environments. △ Less

Submitted 28 April, 2020; originally announced April 2020.

Comments: 19 pages, 21 Figures, accepeted for publication in Astronomy & Astrophysics

Gravitoviscous protoplanetary disks with a dust component. III. Evolution of gas, dust, and pebbles

Authors: Vardan G. Elbakyan, Anders Johansen, Michiel Lambrechts, Vitaly Akimkin, Eduard I. Vorobyov

Abstract: We study the dynamics and growth of dust particles in circumstellar disks of different masses that are prone to gravitational instability during the critical first Myr of their evolution. The dust component is made up of two different components: micron-sized dust and grown dust of evolving size. For the dust component, we considered the dust coagulation, fragmentation, momentum exchange with the… ▽ More We study the dynamics and growth of dust particles in circumstellar disks of different masses that are prone to gravitational instability during the critical first Myr of their evolution. The dust component is made up of two different components: micron-sized dust and grown dust of evolving size. For the dust component, we considered the dust coagulation, fragmentation, momentum exchange with the gas, and dust self-gravity. We found that the micron-sized dust particles grow rapidly in the circumstellar disk, reaching a few cm in size in the inner 100 au of the disk during less than 100 kyr after the disk formation, provided that fragmentation velocity is $30\rm~ms^{-1}$. Due to the accretion of micron dust particles from the surrounding envelope, which serves as a micron dust reservoir, the approximately cm-sized dust particles continue to be present in the disk for more than 900 kyr after the disk formation and maintain a dust-to-gas ratio close to 0.01. We show that a strong correlation exists between the gas and pebble fluxes in the disk. We find that radial surface density distribution of pebbles in the disk shows power-law distribution with an index similar to that of the Minimum-mass solar nebula (MMSN) regardless the disk mass. We also show that the gas surface density in our models agrees well with measurements of dust in protoplanetary disks of AS 209, HD 163296, and DoAr 25 systems. Pebbles are formed during the very early stages of protoplanetary disk evolution. They play a crucial role in the planet formation process. Our disc simulations reveal the early onset ($<10^5$ yr) of an inwards-drifting flux of pebble-sized particles that makes up approximately between one hundredth and one tenth of the gas mass flux, which appears consistent with mm-observations of discs. (Abridged) △ Less

Submitted 31 March, 2020; originally announced April 2020.

Comments: 13 pages, 10 figures

Journal ref: A&A 637, A5 (2020)

The origin of tail-like structures around protoplanetary disks

Authors: E. I. Vorobyov, A. M. Skliarevskii, V. G. Elbakyan, M. Takami, H. B. Liu, S. -Y. Liu, E. Akiyama

Abstract: We study the origin of tail-like structures recently detected around the disk of SU Aurigae and several FU~Orionis-type stars. Dynamic protostellar disks featuring ejections of gaseous clumps and quiescent protoplanetary disks experiencing a close encounter with an intruder star were modeled using the numerical hydrodynamics code FEOSAD. Both the gas and dust dynamics were taken into account, incl… ▽ More We study the origin of tail-like structures recently detected around the disk of SU Aurigae and several FU~Orionis-type stars. Dynamic protostellar disks featuring ejections of gaseous clumps and quiescent protoplanetary disks experiencing a close encounter with an intruder star were modeled using the numerical hydrodynamics code FEOSAD. Both the gas and dust dynamics were taken into account, including dust growth and mutual friction between the gas and dust components. Only plane-of-the-disk encounters were considered. Ejected clumps produce a unique type of tail that is characterized by a bow-shock shape. Such tails originate from the supersonic motion of ejected clumps through the dense envelope that often surrounds young gravitationally unstable protostellar disks. The ejected clumps either sit at the head of the tail-like structure or disperse if their mass is insufficient to withstand the head wind of the envelope. On the other hand, close encounters with quiescent protoplanetary disks produce three types of the tail-like structure; we define these as pre-collisional, post-collisional, and spiral tails. These tails can in principle be distinguished from one another by particular features of the gas and dust flow in and around them. We find that the brown-dwarf-mass intruders do not capture circumintruder disks during the encounter, while the subsolar-mass intruders can acquire appreciable circumintruder disks with elevated dust-to-gas ratios, which can ease their observational detection. However, this is true only for prograde collisions; the retrograde intruders fail to collect appreciable amounts of gas or dust from the disk of the target. The predicted mass of dust in the model tail-like structures is higher than what was inferred for similar structures in SU~Aur, FU~Ori, and Z~CMa, making their observational detection feasible. Abridged. △ Less

Submitted 13 February, 2020; originally announced February 2020.

Comments: Accepted for publication in Astronomy & Astrophysics

Evolution of a Viscous Protoplanetary Disk with Convectively Unstable Regions

Authors: Ya. N. Pavlyuchenkov, A. V. Tutukov, L. A. Maksimova, E. I. Vorobyov

Abstract: The role of convection in the gas-dust accretion disk around a young star is studied. The evolution of a Keplerian disk is modeled using the Pringle equation, which describes the time variations of the surface density under the action of turbulent viscosity. The distributions of the density and temperature in the polar directions are computed simultaneously in the approximation that the disk is hy… ▽ More The role of convection in the gas-dust accretion disk around a young star is studied. The evolution of a Keplerian disk is modeled using the Pringle equation, which describes the time variations of the surface density under the action of turbulent viscosity. The distributions of the density and temperature in the polar directions are computed simultaneously in the approximation that the disk is hydrostatically stable. The computations of the vertical structure of the disk take into account heating by stellar radiation, interstellar radiation, and viscous heating. The main factor governing evolution of the disk in this model is the dependence of the viscosity coefficient on the radius of the disk. The computations of this coefficient take into account the background viscosity providing the continuous accretion of the gas and the convective viscosity, which depends on the parameters of the convection at a given radius. The results of computations of the global evolution and morphology of the disk obtained in this approach are presented. It is shown that, in the adopted model, the accretion has burst-like character: after the inner part of the disk ( R < 3 AU) is filled with matter, this material is relatively fast discharged onto the star, after which the process is repeated. Our results may be useful for explaining the activity of young FU Ori and EX Lup objects. It is concluded that convection may be one of the mechanisms responsible for the non-steady pattern of accretion in protostellar disks. △ Less

Submitted 18 December, 2019; originally announced December 2019.

Comments: Accepted for publication in Astronomy Reports

Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) -- Full Data Release

Authors: Ian W. Stephens, Tyler L. Bourke, Michael M. Dunham, Philip C. Myers, Riwaj Pokhrel, John J. Tobin, Héctor G. Arce, Sarah I. Sadavoy, Eduard I. Vorobyov, Jaime E. Pineda, Stella S. R. Offner, Katherine I. Lee, Lars E. Kristensen, Jes K. Jørgensen, Mark A. Gurwell, Alyssa A. Goodman

Abstract: We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1… ▽ More We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1$^{\prime\prime}$ ($\sim$300 au). The protostars were observed with the 1.3 mm and 850 $μ$m receivers simultaneously to detect continuum at both wavelengths and molecular line emission from CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), N$_2$D$^+$(3-2), CO(3-2), HCO$^+$(4-3), and H$^{13}$CO$^+$(4-3). Some of the observations also used the SMA's recently upgraded correlator, SWARM, whose broader bandwidth allowed for several more spectral lines to be observed (e.g., SO, H$_2$CO, DCO$^+$, DCN, CS, CN). Of the main continuum and spectral tracers observed, 84% of the images and cubes had emission detected. The median C$^{18}$O(2-1) linewidth is $\sim$1.0 km s$^{-1}$, which is slightly higher than those measured with single-dish telescopes at scales of 3000-20000 au. Of the 74 targets, six are suggested to be first hydrostatic core candidates, and we suggest that L1451-mm is the best candidate. We question a previous continuum detection toward L1448 IRS2E. In the SVS13 system, SVS13A certainly appears to be the most evolved source, while SVS13C appears to be hotter and more evolved than SVS13B. The MASSES survey is the largest publicly available interferometric continuum and spectral line protostellar survey to date, and is largely unbiased as it only targets protostars in Perseus. All visibility ($uv$) data and imaged data are publicly available at https://dataverse.harvard.edu/dataverse/full_MASSES/. △ Less

Submitted 19 November, 2019; originally announced November 2019.

Comments: Accepted to ApJS

An ALMA Study of the FU-Ori Type Object V900 Mon: Implications for the Progenitor

Authors: Michihiro Takami, Tsu-Sheng Chen, Hauyu Baobab Liu, Naomi Hirano, Agnes Kospal, Peter Abraham, Eduard I. Vorobyov, Fernando Cruz-Saenz de Miera, Timea Csengeri, Joel Green, Michiel Hogerheijde, Tien-Hao Hsieh, Jennifer L. Karr, Ruobing Dong, Alfonso Trejo, Lei Chen

Abstract: We present ALMA observations of 12CO, 13CO, and C18O J=2--1 lines and the 230 GHz continuum for the FU Ori-type object (FUor) V900 Mon (d~1.5 kpc), for which the accretion burst was triggered between 1953 and 2009. We identified CO emission associated with a molecular bipolar outflow extending up to a ~10^4 au scale and a rotating molecular envelope extending over >10^4 au. The interaction with th… ▽ More We present ALMA observations of 12CO, 13CO, and C18O J=2--1 lines and the 230 GHz continuum for the FU Ori-type object (FUor) V900 Mon (d~1.5 kpc), for which the accretion burst was triggered between 1953 and 2009. We identified CO emission associated with a molecular bipolar outflow extending up to a ~10^4 au scale and a rotating molecular envelope extending over >10^4 au. The interaction with the hot energetic FUor wind, which was observed using optical spectroscopy, appears limited to a region within ~400 au of the star. The envelope mass and the collimation of the extended CO outflow suggest that the progenitor of this FUor is a low-mass Class I young stellar object (YSO). These parameters for V900 Mon, another FUor, and a few FUor-like stars are consistent with the idea that FUor outbursts are associated with normal YSOs. The continuum emission is marginally resolved in our observations with a 0."2x0."15 (~300x225 au) beam, and a Gaussian model provides a deconvolved FWHM of ~90 au. The emission is presumably associated with a dusty circumstellar disk, plus a possible contribution from a wind or a wind cavity close to the star. The warm compact nature of the disk continuum emission could be explained with viscous heating of the disk, while gravitational fragmentation in the outer disk and/or a combination of grain growth and their inward drift may also contribute to its compact nature. △ Less

Submitted 12 September, 2019; originally announced September 2019.

Comments: 21 pages, 11 figures, accepted by ApJ

Gravitoviscous protoplanetary disks with a dust component. II. Spatial distribution and growth of dust in a clumpy disk

Authors: Eduard I. Vorobyov, Vardan G. Elbakyan

Abstract: Spatial distribution and growth of dust in a clumpy protoplanetary disk subject to vigorous gravitational instability and fragmentation is studied numerically with sub-au resolution using the FEOSAD code. Hydrodynamics equations describing the evolution of self-gravitating and viscous protoplanetary disks in the thin-disk limit were modified to include a dust component consisting of two parts: sub… ▽ More Spatial distribution and growth of dust in a clumpy protoplanetary disk subject to vigorous gravitational instability and fragmentation is studied numerically with sub-au resolution using the FEOSAD code. Hydrodynamics equations describing the evolution of self-gravitating and viscous protoplanetary disks in the thin-disk limit were modified to include a dust component consisting of two parts: sub-micron-sized dust and grown dust with a variable maximum radius. The conversion of small to grown dust, dust growth, friction of dust with gas, and dust self-gravity were also considered. We found that the disk appearance is notably time-variable with spiral arms, dusty rings, and clumps, constantly forming, evolving, and decaying. As a consequence, the total dust-to-gas mass ratio is highly non-homogeneous throughout the disk extent, showing order-of-magnitude local deviations from the canonical 1:100 value. Gravitationally bound clumps formed through gravitational fragmentation have a velocity pattern that deviates notably from the Keplerian rotation. Small dust is efficiently converted into grown dust in the clump interiors, reaching a maximum radius of several decimeters. Concurrently, grown dust drifts towards the clump center forming a massive compact central condensation (70-100 $M_\oplus$). We argue that protoplanets may form in the interiors of inward migrating clumps before they disperse through the action of tidal torques. We foresee the formation of protoplanets at orbital distances of several tens of au with initial masses of gas and dust in the protoplanetary seed in the (0.25-1.6) $M_{\rm Jup}$ and (1.0-5.5) $M_\oplus$ limits, respectively. The final masses of gas and dust in the protoplanets may however be much higher due to accretion from surrounding massive metal-rich disks/envelopes. △ Less

Submitted 3 September, 2019; v1 submitted 28 August, 2019; originally announced August 2019.

Comments: 13 pages, 10 figures; The abstract is abridged to meet ArXiv size limit

Journal ref: A&A 631, A1 (2019)

Diagnosing 0.1-10 au scale morphology of the FU\,Ori disk using ALMA and VLTI/GRAVITY

Authors: Hauyu Baobab Liu, Antoine Mérand, Joel D. Green, Sebastián Pérez, Antonio S. Hales, Yao-Lun Yang, Michael M. Dunham, Yasuhiro Hasegawa, Thomas Henning, Roberto Galván-Madrid, Ágnes Kóspál, Michihiro Takami, Eduard I. Vorobyov, Zhaohuan Zhu

Abstract: We report new Atacama Large Millimeter/submillimeter Array Band 3 (86-100 GHz; $\sim$80 mas angular resolution) and Band 4 (146-160 GHz; $\sim$50 mas angular resolution) observations of the dust continuum emission towards the archetypal and ongoing accretion burst young stellar object FU Ori, which simultaneously covered its companion, FU Ori S. In addition, we present near-infrared (2-2.45 $μ$m)… ▽ More We report new Atacama Large Millimeter/submillimeter Array Band 3 (86-100 GHz; $\sim$80 mas angular resolution) and Band 4 (146-160 GHz; $\sim$50 mas angular resolution) observations of the dust continuum emission towards the archetypal and ongoing accretion burst young stellar object FU Ori, which simultaneously covered its companion, FU Ori S. In addition, we present near-infrared (2-2.45 $μ$m) observations of FU Ori taken with the General Relativity Analysis via VLT InTerferometrY (GRAVITY; $\sim$1 mas angular resolution) instrument on the Very Large Telescope Interferometer (VLTI). We find that the emission in both FU Ori and FU Ori S at (sub)millimeter and near infrared bands is dominated by structures inward of $\sim$10 au radii. We detected closure phases close to zero from FU Ori with VLTI/GRAVITY, which indicate the source is approximately centrally symmetric and therefore is likely viewed nearly face-on. Our simple model to fit the GRAVITY data shows that the inner 0.4 au radii of the FU Ori disk has a triangular spectral shape at 2-2.45 $μ$m, which is consistent with the H$_{2}$O and CO absorption features in a $\dot{M}\sim$10$^{-4}$ $M_{\odot}\,yr^{-1}$, viscously heated accretion disk. At larger ($\sim$0.4-10 au) radii, our analysis shows that viscous heating may also explain the observed (sub)millimeter and centimeter spectral energy distribution when we assume a constant, $\sim$10$^{-4}$ $M_{\odot}\,yr^{-1}$ mass inflow rate in this region. This explains how the inner 0.4 au disk is replenished with mass at a modest rate, such that it neither depletes nor accumulates significant masses over its short dynamic timescale. Finally, we tentatively detect evidence of vertical dust settling in the inner 10 au of the FU Ori disk, but confirmation requires more complete spectral sampling in the centimeter bands. △ Less

Submitted 8 August, 2019; originally announced August 2019.

Comments: 16 pages, 9 figures; Accepted to ApJ

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

Authors: D. M. -A. Meyer, A. Kreplin, S. Kraus, E. I. Vorobyov, L. Haemmerle, J. Eisloeffel

Abstract: 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 evolutio… ▽ More 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. △ Less

Submitted 14 September, 2019; v1 submitted 5 June, 2019; originally announced June 2019.

Comments: Published at MNRAS, link: http://adsabs.harvard.edu/abs/2019MNRAS.487.4473M Previous version (2019arXiv190602015M) was identical but had an incomplete title

Gravitoviscous protoplanetary disks with a dust component. I. The importance of the inner sub-au region

Authors: Eduard I. Vorobyov, Aleksandr M. Skliarevskii, Vardan G. Elbakyan, Yaroslav Pavlyuchenkov, Vitaly Akimkin, Manuel Guedel

Abstract: The central region of a circumstellar disk is difficult to resolve in global numerical simulations of collapsing cloud cores, but its effect on the evolution of the entire disk can be significant. We use numerical hydrodynamics simulations to model the long-term evolution of self-gravitating and viscous circumstellar disks in the thin-disk limit. Simulations start from the gravitational collapse o… ▽ More The central region of a circumstellar disk is difficult to resolve in global numerical simulations of collapsing cloud cores, but its effect on the evolution of the entire disk can be significant. We use numerical hydrodynamics simulations to model the long-term evolution of self-gravitating and viscous circumstellar disks in the thin-disk limit. Simulations start from the gravitational collapse of prestellar cores of 0.5--1.0~$M_\odot$ and both gaseous and dusty subsystems were considered, including a model for dust growth. The inner unresolved 1.0 au of the disk is replaced with a central "smart" cell (CSC) -- a simplified model that simulates physical processes that may occur in this region. We found that the mass transport rate through the CSC has an appreciable effect on the evolution of the entire disk. Models with slow mass transport form more massive and warmer disks and they are more susceptible to gravitational instability and fragmentation, including a newly identified episodic mode of disk fragmentation in the T Tauri phase of disk evolution. Models with slow mass transport through the CSC feature episodic accretion and luminosity bursts in the early evolution, while models with fast transport are characterized by a steadily declining accretion rate with low-amplitude flickering. Dust grows to a larger, decimeter size in the slow transport models and efficiently drifts in the CSC, where it accumulates reaching the limit when streaming instability becomes operational. We argue that gravitational instability, together with streaming instability likely operating in the inner disk regions, constitute two concurrent planet-forming mechanisms, which may explain the observed diversity of exoplanetary orbits (Abridged). △ Less

Submitted 16 August, 2019; v1 submitted 27 May, 2019; originally announced May 2019.

Comments: Accepted for publication in Astronomy \& Astrophysics

A Tail Structure Associated with Protoplanetary Disk around SU Aurigae

Authors: Eiji Akiyama, Eduard I. Vorobyov, Hauyu Baobabu Liu, Ruobing Dong, Jerome de Leon, Sheng-Yuan Liu, Motohide Tamura

Abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the CO ($J$=2--1) line emission from the protoplanetary disk around T-Tauri star SU Aurigae (hereafter SU Aur). Previous observations in optical and near infrared wavelengths find a unique structure in SU Aur. One of the highlights of the observational results is that an extended tail-like structure is associated with t… ▽ More We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the CO ($J$=2--1) line emission from the protoplanetary disk around T-Tauri star SU Aurigae (hereafter SU Aur). Previous observations in optical and near infrared wavelengths find a unique structure in SU Aur. One of the highlights of the observational results is that an extended tail-like structure is associated with the disk, indicating mass transfer from or into the disk. Here we report the discovery of the counterpart of the tail-like structure in CO gas extending more than 1000 au long. Based on geometric and kinematic perspectives, both of the disk and the tail-like structure components physically connect to each other. Several theoretical studies predicted the observed tail-like structure via the following possible scenarios, 1) a gaseous stream from the molecular cloud remnant, 2) collision with a (sub)stellar intruder or a gaseous blob from the ambient cloud, and 3) ejection of a planetary or brown dwarf mass object due to gravitational instability via multi-body gravitational interaction. Since the tail-like structures associated with the SU Aur disk is a new example following RW Aurigae, some disks may experience the internal or external interaction and drastically lose mass during disk evolution. △ Less

Submitted 26 February, 2019; originally announced February 2019.

Comments: 10 pages, 5 figures, accepted by AJ

Luminosity outburst chemistry in protoplanetary discs: going beyond standard tracers

Authors: Dmitri S. Wiebe, Tamara S. Molyarova, Vitaly V. Akimkin, Eduard I. Vorobyov, Dmitry A. Semenov

Abstract: The chemical influence of luminosity outbursts on the environments of young solar-type stars is explored. Species are categorised into several types according to their response to the outburst. The first and second types imply chemical changes only during the outburst (with slightly different behaviours). These response types are mostly observed close to the star and are caused by icy mantle evapo… ▽ More The chemical influence of luminosity outbursts on the environments of young solar-type stars is explored. Species are categorised into several types according to their response to the outburst. The first and second types imply chemical changes only during the outburst (with slightly different behaviours). These response types are mostly observed close to the star and are caused by icy mantle evaporation. However, mantles recover after the outburst almost immediately. A notable exception is benzene ice, which is accumulated on dust surfaces during and after the outburst, so that its abundance exceeds the pre-outburst level by orders of magnitude. The third type of response is mostly seen at the disc periphery and implies alteration of abundances during the outburst and preservation of these `abnormal' abundances for centuries. This behaviour is typical of organic compounds, like HCOOCH$_3$, CH$_3$CN, CH$_2$CO. Their presence in the dark disc regions can be a manifestation of the past outburst. CO and CO$_2$ only trace past outbursts at the remote disc regions. The outburst changes the C/O ratio, but it quickly returns to the pre-outburst value almost everywhere in the disc. An important factor determining the sensitivity of molecular composition to the outburst is the dust size distribution. The duration of the pre-outburst stage and of the outburst itself influence the chemical effects, if the burst duration is shorter than 50 yr and the duration of the quiescent phase between the bursts is shorter than 100 kyr. △ Less

Submitted 20 February, 2019; originally announced February 2019.

Comments: Accepted to MNRAS

Episodic excursions of low-mass protostars on the Hertzsprung-Russell diagram

Authors: Vardan G. Elbakyan, Eduard I. Vorobyov, Christian Rab, D. M. -A. Meyer, Manuel Güdel, Takashi Hosokawa, Harold Yorke

Abstract: Following our recent work devoted to the effect of accretion on the pre-main-sequence evolution of low-mass stars, we perform a detailed analysis of episodic excursions of low-mass protostars in the Hertzsprung-Russell (H-R) diagram triggered by strong mass accretion bursts typical of FU Orionis-type objects (FUors). These excursions reveal themselves as sharp increases in the stellar total lumino… ▽ More Following our recent work devoted to the effect of accretion on the pre-main-sequence evolution of low-mass stars, we perform a detailed analysis of episodic excursions of low-mass protostars in the Hertzsprung-Russell (H-R) diagram triggered by strong mass accretion bursts typical of FU Orionis-type objects (FUors). These excursions reveal themselves as sharp increases in the stellar total luminosity and/or effective temperature of the protostar and can last from hundreds to a few thousands of years, depending on the burst strength and characteristics of the protostar. During the excursions, low-mass protostars occupy the same part of the H-R diagram as young intermediate-mass protostars in the quiescent phase of accretion. Moreover, the time spent by low-mass protostars in these regions is on average a factor of several longer than that spent by the intermediate-mass stars in quiescence. During the excursions, low-mass protostars pass close to the position of most known FUors in the H-R diagram, but owing to intrinsic ambiguity the model stellar evolutionary tracks are unreliable in determining the FUor properties. We find that the photospheric luminosity in the outburst state may dominate the accretion luminosity already after a few years after the onset of the outburst, meaning that the mass accretion rates of known FUors inferred from the bolometric luminosity may be systematically overestimated, especially in the fading phase. △ Less

Submitted 28 December, 2018; originally announced December 2018.

Comments: 15 pages, 12 figures

On the episodic excursions of massive protostars in the Hertzsprung-Russell diagram

Authors: D. M. A. Meyer, L. Haemmerlé, E. I. Vorobyov

Abstract: Massive protostars grow and evolve under the effect of rapid accretion of circumstellar gas and dust, falling at high rates ($\ge 10^{-4}$-$10^{-3}\, \rm M_{\odot}\, \rm yr^{-1}$). This mass infall has been shown, both numerically and observationally, to be episodically interspersed by accretion of dense gaseous clumps migrating through the circumstellar disc to the protostellar surface, causing s… ▽ More Massive protostars grow and evolve under the effect of rapid accretion of circumstellar gas and dust, falling at high rates ($\ge 10^{-4}$-$10^{-3}\, \rm M_{\odot}\, \rm yr^{-1}$). This mass infall has been shown, both numerically and observationally, to be episodically interspersed by accretion of dense gaseous clumps migrating through the circumstellar disc to the protostellar surface, causing sudden accretion and luminous bursts. Using numerical gravito-radiation-hydrodynamics and stellar evolution calculations, we demonstrate that, in addition to the known bloating of massive protostars, variable episodic accretion further influences their evolutionary tracks of massive young stellar objects (MYSOs). For each accretion-driven flare, they experience rapid excursions toward more luminous, but colder regions of the Hertzsprung-Russell diagram. During these excursions, which can occur up to the end of the pre-main-sequence evolution, the photosphere of massive protostars can episodically release much less energetic photons and MYSOs surreptitiously adopt the same spectral type as evolved massive (supergiants) stars. Each of these evolutionary loop brings the young high-mass stars close to the forbidden Hayashi region and might make their surrounding H II regions occasionally fainter, before they recover their quiescent, pre-burst surface properties. We interpret such cold, intermittent pre-main-sequence stellar evolutionary excursions and the dipping variability of HII regions as the signature of the presence of a fragmenting circumstellar accretion disc surrounding the MYSOs. We conjecture that this mechanism might equivalently affect young stars in the intermediate-mass regime. △ Less

Submitted 27 December, 2018; originally announced December 2018.

Comments: Accepted to MNRAS

Development and application of fast methods for computing momentum transfer between gas and dust in supercomputer simulation of planet formation

Authors: Olga P. Stoyanovskaya, Vitaly V. Akimkin, Eduard I. Vorobyov, Tatiana A. Glushko, Yaroslav N. Pavlyuchenkov, Valeriy N. Snytnikov, Nikolay V. Snytnikov

Abstract: Circumstellar discs, from which planetary systems are formed, consist of gas, dust and solids. Simulations of self-consistent dynamics of gas, dust and solids in circumstellar discs is a challenging problem. In the paper we present fast algorithms for computing the drag force (momentum transfer) between solid phase and gas. These algorithms (a) are universal and applicable to dust and solids with… ▽ More Circumstellar discs, from which planetary systems are formed, consist of gas, dust and solids. Simulations of self-consistent dynamics of gas, dust and solids in circumstellar discs is a challenging problem. In the paper we present fast algorithms for computing the drag force (momentum transfer) between solid phase and gas. These algorithms (a) are universal and applicable to dust and solids with any sizes smaller than the mean free path of gas molecules, (b) can be used to calculate the momentum transfer between dust and gas instead of one-way effect, as it is done in many models, (c) can perform simulations, without a loss in accuracy, with the time step determined by gas-dynamic parameters rather than by drag force, and (d) are compatible with the widely used parallel algorithms for solving 3D equations of gas dynamics, hydrodynamic equations for dust, and the collisionless Boltzmann equation for large bodies. Preliminary results of supercomputer simulation of the gas-dust disc dynamics within the developed approach are reported. △ Less

Submitted 15 November, 2018; originally announced November 2018.

Comments: 14 pages, 2 figures

Journal ref: Journal of Physics: Conf. Series 1103 (2018) 012008

Science with an ngVLA: Resolving the Radio Complexity of EXor and FUor-type Systems with the ngVLA

Authors: Jacob Aaron White, Marc Audard, Péter Ábrahám, Lucas Cieza, Fernando Cruz-Sáenz de Miera, Michael M. Dunham, Joel D. Green, Manuel Güdel, Nicolas Grosso, Antonio Hales, Lee Hartmann, Kundan Kadam, Joel H. Kastner, Ágnes Kóspál, Sebastian Perez, Andreas Postel, Dary Ruiz-Rodriguez, Christian Rab, Eduard I. Vorobyov, Zhaohuan Zhu

Abstract: Episodic accretion may be a common occurrence in the evolution of young pre-main sequence stars and has important implications for our understanding of star and planet formation. Many fundamental aspects of what drives the accretion physics, however, are still unknown. The ngVLA will be a key tool in understanding the nature of these events. The high spatial resolution, broad spectral coverage, an… ▽ More Episodic accretion may be a common occurrence in the evolution of young pre-main sequence stars and has important implications for our understanding of star and planet formation. Many fundamental aspects of what drives the accretion physics, however, are still unknown. The ngVLA will be a key tool in understanding the nature of these events. The high spatial resolution, broad spectral coverage, and unprecedented sensitivity will allow for the detailed analysis of outburst systems. The proposed frequency range of the ngVLA allows for observations of the gas, dust, and non-thermal emission from the star and disk. △ Less

Submitted 15 October, 2018; originally announced October 2018.

Comments: 8 pages, 1 figure, To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA)

Analysis of Methods for Computing the Trajectories of Dust Particles in a Gas-Dust Circumstellar Disk

Authors: Olga P. Stoyanovskaya, Valeriy N. Snytnikov, Eduard I. Vorobyov

Abstract: A systematic analysis of methods for computing the trajectories of solid-phase particles applied in modern astrophysics codes designed for modeling gas-dust circumstellar disks has been carried out for the first time. The motion of grains whose velocities are determined mainly by the gas drag, that is, for which the stopping time or relaxation time for the velocity of the dust to the velocity of t… ▽ More A systematic analysis of methods for computing the trajectories of solid-phase particles applied in modern astrophysics codes designed for modeling gas-dust circumstellar disks has been carried out for the first time. The motion of grains whose velocities are determined mainly by the gas drag, that is, for which the stopping time or relaxation time for the velocity of the dust to the velocity of the gas t_{stop} is less than or comparable to the rotation period, are considered. The methods are analyzed from the point of view of their stability for computing the motions of small bodies, including dust grains less than 1 micron, which are strongly coupled to the gas. Two test problems are with analytical solutions. Fast first order accurate methods that make it possible to avoid additional restriction on the time step due to gas drag in computations of the motion of grains of any size are presented. For the conditions of a circumstellar disks, the error in the velocity computations obtained when using some stable methods become unacceptably large when the time step size τ>t_{stop}. For the radial migration of bodies that exhibit drifts along nearly Keplerian orbits, an asymptotic approximation, sometimes called the short friction time approximation or drift flux model, gives a relative error for the radial velocity equals to St^2 where St is the Stokes number, the ratio of the stopping time of the body to some fraction of the rotation period (dynamical time scale) in the disk. △ Less

Submitted 4 September, 2018; originally announced September 2018.

Comments: 21 pages, 8 figures

Journal ref: Astronomy Reports, December 2017, Volume 61, Issue 12, pp 1044-1060

Analysis of Numerical Algorithms for Computing Rapid Momentum Transfers between the Gas and Dust in Simulations of Circumstellar Disks

Authors: O. P. Stoyanovskaya, E. I. Vorobyov, V. N. Snytnikov

Abstract: Approaches used in modern numerical simulations of the dynamics of dust and gas in circumstellar disks are tested. The gas and dust are treated like interpenetrating continuous media that can exchange momentum. A stiff coupling between the gas and dust phases is typical for such disks, with the dust stopping time much less than the characteristic dynamical time scale. This imposes high demands on… ▽ More Approaches used in modern numerical simulations of the dynamics of dust and gas in circumstellar disks are tested. The gas and dust are treated like interpenetrating continuous media that can exchange momentum. A stiff coupling between the gas and dust phases is typical for such disks, with the dust stopping time much less than the characteristic dynamical time scale. This imposes high demands on the methods used to simulate the dust dynamics. A grid, piecewise-parabolic method is used as the basic algorithm for solving the gas-dynamical equations. Numerical solutions obtained using various methods to compute the momentum exchanges are presented for the case of monodisperse dust. Numerical solutions are obtained for shock tube problem and the propagation of sound waves in a gas-dust medium. The studied methods are compared in terms of their ability to model media with (a) an arbitrary (short or long) dust stopping time, and (b) an arbitrary dust concentration in the gas (varying the dust to gas mass ratio from 0.01 to 1). A method for computing the momentum exchange with infinite-order accuracy in time is identified, which makes it possible to satisfy the conditions (a) and (b) with minimal computational costs. A first-order method that shows similar results in the test computations is also presented. It is shown that the proposed first-order method for monodisperse dust can be extended to a regime when the dust is polydisperse; i.e., a regime represented by several fractions with different stopping times. Formulas for computing the gas and dust velocities for polydisperse dust with each fraction exchanging momentum with the gas are presented. △ Less

Submitted 8 August, 2018; originally announced August 2018.

Comments: 13 pages, 4 figures

Journal ref: Astronomy Reports, 2018, Vol.62, PP 455-468

Near-Infrared High-Resolution Imaging Polarimetry of FU Ori-Type Objects: Towards A Unified Scheme for Low-Mass Protostellar Evolution

Authors: Michihiro Takami, Guangwei Fu, Hauyu Baobab Liu, Jennifer L. Karr, Jun Hashimoto, Tomoyuki Kudo, Eduard I. Vorobyov, Ágnes Kóspál, Peter Scicluna, Ruobing Dong, Motohide Tamura, Tae-Soo Pyo, Misato Fukagawa, Toru Tsuribe, Michael M. Dunham, Thomas Henning, Jerome de Leon

Abstract: We present near-IR imaging polarimetry of five classical FU Ori-type objects (FU Ori, V1057 Cyg, V1515 Cyg, V1735 Cyg, Z CMa) with a $\sim$0\farcs1 resolution observed using HiCIAO+AO188 at Subaru Telescope. We observed scattered light associated with circumstellar dust around four of them (i.e., all but V1515 Cyg). Their polarized intensity distribution shows a variety of morphologies with arms,… ▽ More We present near-IR imaging polarimetry of five classical FU Ori-type objects (FU Ori, V1057 Cyg, V1515 Cyg, V1735 Cyg, Z CMa) with a $\sim$0\farcs1 resolution observed using HiCIAO+AO188 at Subaru Telescope. We observed scattered light associated with circumstellar dust around four of them (i.e., all but V1515 Cyg). Their polarized intensity distribution shows a variety of morphologies with arms, tails or streams, spikes and fragmented distributions, many of which were reported in our previous paper. The morphologies of these reflection nebulae significantly differ from many other normal young stellar objects (Class I-II objects). These structures are attributed to gravitationally unstable disks, trails of clump ejections, dust blown by a wind or a jet, and a stellar companion. We can consistently explain our results with the scenario that their accretion outbursts (FUor outbursts) are triggered by gravitationally fragmenting disks, and with the hypothesis that many low-mass young stellar objects experience such outbursts. △ Less

Submitted 22 August, 2018; v1 submitted 10 July, 2018; originally announced July 2018.

Comments: 21 pages, 10 figures, accepted by ApJ; resubmitted to astro-ph on August 22, 2018

Gravitational fragmentation and formation of giant protoplanets on tens-of-au orbits

Authors: Eduard I. Vorobyov, Vardan Elbakyan

Abstract: Migration of dense gaseous clumps that form in young protostellar disks via gravitational fragmentation is investigated to determine the likelihood of giant planet formation. High-resolution numerical hydrodynamics simulations in the thin-disk limit are employed to compute the formation and long-term evolution of a gravitationally unstable protostellar disk around a solar-mass star. We show that g… ▽ More Migration of dense gaseous clumps that form in young protostellar disks via gravitational fragmentation is investigated to determine the likelihood of giant planet formation. High-resolution numerical hydrodynamics simulations in the thin-disk limit are employed to compute the formation and long-term evolution of a gravitationally unstable protostellar disk around a solar-mass star. We show that gaseous clumps that form in the outer regions of the disk (>100 AU) through disk fragmentation are often perturbed by other clumps or disk structures, such as spiral arms, and migrate toward the central star on timescales from a few 10^3 to few 10^4 yr. The migration timescale is slowest when stellar motion in response to the disk gravity is considered. When approaching the star, the clumps first gain mass (up to several tens of M_Jup), but then quickly lose most of their diffuse envelopes through tidal torques. Part of the clump envelope can be accreted on the central star causing an FU-Ori-type accretion and luminosity outburst. The tidal mass loss helps the clumps to significantly slow down or even halt their inward migration at a distance of a few tens of AU from the protostar. The resulting clumps are heavily truncated both in mass and size compared to their wider-orbit counterparts, keeping only a dense and hot nucleus. During the inward migration, the temperature in the clump interiors may exceed the molecular hydrogen dissociation limit (2000 K) and the central region of the clump can collapse into a gas giant protoplanet. We argue that FU-Orionis-type luminosity outbursts may be the end product of disk fragmentation and clump inward migration, ushering the formation of giant protoplanets on tens-of-au orbits in systems such as HR~8799. △ Less

Submitted 20 June, 2018; originally announced June 2018.

Comments: 16 pages, 16 figures, accepted for publication in Astronomy & Astrophysics

Mass Assembly of Stellar Systems and their Evolution with the SMA - 1.3 mm Subcompact Data Release

Authors: Ian W. Stephens, Michael M. Dunham, Philip C. Myers, Riwaj Pokhrel, Tyler L. Bourke, Eduard I. Vorobyov, John J. Tobin, Sarah I. Sadavoy, Jaime E. Pineda, Stella S. R. Offner, Katherine I. Lee, Lars E. Kristensen, Jes K. Jørgensen, Alyssa A. Goodman, Héctor G. Arce, Mark Gurwell

Abstract: We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we ca… ▽ More We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we can characterize the evolution of protostars. This paper releases the MASSES 1.3 mm data from the subcompact configuration ($\sim$4$^{\prime\prime}$ or $\sim$1000 au resolution), which is the SMA's most compact array configuration. We release both $uv$ visibility data and imaged data for the spectral lines CO(2-1),$^{13}$CO(2-1), C$^{18}$O(2-1), and N$_2$D$^+$(3-2), as well as for the 1.3 mm continuum. We identify the tracers that are detected toward each source. We also show example images of continuum and CO(2-1) outflows, analyze C$^{18}$O(2-1) spectra, and present data from the SVS 13 star-forming region. The calculated envelope masses from the continuum show a decreasing trend with bolometric temperature (a proxy for age). Typical C$^{18}$O(2-1) linewidths are 1.45 km s$^{-1}$, which is higher than the C$^{18}$O linewidths detected toward Perseus filaments and cores. We find that N$_2$D$^+$(3-2) is significantly more likely to be detected toward younger protostars. We show that the protostars in SVS 13 are contained within filamentary structures as traced by C$^{18}$O(2-1) and N$_2$D$^+$(3-2). We also present the locations of SVS 13A's high velocity (absolute line-of-sight velocities $>$150 km s$^{-1}$) red and blue outflow components. Data can be downloaded from https://dataverse.harvard.edu/dataverse/MASSES . △ Less

Submitted 19 June, 2018; originally announced June 2018.

Comments: Accepted ApJS

A 1.3 mm SMA Survey of 29 Variable Young Stellar Objects

Authors: Hauyu Baobab Liu, Michael M. Dunham, Ilaria Pascucci, Tyler L. Bourke, Naomi Hirano, Steven Longmore, Sean Andrews, Carlos Carrasco-González, Jan Forbrich, Roberto Galván-Madrid, Josep M. Girart, Joel D. Green, Carmen Juárez, Ágnes Kóspál, Carlo F. Manara, Aina Palau, Michihiro Takami, Leonardo Testi, Eduard I. Vorobyov

Abstract: Young stellar objects (YSOs) may undergo periods of active accretion (outbursts), during which the protostellar accretion rate is temporarily enhanced by a few orders of magnitude. Whether or not these accretion outburst YSOs possess similar dust/gas reservoirs to each other, and whether or not their dust/gas reservoirs are similar as quiescent YSOs, are issues not yet clarified. The aim of this w… ▽ More Young stellar objects (YSOs) may undergo periods of active accretion (outbursts), during which the protostellar accretion rate is temporarily enhanced by a few orders of magnitude. Whether or not these accretion outburst YSOs possess similar dust/gas reservoirs to each other, and whether or not their dust/gas reservoirs are similar as quiescent YSOs, are issues not yet clarified. The aim of this work is to characterize the millimeter thermal dust emission properties of a statistically significant sample of long and short duration accretion outburst YSOs (i.e., FUors and EXors) and the spectroscopically identified candidates of accretion outbursting YSOs (i.e., FUor-like objects). We have carried out extensive Submillimeter Array (SMA) observations mostly at $\sim$225 GHz (1.33 mm) and $\sim$272 GHz (1.10 mm), from 2008 to 2017. We covered accretion outburst YSOs located at $<$1 kpc distances from the solar system. We analyze all the existing SMA data of such objects, both published and unpublished, in a coherent way to present a millimeter interferometric database of 29 objects. We obtained 21 detections at $>$3-$σ$ significance. Detected sources except for the two cases of V883 Ori and NGC 2071 MM3 were observed with $\sim$1$"$ angular resolution. Overall our observed targets show a systematically higher millimeter luminosity distribution than those of the $M_{*}>$0.3 $M_{\odot}$ Class II YSOs in the nearby ($\lesssim$400 pc) low-mass star-forming molecular clouds (e.g., Taurus, Lupus, Upp Scorpio, and Chameleon I). In addition, at 1 mm our observed confirmed binaries or triple-system sources are systematically fainter than the rest of the sources even though their 1 mm fluxes are broadly distributed. We may have detected $\sim$30-60\% millimeter flux variability from V2494 Cyg and V2495 Cyg, from the observations separated by $\sim$1 year. △ Less

Submitted 24 October, 2017; originally announced October 2017.

Comments: 23 pages, 17 figures, accepted to A&A

Journal ref: A&A 612, A54 (2018)

Accretion outbursts in massive star formation

Authors: D. M. -A. Meyer, E. I. Vorobyov, R. Kuiper, W. Kley

Abstract: Using the HPC ressources of the state of Baden-Württemberg, we modelled for the first time the luminous burst from a young massive star by accretion of material from its close environment. We found that the surroundings of young massive stars are shaped as a clumpy disk whose fragments provoke outbursts once they fall onto the protostar and concluded that similar strong luminous events observed in… ▽ More Using the HPC ressources of the state of Baden-Württemberg, we modelled for the first time the luminous burst from a young massive star by accretion of material from its close environment. We found that the surroundings of young massive stars are shaped as a clumpy disk whose fragments provoke outbursts once they fall onto the protostar and concluded that similar strong luminous events observed in high-mass star forming regions may be a signature of the presence of such disks. △ Less

Submitted 6 October, 2017; originally announced October 2017.

Comments: Proceedings of the 3rd bwHPC-Symposium, Heidelberg (2016), 2 pages, 2 figures

Journal ref: ISBN 978-3-946531-70-8 - published on 02.10.2017