Analytical continuation of stability of periodic orbits in the restricted three-body problem

A supplement to the theory of analytical continuation of circular orbits in the restricted three-body problem is presented. The first order stability is given analytically to the first power of mass parameter . The theory of the Kirkwood gaps is discussed from this point of view. The stability limit which should determine the size of accretion discs in binaries is found to be in good agreement with earlier numerical experiments for < 1/2.     

Hydrodynamics of accretion columns

A detailed understanding of how the infalling matter in accretion columns is decelerated is essential for the calculation of the emitted radiation. On neutron stars, the deceleration takes place mainly by the interaction of the plasma with radiation, at least for the high-luminosity sources. We report on our two-dimensional calculations of the hydrodynamic flow in such accretion columns. The radiation transport is treated in the diffusion approximation, and we are looking for a stationary solution for the velocity field. The dependence of the results on physical parameters, especially on the accretion rate is discussed. Due to the non-linearity of the problem it turns out that only in certain parameter ranges stationary solutions seem to exist. For accretion rates higher than a critical value there are no stationary accretion flows. This leads us to the conclusion that a time-dependent picture for the accretion is unavoidable.Paper presented at the IAU Colloquium No. 93 on Cataclysmic Variables. Recent Multi-Frequency Observations and Theoretical Developments, held at Dr. Remeis-Sternwarte Bamberg, F.R.G., 16–19 June, 1986.     

Two-dimensional numerical models of the boundary layer of accretion disks in cataclysmic variables

Nearly half of the total available accretion energy can be released in the boundary layer (BL) if the accreting object is slowly rotating. The spectral distribution of the emitted radiation depends crucially on the internal structure of the BL. Up to now no detailed models concerning the BL exist. We have developed an explicit two-dimensional numerical method written for axisymmetric accretion flows including viscosity effects and angular momentum. We display our first models concerning the BL structure incorporating variation of the stellar rotation and of the fraction of the released energy. The first results show a strong dependence of the BL structure on the local rate of cooling and on the rotation of the primary.Paper presented at the IAU Colloquium No. 93 on Cataclysmic Variables. Recent Multi-Frequency Observations and Theoretical Developments, held at Dr. Remeis-Sternwarte Bamberg, F.R.G., 16–19 June, 1986.     

The effect of wind on magnetized advection dominated accretion flows

We study the effects of winds on advection dominated accretion flows in the presence of a global magnetic field under a self-similar treatment. The disk gas is assumed to be isothermal. For a steady state structure of such accretion flows a set of self similar solutions are presented. We consider the wind in a general magnetic field with three components (r,φ,z) in advection-dominated accretion flows. The mass-accretion rate $\dot{M}$ decreases with radius r as $\dot{M}\propto r^{s+1/2}$ , where s is an arbitrary constant. We will see, by increasing the wind parameter s, radial and rotational velocities increase.     

Radiation hydrodynamics of high-luminosity accretion into black holes

To extend Shapiro's (1973a, b) calculations of black hole accretion to the regimes of interstellar gas densities and of black hole masses for which emergent luminosities are expected to be high, the radiation hydrodynamics of spherically symmetric gas flows in static isotropic metrics is discussed. Since for the more luminous objects the optical depth of the accretion volume becomes large, particular attention has to be paid to radiative transfer through non-Euclidean spaces, and a method for solving the full transfer problem is presented. The method is applied to accretion into black holes of mass between 10M and 10⁵ M , under the conservative assumption that all other heat sources, like dissipation of magnetic or turbulent energy, can be neglected in comparison to the compressional work term,p dV. In the interstellar gas parameter range of interest, the radiation field is then dominated by emission and absorption of synchrotron radiation from inner zones of the flow. Temperature stratifications, luminosities and emergent spectra resulting from these processes are calculated.     

An analysis of active close binaries (CB) based on photometric measurements

