Source author record

Eduard I. Vorobyov

Eduard I. Vorobyov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

35works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2026arXiv

On the accuracy of mass and size measurements of young protoplanetary disks

Knowing the masses and sizes of protoplanetary disks is of fundamental importance for the contemporary theories of planet formation. However, their measurements are associated with large uncertainties. In this proof of concept study, we focus on the very early stages of disk evolution, concurrent with the formation of the protostellar seed, because it is then that the initial conditions for subsequent planet formation are likely established. Using three-dimensional hydrodynamic simulations of a protoplanetary disk followed by radiation transfer postprocessing, we constructed synthetic disk images at millimeter wavelengths. We then calculated the synthetic disk radii and masses using an algorithm that is often applied to observations of protoplanetary disks with ALMA, and compared the resulting values with the actual disk mass and size derived directly from hydrodynamic modeling. We paid specific attention to the effects of dust growth on the discrepancy between synthetic and intrinsic disk masses and radii. We find that the dust mass is likely underestimated in Band 6 by factors of 1.4-4.2 when Ossenkopf & Henning opacities and typical dust temperatures are used, but the discrepancy reduces in Band~3, where the dust mass can be even overestimated. Dust growth affects both disk mass and size estimates via the dust-size-dependent opacity, and extremely low values of dust temperature (~ several Kelvin) are required to recover the intrinsic dust mass when dust has grown to mm-sized grains and its opacity has increased. Dust mass estimates are weakly sensitive to the distance to the source, while disk radii may be seriously affected. We conclude that the accuracy of measuring the dust mass and disk radius during the formation of a protoplanetary disk also depends on the progress in dust growth. (Abridged)

preprint2022arXiv

Protostellar-disc fragmentation across all metallicities

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}$.

preprint2021arXiv

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

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.

preprint2021arXiv

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

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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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)

preprint2020arXiv

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

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.

preprint2019arXiv

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

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/.

preprint2016arXiv

An alternative model for the origin of gaps in circumstellar disks

Motivated by recent observational and numerical studies suggesting that collapsing protostellar cores may be replenished from the local environment, we explore the evolution of protostellar cores submerged in the external counter-rotating environment. These models predict the formation of counter-rotating disks with a deep gap in the gas surface density separating the inner disk (corotating with the star) and the outer counter-rotating disk. The properties of these gaps are compared to those of planet-bearing gaps that form in disks hosting giant planets. We employ numerical hydrodynamics simulations of collapsing cores that are replenished from the local counter-rotating environment, as well as numerical hydrodynamic simulations of isolated disks hosting giant planets, to derive the properties of the gaps that form in both cases. Our numerical simulations demonstrate that counter-rotating disks can form for a wide range of mass and angular momentum available in the local environment. The gap that separates both disks has a depletion factor smaller than 1%, can be located at a distance from ten to over a hundred AU from the star, and can propagate inward with velocity ranging from 1 AU/Myr to >100 AU/Myr. Unlike our previous conclusion, the gap can therefore be a long-lived phenomenon, comparable in some cases to the lifetime of the disk itself. For a proper choice of the planetary mass, the viscous α-parameter and the disk mass, the planet-bearing gaps and the gaps in counter-rotating disks may show a remarkable similarity in the gas density profile and depletion factor, which may complicate their observational differentiation.

preprint2016arXiv

Circumstellar Disks of the Most Vigorously Accreting Young Stars

Young stellar objects (YSOs) may not accumulate their mass steadily, as was previously thought, but in a series of violent events manifesting themselves as sharp stellar brightening. These events can be caused by fragmentation due to gravitational instabilities in massive gaseous disks surrounding young stars, followed by migration of dense gaseous clumps onto the star. We report our high angular resolution, coronagraphic near-infrared polarization imaging observations using the High Contrast Instrument for the Subaru Next Generation Adaptive Optics (HiCIAO) of the Subaru 8.2 m Telescope, towards four YSOs which are undergoing luminous accretion outbursts. The obtained infrared images have verified the presence of several hundred AUs scale arms and arcs surrounding these YSOs. In addition, our hydrodynamics simulations and radiative transfer models further demonstrate that these observed structures can indeed be explained by strong gravitational instabilities occurring at the beginning of the disk formation phase. The effect of those tempestuous episodes of disk evolution on star and planet formation remains to be understood.

preprint2016arXiv

Ejection of gaseous clumps from gravitationally unstable protostellar disks

