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

Aims. 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. Methods. We solved the hydrodynamics equations for a self-gravitating and viscous circumstellar disk in a thin-disk limit using the FEOSAD numerical hydrodynamics code. 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. Results. 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 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 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. Conclusions. 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 (<105 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. Such a pebble flux would allow for the formation of planetesimals by streaming instability and the early growth of embryos by pebble accretion. We conclude that unlike the more common studies of isolated steady-state protoplanetary disks, more sophisticated global numerical simulations of circumstellar disk formation and evolution, including the pebble formation from the micron dust particles, are needed for performing realistic planet formation studies. © ESO 2020.We thank the anonymous referee for a insightful report, which helped to improve this paper. Research was financially supported by the Ministry of Science and Higher Education of the Russian Federation (State assignment in the field of scientific activity, Southern Federal University, 2020). V.G.E. acknowledges the Swedish Institute for a visitor grant allowing to conduct research at Lund University. A.J. was supported by the Swedish Research Council (grant 2018-04867), the Knut and Alice Wallenberg Foundation (grant 2012.0150) and the European Research Council (ERC Consolidator Grant 724687-PLANETESYS). M.L. was supported by the Knut and Alice Wallenberg Foundation (grant 2012.0150). V.A. was supported by RFBR grant 18-52-52006

Accretion Bursts in Low-metallicity Protostellar Disks

Aims. The early evolution of protostellar disks with metallicities in the Zâ =â 1.0â-â 0.01â Z⊠range was studied with a particular emphasis on the strength of gravitational instability and the nature of protostellar accretion in low-metallicity systems. Methods. Numerical hydrodynamics simulations in the thin-disk limit were employed that feature separate gas and dust temperatures, and disk mass-loading from the infalling parent 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 the initial conditions. Results. 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⊠models than in their higher metallicity counterparts thanks to the higher rates of mass infall caused by higher gas temperatures (which decouple from lower dust temperatures) in the inner parts of collapsing cores. Protostellar accretion rates are highly variable in the low-metallicity models reflecting the highly dynamic 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⊠models than in the Zâ =â 0.01â Z⊠case. Conclusions. 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. © ESO 2020.Acknowledgements. We are thankful to the anonymous referee for useful comments that helped to improve the manuscript. E. I. V. and M. G. acknowledge support from the Austrian Science Fund (FWF) under research grant P31635-N27. V. G. E. acknowledges the Swedish Institute for a visitor Grant allowing to conduct research at Lund University. The simulations were performed on the Vienna Scientific Cluster

Accretion Bursts in High-Mass Protostars: A New Test Bed for Models of Episodic Accretion

Aims. 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 ∗ ≥ ∼10-4 M· yr-1. Recent observations of high-mass young stellar objects (HMYSOs) with masses M≥ 10 M· uncovered outbursts with accretion rates exceeding ∗ ≥ ∼10-3 M· yr-1. Here, we examine which scenarios proposed in the literature so far to explain accretion bursts of LMYSOs can also apply to the episodic accretion in HMYSOs. Methods. We utilise 1D time-dependent models of protoplanetary discs around HMYSOs to study burst properties. Results. We find that discs around HMYSOs are much hotter than those around their low-mass cousins. As a result, a much more extended region of the disc is prone to the thermal hydrogen ionisation and magnetorotational 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 of an ∗ ≥ and are one to several orders of magnitude too long compared to those observed from HMYSOs to date. 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 Rp 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. Bursts produced by disruptions of planets or more massive objects are not expected to be periodic phenomena, although multiple bursts per protostar are possible. Conclusions. Observations and modelling of episodic accretion bursts across a wide range of young stellar object (YSO) masses is a new promising avenue to break the degeneracy between models of episodic accretion in YSOs. © 2021 ESO.Acknowledgements. We thank the anonymous referee for an insightful report, which helped to improve this Letter. We thank Andrey Sobolev for useful discussions. V. E. and S. N. acknowledge the funding from the UK Science and Technologies Facilities Council, grant No. ST/S000453/1. This work made use of the DiRAC Data Intensive service at Leicester, operated by the University of Leicester IT Services, which forms part of the STFC DiRAC HPC Facility (www.dirac.ac.uk). V. E. also made use of funds from the Ministry of Science and Higher Education of the Russian Federation (State assignment in the field of scientific activity, Southern Federal University, 2020). E. I. Vorobyov acknowledges support from the Russian Science Foundation grant 18-12-00193. A. C. G. has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 743029)

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

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-200 Modot) and ratio of rotational-to-gravitational energies (β = 0.005-0.33). An increase in Mc and/or β 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 β ≤ 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 disc-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 percent 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 ∼ 103 104 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. © 2020 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society

The origin of tail-like structures around protoplanetary disks

