-
A resolved rotating disk wind from a young T Tauri star in the Bok globule CB26
Authors:
R. Launhardt,
Ya. N. Pavlyuchenkov,
V. V. Akimkin,
A. Dutrey,
F. Gueth,
S. Guilloteau,
Th. Henning,
V. Pietu,
K. Schreyer,
D. Semenov,
B. Stecklum,
T. L. Bourke
Abstract:
The disk-outflow connection plays a key role in extracting excess angular momentum from a forming protostar. We have previously reported the discovery of a small molecular outflow from the edge-on T Tauri star in the Bok globule CB26 that shows a peculiar velocity pattern, reminiscent of an outflow that corotates with the disk. We report new, high-resolution mm-interferometric observations of CB26…
▽ More
The disk-outflow connection plays a key role in extracting excess angular momentum from a forming protostar. We have previously reported the discovery of a small molecular outflow from the edge-on T Tauri star in the Bok globule CB26 that shows a peculiar velocity pattern, reminiscent of an outflow that corotates with the disk. We report new, high-resolution mm-interferometric observations of CB26 with the aim of revealing the morphology and kinematics of the outflow at the disk-outflow interface. The IRAM PdBI was used to observe CO(2-1) at 1.3mm with a resolution of 0.5". Using a physical model of the disk, which was derived from the dust emission, we employed chemo-dynamical modeling combined with line radiative transfer to constrain kinematic parameters and to construct a model of the CO emission from the disk that allowed us to separate the emission of the disk from that of the outflow. Our observations confirm the disk-wind nature of the rotating molecular outflow from CB26. The new high-resolution data reveal an X-shaped morphology of the CO emission close to the disk, and vertical streaks extending from the disk surface with a small half-opening angle of ~7deg, which can be traced out to vertical heights of ~500au. We interpret this emission as the combination of the disk atmosphere and a well-collimated disk wind, which we trace down to vertical heights of 40au, where it is launched from the surface of the flared disk at radii of 20-45au. The observed CO outflow has a total momentum flux of 1e-5 Msun km/s/yr, which is nearly three orders of magnitude larger than the maximum thrust that can be provided by the luminosity of the central star. We conclude that photoevaporation cannot be the main driving mechanism for this outflow, but it must be predominantly an MHD disk wind. It is thus far the best-resolved rotating disk wind observed to be launched from a circumstellar disk.
△ Less
Submitted 4 September, 2023;
originally announced September 2023.
-
Effect of Dust Evaporation and Thermal Instability on Temperature Distribution in a Protoplanetary Disk
Authors:
Ya. N. Pavlyuchenkov,
V. V. Akimkin,
A. P. Topchieva,
E. I. Vorobyov
Abstract:
The thermal instability of accretion disks is widely used to explain the activity of cataclysmic variables, but its development in protoplanetary disks has been studied in less detail. We present a semi-analytical stationary model for calculating the midplane temperature of a gas and dust disk around a young star. The model takes into account gas and dust opacities, as well as the evaporation of d…
▽ More
The thermal instability of accretion disks is widely used to explain the activity of cataclysmic variables, but its development in protoplanetary disks has been studied in less detail. We present a semi-analytical stationary model for calculating the midplane temperature of a gas and dust disk around a young star. The model takes into account gas and dust opacities, as well as the evaporation of dust at temperatures above 1000 K. Using this model, we calculate the midplane temperature distributions of the disk under various assumptions about the source of opacity and the presence of dust. We show that when all considered processes are taken into account, the heat balance equation in the region r<1 au has multiple temperature solutions. Thus, the conditions for thermal instability are met in this region. To illustrate the possible influence of instability on the accretion state in a protoplanetary disk, we consider a viscous disk model with alpha parameterization of turbulent viscosity. We show that in such a model the disk evolution is non-stationary, with alternating phases of accumulation of matter in the inner disk and its rapid accretion onto the star, leading to an episodic accretion pattern. These results indicate that this instability needs to be taken into account in evolutionary models of protoplanetary disks.
△ Less
Submitted 28 July, 2023;
originally announced July 2023.
-
Turbulent convection in protoplanetary discs and its role in angular momentum transfer
Authors:
E. P. Kurbatov,
Ya. N. Pavlyuchenkov
Abstract:
We present a model for the transport of anisotropic turbulence in an accretion disc. The model uses the Reynolds stress tensor approach in the mean field approximation. To study the role of convection in a protoplanetary disc, we combine the turbulence model with a radiative transfer calculation, and also include convection using the mixing length approximation. We find that the turbulence generat…
▽ More
We present a model for the transport of anisotropic turbulence in an accretion disc. The model uses the Reynolds stress tensor approach in the mean field approximation. To study the role of convection in a protoplanetary disc, we combine the turbulence model with a radiative transfer calculation, and also include convection using the mixing length approximation. We find that the turbulence generated by convection causes the angular momentum of the accretion disc to be directed outwards. We also confirm the conclusions of other authors that turbulent convection is unable to provide the observed disc accretion rates as well as a heat source sufficient for the convection to be self-sustaining. The reasons for the latter are the strong anisotropy of the turbulence together with the low efficiency of the energy transfer from the background velocity shear to the turbulent stress tensor.
△ Less
Submitted 2 November, 2023; v1 submitted 23 January, 2023;
originally announced January 2023.
-
The shocked molecular layer in RCW 120
Authors:
M. S. Kirsanova,
Ya. N. Pavlyuchenkov,
A. O. H. Olofsson,
D. A. Semenov,
A. F. Punanova
Abstract:
Expansion of wind-blown bubbles or HII regions lead to formation of shocks in the interstellar medium, which compress surrounding gas into dense layers. We made spatially and velocity-resolved observations of the RCW~120 PDR and nearby molecular gas with CO(6-5) and 13CO(6-5) lines and distinguished a bright CO-emitting layer, which we related with the dense shocked molecular gas moving away from…
▽ More
Expansion of wind-blown bubbles or HII regions lead to formation of shocks in the interstellar medium, which compress surrounding gas into dense layers. We made spatially and velocity-resolved observations of the RCW~120 PDR and nearby molecular gas with CO(6-5) and 13CO(6-5) lines and distinguished a bright CO-emitting layer, which we related with the dense shocked molecular gas moving away from the ionizing star due to expansion of HII region. Simulating gas density and temperature, as well as brightness of several CO and C+ emission lines from the PDR, we found reasonable agreement with the observed values. Analysing gas kinematics, we revealed the large-scale shocked PDR and also several dense environments of embedded protostars and outflows. We observe the shocked layer as the most regular structure in the CO(6-5) map and in the velocity space, when the gas around YSOs is dispersed by the outflows.
△ Less
Submitted 16 December, 2022;
originally announced December 2022.
-
Simulation of Thermal Surface Waves in a Protoplanetary Disk in a Two-Dimensional Approximation
Authors:
Ya. N. Pavlyuchenkov,
L. A. Maksimova,
V. V. Akimkin
Abstract:
Theoretical models predict that the obscuration of stellar radiation by irregularities on the surface of a protoplanetary disk can cause self-generating waves traveling towards the star. However, this process is traditionally simulated using the 1+1D approach, the key approximations of which - vertical hydrostatic equilibrium of the disk and vertical diffusion of IR radiation - can distort the pic…
▽ More
Theoretical models predict that the obscuration of stellar radiation by irregularities on the surface of a protoplanetary disk can cause self-generating waves traveling towards the star. However, this process is traditionally simulated using the 1+1D approach, the key approximations of which - vertical hydrostatic equilibrium of the disk and vertical diffusion of IR radiation - can distort the picture. This article presents a two-dimensional radiative hydrodynamic model of the evolution of an axially symmetric gas and dust disk. Within this model, but using simplified assumptions from 1+1D models, we have reproduced the spontaneous generation and propagation of thermal surface waves. The key conclusion of our work is that taking into account two-dimensional hydrodynamics and diffusion of IR radiation suppresses the spontaneous generation and development of thermal waves observed in the 1+1D approximation. The search for the possibility of the existence of surface thermal waves should be continued by studying the problem for various parameters of protoplanetary disks.
△ Less
Submitted 9 November, 2022;
originally announced November 2022.
-
Simulations of thermal surface waves in a protoplanetary disk using 1+1D approximation
Authors:
Ya. N. Pavlyuchenkov,
L. A. Maksimova,
V. V. Akimkin
Abstract:
Heating by the central star is one of the key factors determining the physical structure of protoplanetary disks. Due to the large optical thickness in the radial direction, disk midplane regions are heated by the infrared radiation from the disk surface (atmosphere), which in turn is directly heated by the star. It was previously shown that interception of the stellar radiation by inhomogeneities…
▽ More
Heating by the central star is one of the key factors determining the physical structure of protoplanetary disks. Due to the large optical thickness in the radial direction, disk midplane regions are heated by the infrared radiation from the disk surface (atmosphere), which in turn is directly heated by the star. It was previously shown that interception of the stellar radiation by inhomogeneities on the disk surface can cause perturbations that propagate towards the star. In this work, we investigate the occurrence of such waves within a detailed 1+1D numerical model of the protoplanetary disk. We confirm the previous findings that in the disk, that is optically thick to its own radiation, the surface perturbations indeed occur and propagate towards the star. However, contrary to some analytical predictions, the thermal waves in sufficiently massive disks affect only the upper layers without significant fluctuations of temperature in the midplane. Our results indicate the need to study this instability within more consistent hydrodynamic models.
△ Less
Submitted 13 March, 2022;
originally announced March 2022.
-
Infrared appearance of wind-blown bubbles around young massive stars
Authors:
Maria S. Kirsanova,
Yaroslav N. Pavlyuchenkov
Abstract:
Thousands of ring-like bubbles appear on infrared images of the Galaxy plane. Most of these infrared bubbles form during expansion of HII regions around massive stars. However, the physical effects that determine their morphology are still under debate. Namely, the absence of the infrared emission toward the centres of the bubbles can be explained by pushing the dust grains by stellar radiation pr…
▽ More
Thousands of ring-like bubbles appear on infrared images of the Galaxy plane. Most of these infrared bubbles form during expansion of HII regions around massive stars. However, the physical effects that determine their morphology are still under debate. Namely, the absence of the infrared emission toward the centres of the bubbles can be explained by pushing the dust grains by stellar radiation pressure. At the same time, small graphite grains and PAHs are not strongly affected by the radiation pressure and must be removed by another process. Stellar ultraviolet emission can destroy the smallest PAHs but the photodestruction is ineffective for the large PAHs. Meanwhile, the stellar wind can evacuate all types of grains from HII regions. In the frame of our chemo-dynamical model we vary parameters of the stellar wind and illustrate their influence on the morphology and synthetic infrared images of the bubbles.
△ Less
Submitted 11 March, 2022;
originally announced March 2022.
-
Comparative Analysis of the Model for Exoplanet Atmosphere Outflow
Authors:
P. B. Isakova,
Ya. N. Pavlyuchenkov,
E. S. Kalinicheva,
V. I. Shematovich
Abstract:
Modeling the outflow of planetary atmospheres is important for understanding the evolution of exoplanet systems and for interpreting their observations. Modern theoretical models of exoplanet atmospheres become increasingly detailed and multicomponent, and this makes difficulties for engaging new researchers in the scope. Here, for the first time, we present the results of testing the gas-dynamic…
▽ More
Modeling the outflow of planetary atmospheres is important for understanding the evolution of exoplanet systems and for interpreting their observations. Modern theoretical models of exoplanet atmospheres become increasingly detailed and multicomponent, and this makes difficulties for engaging new researchers in the scope. Here, for the first time, we present the results of testing the gas-dynamic method incorporated in our aeronomic model, which has been proposed earlier. Undertaken tests support the correctness of the method and validate its applicability. For modeling the planetary wind, we propose a new hydrodynamic model equipped with a phenomenological function of heating by stellar UV radiation. The general flow in this model well agrees with results obtained in more detailed aeronomic models. The proposed model can be used for both methodical purposes and testing the gas-dynamic modules of self-consistent chemical-dynamic models of the planetary wind.
△ Less
Submitted 29 May, 2021;
originally announced May 2021.
-
Restoration of the Parameters of a Gas-Dust Disk Based on Its Synthetic Images
Authors:
A. M. Skliarevskii,
Ya. N. Pavlyuchenkov,
E. I. Vorobyov
Abstract:
The topic of the present study is combining a dynamic model of a protoplanetary disk with the computations of radiation transfer for obtaining synthetic spectra and disk images suitable for immediate comparison of the model with observations. Evolution of the disk was computed using the FEOSAD hydrodynamic model, which includes a self-consistent calculation of the dynamics of dust and gas in the 2…
▽ More
The topic of the present study is combining a dynamic model of a protoplanetary disk with the computations of radiation transfer for obtaining synthetic spectra and disk images suitable for immediate comparison of the model with observations. Evolution of the disk was computed using the FEOSAD hydrodynamic model, which includes a self-consistent calculation of the dynamics of dust and gas in the 2D thin disk approximation. Radiation transfer was simulated by the open code RADMC-3D. Three phases of disk evolution were considered: a young gravitationally unstable disk, a disk during an accretion luminosity burst, and an evolved disk. For these stages, the influence of various processes upon the disk's thermal structure was analyzed, as well as the differences between the temperatures obtained in the initial dynamic model and in the model with a detailed calculation of the radiation transfer. It is shown that viscous heating in the inner regions and adiabatic heating in the disk spirals can be important sources of heating. On the basis of the calculated spectral energy distributions, using SED-fitter software package used for the observations, physical parameters of the model disks were reconstructed. A significant spread between reconstructed parameters and initial characteristics of the disk indicates verification necessity of the models within the framework of spatially resolved observations of disks in the different spectral ranges
△ Less
Submitted 21 April, 2021;
originally announced April 2021.
-
Simulations of Dynamical Gas-Dust Circumstellar Disks: Going Beyond the Epstein Regime
Authors:
Olga P. Stoyanovskaya,
Fedor A. Okladnikov,
Eduard I. Vorobyov,
Yaroslav N. Pavlyuchenkov,
Vitaliy V. Akimkin
Abstract:
In circumstellar disks, the size of dust particles varies from submicron to several centimeters, while planetesimals have sizes of hundreds of kilometers. Therefore, various regimes for the aerodynamic drag between solid bodies and gas can be realized in these disks, depending on the grain sizes and velocities: Epstein, Stokes, and Newton, as well as transitional regimes between them. For small bo…
▽ More
In circumstellar disks, the size of dust particles varies from submicron to several centimeters, while planetesimals have sizes of hundreds of kilometers. Therefore, various regimes for the aerodynamic drag between solid bodies and gas can be realized in these disks, depending on the grain sizes and velocities: Epstein, Stokes, and Newton, as well as transitional regimes between them. For small bodies moving in the Epstein regime, the time required to establish the constant relative velocity between the gas and bodies can be much less than the dynamical time scale for the problem - the time for the rotation of the disk about the central body. In addition, the dust may be concentrated in individual regions of the disk, making it necessary to take into account the transfer of momentum between the dust and gas. It is shown that, for a system of equations for gas and monodisperse dust, a semi-implicit first-order approximation scheme in time in which the interphase interaction is calculated implicitly, while other forces, such as the pressure gradient and gravity are calculated explicitly, is suitable for stiff problems with intense interphase interactions and for computations of the drag in non-linear regimes. The piece-wise drag coefficient widely used in astrophysical simulations has a discontinuity at some values of the Mach and Knudsen numbers that are realized in a circumstellar disk. A continuous drag coefficient is presented, which corresponds to experimental dependences obtained for various drag regimes.
△ Less
Submitted 17 February, 2021;
originally announced February 2021.
-
Evolution of a Viscous Protoplanetary Disk with Convectively Unstable Regions. II. Accretion Regimes and Long-Term Dynamics
Authors:
L. A. Maksimova,
Ya. N. Pavlyuchenkov,
A. V. Tutukov
Abstract:
In this article, we proceed to study convection as a possible factor of episodic accretion in protoplanetary disks. Within the model presented in Article~I, the accretion history is analyzed at different rates and areas of matter inflow from the envelope onto the disk. It is shown that the burst-like regime occurs in a wide range of parameters. The long-term evolution of the disk is also modeled,…
▽ More
In this article, we proceed to study convection as a possible factor of episodic accretion in protoplanetary disks. Within the model presented in Article~I, the accretion history is analyzed at different rates and areas of matter inflow from the envelope onto the disk. It is shown that the burst-like regime occurs in a wide range of parameters. The long-term evolution of the disk is also modeled, including the decreasing-with-time matter inflow from the envelope. It is demonstrated that the disk becomes convectively unstable and maintains burst-like accretion onto the star for several million years. Meanwhile, the instability expands to an area of several tens of astronomical units and gradually decreases with time. It is also shown that at early stages in the disk evolution, conditions arise for gravitational instability in the outer parts of the disk and for dust evaporation in the convectively unstable inner regions. The general conclusion of the study is that convection can serve as one of the mechanisms of episodic accretion in protostellar disks, but this conclusion needs to be verified using more consistent hydrodynamic models.
△ Less
Submitted 16 September, 2020;
originally announced September 2020.
-
The PDR structure and kinematics around the compact HII regions S235A and S235C with [CII], [13CII], [OI] and HCO+ line profiles
Authors:
M. S. Kirsanova,
V. Ossenkopf-Okada,
L. D. Anderson,
P. A. Boley,
J. H. Bieging,
Ya. N. Pavlyuchenkov,
M. Luisi,
N. Schneider,
M. Andersen,
M. R. Samal,
A. M. Sobolev,
C. Buchbender,
R. Aladro,
Y. Okada
Abstract:
The aim of the present work is to study structure and gas kinematics in the photodissociation regions (PDRs) around the compact HII regions S235A and S235C. We observe the [CII], [13CII] and [OI] line emission, using SOFIA/upGREAT and complement them by data of HCO+ and CO. We use the [13CII] line to measure the optical depth of the [CII] emission, and find that the [CII] line profiles are influen…
▽ More
The aim of the present work is to study structure and gas kinematics in the photodissociation regions (PDRs) around the compact HII regions S235A and S235C. We observe the [CII], [13CII] and [OI] line emission, using SOFIA/upGREAT and complement them by data of HCO+ and CO. We use the [13CII] line to measure the optical depth of the [CII] emission, and find that the [CII] line profiles are influenced by self-absorption, while the [13CII] line remains unaffected by these effects. Hence, for dense PDRs, [13CII] emission is a better tracer of gas kinematics. The optical depth of the [CII] line is up to 10 in S235A. We find an expanding motion of the [CII]-emitting layer of the PDRs into the front molecular layer in both regions. Comparison of the gas and dust columns shows that gas components visible neither in the [CII] nor in low-J CO lines may contribute to the total column across S235A. We test whether the observed properties of the PDRs match the predictions of spherical models of expanding HII region + PDR + molecular cloud. Integrated intensities of the [13CII], [CII] and [OI] lines are well-represented by the model, but the models do not reproduce the double-peaked [CII] line profiles due to an insufficient column density of C+. The model predicts that the [OI] line could be a more reliable tracer of gas kinematics, but the foreground self-absorbing material does not allow using it in the considered regions.
△ Less
Submitted 30 July, 2020;
originally announced July 2020.
-
Impact of PAH photodissociation on the formation of small hydrocarbons in the Orion Bar and the Horsehead PDRs
Authors:
M. S. Murga,
M. S. Kirsanova,
A. I. Vasyunin,
Ya. N. Pavlyuchenkov
Abstract:
We study whether polycyclic aromatic hydrocarbons (PAHs) can be a weighty source of small hydrocarbons in photo-dissociation regions (PDRs). We modeled the evolution of 20 specific PAH molecules in terms of dehydrogenation and destruction of the carbon skeleton under the physical conditions of two well-studied PDRs, the Orion Bar and the Horsehead nebula which represent prototypical examples of PD…
▽ More
We study whether polycyclic aromatic hydrocarbons (PAHs) can be a weighty source of small hydrocarbons in photo-dissociation regions (PDRs). We modeled the evolution of 20 specific PAH molecules in terms of dehydrogenation and destruction of the carbon skeleton under the physical conditions of two well-studied PDRs, the Orion Bar and the Horsehead nebula which represent prototypical examples of PDRs irradiated by "high" and "low" ultraviolet radiation field. PAHs are described as microcanonical systems. The acetylene molecule is considered as the main carbonaceous fragment of the PAH dissociation as it follows from laboratory experiments and theory. We estimated the rates of acetylene production in gas phase chemical reactions and compared them with the rates of the acetylene production through the PAH dissociation. It is found that the latter rates can be higher than the former rates in the Orion Bar at $A_{\rm V}<1$ and also at $A_{\rm V}>3.5$. In the Horsehead nebula, the chemical reactions provide more acetylene than the PAH dissociation. The produced acetylene participate in the reactions of the formation of small hydrocarbons (C$_2$H, C$_3$H, C$_3$H$^{+}$, C$_3$H$_2$, C$_4$H). Acetylene production via the PAH destruction may increase the abundances of small hydrocarbons produced in gas phase chemical reactions in the Orion Bar only at $A_{\rm V}>3.5$. In the Horsehead nebula, the contribution of PAHs to the abundances of the small hydrocarbons is negligible. We conclude that the PAHs are not a major source of small hydrocarbons in both PDRs except some locations in the Orion Bar.
△ Less
Submitted 13 July, 2020;
originally announced July 2020.
-
Evolution of a Viscous Protoplanetary Disk with Convectively Unstable Regions
Authors:
Ya. N. Pavlyuchenkov,
A. V. Tutukov,
L. A. Maksimova,
E. I. Vorobyov
Abstract:
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 hy…
▽ More
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.
△ Less
Submitted 18 December, 2019;
originally announced December 2019.
-
Molecular envelope around the HII region RCW 120
Authors:
M. S. Kirsanova,
Ya N. Pavlyuchenkov,
D. S. Wiebe,
P. A. Boley,
S. V. Salii,
S. V. Kalenskii,
A. M. Sobolev,
L. D Anderson
Abstract:
The H II region RCW120 is a well-known object, which is often considered as a target to verify theoretical models of gas and dust dynamics in the interstellar medium. However, the exact geometry of RCW120 is still a matter of debate. In this work, we analyse observational data on molecular emission in RCW120 and show that 13CO(2-1) and C18O(2-1) lines are fitted by a 2D model representing a ring-l…
▽ More
The H II region RCW120 is a well-known object, which is often considered as a target to verify theoretical models of gas and dust dynamics in the interstellar medium. However, the exact geometry of RCW120 is still a matter of debate. In this work, we analyse observational data on molecular emission in RCW120 and show that 13CO(2-1) and C18O(2-1) lines are fitted by a 2D model representing a ring-like face-on structure. The changing of the C18O(3-2) line profile from double-peaked to single-peaked from the dense molecular Condensation 1 might be a signature of stalled expansion in this direction. In order to explain a self-absorption dip of the 13CO(2-1) and 13CO(3-2) lines, we suggest that RCW120 is surrounded by a diffuse molecular cloud, and find confirmation of this cloud on a map of interstellar extinction. Optically thick 13CO(2-1) emission and the infrared 8 um PAH band form a neutral envelope of the H II region resembling a ring, while the envelope breaks into separate clumps on images made with optically thin C18O(2-1) line and far-infrared dust emission.
△ Less
Submitted 14 August, 2019;
originally announced August 2019.
-
Development and application of fast methods for computing momentum transfer between gas and dust in supercomputer simulation of planet formation
Authors:
Olga P. Stoyanovskaya,
Vitaly V. Akimkin,
Eduard I. Vorobyov,
Tatiana A. Glushko,
Yaroslav N. Pavlyuchenkov,
Valeriy N. Snytnikov,
Nikolay V. Snytnikov
Abstract:
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…
▽ More
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.
△ Less
Submitted 15 November, 2018;
originally announced November 2018.
-
Survival of a planet in short-period Neptunian desert under effect of photo-evaporation
Authors:
Dmitry E. Ionov,
Yaroslav N. Pavlyuchenkov,
Valery I. Shematovich
Abstract:
Despite the identification of a great number of Jupiter-like and Earth-like planets at close-in orbits, the number of "hot Neptunes" - the planets with 0.6-18 times of Neptune mass and orbital periods less than 3 days - turned out to be very small. The corresponding region in the mass-period distribution was assigned as the "short-period Neptunian desert". The common explanation of this fact is th…
▽ More
Despite the identification of a great number of Jupiter-like and Earth-like planets at close-in orbits, the number of "hot Neptunes" - the planets with 0.6-18 times of Neptune mass and orbital periods less than 3 days - turned out to be very small. The corresponding region in the mass-period distribution was assigned as the "short-period Neptunian desert". The common explanation of this fact is that the gaseous planet with few Neptune masses would not survive in the vicinity of host star due to intensive atmosphere outflow induced by heating from stellar radiation. To check this hypothesis we performed numerical simulations of atmosphere dynamics for a hot Neptune. We adopt the previously developed self-consistent 1D model of hydrogen-helium atmosphere with suprathermal electrons accounted. The mass-loss rates as a function of orbital distances and stellar ages are presented. We conclude that the desert of short-period Neptunes could not be entirely explained by evaporation of planet atmosphere caused by the radiation from a host star. For the less massive Neptune-like planet, the estimated upper limits of the mass loss may be consistent with the photo-evaporation scenario, while the heavier Neptune-like planets could not lose the significant mass through this mechanism. We also found the significant differences between our numerical results and widely used approximate estimates of the mass loss.
△ Less
Submitted 22 May, 2018; v1 submitted 12 March, 2018;
originally announced March 2018.
-
Influence of Photoelectrons on the Structure and Dynamics of the Upper Atmosphere of a Hot Jupiter
Authors:
D. E. Ionov,
V. I. Shematovich,
Ya. N. Pavlyuchenkov
Abstract:
A self-consistent, aeronomic model of the upper atmosphere of a "hot Jupiter" including reactions involving suprathermal photoelectrons is presented. This model is used to compute the height profiles of the gas density, velocity, and temperature in the atmosphere of the exoplanet HD 209458b. It is shown that including suprathermal electrons when computing the heating and cooling functions reduces…
▽ More
A self-consistent, aeronomic model of the upper atmosphere of a "hot Jupiter" including reactions involving suprathermal photoelectrons is presented. This model is used to compute the height profiles of the gas density, velocity, and temperature in the atmosphere of the exoplanet HD 209458b. It is shown that including suprathermal electrons when computing the heating and cooling functions reduces the mass loss rate of the atmosphere by a factor of five.
△ Less
Submitted 7 June, 2017;
originally announced June 2017.
-
Improving the thin-disk models of circumstellar disk evolution. The 2+1-dimensional model
Authors:
Eduard I. Vorobyov,
Yaroslav N. Pavlyuchenkov
Abstract:
Circumstellar disks of gas and dust are naturally formed from contracting pre-stellar molecular cores during the star formation process. To study various dynamical and chemical processes that take place in circumstellar disks prior to their dissipation and transition to debris disks, the appropriate numerical models capable of studying the long-term disk chemodynamical evolution are required. We p…
▽ More
Circumstellar disks of gas and dust are naturally formed from contracting pre-stellar molecular cores during the star formation process. To study various dynamical and chemical processes that take place in circumstellar disks prior to their dissipation and transition to debris disks, the appropriate numerical models capable of studying the long-term disk chemodynamical evolution are required. We present a new 2+1-dimensional numerical hydrodynamics model of circumstellar disk evolution, in which the thin-disk model is complemented with the procedure for calculating the vertical distributions of gas volume density and temperature in the disk. The reconstruction of the disk vertical structure is performed at every time step via the solution of the time-dependent radiative transfer equations coupled to the equation of the vertical hydrostatic equilibrium. We perform a detailed comparison between circumstellar disks produced with our previous 2D model and with the improved 2+1D approach. The structure and evolution of resulting disks, including the differences in temperatures, densities, disk masses and protostellar accretion rates, are discussed in detail. The new 2+1D model yields systematically colder disks, while the in-falling parental clouds are warmer. Both effects act to increase the strength of disk gravitational instability and, as a result, the number of gravitationally bound fragments that form in the disk via gravitational fragmentation as compared to the purely 2D thin-disk simulations with a simplified thermal balance calculation.
△ Less
Submitted 1 June, 2017;
originally announced June 2017.
-
Dust dynamics and evolution in HII regions - II. Effects of dynamical coupling between dust and gas
Authors:
V. V. Akimkin,
M. S. Kirsanova,
Ya. N. Pavlyuchenkov,
D. S. Wiebe
Abstract:
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.…
▽ More
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.
△ Less
Submitted 30 April, 2017;
originally announced May 2017.
-
On the technique for the recovery of the spectrum of turbulence in astrophysical disks
Authors:
D. V. Bisikalo,
E. P. Kurbatov,
Ya. N. Pavlyuchenkov,
A. G. Zhilkin,
P. V. Kaygorodov
Abstract:
We present a method that can be used to recover the spectrum of turbulence from observations of optically thin emission lines formed in astrophysical disks. Within this method we analyze how line intensity fluctuations depend on the angular resolution of the instrument, used for the observations. The method allows us to restore the slope of the power spectrum of velocity turbulent pulsations and e…
▽ More
We present a method that can be used to recover the spectrum of turbulence from observations of optically thin emission lines formed in astrophysical disks. Within this method we analyze how line intensity fluctuations depend on the angular resolution of the instrument, used for the observations. The method allows us to restore the slope of the power spectrum of velocity turbulent pulsations and estimate the upper boundary of the turbulence scale.
△ Less
Submitted 9 March, 2016;
originally announced March 2016.
-
Dust dynamics and evolution in expanding HII regions. I. Radiative drift of neutral and charged grains
Authors:
V. V. Akimkin,
M. S. Kirsanova,
Ya. N. Pavlyuchenkov,
D. S. Wiebe
Abstract:
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 smal…
▽ More
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.
△ Less
Submitted 24 February, 2015;
originally announced February 2015.
-
Thermal structure of a protostellar envelope
Authors:
Ya. N. Pavlyuchenkov,
A. G. Zhilkin,
E. I. Vorobyov,
A. M. Fateeva
Abstract:
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…
▽ More
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.
△ Less
Submitted 17 February, 2015;
originally announced February 2015.
-
Gas-phase CO depletion and N2H+ abundances in starless cores
Authors:
N. Lippok,
R. Launhardt,
D. Semenov,
A. M. Stutz,
Z. Balog,
Th. Henning,
O. Krause,
H. Linz,
M. Nielbock,
Ya. N. Pavlyuchenkov,
M. Schmalzl,
A. Schmiedeke,
H. J. Bieging
Abstract:
Seven isolated, nearby low-mass starless molecular cloud cores have been observed as part of the Herschel key program Earliest Phases of Star formation (EPoS). By applying a ray-tracing technique to the obtained continuum emission and complementary (sub)mm emission maps, we derive the physical structure (density, dust temperature) of these cloud cores. We present observations of the 12CO, 13CO, an…
▽ More
Seven isolated, nearby low-mass starless molecular cloud cores have been observed as part of the Herschel key program Earliest Phases of Star formation (EPoS). By applying a ray-tracing technique to the obtained continuum emission and complementary (sub)mm emission maps, we derive the physical structure (density, dust temperature) of these cloud cores. We present observations of the 12CO, 13CO, and C18O (2-1) and N2H+ (1-0) transitions towards the same cores. Based on the density and temperature profiles, we apply time-dependent chemical and line-radiative transfer modeling and compare the modeled to the observed molecular emission profiles. CO is frozen onto the grains in the center of all cores in our sample. The level of CO depletion increases with hydrogen density and ranges from 46% up to more than 95% in the core centers in the core centers in the three cores with the highest hydrogen density. The average hydrogen density at which 50% of CO is frozen onto the grains is 1.1+-0.4 10^5 cm^-3. At about this density, the cores typically have the highest relative abundance of N2H+. The cores with higher central densities show depletion of N2H+ at levels of 13% to 55%. The chemical ages for the individual species are on average 2+-1 10^5 yr for 13CO, 6+-3 10^4 yr for C18O, and 9+-2 10^4 yr for N2H+. Chemical modeling indirectly suggests that the gas and dust temperatures decouple in the envelopes and that the dust grains are not yet significantly coagulated. We observationally confirm chemical models of CO-freezeout and nitrogen chemistry. We find clear correlations between the hydrogen density and CO depletion and the emergence of N2H+. The chemical ages indicate a core lifetime of less than 1 Myr.
△ Less
Submitted 12 October, 2013; v1 submitted 13 August, 2013;
originally announced August 2013.
-
Polycyclic aromatic hydrocarbons in spatially resolved extragalactic star forming complexes
Authors:
M. S. Khramtsova,
D. S. Wiebe,
P. A. Boley,
Ya. N. Pavlyuchenkov
Abstract:
The abundance of polycyclic aromatic hydrocarbons (PAHs) in low- and high-metallicity galaxies has been widely discussed since the time when detailed infrared data for extragalactic objects were first obtained. On the scales of entire galaxies, a smaller PAH abundance in lower-metallicity galaxies is often observed. We study this relationship for star-forming regions in nearby galaxies, for a samp…
▽ More
The abundance of polycyclic aromatic hydrocarbons (PAHs) in low- and high-metallicity galaxies has been widely discussed since the time when detailed infrared data for extragalactic objects were first obtained. On the scales of entire galaxies, a smaller PAH abundance in lower-metallicity galaxies is often observed. We study this relationship for star-forming regions in nearby galaxies, for a sample containing more than 200 HII complexes, using spatially-resolved observations from the Herschel Space Observatory and Spitzer Space Telescope. We use a model for the dust emission to estimate the physical parameters (PAH abundance, metallicity, ultraviolet radiation field, etc.) of these complexes. The same correlation of PAH abundance with metallicity, as seen for entire galaxies, is apparently preserved at smaller scales, at least when the Kobulnicky & Kewley metallicity calibration is used. We discuss possible reasons for this correlation, noting that traces of less-effective PAH formation in low-metallicity AGB stars should be smeared out by radial mixing in galactic disks. Effective destruction by the harder and more intensive ultraviolet field in low-metallicity environments is qualitatively consistent with our data, as the ultraviolet field intensity, derived from the infrared photometry, is indeed smaller in HII complexes with lower metallicity.
△ Less
Submitted 20 February, 2013;
originally announced February 2013.
-
Stochastic grain heating and mid-infrared emission in protostellar cores
Authors:
Ya. N. Pavlyuchenkov,
D. S. Wiebe,
V. V. Akimkin,
M. S. Khramtsova,
Th. Henning
Abstract:
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μ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…
▽ More
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μ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μ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.
△ Less
Submitted 3 January, 2012;
originally announced January 2012.
-
Determining the Parameters of Massive Protostellar Clouds via Radiative Transfer Modeling
Authors:
Ya. N. Pavlyuchenkov,
D. S. Wiebe,
A. M. Fateeva,
T. S. Vasyunina
Abstract:
A one-dimensional method for reconstructing the structure of prestellar and protostellar clouds is presented. The method is based on radiative transfer computations and a comparison of theoretical and observed intensity distributions at both millimeter and infrared wavelengths. The radiative transfer of dust emission is modeled for specified parameters of the density distribution, central star, an…
▽ More
A one-dimensional method for reconstructing the structure of prestellar and protostellar clouds is presented. The method is based on radiative transfer computations and a comparison of theoretical and observed intensity distributions at both millimeter and infrared wavelengths. The radiative transfer of dust emission is modeled for specified parameters of the density distribution, central star, and external background, and the theoretical distribution of the dust temperature inside the cloud is determined. The intensity distributions at millimeter and IR wavelengths are computed and quantitatively compared with observational data. The best-fit model parameters are determined using a genetic minimization algorithm, which makes it possible to reveal the ranges of parameter degeneracy as well. The method is illustrated by modeling the structure of the two infrared dark clouds IRDC-320.27+029 (P2) and IRDC-321.73+005 (P2). The derived density and temperature distributions can be used to model the chemical structure and spectral maps in molecular lines.
△ Less
Submitted 14 November, 2010;
originally announced November 2010.
-
Modeling of Protostellar Clouds and their Observational Properties
Authors:
A. G. Zhilkin,
Ya. N. Pavlyuchenkov,
S. N. Zamozdra
Abstract:
A physical model and two-dimensional numerical method for computing the evolution and spectra of protostellar clouds are described. The physical model is based on a system of magneto-gasdynamical equations, including ohmic and ambipolar diffusion, and a scheme for calculating the thermal and ionization structure of a cloud. The dust and gas temperatures are determined during the calculations of…
▽ More
A physical model and two-dimensional numerical method for computing the evolution and spectra of protostellar clouds are described. The physical model is based on a system of magneto-gasdynamical equations, including ohmic and ambipolar diffusion, and a scheme for calculating the thermal and ionization structure of a cloud. The dust and gas temperatures are determined during the calculations of the thermal structure of the cloud. The results of computing the dynamical and thermal structure of the cloud are used to model the radiative transfer in continuum and in molecular lines. We presented the results for clouds in hydrostatic and thermal equilibrium. The evolution of a rotating magnetic protostellar cloud starting from a quasi-static state is also considered. Spectral maps for optically thick lines of linear molecules are analyzed. We have shown that the influence of the magnetic field and rotation can lead to a redistribution of angular momentum in the cloud and the formation of a characteristic rotational velocity structure. As a result, the distribution of the velocity centroid of the molecular lines can acquire an hourglass shape. We plan to use the developed program package together with a model for the chemical evolution to interpret and model observed starless and protostellar cores.
△ Less
Submitted 25 June, 2009;
originally announced June 2009.
-
Problems of Star Formation Theory and Prospects of Submillimeter Observations
Authors:
D. Z. Wiebe,
M. S. Kirsanova,
B. M. Shustov,
Ya. N. Pavlyuchenkov
Abstract:
We consider current state of star formation theory and requirements to observations in millimeter and submillimeter ranges which are necessary for resolution of the most actual problems of the physics of star formation. Two key features of star-forming regions which define observational requirements to their studies, are relatively low energy of processes that take place there and smallness of c…
▽ More
We consider current state of star formation theory and requirements to observations in millimeter and submillimeter ranges which are necessary for resolution of the most actual problems of the physics of star formation. Two key features of star-forming regions which define observational requirements to their studies, are relatively low energy of processes that take place there and smallness of corresponding spatial scales. This is especially true for the objects in the latest stages of ``pre-stellar'' evolution, that is, hot cores, hyper- and ultracompact HII regions, and protoplanetary disks. Angular resolution, sensitivity, and spectral coverage in existing projects of ground-based and space telescopes of submillimeter and millimeter range are not completely adequate to necessary requirements. To obtain detailed information on star-forming regions as well as on individual protostars it is necessary to employ a space-based interferometer.
△ Less
Submitted 24 July, 2008;
originally announced July 2008.