A possible mechanism for the formation of tilted disks in intermediate polars

Using 3D gas dynamics, we numerically simulate accretion-disk formation in typical cataclysmic variable intermediate polars with dipolar magnetic fields (B _a = 10⁵?5 × 10⁵ G) and misaligned white-dwarf magnetic and rotation axes. Our simulations confirm that a significant misalignment of the axes results in a significant misalignment of the disk to the orbital plane. However, over time, this disk tilt disappears: early in the simulation, the initial particle positions in the rarefied tilted disk are governed solely by the magnetic field of the white dwarf. Due to the increasing disk mass and hence increasing disk gas pressure, the tilted disk eventually becomes decoupled from the magnetic field. The tidal action of the donor leads to a retrograde (i.e., nodal) precession of the tilted disk’s streamlines, and the disk becomes twisted. When the disk tilt is greater than 4°, the incoming gas stream no longer strikes the disk rim (i.e., bright shocked region). Matter is now transported over and under the disk rim to the inner regions of the disk. Over time, the increased mass of inner parts of the disk due to the action of the colinear gas stream returns the inner-disk regions to a colinear configuration. Meanwhile, the outer regions of the tilted, twisted disk become warped. Our simulations suggest that the lifetime of an intermediate polar’s tilted disk could be several tens to thousands of orbital periods.     

Observational Manifestations of the Precessional Wave in the Accretion Disk of a Cataclysmic Variable Star

Effects due to the interaction of the steam from the inner Lagrangian point with the accretion disk in a cataclysmic variable star are considered. The results of three-dimensional gas-dynamical numerical simulations confirm that the disk thickness in the vicinity of the interaction with the stream is minimum when the component-mass ratio is 0.6. As a consequence, some of the matter from the stream does not collide with the outer edge of the accretion disk, and continues its motion unperturbed toward the accretor. This part of the stream subsequent interacts (collides) with a thickening of the accretion disk due to the presence of a precessional wave in the disk, leading to the appearance of an additional zone of heating at the disk surface. This additional zone of enhanced luminosity (hot spot) is a direct observational manifestation of the precessional wave in the accretion disk.     

Structure of the accretion disk in the cataclysmic variable UX UMa in the transition to its quiescent state

Our analysis of BV RI light curves for the cataclysmic variable UX UMa obtained in intermediate activity states, in the transition between the active and quiescent states of the system on March 12, 1997 and May 3, 2000, shows that the shapes of these light curves cannot be adequately described using the standard hot-spot model. A model with a “hot line” near the edge of the disk and a two-armed spiral structure on the disk surface reproduces much better out-of-eclipse variations in the light curves. The presence of an extended hot line can explain the anomalous shape of the I light curve on March 12, 1997. The decrease in the observed luminosity of the system between March 12, 1997 and May 3, 2000 could be due to a decrease in the disk luminosity by a factor of 2–2.5; the higher disk luminosity on the earlier date is associated with appreciable deviations of the radial temperature distribution of the disk material from that for the standard model. The phases and depths of dips in the out-of-eclipse sections of the UX UMa light curves are due primarily to the parameters of the complex shape of the accretion disk, which has a spiral structure located mainly near its outer edge. The contribution of the spiral arms in the V filter reaches 20–50% of the total disk radiation. The crest of the first spiral wave in our model maintains its approximate position in azimuth; this structure could represent a bulge in a halo at the outer edge of the disk near orbital phases φ ～ 0.7, in the direction of the continuation of the extended shock in the disk itself. The position of the crest of the second spiral arm changes with time. This structure may represent a one-armed spiral wave near the apastron of the weakly elliptical disk. Finally, the observations testify to the presence of another spiral arm that is les clearly manifest in terms of both its luminosity and its height above the unperturbed disk surface. Thus, in an intermediate activity state of UX UMa, the surface of the accretion disk is distorted by the action of a two-armed spiral structure in the outer regions of the disk, which is asymmetric in both its luminosity and dimensions, and a bulge at the disk edge in the region of its interaction with the inflow to the disk.     

Mathematical modeling of the structure of an accretion disk in a stellar binary system

