Source author record

E. I. Vorobyov

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

16works
5topics
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

16 published item(s)

preprint2022arXiv

The influence of metallicity on a combined stellar and disk evolution

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

preprint2020arXiv

Accretion bursts in low-metallicity protostellar disks

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.

preprint2020arXiv

The origin of tail-like structures around protoplanetary disks

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.

preprint2019arXiv

Evolution of a Viscous Protoplanetary Disk with Convectively Unstable Regions

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

preprint2016arXiv

On the existence of accretion-driven bursts in massive star formation

Accretion-driven luminosity outbursts are a vivid manifestation of variable mass accretion onto protostars. They are known as the so-called FU Orionis phenomenon in the context of low-mass protostars. More recently, this process has been found in models of primordial star formation. Using numerical radiation hydrodynamics simulations, we stress that present-day forming massive stars also experience variable accretion and show that this process is accompanied by luminous outbursts induced by the episodic accretion of gaseous clumps falling from the circumstellar disk onto the protostar. Consequently, the process of accretion-induced luminous flares is also conceivable in the high-mass regime of star formation and we propose to regard this phenomenon as a general mechanism that can affect protostars regardless of their mass and/or the chemical properties of the parent environment in which they form. In addition to the commonness of accretion-driven outbursts in the star formation machinery, we conjecture that luminous flares from regions hosting forming high-mass star may be an observational implication of the fragmentation of their accretion disks.

preprint2016arXiv

Self-consistent evolution of accreting low-mass stars and brown dwarfs

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.

preprint2015arXiv

Thermal structure of a protostellar envelope

A numerical hydrodynamical model for the evolution of spherically symmetric collapsing clouds, designed for the calculation of the thermal structure of these objects in both the prestellar and protostellar stages of their evolution, is presented. Distinctive features of the model include the possibility of independently describing the temperatures of the gas and dust, which is extremely important when calculating the thermal structure of prestellar and protostellar clouds, and the account of the radiation flux from the central protostar. This model is used to compare the theoretical density and temperature distributions with observations for nearby sites of star formation obtained with the Herschel Space Observatory. Application of the diffusion approximation with a flux limiter describes well the radial density and temperature distributions in protostellar clouds. However, significant differences between the model and observational density profiles were found for prestellar stages, suggesting the presence of appreciable deviations from equilibrium in the prestellar clouds. An approximate method for calculating the thermal structure of a cloud based on the adaptive $τ$-approximation is presented. Application of the $τ$-approximation yields good agreement with the diffusion approximation for the prestellar phase, but produces appreciable discrepancies for the protostellar phase, when the thermal structure of the accreting envelope is determined by the radiation of the protostar.

preprint2012arXiv

Evolution of the First Supernovae in Protogalaxies: Dynamics of Mixing of Heavy Elements

The paper considers the evolution of the supernova envelopes produced by Population III stars with masses of $M_*\sim 25-200 M_\odot$ located in non-rotating protogalaxies with masses of $M\sim 10^7 M_\odot$ at redshifts $z=12$, with dark-matter density profiles in the form of modified isothermal spheres. The supernova explosion occurs in the ionization zone formed by a single parent star. The properties of the distribution of heavy elements (metals) produced by the parent star are investigated, as well as the efficiency with which they are mixed with the primordial gas in the supernova envelope. In supernovae with high energies ($E\simgt 5\times 10^{52}$ erg), an appreciable fraction of the gas can be ejected from the protogalaxy, but nearly all the heavy elements remain in the protogalaxy. In explosions with lower energies ($E\simlt 3\times 10^{52}$ erg), essentially no gas and heavy elements are lost from the protogalaxy: during the first one to threemillion years, the gas and heavy elements are actively carried from the central region of the protogalaxy ($r\sim 0.1 r_{vir}$, where $r_{vir}$ is the virial radius of the protogalaxy), but an appreciable fraction of the mass of metals subsequently returns when the hot cavity cools and the envelope collapses. Supernovae with high energies ($E\simgt 5\times 10^{52}$ erg) are characterized by a very low efficiency of mixing of metals; their heavy elements are located in the small volume occupied by the disrupted envelope (in a volume comparable with that of the entire envelope), with most of the metals remaining inside the hot, rarified cavity of the envelope. (abridged)

preprint2012arXiv

Instabilities in the Ionization Zones Around the First Stars