We investigate the dynamics of gaseous clumps formed via gravitational fragmentation in young protostellar disks, focusing on the fragments that are ejected from the disk via many-body gravitational interaction. Numerical hydrodynamics simulations were employed to study the evolution of young protostellar disks formed from the collapse of rotating pre-stellar cores with mass in the 1.1-1.6 M_sun range. Protostellar disks formed in our models undergo gravitational fragmentation driven by continuing mass loading from parental collapsing cores. A few fragments can be ejected from the disk during the early evolution, but the low-mass fragments (< 15~M_Jup) disperse creating spectacular bow-type structures while passing through the disk and collapsing core. The least massive fragment that survived the ejection (21 M_Jup) straddles the planetary-mass limit, while the most massive ejected fragments (145 M_Jup) can break up into several pieces, leading to the ejection of wide separation binary clumps in the brown-dwarf mass range. About half of the ejected fragments are gravitationally bound, the majority is supported by rotation against gravity, and all fragments have the specific angular momentum that is much higher than that expected for brown dwarfs. We found that the internal structure of the ejected fragments is distinct from what would be expected for gravitationally contracting clumps formed via cloud core fragmentation, which can help to differentiate their origin. The ejection of fragments is an important process inherent to massive protostellar disks, which produces freely-floating pre-brown dwarf cores, regulates the disk and stellar masses, and potentially enriches the intracluster medium with processed dust and complex organics.

preprint2016arXiv

Misalignment of Outflow Axes in the Proto-Multiple Systems in Perseus

We investigate the alignment between outflow axes in nine of the youngest binary/multiple systems in the Perseus Molecular Cloud. These systems have typical member spacing larger than 1000 AU. For outflow identification, we use 12CO(2-1) and 12CO(3-2) data from a large survey with the Sub- millimeter Array: Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES). The distribution of outflow orientations in the binary pairs is consistent with random or preferentially anti-aligned distributions, demonstrating that these outflows are misaligned. This result suggests that these systems are possibly formed in environments where the distribution of angular momentum is complex and disordered, and these systems do not come from the same co-rotating structures or from an initial cloud with aligned vectors of angular momentum.

preprint2015arXiv

Absence of Significant Cool Disks in Young Stellar Objects Exhibiting Repetitive Optical Outbursts

We report Submillimeter Array (SMA) 1.3 mm high angular resolution observations towards the four EXor type outbursting young stellar objects (YSOs) VY Tau, V1118 Ori, V1143 Ori, and NY Ori. The data mostly show low dust masses $M_{dust}$ in the associated circumstellar disks. Among the sources, NY Ori possesses a relatively massive disk with $M_{dust} \sim 9 \times 10^{-4}$ $M_{\odot}$. V1118 Ori has a marginal detection equivalent to $M_{dust} \sim 6 \times 10^{-5}$ $M_{\odot}$. V1143 Ori has a non-detection also equivalent to $M_{dust} < 6 \times 10^{-5}$ $M_{\odot}$. For the nearest source VY Tau, we get a surprising non-detection which provides a stringent upper limit $M_{dust} < 6 \times 10^{-6}$ $M_{\odot}$. We interpret our findings as suggesting that the gas and dust reservoirs that feed the short duration, repetitive optical outbursts seen in some EXors may be limited to the small scale, innermost region of their circumstellar disks. This hot dust may have escaped our detection limits. Follow-up, more sensitive millimeter observations are needed to improve our understanding of the triggering mechanisms of EXor type outbursts.

preprint2015arXiv

Mass Assembly of Stellar Systems and Their Evolution with the SMA (MASSES). Multiplicity and the Physical Environment in L1448N

We present continuum and molecular line observations at 230 GHz and 345 GHz from the Sub-millimeter Array (SMA) toward three protostars in the Perseus L1448N region. The data are from the large project "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). Three dust continuum sources, Source B, Source NW, and Source A, are detected at both frequencies. These sources have corresponding emission peaks in C18O (J=2-1), 13CO (J=2-1), and HCO+ (J=4-3), and have offsets with N2D+ (J=3-2) peaks. High angular resolution data from a complimentary continuum survey with the Karl G. Jansky Very Large Array show that Source B is associated with three 8 mm continuum objects, Source NW with two, and Source A remains single. These results suggest that multiplicity in L1448N exists at different spatial scales from a few thousand AU to < 100 AU. Velocity gradients in each source obtained from two-dimensional fits to the SMA C18O emission are found to be perpendicular to within 20 degrees of the outflow directions as revealed by 12CO (J=2-1). We have observed that Sources B and NW with multiplicity have higher densities than Source A without multiplicity. This suggests that thermal Jeans fragmentation can be relevant in the fragmentation process. However, we have not observed a difference in the ratio between rotational and gravitational energy between sources with and without multiplicity. We also have not observed a trend between non-thermal velocity dispersions and the level of fragmentation. Our study has provided the first direct and comprehensive comparison between multiplicity and core properties in low-mass protostars, although based on small number statistics.

preprint2015arXiv

Stellar hydrodynamical modeling of dwarf galaxies: simulation methodology, tests, and first results

