-
Multi-wavelength detection of an ongoing FUOr-type outburst on a low-mass YSO
Authors:
Zhen Guo,
P. W. Lucas,
R. G. Kurtev,
J. Borissova,
V. Elbakyan,
C. Morris,
A. Bayo,
L. Smith,
A. Caratti o Garatti,
C. Contreras Peña,
D. Minniti,
J. Jose,
M. Ashraf,
J. Alonso-García,
N. Miller,
H. D. S. Muthu
Abstract:
During the pre-main-sequence evolution, Young Stellar Objects (YSOs) assemble most of their mass during the episodic accretion process. The rarely seen FUOr-type events (FUOrs) are valuable laboratories to investigate the outbursting nature of YSOs. Here, we present multi-wavelength detection of a high-amplitude eruptive source in the young open cluster VdBH 221 with an ongoing outburst, including…
▽ More
During the pre-main-sequence evolution, Young Stellar Objects (YSOs) assemble most of their mass during the episodic accretion process. The rarely seen FUOr-type events (FUOrs) are valuable laboratories to investigate the outbursting nature of YSOs. Here, we present multi-wavelength detection of a high-amplitude eruptive source in the young open cluster VdBH 221 with an ongoing outburst, including optical to mid-infrared time series and near-infrared spectra. The initial outburst has an exceptional amplitude of $>$6.3 mag in Gaia and 4.6 mag in $K_s$, with a peak luminosity up to 16 $L_{\odot}$ and a peak mass accretion rate of 1.4 $\times$ 10$^{-5}$ $M_\odot$ yr$^{-1}$. The optical to infrared spectral energy distribution (SED) of this object is consistent with a low-mass star (0.2$M_\odot$) with a modest extinction ($A_V < 2$ mag). A 100-d delay between optical and infrared rising stages is detected, suggesting an outside-in origin of the instability. The spectroscopic features of this object reveal a self-luminous accretion disc, very similar to FU Orionis, with a low line-of-sight extinction. Most recently, there has been a gradual increase in brightness throughout the wavelength range, possibly suggesting an enhancement of the mass accretion rate.
△ Less
Submitted 25 January, 2024;
originally announced January 2024.
-
Dust growth and pebble formation in the initial stages of protoplanetary disk evolution
Authors:
Eduard Vorobyov,
Igor Kulikov,
Vardan Elbakyan,
James McKevitt,
Manuel Guedel
Abstract:
Aims. The initial stages of planet formation may start concurrently with the formation of a gas-dust protoplanetary disk. This makes the study of the earliest stages of protoplanetary disk formation crucially important. Here we focus on dust growth and pebble formation in a protoplanetary disk that is still accreting from a parental cloud core. Methods. We have developed an original three-dimensio…
▽ More
Aims. The initial stages of planet formation may start concurrently with the formation of a gas-dust protoplanetary disk. This makes the study of the earliest stages of protoplanetary disk formation crucially important. Here we focus on dust growth and pebble formation in a protoplanetary disk that is still accreting from a parental cloud core. Methods. We have developed an original three-dimensional numerical hydrodynamics code, which computes the collapse of rotating clouds and disk formation on nested meshes using a novel hybrid Coarray Fortran-OpenMP approach for distributed and shared memory parallelization. Dust dynamics and growth are also included in the simulations. Results. We found that the dust growth from $\sim 1~μ$m to 1-10~mm already occurs in the initial few thousand years of disk evolution but the Stokes number hardly exceeds 0.1 because of higher disk densities and temperatures compared to the minimum mass Solar nebular. The ratio of the dust-to-gas vertical scale heights remains rather modest, 0.2--0.5, which may be explained by the perturbing action of spiral arms that develop in the disk soon after its formation. The dust-to-gas mass ratio in the disk midplane is highly nonhomogeneous throughout the disk extent and is in general enhanced by a factor of several compared to the fiducial 1:100 value. Low St hinders strong dust accumulation in the spiral arms compared to the rest of the disk and the nonsteady nature of the spirals is also an obstacle. The spatial distribution of pebbles in the disk midplane exhibits a highly nonhomogeneous and patchy character. The total mass of pebbles in the disk increases with time and reaches a few tens of Earth masses after a few tens of thousand years of disk evolution. Abridged.
△ Less
Submitted 4 January, 2024;
originally announced January 2024.
-
On the origin of accretion bursts in FUORs
Authors:
Sergei Nayakshin,
Vardan Elbakyan
Abstract:
Accretion luminosity of young star FU Ori increased from undetectable levels to hundreds of Solar luminosities in 1937 and remains nearly as high at the present time. In a recent paper we showed how Extreme Evaporation (EE) of a young gas giant planet that migrated to a 10 day orbit around the star may power FU Ori. However, our model assumed a power-law mass-radius relation for the evaporating pl…
▽ More
Accretion luminosity of young star FU Ori increased from undetectable levels to hundreds of Solar luminosities in 1937 and remains nearly as high at the present time. In a recent paper we showed how Extreme Evaporation (EE) of a young gas giant planet that migrated to a 10 day orbit around the star may power FU Ori. However, our model assumed a power-law mass-radius relation for the evaporating planet. Here we employ a stellar evolution code to model mass losing planets. We find that adiabatic planets expand rapidly, which results in runaway FUOR bursts. Super-adiabatic planets contract while losing mass; their outbursts are dimming with time. Long steadily declining bursts such as FU Ori require relatively fine tuned internal planetary structure, which may be rare. More commonly we find that super-adiabatic planets contract too rapidly and their EE falters, leading to FUOR burst stutter. This stutter allows a single planet to produce many short repeating bursts, which may be relevant to bursts observed in V346 Nor, V899, V1647. We compute broad band spectra of our best fitting scenario for FU Ori. Since the outburst is triggered behind the planet location, the mid-IR emission rises many months before the optical, similar to bursts in Gaia-17bpi and Gaia-18dvy. We show that in outbursts powered by the classic thermal instability, mid-IR lags the optical, whereas the dead zone activation models predict mid-IR light precede the optical burst by many years to decades. We comment on the stellar flyby scenario for FU Ori.
△ Less
Submitted 21 September, 2023;
originally announced September 2023.
-
Distinguishing Magnetized Disc Winds from Turbulent Viscosity through Substructure Morphology in Planet-forming Discs
Authors:
Yinhao Wu,
Yi-Xian Chen,
Haochang Jiang,
Ruobing Dong,
Enrique Macías,
Min-Kai Lin,
Giovanni P. Rosotti,
Vardan Elbakyan
Abstract:
The traditional paradigm of viscosity-dominated evolution of protoplanetary discs has been recently challenged by magnetized disc winds. However, distinguishing wind-driven and turbulence-driven accretion through observations has been difficult. In this study, we present a novel approach to identifying their separate contribution to angular momentum transport by studying the gap and ring morpholog…
▽ More
The traditional paradigm of viscosity-dominated evolution of protoplanetary discs has been recently challenged by magnetized disc winds. However, distinguishing wind-driven and turbulence-driven accretion through observations has been difficult. In this study, we present a novel approach to identifying their separate contribution to angular momentum transport by studying the gap and ring morphology of planet-forming discs in the ALMA continuum. We model the gap-opening process of planets in discs with both viscous evolution and wind-driven accretion by 2D multi-fluid hydrodynamical simulations. Our results show that gap-opening planets in wind-driven accreting discs generate characteristic substructures that differ from those in purely viscous discs. Specifically, we demonstrate that discs, where wind-driven accretion dominates the production of substructures, exhibit significant asymmetries. Based on the diverse outputs of mock images in the ALMA continuum, we roughly divide the planet-induced features into four regimes (moderate-viscosity dominated, moderate-wind dominated, strong-wind dominated, inviscid). The classification of these regimes sets up a potential method to constrain the strength of magnetized disc wind and viscosity based on the observed gap and ring morphology. We discuss the asymmetry feature in our mock images and its potential manifestation in ALMA observations.
△ Less
Submitted 27 May, 2023; v1 submitted 14 May, 2023;
originally announced May 2023.
-
Extreme evaporation of planets in hot thermally unstable protoplanetary discs: the case of FU Ori
Authors:
Sergei Nayakshin,
James E. Owen,
Vardan Elbakyan
Abstract:
Disc accretion rate onto low mass protostar FU Ori suddenly increased hundreds of times 85 years ago and remains elevated to this day. We show that the sum of historic and recent observations challenges existing FU Ori models. We build a theory of a new process, Extreme Evaporation (EE) of young gas giant planets in discs with midplane temperatures exceeding 30, 000 K. Such temperatures are reache…
▽ More
Disc accretion rate onto low mass protostar FU Ori suddenly increased hundreds of times 85 years ago and remains elevated to this day. We show that the sum of historic and recent observations challenges existing FU Ori models. We build a theory of a new process, Extreme Evaporation (EE) of young gas giant planets in discs with midplane temperatures exceeding 30, 000 K. Such temperatures are reached in the inner 0.1 AU during thermal instability bursts. In our 1D time-dependent code the disc and an embedded planet interact through gravity, heat, and mass exchange. We use disc viscosity constrained by simulations and observations of dwarf novae instabilities, and we constrain planet properties with a stellar evolution code. We show that dusty gas giants born in the outer self-gravitating disc reach the innermost disc in a $\sim$ 10,000 years with radius of $\sim 10 R_J$. We show that their EE rates are $\sim 10^{-5}$ Msun/yr; if this exceeds the background disc accretion activity then the system enters a planet-sourced mode. Like a stellar secondary in mass-transferring binaries, the planet becomes the dominant source of matter for the star, albeit for $\sim$ O(100) years. We find that a $\sim$ 6 Jupiter mass planet evaporating in a disc fed at a time-averaged rate of $\sim 10^{-6}$ Msun/yr appears to explain all that we currently know about FU Ori accretion outburst. More massive planets and/or planets in older less massive discs do not experience EE process. Future FUOR modelling may constrain planet internal structure and evolution of the earliest discs.
△ Less
Submitted 5 May, 2023;
originally announced May 2023.
-
Computing the gravitational potential on nested meshes using the convolution method
Authors:
Eduard Vorobyov,
James McKevitt,
Igor Kulikov,
Vardan Elbakyan
Abstract:
Aims. Our aim is to derive a fast and accurate method for computing the gravitational potential of astrophysical objects with high contrasts in density, for which nested or adaptive meshes are required. Methods. We present an extension of the convolution method for computing the gravitational potential to the nested Cartesian grids. The method makes use of the convolution theorem to compute the gr…
▽ More
Aims. Our aim is to derive a fast and accurate method for computing the gravitational potential of astrophysical objects with high contrasts in density, for which nested or adaptive meshes are required. Methods. We present an extension of the convolution method for computing the gravitational potential to the nested Cartesian grids. The method makes use of the convolution theorem to compute the gravitational potential using its integral form. Results. A comparison of our method with the iterative outside-in conjugate gradient and generalized minimal residual methods for solving the Poisson equation using nonspherically symmetric density configurations has shown a comparable performance in terms of the errors relative to the analytic solutions. However, the convolution method is characterized by several advantages and outperforms the considered iterative methods by factors 10--200 in terms of the runtime, especially when graphics processor units are utilized. The convolution method also shows an overall second-order convergence, except for the errors at the grid interfaces where the convergence is linear. Conclusions. High computational speed and ease in implementation can make the convolution method a preferred choice when using a large number of nested grids. The convolution method, however, becomes more computationally costly if the dipole moments of tightly spaced gravitating objects are to be considered at coarser grids.
△ Less
Submitted 17 January, 2023;
originally announced January 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.
-
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.
-
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.
-
Gap opening by planets in discs with magnetised winds
Authors:
Vardan Elbakyan,
Yinhao Wu,
Sergei Nayakshin,
Giovanni Rosotti
Abstract:
Planets open deep gaps in protoplanetary discs when their mass exceeds a gap opening mass, $M_{\rm gap}$. We use one- and two-dimensional simulations to study planet gap opening in discs with angular momentum transport powered by MHD disc winds. We parameterise the efficiency of the MHD disc wind angular momentum transport through a dimensionless parameter $α_{\rm dw}$, which is an analogue to the…
▽ More
Planets open deep gaps in protoplanetary discs when their mass exceeds a gap opening mass, $M_{\rm gap}$. We use one- and two-dimensional simulations to study planet gap opening in discs with angular momentum transport powered by MHD disc winds. We parameterise the efficiency of the MHD disc wind angular momentum transport through a dimensionless parameter $α_{\rm dw}$, which is an analogue to the turbulent viscosity $α_{\rm v}$. We find that magnetised winds are much less efficient in counteracting planet tidal torques than turbulence is. For discs with astrophysically realistic values of $α_{\rm dw}$, $M_{\rm gap}$ is always determined by the residual disc turbulence, and is a factor of a few to ten smaller than usually obtained for viscous discs. We introduce a gap opening criterion applicable for any values of $α_{\rm v}$ and $α_{\rm dw}$ that may be useful for planet formation population synthesis. We show that in discs powered by magnetised winds, growing planets detach from the disc at planet masses below $\sim 0.1M_{\rm Jup}$ inside 10 AU. This promotes formation of super-Earth planets rather than gas giants in this region, in particular precluding formation of hot jupiters in situ. On larger scales, ALMA gap opening planet candidates may be less massive than currently believed. Future high-resolution observations with instruments such as the extended ALMA, ngVLA, and SKA are likely to show abundant narrow annular features at $R < 10$ AU due to ubiquitous super-Earth planets.
△ Less
Submitted 23 June, 2022;
originally announced June 2022.
-
ALMA constraints on assembly of Core Accretion planets
Authors:
Sergei Nayakshin,
Vardan Elbakyan,
Giovanni Rosotti
Abstract:
Resolved dust continuum and CO line ALMA imaging, and in some cases detection of H$α$ emission, hint that young massive planets are abundant at wide separations in protoplanetary discs. Here, we show how these observations can probe the runaway phase of planetary growth in the Core Accretion theory. Planets in this phase have the right range of masses to account for the predominantly moderate cont…
▽ More
Resolved dust continuum and CO line ALMA imaging, and in some cases detection of H$α$ emission, hint that young massive planets are abundant at wide separations in protoplanetary discs. Here, we show how these observations can probe the runaway phase of planetary growth in the Core Accretion theory. Planets in this phase have the right range of masses to account for the predominantly moderate contrast gaps and rings seen in ALMA observations. However, we find that these planets gain mass and migrate inward very rapidly. As a result, the phase when they could produce gaps with properties similar to those observed is very short, i.e., $t_{\rm gap} \lesssim 0.1$~Myr, independently of the disc viscosity parameter. This would require many tens to hundreds of gas giant planets to be born per ALMA system, violating the available mass budget of solids in realistic discs. This also predicts preponderance of discs with very wide gaps or complete inner disc holes, which is not observed. We show that suppression of both planet accretion and migration by a factor of at least ten is a possible solution to these serious problems. Future population synthesis models of planet formation should aim to address both exoplanetary data of older discless planetary systems and ALMA discs with embedded planets in one framework.
△ Less
Submitted 11 April, 2022;
originally announced April 2022.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Knotty protostellar jets as a signature of episodic protostellar accretion?
Authors:
Eduard Vorobyov,
Vardan Elbakyan,
Adele Plunkett,
Michael Dunham,
Marc Audard,
Manuel Guedel,
Odysseas Dionatos
Abstract:
We aim at studying the causal link between the knotty jet structure in CARMA 7, a young Class 0 protostar in the Serpens South cluster, and episodic accretion in young protostellar disks. We used numerical hydrodynamics simulations to derive the protostellar accretion history in gravitationally unstable disks around solar-mass protostars. We compared the time spacing between luminosity bursts Δτ_m…
▽ More
We aim at studying the causal link between the knotty jet structure in CARMA 7, a young Class 0 protostar in the Serpens South cluster, and episodic accretion in young protostellar disks. We used numerical hydrodynamics simulations to derive the protostellar accretion history in gravitationally unstable disks around solar-mass protostars. We compared the time spacing between luminosity bursts Δτ_mod, caused by dense clumps spiralling on the protostar, with the differences of dynamical timescales between the knots Δτ_obs in CARMA 7. We found that the time spacing between the bursts have a bi-modal distribution caused by isolated and clustered luminosity bursts. The former are characterized by long quiescent periods between the bursts with Δτ_mod = a few * (10^3-10^4) yr, whereas the latter occur in small groups with time spacing between the bursts Δτ_mod= a few * (10-10^2) yr. For the clustered bursts, the distribution of Δτ_mod in our models can be fit reasonably well to the distribution of Δτ_obs in the protostellar jet of CARMA 7, if a certain correction for the (yet unknown) inclination angle with respect to the line of sight is applied. The K-S test on the model and observational data sets suggests the best-fit values for the inclination angles of 55-80 deg., which become narrower (75-80 deg.) if only strong luminosity bursts are considered. The dynamical timescales of the knots in the jet of CARMA 7 are too short for a meaningful comparison with the long time spacings between isolated bursts in our models. The exact sequences of time spacings between the luminosity bursts in our models and knots in the jet of CARMA 7 were found difficult to match. (abridged)
△ Less
Submitted 20 January, 2018;
originally announced January 2018.
-
Formation of freely floating sub-stellar objects via close encounters
Authors:
Eduard I. Vorobyov,
Maria Steinrueck,
Vardan Elbakyan,
Manuel Guedel
Abstract:
We numerically studied close encounters between a young stellar system hosting a massive, gravitationally fragmenting disk and an intruder diskless star with the purpose to determine the evolution of fragments that have formed in the disk prior to the encounter. Numerical hydrodynamics simulations in the non-inertial frame of reference of the host star were employed to simulate the prograde and re…
▽ More
We numerically studied close encounters between a young stellar system hosting a massive, gravitationally fragmenting disk and an intruder diskless star with the purpose to determine the evolution of fragments that have formed in the disk prior to the encounter. Numerical hydrodynamics simulations in the non-inertial frame of reference of the host star were employed to simulate the prograde and retrograde co-planar encounters. The initial configuration of the target system (star plus disk) was obtained via a separate numerical simulation featuring the gravitational collapse of a solar-mass pre-stellar core. We found that close encounters can lead to the ejection of fragments that have formed in the disk of the target prior to collision. In particular, prograde encounters are more efficient in ejecting the fragments than the retrograde encounters. The masses of ejected fragments are in the brown-dwarf mass regime. They also carry away an appreciable amount of gas in their gravitational radius of influence, implying that these objects may possess extended disks or envelopes, as also suggested by Thies et al. (2015). Close encounters can also lead to the ejection of entire spiral arms, followed by fragmentation and formation of freely-floating objects straddling the planetary mass limit. However, numerical simulations with a higher resolution are needed to confirm this finding.
△ Less
Submitted 23 August, 2017;
originally announced August 2017.
-
The effect of accretion on the pre-main-sequence evolution of low-mass stars and brown dwarfs
Authors:
Eduard Vorobyov,
Vardan Elbakyan,
Takashi Hosokawa,
Yuya Sakurai,
Manuel Guedel,
Harold Yorke
Abstract:
The pre-main-sequence evolution of low-mass stars and brown dwarfs is studied numerically starting from the formation of a protostellar/proto-brown dwarf seed and taking into account the mass accretion onto the central object during the initial several Myr of evolution. The stellar evolution was computed using the STELLAR evolution code developed by Yorke & Bodenheimer with recent modifications by…
▽ More
The pre-main-sequence evolution of low-mass stars and brown dwarfs is studied numerically starting from the formation of a protostellar/proto-brown dwarf seed and taking into account the mass accretion onto the central object during the initial several Myr of evolution. The stellar evolution was computed using the STELLAR evolution code developed by Yorke & Bodenheimer with recent modifications by Hosokawa et al. The mass accretion rates were taken from numerical hydrodynamics models of Vorobyov & Basu computing the circumstellar disk evolution starting from the gravitational collapse of pre-stellar cloud cores of various mass and angular momentum. The resulting stellar evolution tracks were compared with the isochrones and isomasses calculated using non-accreting models. We find that mass accretion in the initial several Myr of protostellar evolution can have a strong effect on the subsequent evolution of young stars and brown dwarfs. The disagreement between accreting and non-accreting models in terms of the total stellar luminosity L_st, stellar radius R_st and effective temperature T_eff depends on the thermal efficiency of accretion, i.e., on the fraction of accretion energy absorbed by the central object. The largest mismatch is found for the cold accretion case, in which essentially all accretion energy is radiated away. The relative deviations in L_st and R_st in this case can reach 50% for 1.0-Myr-old objects and remain notable even for 10-Myr-old objects. In the hot and hybrid accretion cases, in which a constant fraction of accretion energy is absorbed, the disagreement between accreting and non-accreting models becomes less pronounced, but still remains notable for 1.0-Myr-old objects. These disagreements may lead to the wrong age estimate for objects of (sub-)solar mass when using the isochrones based on non-accreting models (abridged).
△ Less
Submitted 1 June, 2017;
originally announced June 2017.
-
The chemistry of episodic accretion in embedded objects. 2D radiation thermo-chemical models of the post-burst phase
Authors:
Ch. Rab,
V. Elbakyan,
E. Vorobyov,
M. Güdel,
O. Dionatos,
M. Audard,
I. Kamp,
W. -F. Thi,
P. Woitke,
A. Postel
Abstract:
Episodic accretion is an important process in the evolution of young stars and their environment. The observed strong luminosity bursts of young stellar objects likely have a long lasting impact on the chemical evolution of the disk and envelope structure. We want to investigate observational signatures of the chemical evolution in the post-burst phase for embedded sources. With such signatures it…
▽ More
Episodic accretion is an important process in the evolution of young stars and their environment. The observed strong luminosity bursts of young stellar objects likely have a long lasting impact on the chemical evolution of the disk and envelope structure. We want to investigate observational signatures of the chemical evolution in the post-burst phase for embedded sources. With such signatures it is possible to identify targets that experienced a recent luminosity burst. We present a new model for episodic accretion chemistry based on the 2D, radiation thermo-chemical disk code ProDiMo. We have extended ProDiMo with a proper treatment for envelope structures. For a representative Class I model, we calculated the chemical abundances in the post-burst phase and produced synthetic observables like intensity maps and radial profiles. During a burst many chemical species, like CO, sublimate from the dust surfaces. As the burst ends they freeze out again (post-burst phase). This freeze-out happens from inside-out due to the radial density gradient in the disk and envelope structure. This inside-out freeze-out produces clear observational signatures in spectral line emission, like rings and distinct features in the slope of radial intensity profiles. We fitted synthetic C18O J=2-1 observations with single and two component fits and find that post-burst images are much better matched by the latter. Comparing the quality of such fits allows identification of post-burst targets in a model-independent way. Our models confirm that it is possible to identify post-burst objects from spatially resolved CO observations. However, to derive proper statistics, like frequencies of bursts, from observations it is important to consider aspects like the inclination and structure of the target and also dust properties as those have a significant impact on the freeze-out timescale.
△ Less
Submitted 1 June, 2017; v1 submitted 10 May, 2017;
originally announced May 2017.
-
The nature of very low luminosity objects (VeLLOs)
Authors:
Eduard I. Vorobyov,
Vardan Elbakyan,
Michael Dunham,
Manuel Guedel
Abstract:
The nature of very low luminosity objects with the internal luminosity L_obj <= 0.1 Lsun is investigated by means of numerical modeling coupling the core collapse simulations with the stellar evolution calculations. The gravitational collapse of a large sample of model cores in the mass range 0.1-2.0 Msun is investigated. Numerical simulations were started at the pre-stellar phase and terminated a…
▽ More
The nature of very low luminosity objects with the internal luminosity L_obj <= 0.1 Lsun is investigated by means of numerical modeling coupling the core collapse simulations with the stellar evolution calculations. The gravitational collapse of a large sample of model cores in the mass range 0.1-2.0 Msun is investigated. Numerical simulations were started at the pre-stellar phase and terminated at the end of the embedded phase when 90% of the initial core mass had been accreted onto the forming protostar plus disk system. The disk formation and evolution was studied using numerical hydrodynamics simulations, while the formation and evolution of the central star was calculated using a stellar evolution code. Three scenarios for mass accretion from the disk onto the star were considered: hybrid accretion in which a fraction of accreted energy absorbed by the protostar depends on the accretion rate, hot accretion wherein a fraction of accreted energy is constant, and cold accretion wherein all accretion energy is radiated away. Our conclusions on the nature of VeLLOs depend crucially on the character of protostellar accretion. In the hybrid accretion scenario, most VeLLOs (90.6%) are expected to be the first hydrostatic cores (FHSCs) and only a small fraction (9.4%) are true protostars. In the hot accretion scenario, all VeLLOs are FHSCs due to overly high photospheric luminosity of protostars. In the cold accretion scenario, on the contrary, the majority of VeLLOs belong to the Class I phase of stellar evolution. The reason is that the stellar photospheric luminosity, which sets the floor for the total internal luminosity of a young star, is lower in cold accretion, thus enabling more VeLLOs in the protostellar stage. VeLLOs are relatively rare objects occupying 7%-11% of the total duration of the embedded phase and their masses do not exceed 0.3 Msun. (abridged).
△ Less
Submitted 12 December, 2016;
originally announced December 2016.
-
Self-consistent evolution of accreting low-mass stars and brown dwarfs
Authors:
I. Baraffe,
V. G. Elbakyan,
E. I. Vorobyov,
G. Chabrier
Abstract:
We present self-consistent calculations coupling numerical hydrodynamics simulations of collapsing pre-stellar cores and stellar evolution models of accreting objects. We analyse the main impact of consistent accretion history on the evolution and lithium depletion of young low-mass stars and brown dwarfs. These consistent models confirm the generation of a luminosity spread in Herzsprung-Russell…
▽ More
We present self-consistent calculations coupling numerical hydrodynamics simulations of collapsing pre-stellar cores and stellar evolution models of accreting objects. We analyse the main impact of consistent accretion history on the evolution and lithium depletion of young low-mass stars and brown dwarfs. These consistent models confirm the generation of a luminosity spread in Herzsprung-Russell diagrams at ages $\sim$ 1-10 Myr. They also confirm that early accretion can produce objects with abnormal Li depletion, as found in a previous study that was based on arbitrary accretion rates. The results strengthen that objects with anomalously high level of Li depletion in young clusters should be extremely rare. We also find that early phases of burst accretion can produce coeval models of similar mass with a range of different Li surface abundances, and in particular with Li-excess compared to the predictions of non-accreting counterparts. This result is due to a subtle competition between the effect of burst accretion and its impact on the central stellar temperature, the growth of the stellar radiative core and the accretion of fresh Li from the accretion disk. Only consistent models could reveal such a subtle combination of effects. This new result could explain the recent, puzzling observations of Li-excess of fast rotators in the young cluster NGC 2264. Present self-consistent accreting models are available in electronic form.
△ Less
Submitted 26 August, 2016;
originally announced August 2016.