Numerical simulations of gas-dynamical processes taking place in the accretion disk of a stellar binary system are presented. The initial state of the disk is an equilibrium gaseous configuration. Mechanisms for the development of spiral waves and associated variations in the angular momentum of the gas are considered. The influence of the ratio of the binary-component masses and the initial disk configuration are investigated. It is concluded that the existence of a steady-state disk is impossible without a flow of gas from the donor star.     

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.     

Three-dimensional hydrodynamical modeling of accretion-disk formation in microquasars. Cen X-3

We have carried out three-dimensional hydrodynamical modeling of the formation of an accretion disk around a compact object due to radiative wind of a massive donor in a close binary system. The massive X-ray binary Cen X-3, which has a precessing accretion disk and may possess relativistic jets, is considered as an example. The computations show that, when the action of the central compact object on the formation of the wind is taken into account, the radiative wind forms an accretion disk with a radius of 0.16 (in units of the orbital separation), which accretes at a rate close to 1 × 10^?8 M _⊙/yr. In this model, the disk is spherically symmetrical and geometrically thick, with a tunnel going from the accretor to the upper layers of the disk along the accretor’s rotational axis at the disk center. The number density of the gas in the tunnel is five orders of magnitude lower than in the disk. The wind-disk interaction at the outer boundary of the disk produces a strong shock (wind-disk shock) directed toward the donor. The black-body emission of the disk and tunnel is nonstationary, and resembles the outbursts observed in Cen X-3. An analysis of the location of the region of nonstationary emission suggests that the outbursts occur in the wind-disk shock.     

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.     

The possible origin of the mini-spiral in the nuclear disk of the Galaxy as a result of the super-reflection instability

It is shown that the mini-spiral in the nuclear disk of the Galaxy could arise as a result of super-reflection instability in the disk. This instability develops at a break in the profile of the disk rotation. The structure of the perturbation obtained via numerical modeling is visually similar to the observed mini-spiral.     

Formation of ring structures in galactic disks during close passages of galaxies

The formation of ring structures in galactic disks is investigated. It is shown that, in addition to the known mechanism of forming rings in “head-on” collisions between galaxies, ring structures can be formed during close passages of galaxies if the perturbing galaxy moves in a plane close to the equatorial plane of the perturbed disk galaxy, opposite to the direction of rotation of the disk. Numerical simulations of the formation of structures in the disk of a massive galaxy undergoing a passage with another galaxy are considered. The results of these cmputations show the formation of pronounced ring structures in the galactic disk when the initial inclination of the trajectory of the perturbing galaxy to the equatorial plane of the perturbed galaxy is no more than ~25°. However, the probability of close passages of galaxies with these parameters is small, as is the probability of head-on collisions. The characteristic time scale for the existence of pronounced rings is of order the dynamical time scale at the edge of the galaxy, 200–300 million years, close to the corresponding time for head-on collisions. The evolution of the rings has the same character in both cases: they gradually expand and move toward the periphery of the galaxy. The results of these simulations can also be applied to a close passage of one star by another star with a protoplanetary disk. According to the computation results, the characteristic time scale for the existence of pronounced rings in such a protoplanetary disk depends mainly on the size of the disk; this time scale can reach several tens of thousands of years for a disk radius of about 1000 AU. The formation of ring structures in such a disk could influence the formation and evolution of planetesimals, and possibly the character of the formation of planets and the distribution of their orbital semi-major axes.     

The structure of matter flows in semi-detached binaries after the termination of mass transfer

We present the results of three-dimensional gas-dynamical simulations of matter flows in semi-detached binaries after termination of the mass transfer between the components of the system. The structure of the residual accretion disk is studied. When the mass transfer has ended, the quasi-elliptical disk becomes circular and its structure changes: tidal interactions result in the formation of a second arm in the spiral shock wave. In addition, a condensation (blob) moving through the disk with variable velocity is formed. The blob is maintained by interactions with the arms of the spiral shock and exists essentially over the entire lifetime of the disk. We also show that, for a viscosity corresponding to α～0.01 (typical for observed accretion disks), the lifetime of the residual disk is about 50 orbital periods.     

Three-dimensional hydrodynamical modeling of the formation of the accretion disk in the SS 433 binary system