Cosmological simulations still lack numerical resolution or physical processes to simulate dwarf galaxies in sufficient details. Accurate numerical simulations of individual dwarf galaxies are thus still in demand. We aim at (i) studying in detail the coupling between stars and gas in a galaxy, exploiting the so-called stellar hydrodynamical approach, and (ii) studying the chemo-dynamical evolution of individual galaxies starting from self-consistently calculated initial gas distributions. We present a novel chemo-dynamical code in which the dynamics of gas is computed using the usual hydrodynamics equations, while the dynamics of stars is described by the stellar hydrodynamics approach, which solves for the first three moments of the collisionless Boltzmann equation. The feedback from stellar winds and dying stars is followed in detail. In particular, a novel and detailed approach has been developed to trace the aging of various stellar populations, which enables an accurate calculation of the stellar feedback depending on the stellar age. We build initial equilibrium models of dwarf galaxies that take gas self-gravity into account and present different levels of rotational support. Models with high rotational support develop prominent bipolar outflows; a newly-born stellar population in these models is preferentially concentrated to the galactic midplane. Models with little rotational support blow away a large fraction of the gas and the resulting stellar distribution is extended and diffuse. The stellar dynamics turns out to be a crucial aspect of galaxy evolution. If we artificially suppress stellar dynamics, supernova explosions occur in a medium heated and diluted by the previous activity of stellar winds, thus artificially enhancing the stellar feedback (abridged).

preprint2015arXiv

Supermassive star formation via episodic accretion: protostellar disc instability and radiative feedback efficiency

The formation of SMSs is a potential pathway to seed SMBHs in the early universe. A critical issue for forming SMSs is stellar UV feedback, which may limit the stellar mass growth via accretion. In this paper we study the evolution of an accreting SMS and its UV emissivity under conditions of realistic variable accretion from a self-gravitating circumstellar disc. First we conduct a 2D hydrodynamical simulation to follow the long-term protostellar accretion until the stellar mass exceeds $10^4~M_\odot$. The disc fragments due to gravitational instability, creating a number of small clumps that rapidly migrate inward to fall onto the star. The resulting accretion history is thus highly time-dependent: short episodic accretion bursts are followed by longer, relative quiescent phases. We show that the circumstellar disc for the so-called direct collapse model is more unstable and generates greater variability over shorter timescales than normal Pop III cases. We conduct a post-process stellar evolution calculation using the obtained accretion history. Our results show that, regardless of the strong variability of the accretion rates, the stellar radius monotonically increases with almost constant effective temperature at $T_{\rm eff} \simeq 5000$ K as the stellar mass increases. The resulting UV feedback is too weak to hinder mass accretion due to the low flux of stellar UV photons, thus verifying our implicit assumption of no stellar feedback during the hydrodynamic simulations. The insensitivity of stellar evolution to variable accretion is attributed to the fact that typical timescales of variability, $\lesssim 10^3$ years, are too short to affect the stellar structure. We argue that this evolution will continue until the SMS eventually collapses to produce a massive black hole by the general relativistic instability after the stellar mass reaches $\gtrsim 10^5~M_\odot$.

preprint2015arXiv

Variable protostellar accretion with episodic bursts

We present the latest development of the disk gravitational instability and fragmentation model, originally introduced by us to explain episodic accretion bursts in the early stages of star formation. Using our numerical hydrodynamics model with improved disk thermal balance and star-disk interaction, we computed the evolution of protostellar disks formed from the gravitational collapse of prestellar cores. In agreement with our previous studies, we find that cores of higher initial mass and angular momentum produce disks that are more favorable to gravitational instability and fragmentation, while a higher background irradiation and magnetic fields moderate the disk tendency to fragment. The protostellar accretion in our models is time-variable, thanks to the nonlinear interaction between different spiral modes in the gravitationally unstable disk, and can undergo episodic bursts when fragments migrate onto the star owing to the gravitational interaction with other fragments or spiral arms. Most bursts occur in the partly embedded Class I phase, with a smaller fraction taking place in the deeply embedded Class 0 phase and a few possible bursts in the optically visible Class II phase. The average burst duration and mean luminosity are found to be in good agreement with those inferred from observations of FU-Orionis-type eruptions. The model predicts the existence of two types of bursts: the isolated ones, showing well-defined luminosity peaks separated with prolonged periods (~ 10^4 yr) of quiescent accretion, and clustered ones, demonstrating several bursts occurring one after another during just a few hundred years. Finally, we estimate that 40\%--70\% of the star-forming cores can display bursts after forming a star-disk system.

preprint2014arXiv

Episodic Accretion in Young Stars

