Source author record

Eduard Vorobyov

Eduard 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

14works
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

14 published item(s)

preprint2026arXiv

Solar-like oscillations in accreting pre-main sequence stars: Insights and prospects

We present theoretical predictions for solar-like oscillators in the pre-main sequence phase of stellar evolution. Our pre-main sequence models start from a stellar seed of 0.01 solar masses that gains mass through accretion, offering an alternative description to the classical approach segmented into the Hayashi and Henyey tracks. Evolutionary models are calculated using the \mesa stellar evolution code with a custom accretion routine and pulsation properties are investigated using the \GYRE oscillation code. We present evolutionary tracks and internal structures for accreting pre-main sequence solar-like stars in the mass range from 0.7 to 1.6 solar masses, adopting 35 mass accretion histories previously extracted from two-dimensional magneto-hydrodynamical simulations. Atmospheric parameters of our models constrain characteristic frequencies of pre-main sequence solar-like oscillators to be generally greater than 500 $μ\text{Hz}$. We highlight the imprint of accretion on the buoyancy and Lamb profiles and illustrate the effects on the small- and large frequency separations. We additionally quantify individual frequency differences across the 35 accretion histories at the zero-age main-sequence, showcasing differences no larger than $20 \ μ\text{Hz}$ for an exemplary model. Finally, we discuss the potential of detecting solar-like oscillations in pre-main sequence stars with the upcoming ESA PLATO mission.

preprint2024arXiv

Dust growth and pebble formation in the initial stages of protoplanetary disk evolution

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.

preprint2022arXiv

Primordial Dusty Rings and Episodic Outbursts in Protoplanetary Discs

We investigate the formation and evolution of "primordial" dusty rings occurring in the inner regions of protoplanetary discs, with the help of long-term, coupled dust-gas, magnetohydrodynamic simulations. The simulations are global and start from the collapse phase of the parent cloud core, while the dead zone is calculated via an adaptive $α$ formulation by taking into account the local ionization balance. The evolution of the dusty component includes its growth and back reaction on to the gas. Previously, using simulations with only a gas component, we showed that dynamical rings form at the inner edge of the dead zone. We find that when dust evolution as well as magnetic field evolution in the flux-freezing limit are included, the dusty rings formed are more numerous and span a larger radial extent in the inner disc, while the dead zone is more robust and persists for a much longer time. We show that these dynamical rings concentrate enough dust mass to become streaming unstable, which should result in rapid planetesimal formation even in the embedded phases of the system. The episodic outbursts caused by the magnetorotational instability have significant impact on the evolution of the rings. The outbursts drain the inner disc of grown dust, however, the period between bursts is sufficiently long for the planetesimal growth via streaming instability.The dust mass contained within the rings is large enough to ultimately produce planetary systems with the core accretion scenario. The low mass systems rarely undergo outbursts and thus, the conditions around such stars can be especially conducive for planet formation.

preprint2021arXiv

Eruptive Behavior of Magnetically Layered Protoplanetary Disks in Low Metallicity Environments

A protoplanetary disk typically forms a dead zone near its midplane at the distance of a few au from the central protostar. Accretion through such a magnetically layered disk can be intrinsically unstable and has been associated with episodic outbursts in young stellar objects. We present the first investigation into the effects of low metallicity environment on the structure of the dead zone as well as the resulting outbursting behavior of the protoplanetary disk. We conducted global numerical hydrodynamic simulations of protoplanetary disk formation and evolution in the thin-disk limit. The consequences of metallicity were considered via its effects on the gas and dust opacity of the disk, the thickness of the magnetically active surface layer, and the temperature of the prestellar cloud core. We show that the metal poor disks accumulate much more mass in the innermost regions, as compared to the solar metallicity counterparts. The duration of the outbursting phase also varies with metallicity - the low metallicity disks showed more powerful luminosity eruptions with a shorter burst phase, which was confined mostly to the early, embedded stages of the disk evolution. The lowest metallicity disks with the higher cloud core temperature showed the most significant differences. The occurrence of outbursts was relatively rare in the disks around low mass stars and this was especially true at lowest metallicities. We conclude that the metal content of the disk environment can have profound effects on both the disk structure and evolution in terms of episodic accretion.

