A coupled,dynamical and chemical model for the prestellar core L1544: Comparison of modeled and observed C<Superscript>18</Superscript>O,HCO<Superscript>+</Superscript>, and CS emission spectra

We consider radiative transfer in C¹⁸O, HCO⁺, and CS molecular lines in a spherically symmetrical, coupled, dynamical and chemical model of a prestellar core whose evolution is determined by ambipolar diffusion. Theoretical and observed line profiles are compared for the well-studied core L1544, which may be a collapsing protostellar cloud. We study the relationship between the line shapes and model parameters. The structure of the envelope and kinematic parameters of the cloud are the most important factors determining the shape of the lines. Varying the input model parameters for the radiative transfer—the kinetic temperature and microturbulent velocity—within the limits imposed by observations does not result in any substantial variations of the line profiles. The comparison between the model and observed spectra indicates that L1544 displays a flattened structure, and is viewed at an oblique angle. A two-dimensional model is needed to reproduce this structure.     

UV-controlled physical and chemical structure of protoplanetary disks

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.     

Modeling of protostellar clouds and their observational properties

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-gas-dynamical 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 when calculating the thermal structure of the cloud. The results of computing the dynamical and thermal structure of the cloud are used to model the transfer of continuum and molecular-line radiation in the cloud. Results are presented for clouds in hydrostatic and thermal equilibrium. The evolution of a rotating magnetic protostellar cloud that is compressed starting from a quasi-static equilibrium state is also considered. Spectral maps for optically thick lines of linear molecules are analyzed. 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. It is planned in future to use the developed program package and a model for the chemical evolution of clouds to interpret and model in detail observed starless and protostellar cores.     

Infrared emission and the destruction of dust in HII regions

The generation of infrared (IR) radiation and the observed IR-intensity distribution at wavelengths of 8, 24, and 100 µm in the ionized hydrogen region around a young, massive star is investigated. The evolution of the HII region is treated using a self-consistent chemical-dynamical model in which three dust populations are included—large silicate grains, small graphite grains, and polycyclic, aromatic hydrocarbons (PAHs). A radiative transfer model taking into account stochastic heating of small grains and macromolecules is used to model the IR spectral energy distribution. The computational results are compared with Spitzer and Herschel observations of the RCW 120 nebula. The contributions of collisions with gas particles and the radiation field of the star to stochastic heating of small grains are investigated. It is shown that a model with a homogeneous PAH content cannot reproduce the ring-like IR-intensity distribution at 8 µm. A model in which PAHs are destroyed by ultraviolet radiation of the star, generating region HII, provides a means to explain this intensity distribution. This model is in agreement with observations for realistic characteristic destruction times for the PAHs.     

Models of astrophysical decretion-accretion diffusional disks

The formation of gaseous diffusional accretion-decretion disks is an important stage in the evolution of numerous astronomical objects. Matter is accreted onto the object in the accretion part of these disks, while the angular momentum of the accreted matter is transported from the central region to the periphery in the decretion part. Here, we consider general questions connected with the formation and evolution of diffusive accretion-decretion disks in various astrophysical objects. Such disks can be described using nonstationary diffusion models. The phenomenological parameters of these models are the coefficients in the relations for the characteristic turbulent velocity and mean free path of diffusion elements in the disk. We have developed a numerical technique to compute the disk evolution for a number of models (a massive disk, a disk with continuous accretion, a purely decretion disk). Analytical expressions estimating the basic parameters of accretion-decretion disks are presented. We discuss the relationship between the models considered and the classical α model of an accretion disk.     

Evolution of a Viscous Protoplanetary Disk with Convectively Unstable Regions

The role of convection in the gas-dust accretion disk around a young star is studied. The evolution of a Keplerian disk is modeled using the Pringle equation, which describes the time variations of the surface density under the action of turbulent viscosity. The distributions of the density and temperature in the polar directions are computed simultaneously in the approximation that the disk is hydrostatically stable. The computations of the vertical structure of the disk take into account heating by stellar radiation, interstellar radiation, and viscous heating. The main factor governing evolution of the disk in this model is the dependence of the viscosity coefficient on the radius of the disk. The computations of this coefficient take into account the background viscosity providing the continuous accretion of the gas and the convective viscosity, which depends on the parameters of the convection at a given radius. The results of computations of the global evolution and morphology of the disk obtained in this approach are presented. It is shown that, in the adopted model, the accretion has burst-like character: after the inner part of the disk ($$R < 3$$ AU) is filled with matter, this material is transferred relatively rapidly 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.     

Influence of luminosity bursts on properties of protostellar disks

A (2+ 1)-dimensional numerical model for the formation and evolution of young stellar objects with sub-solar masses is presented. The numerical hydrodynamicall code describing the formation and evolution of a pølar disk in a two-dimensional approximation is supplemented by one-dimensional code for the evolution of the star and an algorithm for establishing the vertical structure of the disk. This code is used to investigate the influence of luminosity bursts with intensities similar to those observed in FU Orionis objects (FUors) on the properties and thermal balance of pølar disks. A model with gravitational instability and fragmentation of the disk, with subsequent migration of the fragments onto the protostar, is used as a basic model for FUors. Typical FUor bursts (L ～ 100 L ^⊙) can appreciably influence the thermal balance of their disks and parent envelopes, leading to an increase in the disk temperature by more than a factor of two. On the other hand, massive fragments in the disk are only weakly perturbed by such bursts, partially due to screening by the disk and partially due to their high temperature brought about by adiabatic heating. Apart from massive fragments, the characteristic thermal time scales are appreciably shorter than the dynamical time scales throughout the radial extent of the disk and envelope; this enables the use of a stationary radiative-transfer equation when determining the vertical structure of the disk.     

Structure of CB 26 protoplanetary disk derived from millimeter dust continuum maps

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 the 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 diameter of approximately 35 AU, in agreement with the conclusions of other studies. The inclination of the disk is 78°, 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.     

A multicomponent model for computing the thermal structure of collapsing protostellar clouds

A model for simulating the thermal and dynamical evolution of protostellar clouds is presented. In the model, the dust and gas temperatures are treated separately, making it possible to more precisely describe the initial stages of the cloud’s gravitational contraction and collapse. The model is fast enough to be applied in hydrodynamical computations, and has a high enough accuracy for the results to be used to compute emission spectra and comparing them with observational data. Two problems are considered as test examples and simple applications: calculation of the structure of clouds in thermal and hydrostatic equilibrium, and modeling the evolution of a protostellar cloud in a spherically symmetric approximation, including the formation of the first hydrostatic core.     

Influence of photoelectrons on the structure and dynamics of the upper atmosphere of a hot Jupiter

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.