In the last twenty years, the topic of episodic accretion has gained significant interest in the star formation community. It is now viewed as a common, though still poorly understood, phenomenon in low-mass star formation. The FU Orionis objects (FUors) are long-studied examples of this phenomenon. FUors are believed to undergo accretion outbursts during which the accretion rate rapidly increases from typically $10^{-7}$ to a few $10^{-4}$ $M_\odot$ yr$^{-1}$, and remains elevated over several decades or more. EXors, a loosely defined class of pre-main sequence stars, exhibit shorter and repetitive outbursts, associated with lower accretion rates. The relationship between the two classes, and their connection to the standard pre-main sequence evolutionary sequence, is an open question: do they represent two distinct classes, are they triggered by the same physical mechanism, and do they occur in the same evolutionary phases? Over the past couple of decades, many theoretical and numerical models have been developed to explain the origin of FUor and EXor outbursts. In parallel, such accretion bursts have been detected at an increasing rate, and as observing techniques improve each individual outburst is studied in increasing detail. We summarize key observations of pre-main sequence star outbursts, and review the latest thinking on outburst triggering mechanisms, the propagation of outbursts from star/disk to disk/jet systems, the relation between classical EXors and FUors, and newly discovered outbursting sources -- all of which shed new light on episodic accretion. We finally highlight some of the most promising directions for this field in the near- and long-term.

preprint2014arXiv

On the reliability of protostellar disc mass measurements and the existence of fragmenting discs

We couple non-magnetic, hydrodynamical simulations of collapsing protostellar cores with radiative transfer evolutionary models to generate synthetic observations. We then use these synthetic observations to investigate the extent to which a simple method for measuring protostellar disc masses used in the literature recovers the intrinsic masses of the discs formed in the simulations. We evaluate the effects of contamination from the surrounding core, partially resolving out the disc, optical depth, fixed assumed dust temperatures, inclination, and the dust opacity law. We show that the combination of these effects can lead to disc mass underestimates by up to factors of 2-3 at millimeter wavelengths and up to an order of magnitude or larger at submillimeter wavelengths. The optically thin portions of protostellar discs are generally cooler in the Class I stage than the Class 0 stage since Class I discs are typically larger and more optically thick, and thus more shielded. The observed disc mass distribution closely resembles the intrinsic distribution if this effect is taken into account, especially at millimeter wavelengths where optical depth effects are minimized. Approximately 50%-70% of protostellar discs observed to date with this method are consistent with the masses of the gravitationally unstable discs formed in the simulations, suggesting that at least some protostellar discs are likely sufficiently massive to fragment. We emphasize key future work needed to confirm these results, including assembling larger, less biased samples, and using molecular line observations to distinguish between rotationally supported, Keplerian discs and magnetically supported pseudodiscs.

preprint2014arXiv

The effect of external environment on the evolution of protostellar disks

Using numerical hydrodynamics simulations we studied the gravitational collapse of pre-stellar cores of sub-solar mass embedded into a low-density external environment. Four models with different magnitude and direction of rotation of the external environment with respect to the central core were studied and compared with an isolated model. We found that the infall of matter from the external environment can significantly alter the disk properties as compared to those seen in the isolated model. Depending on the magnitude and direction of rotation of the external environment, a variety of disks can form including compact (<= 200 AU) ones shrinking in size due to infall of external matter with low angular momentum, as well as extended disks forming due to infall of external matter with high angular momentum. The former are usually stable against gravitational fragmentation, while the latter are prone to fragmentation and formation of stellar systems with sub-stellar/very-low-mass companions. In the case of counterrotating external environment, very compact (< 5 AU) and short-lived (<= a few * 10^5 yr) disks can form when infalling material has low angular momentum. The most interesting case is found for the infall of counterrotating external material with high angular momentum, leading to the formation of counterrotating inner and outer disks separated by a deep gap at a few tens AU. The gap migrates inward due to accretion of the inner disk onto the protostar, turns into a central hole, and finally disappears giving way to the outer strongly gravitationally unstable disk. This model may lead to the emergence of a transient stellar system with sub-stellar/very-low-mass components counterrotating with respect to that of the star.

preprint2014arXiv

The Evolution of Protostars: Insights from Ten Years of Infrared Surveys with Spitzer and Herschel

Stars form from the gravitational collapse of dense molecular cloud cores. In the protostellar phase, mass accretes from the core onto a protostar, likely through an accretion disk, and it is during this phase that the initial masses of stars and the initial conditions for planet formation are set. Over the past decade, new observational capabilities provided by the Spitzer Space Telescope and Herschel Space Observatory have enabled wide-field surveys of entire star-forming clouds with unprecedented sensitivity, resolution, and infrared wavelength coverage. We review resulting advances in the field, focusing both on the observations themselves and the constraints they place on theoretical models of star formation and protostellar evolution. We also emphasize open questions and outline new directions needed to further advance the field.

preprint2014arXiv

The Role for the Inner Disk in Mass Accretion to the Star in the Early Phase of Star Formation