preprint2021arXiv

Evolution of dust in protoplanetary disks of eruptive stars

Luminosity bursts in young FU Orionis-type stars warm up the surrounding disks of gas and dust, thus inflicting changes on their morphological and chemical composition. In this work, we aim at studying the effects that such bursts may have on the spatial distribution of dust grain sizes and the corresponding spectral index in protoplanetary disks. We use the numerical hydrodynamics code FEOSAD, which simulates the co-evolution of gas, dust, and volatiles in a protoplanetary disk, taking dust growth and back reaction on gas into account. The dependence of the maximum dust size on the water ice mantles is explicitly considered. The burst is initialized by increasing the luminosity of the central star to 100-300 L_sun for a time period of 100 yr. The water snowline shifts during the burst to a larger distance, resulting in the drop of the maximum dust size interior to the snowline position because of more efficient fragmentation of bare grains. After the burst, the water snowline shifts quickly back to its preburst location followed by renewed dust growth. The timescale of dust regrowth after the burst depends on the radial distance so that the dust grains at smaller distances reach the preburst values faster than the dust grains at larger distances. As a result, a broad peak in the radial distribution of the spectral index in the millimeter dust emission develops at \approx 10 au, which shifts further out as the disk evolves and dust grains regrow to preburst values at progressively larger distances. This feature is most pronounced in evolved axisymmetric disks rather than in young gravitationally unstable counterparts, although young disks may still be good candidates if gravitational instability is suppressed. Abridged.

preprint2020arXiv

Chemical modelling of FU Ori protoplanetary disks

Luminosity outbursts of the FUOri type stars, which have a magnitude of $\sim100 L_{\odot}$ and last for decades, may affect chemical composition of the surrounding protoplanetary disk. Using astrochemical modeling we analyze the changes induced by the outburst and search for species sensitive to the luminosity rise. Some changes in the disk molecular composition appear not only during the outburst itself but can also retain for decades after the end of the outburst. We analyze main chemical processes responsible for these effects and assess timescales at which chemically inert species return to the pre-outburst abundances.

preprint2020arXiv

Global Protoplanetary Disk Simulations: Dead Zone Formation and FUor Outbursts

We conducted global hydrodynamic simulations of protoplanetary disk evolution with an adaptive Shakura-Sunyaev α prescription to represent the layered disk structure, and starting with the collapse phase of the molecular cloud. With the canonical values of model parameters, self-consistent dead zones formed at the scale of a few au. The instabilities associated with the dead zone and corresponding outbursts, similar to FUor eruptions, were also observed in the simulations.

preprint2020arXiv

Outbursts in Global Protoplanetary Disk Simulations

While accreting through a circumstellar disk, young stellar objects are observed to undergo sudden and powerful accretion events known as FUor or EXor outbursts. Although such episodic accretion is considered to be an integral part of the star formation process, the triggers and mechanisms behind them remain uncertain. We conducted global numerical hydrodynamics simulations of protoplanetary disk formation and evolution in the thin-disk limit, assuming both magnetically layered and fully magnetorotational instability (MRI)-active disk structure. In this paper, we characterize the nature of the outbursts occurring in these simulations. The instability in the dead zone of a typical layered disk results in "MRI outbursts". We explore their progression and their dependence on the layered disk parameters as well as cloud core mass. The simulations of fully MRI-active disks showed an instability analogous to the classical thermal instability. This instability manifested at two temperatures--above approximately 1400 K and 3500 K--due to the steep dependence of Rosseland opacity on the temperature. The origin of these thermally unstable regions is related to the bump in opacity resulting from molecular absorption by water vapor and may be viewed as a novel mechanism behind some of the shorter duration accretion events. Although we demonstrated local thermal instability in the disk, more investigations are needed to confirm that a large-scale global instability will ensue. We conclude that the magnetic structure of a disk, its composition, as well as the stellar mass, can significantly affect the nature of episodic accretion in young stellar objects.