For the purpose of analysing light curves of active CB with an accretion disc being at the evolutional phase of an intensive matter exchange between the components a model for light-curve synthesis has been realized where the attention is given to systems like W Ser not sufficiently examined yet with regard that in them the accretion disc is formed around an ordinary star. In the paper one uses the elements presented in the first paper of this series. The model can successfully describe the essential characteristics of the observed light curves due to existence of an accretion disc and a hot spot, as well as those originated in the temperature distribution along the disc radius. The system components are considered in the framework of the nonsynchronous Roche model and the accretion disc of a constant thickness lies in the orbital plane around the star capturing the matter of the neighbouring component.The primary surrounded by the disc is situated relatively well within the Roche oval and its rotation can be significantly nonsynchronous. Near the Lagrange equilibrium pointL ₁ flows from the secondary (which fills the Roche limit) the gas stream nourishing the disc. In the zone where the stream touches the lateral side of the disc a hot spot is formed.The proposed model enables estimating of the basic orbital and physical parameters of active type W Ser CB (Djuraevi, 1991) on the basis of photometric measurements by applying the inverse-problem method.     

The formation and evolution of low mass protostars

A review is presented of the earliest stages of protostellar evolution. Observations of prestellar cores, which are believed to represent the initial conditions for protostellar collapse, depart significantly from the scale-free density distribution which is usually taken as the starting point for the formation of a low-mass protostar. Pre-stellar cores are observed to have radial density profiles which have flat inner regions, steepening towards their edges. This is seen to qualitatively match the predictions of the Bonnor-Ebert stability criterion for pressure-bounded self-gravitating gas clouds. From these initial conditions, theoretical modelling of cores threaded by magnetic fields predicts that quasi-static evolution by the process of ambipolar diffusion will lead to a significantly different starting point for collapse than the static singular isothermal sphere.This departure from a scale-free density distribution for the initial conditions has recently been shown to produce an ensuing protostellar collapse which has a non-constant accretion rate. Recently published observations of low-mass protostars in the Ophiuchi cluster are demonstrated to be consistent with such a non-constant protostellar mass accretion rate, contrary to the standard protostellar collapse model. Instead, the data appear consistent with an initially high accretion rate, which subsequently decays. The initial phase of high accretion rate is labelled the main accretion phase, during which 50 per cent of the circumstellar envelope mass is accreted in 10 per cent of the total accretion time, and which is equated observationally with Class 0 objects. The subsequent phase with roughly an order of magnitude lower accretion rate is labelled the late accretion phase, during which the remainder of the envelope mass is accreted in the remaining 90 per cent of the total accretion time, at an order of magnitude lower accretion rate, and which is equated observationally with Class I objects. The growth of circumstellar discs begins in the Class 0 stage, and proceeds through the Class I and II stages. Published data of the Taurus star-forming region currently available appear also to be consistent with this scenario.     

Giant planet formation

The accumulation of giant planets involves processes typical for terrestrial planet formation as well as gasdynamic processes that were previously known only in stars. The condensible element cores of the gas-giants grow by solid body accretion while envelope formation is governed by stellar-like equilibria and the dynamic departures thereof. Two hypotheses for forming Uranus/Neptune-type planets — at sufficiently large heliocentric distances while allowing accretion of massive gaseous envelopes, i.e. Jupiter-type planets at intermediate distances — have been worked out in detailed numerical calculations: (1) Hydrostatic gas-accretion models with time-dependent solid body accretion-rates show a slow-down of core-accretion at the appropriate masses of Uranus and Neptune. As a consequence, gas-accretion also stagnates and a window is opened for removing the solar nebula during a time of roughly constant envelope mass. (2) Gasdynamic calculations of envelope accretion for constant planetesimal accretion-rates show a dynamic transition to new envelope equilibria at the so called critical mass. For a wide range of solar nebula conditions the new envelopes have respective masses similar to those of Uranus and Neptune and are more tightly bound to the cores. The transitions occur under lower density conditions typical for the outer parts of the solar nebula, whereas for higher densities, i.e. closer to the Sun, gasdynamic envelope accretion sets in and is able to proceed to Jupiter-masses.     

The visco-resistive stability of arcades in the corona of the Sun