Three-dimensional hydrodynamical modeling of the formation of the accretion disk in the SS 433 binary system is carried out with various types of cooling and numerical grids. These computations show that a thick accretion disk with a height of 0.25–0.30 (in units of the component separation) is formed around the compact object, from a flow with a large radius (0.2–0.3 in the same units) that forms in the vicinity of the inner Lagrangian point. This disk has the form of a flattened torus. The number of orbits of a particle of gas in the disk is 100–150, testifying to a minimal influence of numerical viscosity in these computations. The computations also show that the stream flowing from L₁ is nearly conservative, and spirals in the disk are not formed due to the influence of the donor gravitation.     

The role of close passages of galaxies and the asymmetry of their dark haloes in the formation of their spiral patterns

The influence of close passages of galaxies on the shapes of disk galaxies and the distribution of stars in them is studied for several types of interactions in the framework of the restricted N-body problem. Depending on the conditions adopted, either two spiral density waves or ring structures are formed in the stellar disk of the galaxy. These structures can generate star formation fronts with the corresponding shape, as are observed in disk galaxies. Our calculations can also be applied to study the influence of the passage of a nearby star on a protoplanetary disk. The formation of ring structures there could specify the type of planet formation in the outer regions of the planetary system and the distribution of semimajor axes for the planetary orbits. We use the same model to study the generation and evolution of spiral density waves in the stellar disks of galaxies as a result of the recently found asymmetry of the gravitational potential in the massive dark haloes in disk galaxies. The dipole component of the gravitational field of the halo can continuously permanently generate the spiral structure in disk galaxies.     

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.     

Kinematic properties of the Galactic thick disk from data from the RAVE DR4 catalog

The analysis of the kinematic properties of the Galactic thick disk based on data from modern catalogs of stellar radial velocities and proper motions is presented. A new aspect of new determination of the kinematic characteristics of the thick disk is that the selected objects define this disk’s properties near the plane of symmetry. The velocity dispersion of stars in the Galactic thick disk in the radial direction and the direction of the Galactic rotation have been determined. The stellar-velocity distribution in the direction of the rotation is asymmetric. The parameters of this asymmetry have been determined, and the lag of the rotational velocity of the thick disk relative to objects in the thin disk estimated. The value of this “asymmetric drift,” about 20 km/s, suggests larger spatial scales for the kinematic characteristics in the radial direction for the Galactic thick disk than for the thin 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.     

Estimating the masses of the spherical and disk components of galaxies via numerical simulations

A series of numerical N-body simulations is performed in order to dynamically model the properties of four galaxies (NGC 5603, NGC 3198, NGC 891, and NGC 1566) with known rotation curves, radial disk scales L, and velocity dispersions of old disk stars at various galactocentric distances r. Each model includes a three-dimensional collisionless disk and rigid spherical components, whose relative mass μ was treated as a free parameter that differed from simulation to simulation. The observed disk stellar velocity dispersions were assumed to be equal to or (in the general case) greater than the corresponding line-of-sight projections of the simulated values for the adopted μ after the initially unstable disk is heated and arrives at a steady state. A comparison of the simulated and observed rotational velocities and velocity dispersions provides evidence for “light” disks with μ≥2 in the disk (r<4L).     

A possible mechanism for ionization of the Magellanic stream by disk stars

In our model describing the leakage of ionizing radiation from the Galactic disk into the halo, disk stars can contribute substantially to the ionization of halo objects such as high-velocity clouds and the Magellanic stream. This ionization is produced by a relatively hard radiation field, which can maintain its ionizing effect even at a considerable distance from the plane of the disk.     

Amplification of spiral density waves in gaseous and dust galactic disks as a result of interaction with the halo

The effect of the discrete structure of the halo on local oscillations of a galactic disk is analyzed. Such effects have much in common with dynamical friction. Gaseous and stellar disks are considered; in both cases, some leading spiral density waves are unstable. Bending oscillations of the disk can also be unstable when the disk interacts with the halo.     

The development of large-scale instability in Keplerian stellar accretion disks

Our previous studies of large-scale vortical flows arising in shear flows of stellar accretion disks with Keplerian azimuthal velocity distributions as a result of the development of small perturbations are continued. The development of large-scale instability in an accretion disk is investigated via mathematical modeling. One result obtained is the change of the disk flow structure due to the formation of large vortices. In the limiting case, sufficiently long evolution leads to the formation of several asymmetric spiral structures of the flow of disk matter. The presence of large-scale structures leads to angular-momentum redistribution in the disk.