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)

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.

The collisional evolution of solid material in protoplanetary disks is a crucial step in the formation of planetesimals, comets, and planets. Although dense protoplanetary environments favor fast dust coagulation, there are several factors that limit the straightforward pathway from interstellar micron-size grains to pebble-size aggregates. Apart from the grain bouncing, fragmentation, and fast drift to the central star, a notable limiting factor is the electrostatic repulsion of like-charged grains. In this study we aim at theoretical modeling of the dust coagulation coupled with the dust charging and disk ionization calculations. We show that the electrostatic barrier is a strong restraining factor to the coagulation of micrometer-size dust in dead zones of the disk (where the turbulence is suppressed). While the sustained turbulence helps to overcome the electrostatic barrier, low fractal dimensions of dust aggregates can potentially block their further coagulation even in this case. Coulomb repulsion may keep a significant fraction of small dust in the disk atmosphere and outer regions.

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

We present a numerical tool Shiva designed to simulate the dust destruction in warm neutral, warm ionized, and hot ionized media under the influence of photo-processing, sputtering, and shattering. The tool is designed primarily to study the evolution of hydrogenated amorphous carbons (HACs), but options to simulate polycyclic aromatic hydrocarbons (PAHs), silicate and graphite grains are also implemented. HAC grain photo-processing includes both dehydrogenation and carbon atom loss. Dehydrogenation leads to material transformation from aliphatic to aromatic structure. Simultaneously, some other physical properties (band gap energy, optical properties, etc.) of the material change as well. The Shiva tool allows calculating the time-dependent evolution of the dust size distribution depending on hydrogen, helium, and carbon number densities and ionization state, gas temperature, radiation flux, relative gas-dust and grain-grain velocities. For HAC grains the evolution of band gap energy distribution is also computed. We describe a dust evolution model, on which the tool relies, and present evolutionary time-scales for dust grains of different sizes depending on external conditions. This allows a user to estimate quickly a lifetime of a specific dust grain under relevant conditions. As an example of the tool usage, we demonstrate how grain properties and corresponding infrared spectra evolve in photo-dissociation regions, H II regions, and supernova remnant shocks.

The fact that infrared ring nebulae (IRRNs) are frequently associated with HII regions, provides us the opportunity to study dust at the interface between ionized and neutral gas. In this paper, we analyze the infrared radiation in the range from 8 to 500~$\mu$m in outer and inner parts of 32 IRRNs showing round shape. We aim to determine the morphology of these objects and possible dust evolution processes on the base of the comparison of IR radiation towards the ionized and neutral regions. We calculate six slopes between adjacent wavelengths in their spectral energy distributions to trace the difference in the physical conditions inside and outside ionized regions. Using the data on these 32 objects we show that their morphology is likely 3D spherical rather than 2D plane-like. The slope between 70 and 160~$\mu$m is the most appropriate tracer of the dust temperature in the outer envelope. The larger 8-to-24~$\mu$m intensity ratio is associated with smaller intensities at mid-IR indicating that the PAHs may indeed be generated due to larger grain destruction. These data are important for the subsequent theoretical modeling and determining the dust evolution in HII regions and their envelopes.

Dust is known to drift and grow in protoplanetary discs, which results in dust segregation over the disc extent. Maps of the spectral index $\alpha$ are a common tool for studying the dust content in protoplanetary discs. The analysis of observationally derived maps reveals significant gradients of the spectral index, confirming that dust evolves in the disc, but a more detailed information about the dust redistribution is required to make inferences about the early stages of dust growth. We calculated the spectral index maps based on the results of numerical hydrodynamical simulations using the FEOSAD code, which allows studying a long-term dynamics of a self-gravitating viscous disc populated with coagulating, drifting, and fragmenting dust. Here we demonstrate that values of the spectral index estimated for different wavelength intervals within the far-infrared and radio bands reveal the presence of dust grains of various sizes. Specifically, we show that the disc regions with the maximal spectral index in a specific wavelength interval are the regions with the prevalence of dust grains of a specific size. Thus, a set of spectral index maps derived using different wavelength intervals can be used to recover the dust size-distribution over the disc extent.