The stability of ballooning modes in coronal arcades is studied using linear visco-resistive MHD. Rigid wall conditions are adopted for modelling the photospheric line-tying of the magnetic field. The full Braginskii viscosity stress tensor is used and particular attention is given to the effect of the viscosity coefficient ₃ which was left out of an earlier investigation by Van der Linden, Goossens, and Hood (1987, 1988). The numerical results for shearless arcades show that the coefficient ₃ has a stabilizing effect. However, for realistic values of the equilibrium quantities the stabilizing effect by ₃ can be neglected in comparison with the strong stabilizing effect of the perpendicular viscosity. The effect of magnetic field strength and mode number on stability are determined. In particular it is found that there exists a critical field strength for every mode number such that the mode is stable for weaker fields and unstable for stronger fields.     

Supermassive black hole in an elliptical galaxy: Accretion of a hot gas with a low but finite angular momentum

The accretion of hot slowly rotating gas onto a supermassive black hole is considered. The important case where the velocities of turbulent pulsations at the Bondi radius r _B are low, compared to the speed of sound c _s, is studied. Turbulence is probably responsible for the appearance of random average rotation. Although the angular momentum at r _B is low, it gives rise to the centrifugal barrier at a depth r _c = l ² /GM _BH ? r _B, that hinders supersonic accretion. The numerical solution of the problem of hot gas accretion with finite angular momentum is found taking into account electron thermal conductivity and bremsstrahlung energy losses of two temperature plasma for density and temperature near Bondi radius similar to observed in M87 galaxy. The saturation of the Spitzer thermal conductivity was also taken into account. The parameters of the saturated electron thermal conductivity were chosen similar to the parameters used in the numerical simulations of interaction of the strong laser beam radiation with plasma targets. These parameters are confirmed in the experiments. It is shown that joint action of electron thermal conductivity and free-free radiation leads to the effective cooling of accreting plasma and formation of the subsonic settling of accreting gas above the zone of a centrifugal barrier. A toroidal condensation and a hollow funnel that separates the torus from the black hole emerge near the barrier. The barrier divides the flow into two regions: (1) the settling zone with slow subKeplerian rotation and (2) the zone with rapid supersonic nearly Keplerian rotation. Existence of the centrifugal barrier leads to significant decrease of the accretion rate ? in comparison with the critical Bondi solution for γ = 5/3 for the same values of density and temperature of the hot gas near Bondi radius. Shear instabilities in the torus and related friction cause the gas to spread slowly along spirals in the equatorial plane in two directions.As a result, outer (r > r _c) and inner (r < r _c) disks are formed. The gas enters the immediate neighborhood of the black hole or the zone of the internal ADAF flow along the accretion disk (r < r _c). Since the angular momentum is conserved, the outer disk removes outward an excess of angular momentum along with part of the matter falling into the torus. It is possible, that such outer Keplerian disk was observed by Hubble Space Telescope around the nucleus of the M87 galaxy in the optical emission lines. We discuss shortly the characteristic times during which the accretion of the gas with developed turbulence should lead to the changes in the orientation of the torus, accretion disk and, possibly, of the jet.     

Similarity solution for unsteady accretion flow

Similarity solution for unsteady accretion flow in a gravitational field of a point mass is obtained. Characteristic features of the flow pattern are discussed. It is shown that the shock waves appeared in the accretion flow propagate outward asr _s t ^2/3.     

Radial self-gravity of accretion disc around supermassive black holes

This paper studies the properties of self-gravity of accretion discs around supermassive black holes. With integration on the thin disc configuration, this paper has calculated the radial and vertical components of self-gravity of accretion discs. The discussion mainly concentrates on the radial component, and the results are briefly as follows: for accretion discs around supermassive black holes (M10 ⁸–10¹⁰ M). At the distance where (R/R _g)10⁵–10⁴, the radial component of self-gravity dominates over the central one where the dynamical structure of the accretion discs completely differs from that of Keplerian disc. A turbulence driven by radial self-gravity instability as a kind of energy source is proposed. This paper has two criteria for the comparison of magnitude between the self-gravity of accretion discs and the gravity of the central black hole, from which an analytic estimation for the outer radius of the accretion discs has been derived. The results of this paper may be used to explain the accretion discs of quasars and AGNs.This research was supported by the National Natural Science Foundation of China.     

The accretion of interstellar matter by white dwarfs