We consider the evolution of the ionization zone around Population III stars with $M_*\sim 25-200 M_\odot$ in protogalaxies with $M\sim 10^7 M_\odot$ at redshifts $z = 12$, assuming that the dark matter profile is a modified isothermal sphere. We study the conditions for the growth of instabilities in the ionization zones. The Rayleigh-Taylor and thermal instabilities develop efficiently in the ionization zones around 25-40 $M_\odot$ stars, while this efficiency is lower for stars with $\sim 120 M_\odot$. For more massive stars ($\sim 200 M_\odot$), the flux of ionizing photons is strong enough to considerably reduce the gas density in the ionization zone, and the typical lifetimes of stars ($\sim 2$ Myr) are insufficient for the growth of instabilities. The gas in a protogalaxy with $M\sim 10^7 M_\odot$ with a 200 $M_\odot$ central star is completely ionized by the end of the star's lifetime; in the case of a 120 $M_\odot$ central star, only one-third of the total mass of gas is ionized. Thus, ionizing photons from stars with $M_*\simlt 120 M_\odot$ cannot leave protogalaxies with $M\simgt 10^7 M_\odot$. If the masses of the central stars are 25 and 40 $M_\odot$, the gas in protogalaxies of this mass remains essentially neutral. We discuss the consequences of the evolution of the ionization zones for the propagation of the envelope after the supernova explosions of the stars and the efficiency of enrichment of the intergalactic medium in heavy elements.

preprint2012arXiv

The Burst Mode of Accretion in Primordial Star Formation

We present simulation results for the formation and long-term evolution of a primordial protostellar disk harbored by a first star. Using a 2+1D nonaxisymmetric thin disk numerical simulation, together with a barotropic relation for the gas, we are able to probe ~20 kyr of the disk's evolution. During this time period we observe fragmentation leading to loosely bound gaseous clumps within the disk. These are then torqued inward and accreted onto the growing protostar, giving rise to a burst phenomenon. The luminous feedback produced by this mechanism may have important consequences for the subsequent growth of the protostar.

preprint2009arXiv

The age of blue LSB galaxies

(Abridged). Low metallicities, large gas-to-star mass ratios, and blue colors of most low surface brightness (LSB) galaxies imply that these systems may be younger than their high surface brightness counterparts. We seek to find observational signatures that can help to constrain the age of blue LSB galaxies. We use numerical hydrodynamic modelling to study the long-term (~13 Gyr) dynamical and chemical evolution of blue LSB galaxies adopting a sporadic scenario for star formation. Our models utilize various rates of star formation and different shapes of the initial mass function (IMF). We complement hydrodynamic modelling with population synthesis modelling to produce the integrated B-V colors and Halpha equivalent widths (EW(Ha)). We find that the mean oxygen abundances, B-V colors, EW(Ha), and the radial fluctuations in the oxygen abundance, when considered altogether, can be used to constrain the age of blue LSB galaxies if some independent knowledge of the IMF is available. Our modelling strongly suggests the existence of a minimum age for blue LSB galaxies. Model B-V colors and mean oxygen abundances set a tentative minimum age at 1.5-3.0 Gyr, whereas model EW(Ha) suggest a larger value of order 5-6 Gyr. The latter value may decrease somewhat, if blue LSB galaxies host IMFs with a truncated upper mass limit. We found no firm evidence that the age of blue LSB galaxies is significantly smaller than 13 Gyr.

preprint2009arXiv

The Bimodality of Accretion In T Tauri Stars and Brown Dwarfs

We present numerical solutions of the collapse of prestellar cores that lead to the formation and evolution of circumstellar disks. The disk evolution is then followed for up to three million years. A variety of models of different initial masses and rotation rates allows us to study disk accretion around brown dwarfs and low-mass T Tauri stars, with central object mass $M_st < 0.2 Msun$, as well as intermediate and upper-mass T Tauri stars (0.2 Msun < M_st < 3.0 Msun). Our models include self-gravity and allow for nonaxisymmetric motions. In addition to the self-consistently generated gravitational torques, we introduce an effective turbulent α-viscosity with α= 0.01, which allows us particularly to model accretion in the low-mass regime where disk self-gravity is diminishing. A range of models with observationally-motivated values of the initial ratio of rotational to gravitational energy yields a correlation between mass accretion rate \dot{M} and M_st that is relatively steep, as observed. Additionally, our modeling reveals evidence for a bimodality in the \dot{M}--M_st correlation, with a steeper slope at lower masses and a shallower slope at intermediate and upper masses, as also implied by observations. We show that the neglect of disk self-gravity leads to a much steeper \dot{M}--M_st relation for intermediate and upper-mass T Tauri stars. This demonstrates that an accurate treatment of global self-gravity is essential to understanding observations of circumstellar disks.

preprint2009arXiv

Variable accretion in the embedded phase of stellar evolution

Motivated by the recent detection of a large number of embedded young stellar objects (YSOs) with mass accretion rates that are inconsistent with the predictions of the standard model of inside-out collapse (Shu 1977), we perform a series on numerical hydrodynamic simulations of the gravitational collapse of molecular cloud cores with various initial masses, rotation rates, and sizes. We focus on the early Class I stage of stellar evolution when circumstellar disks are exposed to high rates of mass deposition from infalling envelopes. Our numerical modeling reproduces the large observed spread in accretion rates inferred for embedded YSOs in Perseus, Serpens, and Ophiuchus star forming regions by Enoch et al. (2009), yielding 37%--75% of objects with "sub-Shu" accretion rates \dot{M} \la 10^{-6} Msun/yr and 1%--2% of objects with "super-Shu" accretion rates \dot{M}>10^{-5} Msun/yr. Mass accretion rates in the Class I stage have a log-normal distribution, with its shape controlled by disk viscosity and disk temperature. The spread in $\dot{M}$ is greater in models with lower viscosity and smaller in models with higher viscosity and higher disk temperature, suggesting that gravitational instability may be a dominant cause of the observed diversity in $\dot{M}$ in embedded YSOs. Our modeling predicts a weak dependence between the mean mass accretion rates and stellar masses in the Class I stage, in sharp contrast to the corresponding steep dependence for evolved T Tauri stars and brown dwarfs.