A physical mechanism that drives FU Orionis-type outbursts is reconsidered. We study the effect of inner part of a circumstellar disk covering a region from near the central star to the radius of approximately $5$ AU (hereafter, the inner disk). Using the fluctuated mass accretion rate onto the inner disk $\dot{M}_{\rm out}$, we consider the viscous evolution of the inner disk and the time variability of the mass accretion rate onto the central star $\dot{M}_{\rm in}$ by means of numerical calculation of an unsteady viscous accretion disk in a one-dimensional axisymmetric model. First, we calculate the evolution of the inner disk assuming an oscillating $\dot{M}_{\rm out}$. It is shown that the time variability of $\dot{M}_{\rm in}$ does not coincide with $\dot{M}_{\rm out}$ due to viscous diffusion. Second, we investigate the properties of spontaneous outbursts with temporally constant $\dot{M}_{\rm out}$. Outburst occur only in a limited range of mass accretion rates onto the inner disk $10^{-10}<\dot{M}_{\rm out}< 3\times 10^{-6}~{\rm M}_{\odot} {\rm yr}^{-1}$ due to gravo-magneto limit cycle (GML). Finally, we discuss the case with a combination of episodic $\dot{M}_{\rm out}$ and accretion outbursts cause by the GML in the inner disk. The GML can drive accretion outbursts onto the star even for the case of fluctuating $\dot{M}_{\rm out}$, although fluctuations of $\dot{M}$ decay during transmitting the inner disk inwards. We newly identified two modes of outburst which are spontaneous one and stimulated one. In a stimulated mode of outburst, $\dot{M}_{\rm out}$ does appear directly in $\dot{M}_{\rm in}$ (the latter defining the stellar accretion luminosity). In a spontaneous mode of outburst, $\dot{M}_{\rm out}$ appears as the interval between outbursts.

preprint2013arXiv

Formation of giant planets and brown dwarfs on wide orbits

Context: We studied numerically the formation of giant planet (GP) and brown dwarf (BD) embryos in gravitationally unstable protostellar disks and compared our findings with directly-imaged, wide-orbit (>= 50 AU) companions known to-date. The viability of the disk fragmentation scenario for the formation of wide-orbit companions in protostellar disks around (sub-)solar mass stars was investigated. Methods: We used numerical hydrodynamics simulations of disk formation and evolution with an accurate treatment of disk thermodynamics. The use of the thin-disk limit allowed us to probe the long-term evolution of protostellar disks. We focused on models that produced wide-orbit GP/BD embryos, which opened a gap in the disk and showed radial migration timescales similar to or longer than the typical disk lifetime. Results: While disk fragmentation was seen in the majority of our models, only 6 models out of 60 revealed the formation of quasi-stable, wide-orbit GP/BD embryos. Disk fragmentation produced GP/BD embryos with masses in the 3.5-43 M_J range, covering the whole mass spectrum of directly-imaged, wide-orbit companions to (sub-)solar mass stars. On the other hand, our modelling failed to produce embryos on orbital distances <= 170 AU, whereas several directly-imaged companions were found at smaller orbits down to a few AU. Disk fragmentation also failed to produce wide-orbit companions around stars with mass <= 0.7 Msun, in disagreement with observations. Conclusions: Disk fragmentation is unlikely to explain the whole observed spectrum of wide-orbit companions to (sub-)solar-mass stars and other formation mechanisms, e.g., dynamical scattering of closely-packed companions onto wide orbits, should be invoked to account for companions at orbital distance from a few tens to \approx 150 AU and wide-orbit companions with masses of the host star <= 0.7 Msun. (abridged)

preprint2013arXiv

The Burst Mode of Accretion in Primordial Protostars

We study the formation and long-term evolution of primordial protostellar disks harbored by first stars using numerical hydrodynamics simulations in the thin-disk limit. The initial conditions are specified by pre-stellar cores with distinct mass, angular momentum, and temperature. This allows us to probe several tens of thousand years of the disk's initial evolution, during which we observe multiple episodes of fragmentation leading to the formation of gravitationally bound gaseous clumps within spiral arms. These fragments are torqued inward due to gravitational interaction with the spiral arms on timescales of 10^3 - 10^4 yr and accreted onto the growing protostar, giving rise to accretion and luminosity bursts. The burst phenomenon is fueled by continuing accretion of material falling onto the disk from the collapsing parent core, which replenishes the mass lost by the disk due to accretion, and triggers repetitive episodes of disk fragmentation. We show that the burst phenomenon is expected to occur for a wide spectrum of initial conditions in primordial pre-stellar cores and speculate on how the intense luminosities (~10^7 solar luminosities) produced by this mechanism may have important consequences for the disk evolution and subsequent growth of the protostar.

preprint2013arXiv

The effect of episodic accretion on the phase transition of CO and CO_2 in low-mass star formation