preprint2016arXiv

Signatures of Gravitational Instability in Resolved Images of Protostellar Disks

Protostellar (class 0/I) disks, having masses comparable to those of their nascent host stars, and fed continuously from their natal infalling envelopes, are prone to gravitational instability (GI). Motivated by advances in near-infrared (NIR) adaptive optics imaging and mm-wave interferometry, we explore the observational signatures of GI in disks, using hydrodynamical and Monte Carlo radiative transfer simulations to synthesize NIR scattered light images and mm dust continuum maps. Spiral arms induced by GI, located at disk radii of hundreds of AUs, are local overdensities and have their photospheres displaced to higher altitudes above the disk midplane, arms therefore scatter more NIR light from their central stars than inter-arm regions, and are detectable at distances up to 1 kpc by Gemini/GPI, VLT/SPHERE, and Subaru/HiCIAO/SCExAO. By contrast, collapsed clumps formed by disk fragmentation have such strong local gravitational fields that their scattering photospheres are at lower altitudes, such fragments appear fainter than their surroundings in NIR scattered light. Spiral arms and streamers recently imaged in four FU Ori systems at NIR wavelengths resemble GI-induced structures and support the interpretation that FUors are gravitationally unstable protostellar disks. At mm wavelengths, both spirals and clumps appear brighter in thermal emission than the ambient disk and can be detected by ALMA at distances up to 0.4 kpc with one-hour integration times at ~0.1 arcsec resolution. Collapsed fragments having masses >=1 MJ can be detected by ALMA within ~10 minutes.

preprint2016arXiv

Using the PPML approach for constructing a low-dissipation, operator-splitting scheme for numerical simulations of hydrodynamic flows

An approach for constructing a low-dissipation numerical method is described. The method is based on a combination of the operator-splitting method, Godunov method, and piecewise-parabolic method on the local stencil. Numerical method was tested on a standard suite of hydrodynamic test problems. In addition, the performance of the method is demonstrated on a global test problem showing the development of a spiral structure in a gravitationally unstable gaseous galactic disk.

preprint2013arXiv

Fragmenting protostellar disks: properties and observational signatures

Using numerical hydrodynamic simulations, we study the gravitational fragmentation of an unstable protostellar disc formed during the collapse of a pre-stellar core with a mass of 1.2 M_sun. The forming fragments span a mass range from about a Jupiter mass to very-low-mass protostars and are located at distances from a few tens to a thousand AU, with a dearth of objects at < 100 AU. We explore the possibility of observational detection of the fragments in discs viewed through the outflow cavity at a distance of 250 pc. We demonstrate that one hour of integration time with the Atacama Large Millimeter/sub-millimeter Array (ALMA) is sufficient to detect the fragments with masses as low as 1.5 M_Jup at orbital distances up to 800 AU from the protostar. The ALMA resolution sets the limit on the minimum orbital distance of detectable fragments. For the adopted resolution of our simulated ALMA images of 0.1", the fragments can be detected at distances down to 50 AU. At smaller distances, the fragments usually merge with the central density peak. The likelihood for detecting the fragments reduces significantly for a lower resolution of 0.5". Some of the most massive fragments, regardless of their orbital distance, can produce characteristic peaks at approximately 5 micron and hence their presence can be indirectly inferred from the observed spectral energy distributions of protostars.

preprint2012arXiv

Observed luminosity spread in young clusters and Fu Ori stars: a unified picture