preprint2008arXiv

Disk masses in the embedded and T Tauri phases of stellar evolution

(Abridged). Motivated by a growing concern that masses of circumstellar disks may have been systematically underestimated by conventional observational methods, we present a numerical hydrodynamics study of time-averaged disk masses (<M_d>) around low-mass Class 0, Class I, and Class II objects. Mean disk masses (\overline{M}_d}) are then calculated by weighting the time-averaged disk masses according to the corresponding stellar masses using a power-law weight function with a slope typical for the Kroupa initial mass function of stars. Two distinct types of disks are considered: self-gravitating disks, in which mass and angular momentum are redistributed exclusively by gravitational torques, and viscous disks, in which both the gravitational and viscous torques are at work. We find that self-gravitating disks have mean masses that are slowly increasing along the sequence of stellar evolution phases. More specifically, Class 0/I/II self-gravitating disks have mean masses \overline{M}_d=0.09, 0.10, and 0.12 M_sun, respectively. Viscous disks have similar mean masses (\overline{M}_d=0.10-0.11 M_sun) in the Class 0/I phases but almost a factor of 2 lower mean mass in the Class II phase (\overline{M}_d,CII=0.06 M_sun). In each evolution phase, time-averaged disk masses show a large scatter around the mean value. Our obtained mean disk masses are larger than those recently derived by Andrews & Williams and Brown et al., regardless of the physical mechanisms of mass transport in the disk.

preprint2007arXiv

Propagating Star Formation in the Collisional Ring Galaxy Arp 10

Propagating star formation in a collisional ring galaxy Arp10 is investigated by a complex approach, which includes the broad- and narrow-band photometry, long- slit spectroscopy, and scanning Fabry-Perot spectroscopy. The ionized gas velocity field obtained with best spatial resolution to date indicates a non- isotropic expansion of the outer ring with a maximum velocity 110km/s. Strong vertical and non-circular motions are also seen in the vicinity of the inner ring. Our kinematic data suggest that Arp10 has a small inclination i=22\degr and high total (luminous plus dark matter) mass of about $10^{12} M_{\odot}$ within a 50 kpc radius.The abundance of oxygen 12 + log(O/H) in both star- forming rings is about 8.6.The analysis of spectral indices provides an estimate on the propagation velocities of both rings and metallicity of the pre-collision stellar population.A small "knot" near the nucleus of Arp10,which was previously suspected as a possible candidate for collision, is now unambiguously identified as the "intruder" of at least 1/4 of the total mass of Arp~10.We use a simplified two-dimensional hydrodynamic modeling of galaxy collisions to test a collisional origin of Arp10. We confirm that the sizes of the inner and outer rings, maximum expansion velocity of the outer ring, and radial profile of the gas circular velocity can be reproduced by a near-central collision with the intruder galaxy, which occurred approximately 85Myr ago. We acknowledge that an apparent crescent- shaped distribution of H$α$ emission in the outer ring is caused by a star formation threshold in the gas disk of Arp10.

preprint2003arXiv

Propagation velocities of gas rings in collisional ring galaxies

The propagation velocity of the first gas ring in collisional ring galaxies, i.e. the velocity at which the maximum in the radial gas density profile propagates radially in the galactic disk, is usually inferred from the radial expansion velocity of gas in the first ring. Our numerical hydrodynamics modeling of ring galaxy formation however shows that the maximum radial expansion velocity of gas in the first ring ($v_{gas}$) is invariably below the propagation velocity of the first gas ring itself ($v_{ring}$). Modeling of the Cartwheel galaxy indicates that the outer ring is currently propagating at $v_{ring} \approx$ 100 km/s, while the maximum radial expansion velocity of gas in the outer ring is currently $v_{gas} \approx$ 65 km/s. Modeling of the radial B-V/V-K color gradients of the Cartwheel ring galaxy also indicates that the outer ring is propagating at $v_{ring} \ge $ 90 km/s. We show that a combined effect of inclination, finite thickness, and warping of the Cartwheel's disk might be responsible for the lack of angular difference in the peak positions found for the azimuthally averaged $Hα$, K and B surface brightness profiles of the Cartwheel's outer ring. Indeed, the radial $Hα$ surface brightness profiles obtained along the Cartwheel's major axis, where effects of inclination and finite thickness are minimized, do peak exterior to those at K- and B-bands. The angular difference in peak positions implies $v_{ring}$ = 110 km/s, which is in agreement with the model predictions. We briefly discuss the utility of radio continuum emission and spectral line equivalent widths for determining the propagation velocity of gas rings in collisional ring galaxies.