We study the evaporation and condensation of CO and CO_2 during the embedded stages of low-mass star formation by using numerical simulations. We focus on the effect of luminosity bursts, similar in magnitude to FUors and EXors, on the gas-phase abundance of CO and CO_2 in the protostellar disk and infalling envelope. The evolution of a young protostar and its environment is followed based on hydrodynamical models using the thin-disk approximation, coupled with a stellar evolution code and phase transformations of CO and CO_2. The accretion and associated luminosity bursts in our model are caused by disk gravitational fragmentation followed by quick migration of the fragments onto the forming protostar. We found that bursts with luminosity on the order of 100-200 L_sun can evaporate CO ices in part of the envelope. The typical freeze-out time of the gas-phase CO onto dust grains in the envelope (a few kyr) is much longer than the burst duration (100-200 yr). This results in an increased abundance of the gas-phase CO in the envelope long after the system has returned into a quiescent stage. In contrast, luminosity bursts can evaporate CO_2 ices only in the disk, where the freeze-out time of the gas-phase CO_2 is comparable to the burst duration. We thus confirm that luminosity bursts can leave long-lasting traces in the abundance of gas-phase CO in the infalling envelope, enabling the detection of recent bursts as suggested by previous semi-analytical studies.

preprint2012arXiv

A Hybrid Scenario for the Formation of Brown Dwarfs and Very Low Mass Stars

We present a calculation of protostellar disk formation and evolution in which gaseous clumps (essentially, the first Larson cores formed via disk fragmentation) are ejected from the disk during the early stage of evolution. This is a universal process related to the phenomenon of ejection in multiple systems of point masses. However, it occurs in our model entirely due to the interaction of compact, gravitationally-bound gaseous clumps and is free from the smoothing-length uncertainty that is characteristic of models using sink particles. Clumps that survive ejection span a mass range of 0.08--0.35 $M_\odot$, and have ejection velocities $0.8 \pm 0.35$ km s$^{-1}$, which are several times greater than the escape speed. We suggest that, upon contraction, these clumps can form substellar or low-mass stellar objects with notable disks, or even close-separation very-low-mass binaries. In this hybrid scenario, allowing for ejection of clumps rather than finished protostars/proto--brown-dwarfs, disk formation and the low velocity dispersion of low-mass objects are naturally explained, while it is also consistent with the observation of isolated low-mass clumps that are ejection products. We conclude that clump ejection and the formation of isolated low mass stellar and substellar objects is a common occurrence, with important implications for understanding the initial mass function, the brown dwarf desert, and the formation of stars in all environments and epochs.

preprint2012arXiv

Collisionless Stellar Hydrodynamics as an Efficient Alternative to N-body Methods

For simulations that deal only with dark matter or stellar systems, the conventional N-body technique is fast, memory efficient, and relatively simple to implement. However when including the effects of gas physics, mesh codes are at a distinct disadvantage compared to SPH. Whilst implementing the N-body approach into SPH codes is fairly trivial, the particle-mesh technique used in mesh codes to couple collisionless stars and dark matter to the gas on the mesh, has a series of significant scientific and technical limitations. These include spurious entropy generation resulting from discreteness effects, poor load balancing and increased communication overhead which spoil the excellent scaling in massively parallel grid codes. We propose the use of the collisionless Boltzmann moment equations as a means to model collisionless material as a fluid on the mesh, implementing it into the massively parallel FLASH AMR code. This approach, which we term "collisionless stellar hydrodynamics" enables us to do away with the particle-mesh approach. Since the parallelisation scheme is identical to that used for the hydrodynamics, it preserves the excellent scaling of the FLASH code already demonstrated on peta-flop machines. We find the classic hydrodynamic equations and Boltzmann moment equations can be reconciled under specific conditions, allowing us to generate analytic solutions for collisionless systems using conventional test problems. We confirm the validity of our approach using a suite of demanding test problems, including the use of a modified Sod shock test. We conclude by demonstrating the ability of our code to model complex phenomena by simulating the evolution of a spiral galaxy whose properties agree with those predicted by swing amplification theory. (Abridged)

preprint2012arXiv

Self-gravitating equilibrium models of dwarf galaxies and the minimum mass for star formation