The chemical influence of luminosity outbursts on the environments of young solar-type stars is explored. Species are categorised into several types according to their response to the outburst. The first and second types imply chemical changes only during the outburst (with slightly different behaviours). These response types are mostly observed close to the star and are caused by icy mantle evaporation. However, mantles recover after the outburst almost immediately. A notable exception is benzene ice, which is accumulated on dust surfaces during and after the outburst, so that its abundance exceeds the pre-outburst level by orders of magnitude. The third type of response is mostly seen at the disc periphery and implies alteration of abundances during the outburst and preservation of these `abnormal' abundances for centuries. This behaviour is typical of organic compounds, like HCOOCH$_3$, CH$_3$CN, CH$_2$CO. Their presence in the dark disc regions can be a manifestation of the past outburst. CO and CO$_2$ only trace past outbursts at the remote disc regions. The outburst changes the C/O ratio, but it quickly returns to the pre-outburst value almost everywhere in the disc. An important factor determining the sensitivity of molecular composition to the outburst is the dust size distribution. The duration of the pre-outburst stage and of the outburst itself influence the chemical effects, if the burst duration is shorter than 50 yr and the duration of the quiescent phase between the bursts is shorter than 100 kyr.

Circumstellar discs, from which planetary systems are formed, consist of gas, dust and solids. Simulations of self-consistent dynamics of gas, dust and solids in circumstellar discs is a challenging problem. In the paper we present fast algorithms for computing the drag force (momentum transfer) between solid phase and gas. These algorithms (a) are universal and applicable to dust and solids with any sizes smaller than the mean free path of gas molecules, (b) can be used to calculate the momentum transfer between dust and gas instead of one-way effect, as it is done in many models, (c) can perform simulations, without a loss in accuracy, with the time step determined by gas-dynamic parameters rather than by drag force, and (d) are compatible with the widely used parallel algorithms for solving 3D equations of gas dynamics, hydrodynamic equations for dust, and the collisionless Boltzmann equation for large bodies. Preliminary results of supercomputer simulation of the gas-dust disc dynamics within the developed approach are reported.

The FU Ori-type young stellar objects are characterized by a sudden increase in luminosity by 1$-$2 orders of magnitude, followed by slow return to the pre-outburst state on timescales of $\sim$10$-$100 yr. The outburst strongly affects the entire disk, changing its thermal structure and radiation field. In this paper, using a detailed physical-chemical model we study the impact of the FU Ori outburst on the disk chemical inventory. Our main goal is to identify gas-phase molecular tracers of the outburst activity that could be observed after the outburst with modern telescopes such as ALMA and NOEMA. We find that the majority of molecules experience a considerable increase in the total disk gas-phase abundances due to the outburst, mainly due to the sublimation of their ices. Their return to the pre-outburst chemical state takes different amounts of time, from nearly instantaneous to very long. Among the former ones we identify CO, NH$_3$, C$_2$H$_6$, C$_3$H$_4$, etc. Their abundance evolution tightly follows the luminosity curve. For CO the abundance increase does not exceed an order of magnitude, while for other tracers the abundances increase by 2$-$5 orders of magnitude. Other molecules like H$_2$CO and NH$_2$OH have longer retention timescales, remaining in the gas phase for $\sim10-10^3$ yr after the end of the outburst. Thus H$_2$CO could be used as an indicator of the previous outbursts in the post-outburst FU Ori systems. We investigate the corresponding time-dependent chemistry in detail and present the most favorable transitions and ALMA configurations for future observations.

The long-term evolution of a circumstellar disk starting from its formation and ending in the T Tauri phase was simulated numerically with the purpose of studying the evolution of dust in the disk with distinct values of viscous \alpha-parameter and dust fragmentation velocity v_frag. We solved numerical hydrodynamics equations in the thin-disk limit, which are modified to include a dust component consisting of two parts: sub-micron-sized dust and grown dust with a maximum radius a_r. The former is strictly coupled to the gas, while the latter interacts with the gas via friction. The conversion of small to grown dust, dust growth, and dust self-gravity are also considered. We found that the process of dust growth known for the older protoplanetary phase also holds for the embedded phase of disk evolution. The dust growth efficiency depends on the radial distance from the star - a_r is largest in the inner disk and gradually declines with radial distance. In the inner disk, a_r is limited by the dust fragmentation barrier. The process of small-to-grown dust conversion is very fast once the disk is formed. The total mass of grown dust in the disk (beyond 1 AU) reaches tens or even hundreds of Earth masses already in the embedded phase of star formation and even a greater amount of grown dust drifts in the inner, unresolved 1 AU of the disk. Dust does not usually grow to radii greater than a few cm. A notable exception are models with \alpha <= 10^{-3}, in which case a zone with reduced mass transport develops in the inner disk and dust can grow to meter-sized boulders in the inner 10 AU. Grown dust drifts inward and accumulates in the inner disk regions. This effect is most pronounced in the \alpha <= 10^{-3} models where several hundreds of Earth masses can be accumulated in a narrow region of several AU from the star by the end of embedded phase. (abridged).

Protoplanetary disk mass is a key parameter controlling the process of planetary system formation. CO molecular emission is often used as a tracer of gas mass in the disk. In this study we consider the ability of CO to trace the gas mass over a wide range of disk structural parameters and search for chemical species that could possibly be used as alternative mass tracers to CO. Specifically, we apply detailed astrochemical modeling to a large set of models of protoplanetary disks around low-mass stars, to select molecules with abundances correlated with the disk mass and being relatively insensitive to other disk properties. We do not consider sophisticated dust evolution models, restricting ourselves with the standard astrochemical assumption of $0.1~\mu $m dust. We find that CO is indeed the best molecular tracer for total gas mass, despite the fact that it is not the main carbon carrier, provided reasonable assumptions about CO abundance in the disk are used. Typically, chemical reprocessing lowers the abundance of CO by a factor of 3, compared to the case of photo-dissociation and freeze-out as the only ways of CO depletion. On average only 13% C-atoms reside in gas-phase CO, albeit with variations from 2 to 30%. CO$_2$, H$_2$O and H$_2$CO can potentially serve as alternative mass tracers, the latter two being only applicable if disk structural parameters are known.

In this paper, we extend the study initiated in PaperI by modelling grain ensemble evolution in a dynamical model of an expanding HII region and checking the effects of momentum transfer from dust to gas. The radiation pressure on the dust, the dust drift, and the lug on the gas by the dust are all important process that should be considered simultaneously to describe the dynamics of HII regions. With accounting for the momentum transfer from the dust to the gas, the expansion time of the HII region is notably reduced (for our model of RCW120, the time to reach the observed radius of the HII region is reduced by a factor of 1.5). Under the common approximation of frozen dust, where there is no relative drift between the dust and gas, the radiation pressure from the ionizing star drives the formation of the very deep gas cavity near the star. Such a cavity is much less pronounced when the dust drift is taken into account. The dust drift leads to the two-peak morphology of the dust density distribution and significantly reduces the dust-to-gas ratio in the ionized region (by a factor of 2 to 10). The dust-to-gas ratio is larger for higher temperatures of the ionizing star since the dust grains have a larger electric charge and are more strongly coupled to the gas.

Ionization-recombination balance in dense interstellar and circumstellar environments is a key factor for a variety of important physical processes, such as chemical reactions, dust charging and coagulation, coupling of the gas with magnetic field and development of instabilities in protoplanetary disks. We determine a critical gas density above which the recombination of electrons and ions on the grain surface dominates over the gas-phase recombination. For this regime, we present a self-consistent analytical model which allows us to exactly calculate abundances of charged species in dusty gas, without making assumptions on the grain charge distribution. To demonstrate the importance of the proposed approach, we check whether the conventional approximation of low grain charges is valid for typical protoplanetary disks, and discuss the implications for dust coagulation and development of the "dead zone" in the disk. The presented model is applicable for arbitrary grain-size distributions and, for given dust properties and conditions of the disk, has only one free parameter - the effective mass of the ions, shown to have a low effect on results. The model can be easily included in numerical simulations following the dust evolution in dense molecular clouds and protoplanetary disks.

The coagulation of dust particles under the conditions in protoplanetary disks is investigated. The study focuses on the repulsive electrostatic barrier against growth of charged dust grains. Taking into account the photoelectric effect leads to the appearance of a layer at intermediate heights where the dust has a close to zero charge, enabling the dust grains to grow efficiently. An increase in the coagulation rate comes about not only due to the lowering of the Coulomb barrier, but also because of the electrostatic attraction between grains of opposite charge due to the non-zero dispersion of the near-zero charge. Depending on the efficiency of mixing in the disk, the acceleration of the evolution of the dust in this layer could be important, both in the quasi-stationary stage of the disk evolution and during its dispersal.

We consider dust drift under the influence of stellar radiation pressure during the pressure-driven expansion of an HII region using the chemo-dynamical model MARION. Dust size distribution is represented by four dust types: conventional polycyclic aromatic hydrocarbons (PAHs), very small grains (VSGs), big grains (BGs) and also intermediate-sized grains (ISGs), which are larger than VSGs and smaller than BGs. The dust is assumed to move at terminal velocity determined locally from the balance between the radiation pressure and gas drag. As Coulomb drag is an important contribution to the overall gas drag, we evaluate a grain charge evolution within the HII region for each dust type. BGs are effectively swept out of the HII region. The spatial distribution of ISGs within the HII region has a double peak structure, with a smaller inner peak and a higher outer peak. PAHs and VSGs are mostly coupled to the gas. The mean charge of PAHs is close to zero, so they can become neutral from time to time because of charge fluctuations. These periods of neutrality occur often enough to cause the removal of PAHs from the very interior of the HII region. For VSGs, the effect of charge fluctuations is less pronounced but still significant. We conclude that accounting for charge dispersion is necessary to describe the dynamics of small grains.

Observations of the circumstellar disk in the Bok globule CB 26 at 110, 230, and 270 GHz are presented together with the results of the simulations and estimates of the disk parameters. These observations were obtained using the SMA, IRAM Plateau de Bure, and OVRO interferometers. The maps have relatively high angular resolutions (0.4-1"), making it possible to study the spatial structure of the gas-dust disk. The disk parameters are reconstructed via a quantitative comparison of observational and theoretical intensity maps. The disk model used to construct the theoretical maps is based on the assumption of hydrostatic and radiative equilibrium in the vertical direction, while the radial surface density profile is described phenomenologically. The system of equations for the transfer of the infrared and ultraviolet radiation is solved in the vertical direction, in order to compute the thermal structure of the disk. The disk best-fit parameters are derived for each map and all the maps simultaneously, using a conjugate gradient method. The degrees of degeneracy of the parameters describing the thermal structure and density distribution of the disk are analyzed in detail. All three maps indicate the presence of an inner dust-free region with a radius of approximately 35 AU, in agreement with the conclusions of other studies. The inclination of the disk is 78 deg, which is smaller than the value adopted in our earlier study of rotating molecular outflows from CB 26. The model does not provide any evidence for the growth of dust particles above a_max=0.02 cm.

We present a self-consistent model of a protoplanetary disk: 'ANDES' ('AccretioN disk with Dust Evolution and Sedimentation'). ANDES is based on a flexible and extendable modular structure that includes 1) a 1+1D frequency-dependent continuum radiative transfer module, 2) a module to calculate the chemical evolution using an extended gas-grain network with UV/X-ray-driven processes surface reactions, 3) a module to calculate the gas thermal energy balance, and 4) a 1+1D module that simulates dust grain evolution. For the first time, grain evolution and time-dependent molecular chemistry are included in a protoplanetary disk model. We find that grain growth and sedimentation of large grains to the disk midplane lead to a dust-depleted atmosphere. Consequently, dust and gas temperatures become higher in the inner disk (R < 50 AU) and lower in the outer disk (R > 50 AU), in comparison with the disk model with pristine dust. The response of disk chemical structure to the dust growth and sedimentation is twofold. First, due to higher transparency a partly UV-shielded molecular layer is shifted closer to the dense midplane. Second, the presence of big grains in the disk midplane delays the freeze-out of volatile gas-phase species such as CO there, while in adjacent upper layers the depletion is still effective. Molecular concentrations and thus column densities of many species are enhanced in the disk model with dust evolution, e.g., CO2, NH2CN, HNO, H2O, HCOOH, HCN, CO. We also show that time-dependent chemistry is important for a proper description of gas thermal balance.

Stochastic heating of small grains is often mentioned as a primary cause of large infrared (IR) fluxes from star-forming galaxies, e.g. at 24\mu m. If the mechanism does work at a galaxy-wide scale, it should show up at smaller scales as well. We calculate temperature probability density distributions within a model protostellar core for four dust components: large silicate and graphite grains, small graphite grains, and polycyclic aromatic hydrocarbon particles. The corresponding spectral energy distributions are calculated and compared with observations of a representative infrared dark cloud core. We show that stochastic heating, induced by the standard interstellar radiation field, cannot explain high mid-IR emission toward the centre of the core. In order to reproduce the observed emission from the core projected centre, in particular, at 24\mu m, we need to increase the ambient radiation field by a factor of about 70. However, the model with enhanced radiation field predicts even higher intensities at the core periphery, giving it a ring-like appearance, that is not observed. We discuss possible implications of this finding and also discuss a role of other non-radiative dust heating processes.

We study details of the UV radiation transfer in a protoplanetary disk, paying attention to the influence of dust growth and sedimentation on the disk density and temperature. Also, we show how the dust evolution affects photoreaction rates of key molecules, like CN and CS.