The problem of the energy released by interstellar matter in its accretion by low-luminosity stars is examined. It is found that nearly all the gravitational energy is available as the energy source for the accreting star. When this source of energy becomes predominant, a luminosity-temperature relation can be derived.Evolutionary consequences for the white dwarfs are discussed.The behaviour of selected classes of white dwarfs is recognized in the theoretical H.R. diagram.The evidence for a separation between field white dwarfs and Hyades cluster's white dwarfs is shortly discussed. For the Hyades white dwarfs a well defined sequence in good agreement with a constant accretion locus is found.A critical discussion is finally given on the efficiency of the accretion mechanism in the Hyades dwarfs.

---

Sommario Viene discusso il problema dell'energia ceduta dalla materia interstellare nei processi di cattura gravitazionale da parte di oggetti stellari di bassa luminosità. Si mostra come le maggior parte di questa energia sia effettivamente ceduta alla struttura stellare.Nei casi in cui l'energia di accrescimento predomini sulle altre sorgenti di energia è possibile ricavare una relazione luminosità-temperatura. Sono discusse possibili conseguenze sull'evoluzione delle nane bianche.E'esaminato il comportamento di alcune classi di nane bianche nel diagramma H.R. teorico, e viene brevemente discussa l'evidenza per una reale separazione delle nane bianche di campo da quelle nell'ammasso delle Iadi. Per queste ultime in particolare si trova una sequenza molto ben definita, in ottimo accordo con un luogo di costante accrescimento.Un esame critico è infine compiuto riguardo la possibile efficienza dell'accrescimento nelle nane delle Iadi.

     

Doppler tomography of accretion regions in algols

The technique of Doppler Tomography has been used to image the accretion regions in five short-period Algols. There is clear evidence of gas flows along the predicted free-fall path of the gas stream as well as asymmetric disk-like structures around the mass gainer. Another source of H emission is associated with the cool magnetically active star.     

The accretion of matter by a collapsing star in the presence of a magnetic field. II. Selfconsistent stationary picture

The stationary two-dimensional magnetohydrodynamic solution for the accretion of the matter without pressure into a gravitating centre of a black hole is obtained. It is assumed that the magnetic field far from the collapsed star is homogeneous and its influence on the flow is negligible. Around the star, at the plane perpendicular to the direction of the magnetic field, the dense quasistationary disc is formed, the structure of which in a large extent is determined by dissipation processes. The structure is calculated for (a) a laminar disc with the Coulomb mechanism of dissipation; and (b) a turbulent disc.The estimations of the parameters of the shock which result from the infall of the matter onto the disc are given. In the last section the numerical estimation and approximate character of the radiation spectrum of the disc and the shock are obtained for two cases of 10M and 10⁵ M . The luminosity of collapsed objects withM=10M appears to be about solar, thus its observation is possibly only at the distances less than 300–1000 pc. The collapsed objects in the Galaxy withM=10⁵ M could constitute very bright sources in spectral regions from optical up to X-ray. The spectra of a laminar and a turbulent disc for 10M black hole are very different. The laminar disc radiates primarily in the ultraviolet. The turbulent disc radiates a large part of its flux in the infrared. Therefore, one cannot exclude the possibility that some of the galactic infrared star-like sources are individual black holes in the accretion state.     

Description of accretion induced outflows from ultra-luminous sources to under-luminous AGNs

We study the energetics of the accretion-induced outflow and then plausible jet around black holes/compact objects using a newly developed disc-outflow coupled model. Inter-connecting dynamics of outflow and accretion essentially upholds the conservation laws. The energetics depend strongly on the viscosity parameter $\alpha$ and the cooling factor f which exhibit several interesting features. The bolometric luminosities of ultra-luminous X-ray binaries (e.g. SS433) and family of highly luminous AGNs and quasars can be reproduced by the model under the super-Eddington accretion flows. Under appropriate conditions, low-luminous AGNs (e.g. Sagittarius $A^1$) also fit reasonably well with the luminosity corresponding to a sub-Eddington accretion flow with $f\to 1$.     

Critical accretion flow of gas onto compact stars