The idea that non steady accretion during the embedded phase of protostar evolution can produce the observed luminosity spread in the Herzsprung-Russell diagram (HRD) of young clusters has recently been called into question. Observations of Fu Ori, for instance, suggest an expansion of the star during strong accretion events whereas the luminosity spread implies a contraction of the accreting objects, decreasing their radiating surface. In this paper, we present a global scenario based on calculations coupling episodic accretion histories derived from numerical simulations of collapsing cloud prestellar cores of various masses and subsequent protostar evolution. Our calculations show that, assuming an initial protostar mass $\mi \sim 1\,\mjup$, typical of the second Larson's core, both the luminosity spread in the HRD and the inferred properties of Fu Ori events (mass, radius, accretion rate) can be explained by this scenario, providing two conditions. First, there must be some variation within the fraction of accretion energy absorbed by the protostar during the accretion process. Second the range of this variation should increase with increasing accretion burst intensity, and thus with the initial core mass and final star mass. The numerical hydrodynamics simulations of collapsing cloud prestellar cores indeed show that the intensity of the accretion bursts correlates with the mass and initial angular momentum of the prestellar core. Massive prestellar cores with high initial angular momentum are found to produce intense bursts characteristic of Fu Ori like events. Our results thus suggest a link between the burst intensities and the fraction of accretion energy absorbed by the protostar, with some threshold in the accretion rate, of the order of $10^{-5}\msolyr$, delimitating the transition from "cold" to "hot" accretion. [Abridged]

preprint2010arXiv

Lifetime of the embedded phase of low-mass star formation and the envelope depletion rates

Motivated by a considerable scatter in the observationally inferred lifetimes of the embedded phase of star formation, we study the duration of the Class 0 and Class I phases in upper-mass brown dwarfs and low-mass stars using numerical hydrodynamics simulations of the gravitational collapse of a large sample of cloud cores. We resolve the formation of a star/disk/envelope system and extend our numerical simulations to the late accretion phase when the envelope is nearly totally depleted of matter. We adopted a classification scheme of Andre et al. and calculate the lifetimes of the Class 0 and Class I phases (τ_C0 and τ_CI, respectively) based on the mass remaining in the envelope. When cloud cores with various rotation rates, masses, and sizes (but identical otherwise) are considered, our modeling reveals a sub-linear correlation between the Class 0 lifetimes and stellar masses in the Class 0 phase with the least-squares fit exponent m=0.8 \pm 0.05. The corresponding correlation between the Class I lifetimes and stellar masses in the Class I is super-linear with m=1.2 \pm 0.05. If a wider sample of cloud cores is considered, which includes possible variations in the initial gas temperature, cloud core truncation radii, density enhancement amplitudes, initial gas density and angular velocity profiles, and magnetic fields, then the corresponding exponents may decrease by as much as 0.3. The duration of the Class I phase is found to be longer than that of the Class~0 phase in most models, with a mean ratio τ_CI / τ_C0 \approx 1.5--2. A notable exception are YSOs that form from cloud cores with large initial density enhancements, in which case τ_C0 may be greater than τ_CI. Moreover, the upper-mass (>= 1.0 Msun) cloud cores with frozen-in magnetic fields and high cloud core rotation rates may have the τ_CI / τ_C0 ratios as large as 3.0--4.0. (Abdridged).

preprint2008arXiv

Stellar Velocity Distribution in Galactic Disks

We present numerical studies of the properties of the stellar velocity distribution in galactic disks which have developed a saturated, two-armed spiral structure. In previous papers we used the Boltzmann moment equations (BME) up to second order for our studies of the velocity structure in self-gravitating stellar disks. A key assumption of our BME approach is the zero-heat flux approximation, i.e. the neglection of third order velocity terms. We tested this assumption by performing test particle simulations for stars in a disk galaxy subject to a rotating spiral perturbation. As a result we corroborated qualitatively the complex velocity structure found in the BME approach. It turned out that an equilibrium configuration in velocity space is only slowly established on a typical timescale of 5 Gyrs or more. Since many dynamical processes in galaxies (like the growth of spirals or bars)act on shorter timescales, pure equilibrium models might not be fully appropriate for a detailed comparison with observations like the local Galactic velocity distribution. Third order velocity moments were typically small and uncorrelated over almost all of the disk with the exception of the 4:1 resonance region (UHR). Near the UHR (normalized) fourth and fifth order velocity moments are still of the same order as the second and third order terms. Thus, at the UHR higher order terms are not negligible.