We construct a series of model galaxies in rotational equilibrium consisting of gas, stars, and a fixed dark matter (DM) halo and study how these equilibrium systems depend on the mass and form of the DM halo, gas temperature, non-thermal and rotation support against gravity, and also on the redshift of galaxy formation. For every model galaxy we find the minimum gas mass M_g^min required to achieve a state in which star formation (SF) is allowed according to contemporary SF criteria. The obtained M_g^min--M_DM relations are compared against the baryon-to-DM mass relation M_b--M_DM inferred from the \LambdaCDM theory and WMAP4 data. Our aim is to construct realistic initial models of dwarf galaxies (DGs), which take into account the gas self-gravity and can be used as a basis to study the dynamical and chemical evolution of DGs. Rotating equilibria are found by solving numerically the steady-state momentum equation for the gas component in the combined gravitational potential of gas, stars, and DM halo using a forward substitution procedure. We find that for a given M_DM the value of M_g^min depends crucially on the gas temperature T_g, gas spin parameter α, degree of non-thermal support σ_eff, and somewhat on the redshift for galaxy formation z_gf. Depending on the actual values of T_g, α, σ_eff, and z_gf, model galaxies may have M_g^min that are either greater or smaller than M_b. Galaxies with M_DM \ga 10^9 M_sun are usually characterized by M_g^min \la M_b, implying that SF in such objects is a natural outcome as the required gas mass is consistent with what is available according to the \LambdaCDM theory. On the other hand, models with M_DM \la 10^9 M_sun are often characterized by M_g^min >> M_b, implying that they need much more gas than available to achieve a state in which SF is allowed. Abridged.

preprint2012arXiv

The Ejection of Low Mass Clumps During Star Formation

Modeling of the self-consistent formation and evolution of disks as a result of prestellar core collapse reveals an intense early phase of recurrent gravitational instability and clump formation. These clumps generally migrate inward due to gravitational interaction with trailing spiral arms, and can be absorbed into the central object. However, in situations of multiple clump formation, gravitational scattering of clumps can result in the ejection of a low mass clump. These clumps can then give rise to free-floating low mass stars, brown dwarfs, or even giant planets. Detailed modeling of this process in the context of present-day star formation reveals that these clumps start out essentially as Larson first cores and grow subsequently by accretion. In the context of Pop III star formation, preliminary indications are that the disk clumps may also be of low mass. This mechanism of clump formation and possible ejection provides a channel for the formation of low mass objects in the first generation of stars.

preprint2012arXiv

The Migrating Embryo Model for Disk Evolution

A new view of disk evolution is emerging from self-consistent numerical simulation modeling of the formation of circumstellar disks from the direct collapse of prestellar cloud cores. This has implications for many aspects of star and planet formation, including the growth of dust and high-temperature processing of materials. A defining result is that the early evolution of a disk is crucially affected by the continuing mass loading from the core envelope, and is driven into recurrent phases of gravitational instability. Nonlinear spiral arms formed during these episodes fragment to form gaseous clumps in the disk. These clumps generally migrate inward due to gravitational torques arising from their interaction with a trailing spiral arm. Occasionally, a clump can open up a gap in the disk and settle into a stable orbit, revealing a direct pathway to the formation of companion stars, brown dwarfs, or giant planets. At other times, when multiple clumps are present, a low mass clump may even be ejected from the system, providing a pathway to the formation of free-floating brown dwarfs and giant planets in addition to low mass stars. Finally, it has been suggested that the inward migration of gaseous clumps can provide the proper conditions for the transport of high-temperature processed solids from the outer disk to the inner disk, and even possibly accelerate the formation of terrestrial planets in the inner disk. All of these features arising from clump formation and migration can be tied together conceptually in a Migrating Embryo model for disk evolution that can complement the well-known Core Accretion model for planet formation.

preprint2011arXiv

Destruction of massive fragments in protostellar disks and crystalline silicate production

We present a mechanism for the crystalline silicate production associated with the formation and subsequent destruction of massive fragments in young protostellar disks. The fragments form in the embedded phase of star formation via disk fragmentation at radial distances \ga 50-100 AU and anneal small amorphous grains in their interior when the gas temperature exceeds the crystallization threshold of ~ 800 K. We demonstrate that fragments that form in the early embedded phase can be destroyed before they either form solid cores or vaporize dust grains, thus releasing the processed crystalline dust into various radial distances from sub-AU to hundred-AU scales. Two possible mechanisms for the destruction of fragments are the tidal disruption and photoevaporation as fragments migrate radially inward and approach the central star and also dispersal by tidal torques exerted by spiral arms. As a result, most of the crystalline dust concentrates to the disk inner regions and spiral arms, which are the likely sites of fragment destruction.

preprint2011arXiv

Embedded protostellar disks around (sub-)solar stars. II. Disk masses, sizes, densities, temperatures and the planet formation perspective

We present basic properties of protostellar disks in the embedded phase of star formation (EPSF), which is difficult to probe observationally using available observational facilities. We use numerical hydrodynamics simulations of cloud core collapse and focus on disks formed around stars in the 0.03-1.0 Msun mass range. Our obtained disk masses scale near-linearly with the stellar mass. The mean and median disk masses in the Class 0 and I phases (M_{d,C0}^{mean}=0.12 Msun, M_{d,C0}^{mdn}=0.09 Msun and M_{d,CI}^{mean}=0.18 Msun, M_{d,CI}^{mdn}=0.15 Msun, respectively) are greater than those inferred from observations by (at least) a factor of 2--3. We demonstrate that this disagreement may (in part) be caused by the optically thick inner regions of protostellar disks, which do not contribute to millimeter dust flux. We find that disk masses and surface densities start to systematically exceed that of the minimum mass solar nebular for objects with stellar mass as low as M_st=0.05-0.1 Msun. Concurrently, disk radii start to grow beyond 100 AU, making gravitational fragmentation in the disk outer regions possible. Large disk masses, surface densities, and sizes suggest that giant planets may start forming as early as in the EPSF, either by means of core accretion (inner disk regions) or direct gravitational instability (outer disk regions), thus breaking a longstanding stereotype that the planet formation process begins in the Class II phase.