The critical accretion flow of gas onto compact stars with mass of 0.6M is investigated by numerical integrations of the time-dependent hydrodynamic equations in the sphericallysymmetric and optically thick case. For the compact stars surrounded by such a dense cloud of gas, the radiation pressure force decelerates the infall gas significantly and free fall regime of the gas is not at all attained. This results in incident low velocities at the standing shock front close to the stellar surface, low temperatures of the gas around the compact stars, and no X-ray in white dwarfs but soft X-rays in neutron stars, respectively. Some applications of the results to the X-ray sources are discussed.     

Accretion of the gaseous envelope of Jupiter around a 5-10 Earth-mass core

New numerical simulations of the formation and evolution of Jupiter are presented. The formation model assumes that first a solid core of several M_⊕ accretes from the planetesimals in the protoplanetary disk, and then the core captures a massive gaseous envelope from the protoplanetary disk. Earlier studies of the core accretion-gas capture model [Pollack, J.B., Hubickyj, O., Bodenheimer, P., Lissauer, J.J., Podolak, M., Greenzweig, Y., 1996. Icarus 124, 62-85] demonstrated that it was possible for Jupiter to accrete with a solid core of 10-30 M_⊕ in a total formation time comparable to the observed lifetime of protoplanetary disks. Recent interior models of Jupiter and Saturn that agree with all observational constraints suggest that Jupiter's core mass is 0-11 M_⊕ and Saturn's is 9-22 M_⊕ [Saumon, G., Guillot, T., 2004. Astrophys. J. 609, 1170-1180]. We have computed simulations of the growth of Jupiter using various values for the opacity produced by grains in the protoplanet's atmosphere and for the initial planetesimal surface density, σinit,Z, in the protoplanetary disk. We also explore the implications of halting the solid accretion at selected core mass values during the protoplanet's growth. Halting planetesimal accretion at low core mass simulates the presence of a competing embryo, and decreasing the atmospheric opacity due to grains emulates the settling and coagulation of grains within the protoplanet's atmosphere. We examine the effects of adjusting these parameters to determine whether or not gas runaway can occur for small mass cores on a reasonable timescale. We compute four series of simulations with the latest version of our code, which contains updated equation of state and opacity tables as well as other improvements. Each series consists of a run without a cutoff in planetesimal accretion, plus up to three runs with a cutoff at a particular core mass. The first series of runs is computed with an atmospheric opacity due to grains (hereafter referred to as ‘grain opacity’) that is 2% of the interstellar value and . Cutoff runs are computed for core masses of 10, 5, and 3 M_⊕. The second series of Jupiter models is computed with the grain opacity at the full interstellar value and . Cutoff runs are computed for core masses of 10 and 5 M_⊕. The third series of runs is computed with the grain opacity at 2% of the interstellar value and . One cutoff run is computed with a core mass of 5 M_⊕. The final series consists of one run, without a cutoff, which is computed with a temperature dependent grain opacity (i.e., 2% of the interstellar value for ramping up to the full interstellar value for ) and . Our results demonstrate that reducing grain opacities results in formation times less than half of those for models computed with full interstellar grain opacity values. The reduction of opacity due to grains in the upper portion of the envelope with has the largest effect on the lowering of the formation time. If the accretion of planetesimals is not cut off prior to the accretion of gas, then decreasing the surface density of planetesimals lowers the final core mass of the protoplanet, but increases the formation timescale considerably. Finally, a core mass cutoff results in a reduction of the time needed for a protoplanet to evolve to the stage of runaway gas accretion, provided the cutoff mass is sufficiently large. The overall results indicate that, with reasonable parameters, it is possible that Jupiter formed at 5 AU via the core accretion process in 1 Myr with a core of 10 M_⊕ or in 5 Myr with a core of 5 M_⊕.     

On disk accretion

A simple but sufficiently accurate method for calculating an accretion disk structure is developed. The detailed analysis of the accretion disk is fulfilled by using this method. The effect of turbulent heat transfer on the disk structure is taken into account along with the effect of radiative transfer. The turbulent heat transfer is shown to play an important role, and may be even predominant in the inner disk region. The dependence of temperatureT and density on thez-coordinate is found. Simple analytic expressions are proposed for the run of the density in all the disk zones. It is shown that the inner disk region is convectively stable. The main parameters of all the zones are derived. The geometry of regular motions is studied; the regular hydrodynamical flows are found to appear directed both toward and outside the central object. A detailed comparison is made with the results of other authors.