Aims. We study the origin of tail-like structures recently detected around the disk of SU Aurigae and several FU Orionis-type stars. Methods. 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. Results. 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 mass of gas in the tail varies in the range 0.85-11.8 MJup, while the total mass of dust lies in the 1.75-30.1 M⊙ range, with the spiral tails featuring the highest masses. The predicted mass of dust in the model tail-like structures is therefore higher than what was inferred for similar structures in SU Aur, FU Ori, and Z CMa, making their observational detection feasible. Conclusions. Tail-like structures around protostellar and protoplanetary disks can be used to infer interesting phenomena such as clump ejection or close encounters. In particular, the bow-shock morphology of the tails could point to clump ejections as a possible formation mechanism. Further numerical and observational studies are needed to better understand the detectability and properties of the tails. © ESO 2020.Work was supported by the Russian Fund for Fundamental Research, Russian-Taiwanese project 19-52-52011 and MoST project 108-2923-M-001-006-MY3. H.B.L. is supported by the Ministry of Science and Technology (MoST) of Taiwan, grant No. 108-2112-M-001-002-MY3. V.G.E. acknowledges the Swedish Institute for a travel grant allowing to visit Lund University. The simulations were performed on the Vienna Scientific Cluster

Burst occurrence in young massive stellar objects

Episodic accretion-driven outbursts are an extreme manifestation of accretion variability. It has been proposed that the development of gravitational instabilities in the proto-circumstellar medium of massive young stellar objects (MYSOs) can lead to such luminous bursts, when clumps of fragmented accretion discs migrate on to the star. We simulate the early evolution of MYSOs formed by the gravitational collapse of rotating 100-M pre-stellar cores and analyse the characteristics of the bursts that accompany their strongly time-variable protostellar light curve episodically. We predict that MYSOs spend ≈103 yr (≈1.7 per cent) of their modelled early 60 kyr experiencing eruptive phases, during which the peak luminosity exceeds the quiescent pre-burst values by factors from 2.5 to more than 40. Throughout these short time periods, they can acquire a substantial fraction (up to ≈50 per cent) of their zero-age main-sequence mass. Our findings show that fainter bursts are more common than brighter ones. We discuss our results in the context of the known bursting MYSOs, e.g. NGC 6334I-MM1 and S255IR-NIRS3, and propose that these monitored bursts are part of a long-time ongoing series of eruptions, which might, in the future, be followed by other luminous flares. © 2018 The Author(s)

Hadron attenuation in deep inelastic lepton-nucleus scattering

We present a detailed theoretical investigation of hadron attenuation in deep inelastic scattering (DIS) off complex nuclei in the kinematic regime of the HERMES experiment. The analysis is carried out in the framework of a probabilistic coupled-channel transport model based on the Boltzmann-Uehling-Uhlenbeck (BUU) equation, which allows for a treatment of the final-state interactions (FSI) beyond simple absorption mechanisms. Furthermore, our event-by-event simulations account for the kinematic cuts of the experiments as well as the geometrical acceptance of the detectors. We calculate the multiplicity ratios of charged hadrons for various nuclear targets relative to deuterium as a function of the photon energy nu, the hadron energy fraction z_h=E_h/nu and the transverse momentum p_T. We also confront our model results on double-hadron attenuation with recent experimental data. Separately, we compare the attenuation of identified hadrons (pi^\pm, \pi^0, K^\pm, p and pbar) on Ne and Kr targets with the data from the HERMES Collaboration and make predictions for a Xe target. At the end we turn towards hadron attenuation on Cu nuclei at EMC energies. Our studies demonstrate that (pre-)hadronic final-state interactions play a dominant role in the kinematic regime of the HERMES experiment while our present approach overestimates the attenuation at EMC energies.Comment: 61 pages, 19 figures, version accepted for publication in Phys. Rev.

Multiwavelength detection of an ongoing FUOr-type outburst on a low-mass YSO

©2024 The Author(s). Published by Oxford University Press on behalf of Royal Astronomical Society. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY), https://creativecommons.org/licenses/by/4.0/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 are valuable laboratories to investigate the outbursting nature of YSOs. Here, we present multiwavelength 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 Ks, with a peak luminosity up to 16 L⊙ and a peak mass accretion rate of 1.4 × 10-5 M⊙ yr-1. The optical to infrared spectral energy distribution of this object is consistent with a low-mass star (0.2 M⊙) with a modest extinction (AV < 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.Peer reviewe

Measurement of the Spin-Dependence of the pbar-p Interaction at the AD-Ring

We propose to use an internal polarized hydrogen storage cell gas target in the AD ring to determine for the first time the two total spin-dependent pbar-p cross sections sigma_1 and sigma_2 at antiproton beam energies in the range from 50 to 450 MeV. The data obtained are of interest by themselves for the general theory of pbar-p interactions since they will provide a first experimental constraint of the spin-spin dependence of the nucleon-antinucleon potential in the energy range of interest. In addition, measurements of the polarization buildup of stored antiprotons are required to define the optimum parameters of a future, dedicated Antiproton Polarizer Ring (APR), intended to feed a double-polarized asymmetric pbar-p collider with polarized antiprotons. Such a machine has recently been proposed by the PAX collaboration for the new Facility for Antiproton and Ion Research (FAIR) at GSI in Darmstadt, Germany. The availability of an intense stored beam of polarized antiprotons will provide access to a wealth of single- and double-spin observables, thereby opening a new window on QCD spin physics.Comment: 51 pages, 23 figures, proposal submitted to the SPS committee of CER