preprint2011arXiv

Resolving the Luminosity Problem in Low-Mass Star Formation

We determine the observational signatures of protostellar cores by coupling two-dimensional radiative transfer calculations with numerical hydrodynamical simulations that predict accretion rates that both decline with time and feature short-term variability and episodic bursts caused by disk gravitational instability and fragmentation. We calculate the radiative transfer of the collapsing cores throughout the full duration of the collapse, using as inputs the core, disk, and protostellar masses, radii, and mass accretion rates predicted by the hydrodynamical simulations. From the resulting spectral energy distributions, we calculate standard observational signatures (bolometric luminosity, bolometric temperature, ratio of bolometric to submillimeter luminosity) to directly compare to observations. We show that the accretion process predicted by these models reproduces the full spread of observed protostars in both Lbol - Tbol and Lbol - core mass space, including very low luminosity objects, provides a reasonable match to the observed protostellar luminosity distribution, and resolves the long-standing luminosity problem. These models predict an embedded phase duration shorter than recent observationally determined estimates (0.12 Myr vs. 0.44 Myr), and a fraction of total time spent in Stage 0 of 23%, consistent with the range of values determined by observations. On average, the models spend 1.3% of their total time in accretion bursts, during which 5.3% of the final stellar mass accretes, with maximum values being 11.8% and 35.5% for the total time and accreted stellar mass, respectively. Time-averaged models that filter out the accretion variability and bursts do not provide as good of a match to the observed luminosity problem, suggesting that the bursts are required.

preprint2010arXiv

Embedded protostellar disks around (sub-)solar protostars. I. Disk structure and evolution

We perform a comparative numerical hydrodynamics study of embedded protostellar disks formed as a result of the gravitational collapse of cloud cores of distinct mass (M_cl=0.2--1.7 M_sun) and ratio of rotational to gravitational energy (β=0.0028--0.023). An increase in M_cl and/or βleads to the formation of protostellar disks that are more susceptible to gravitational instability. Disk fragmentation occurs in most models but its effect is often limited to the very early stage, with the fragments being either dispersed or driven onto the forming star during tens of orbital periods. Only cloud cores with high enough M_cl or βmay eventually form wide-separation binary/multiple systems with low mass ratios and brown dwarf or sub-solar mass companions. It is feasible that such systems may eventually break up, giving birth to rogue brown dwarfs. Protostellar disks of {\it equal} age formed from cloud cores of greater mass (but equal β) are generally denser, hotter, larger, and more massive. On the other hand, protostellar disks formed from cloud cores of higher β(but equal M_cl) are generally thinner and colder but larger and more massive. In all models, the difference between the irradiation temperature and midplane temperature \triangle T is small, except for the innermost regions of young disks, dense fragments, and disk's outer edge where \triangle T is negative and may reach a factor of two or even more. Gravitationally unstable, embedded disks show radial pulsations, the amplitude of which increases along the line of increasing M_cl and βbut tends to diminish as the envelope clears. We find that single stars with a disk-to-star mass ratio of order unity can be formed only from high-βcloud cores, but such massive disks are unstable and quickly fragment into binary/multiple systems.

preprint2010arXiv

Formation and survivability of giant planets on wide orbits

Motivated by the recent discovery of massive planets on wide orbits, we present a mechanism for the formation of such planets via disk fragmentation in the embedded phase of star formation. In this phase, the forming disk intensively accretes matter from the natal cloud core and undergoes several fragmentation episodes. However, most fragments are either destroyed or driven into the innermost regions (and probably onto the star) due to angular momentum exchange with spiral arms, leading to multiple FU-Ori-like bursts and disk expansion. Fragments that are sufficiently massive and form in the late embedded phase (when the disk conditions are less extreme) may open a gap and evolve into giant planets on typical orbits of several tens to several hundreds of AU. For this mechanism to work, the natal cloud core must have sufficient mass and angular momentum to trigger the burst mode and also form extended disks of the order of several hundreds of AU. When mass loading from the natal cloud core diminishes and the main fragmentation phase ends, such extended disks undergo a transient episode of contraction and density increase, during which they may give birth to a last and survivable set of giant planets on wide and relatively stable orbits.