Trans-Neptunian objects (TNOs) with diameter greater than 100 km currently moving in not too eccentric orbits could be formed directly by the contraction of large rarefied condensations. Along with the gravitational influence of planets, gravitational interactions of TNOs played a certain role in their orbital evolution as well. More than 20% of Earth-crossing objects could have come from the trans-Neptunian belt. TNOs and Centaurs (invisible comets mainly beyond Jupiter) could produce an important contribution to the dust content of the interplanetary dust cloud.

Our efficient numerical approach has been applied to modeling the asymmetric circumstellar dust disk around Beta Pictoris as observed with the HST/STIS. We present a new model on the origin of the warping of the Beta Pic disk. We suggest that the observed warp is formed by the gravitational influence of a planet with a mass of about 10 masses of Earth, at a distance of 70 AU, and a small inclination ($\sim 2.5^\circ$) of the planetary orbit to the main dust disk. Results of our modeling are compared with the STIS observations.

(abridged) This review is based on an extensive work done in collaboration with N. Gorkavyi, J. Mather, and T. Taidakova, which aimed at the physical modelling of the interplanetary dust (IPD) cloud in the Solar system, i.e., establishing a link between the observable characteristics of the zodiacal cloud and the dynamical and physical properties of the parent minor bodies. Our computational approach permits with modest computer resources to integrate the trajectories of hundreds of particles and to effectively store up to $10^{10-11}$ positions, which provides a high fidelity 3D distribution of the dust. Our numerical codes account for the major dynamical effects that govern the motion of IPD particles. We have demonstrated the efficacy of our codes by performing the following analyses: (i) simulation of the distribution of Centaurs (comets scattered in their journey from the Kuiper belt inward the Solar system) and revealing the effects of the outer planets in producing `cometary belts'; (ii) detailed inspection of a rich resonant structure found in these belts, which predicts the existence of gaps similar to the Kirkwood gaps in the main asteroid belt; (iii) a preliminary 3-D physical model of the IPD cloud, which includes three dust components -- asteroidal, cometary, and kuiperoidal -- and is consistent with the available data of Pioneers and Voyagers dust detectors; (iv) modeling of the IPD cloud, which provides a zodiacal light distribution in accord, to the order of 1%, with a subset of the COBE/DIRBE observations; and (v) showing that the resonant structure in dusty circumstellar disks of Vega and Epsilon Eridani is a signature of embedded extrasolar planets. Further improvements of our modelling and their importance for astronomy and cosmology are outlined.

Our numerical simulations of the Edgeworth-Kuiper belt objects gravitationally scattered by the four giant planets accounting for mean motion resonances reveal numerous resonant gaps in the distribution of the scattered population.

We apply our recently elaborated, powerful numerical approach to the high-resolution modeling of the structure and emission of circumstellar dust disks, incorporating all relevant physical processes. Specifically, we examine the resonant structure of a dusty disk induced by the presence of one planet. It is shown that the planet, via resonances and gravitational scattering, produces (1) an asymmetric resonant dust belt with one or more clumps intermittent with one or a few off-center cavities; and (2) a central cavity void of dust. These features can serve as indicators of a planet embedded in the circumstellar dust disk and, moreover, can be used to determine its major orbital parameters and even the mass of the planet. The results of our study reveal a remarkable similarity with various types of highly asymmetric circumstellar disks observed with the James Clerk Maxwell Telescope around Epsilon Eridani and Vega. The proposed interpretation of the clumps in those disks as being resonant patterns is testable -- it predicts the asymmetric design around the star to revolve, viz., by 1.2--1.6 deg/yr about Vega and 0.6--0.8 deg/yr about Epsilon Eri.

(Abridged) Using an efficient computational approach, we have reconstructed the structure of the dust cloud in the Solar system between 0.5 and 100 AU produced by the Kuiper belt objects. Our simulations offer a 3-D physical model of the `kuiperoidal' dust cloud based on the distribution of 280 dust particle trajectories produced by 100 known Kuiper belt objects ; the resulting 3-D grid consists of $1.9\times 10^6$ cells containing $1.2\times 10^{11}$ particle positions. The following processes that influence the dust particle dynamics are taken into account: 1) gravitational scattering on the eight planets (neglecting Pluto); 2) planetary resonances; 3) radiation pressure; and 4) the Poynting-Robertson (P-R) and solar wind drags. We find the dust distribution highly non-uniform: there is a minimum in the kuiperoidal dust between Mars and Jupiter, after which both the column and number densities of kuiperoidal dust sharply increase with heliocentric distance between 5 and 10 AU, and then form a plateau between 10 and 50 AU. Between 25 and 45 AU, there is an appreciable concentration of kuiperoidal dust in the form of a broad belt of mostly resonant particles associated with Neptune. In fact, each giant planet possesses its own circumsolar dust belt consisting of both resonant and gravitationally scattered particles. As with the cometary belts simulated in our related papers (Ozernoy, Gorkavyi, & Taidakova 2000a,b), we reveal a complex resonant structure of the dust belts containing rich families of resonant peaks and gaps. An important result is that both the column and number dust density are more or less flat between 10 and 50 AU, which might explain the surprising data obtained by Pioneers 10 & 11 and Voyager that the dust number density remains approximately distance-independent in this region.

We have applied a powerful numerical approach to compute, with a high resolution, the structure of dusty circumstellar disks with embedded planets. We emphasize some testable implications of our simulations which would verify the presence of a planet via thermal emission of one or more dusty clumps which are in mean motion resonances with the planet. In particular, our simulations indicate that Vega may have a massive planet of m~2 m_J (m_J being Jupiter's mass) at a distance of 50-60 AU, and Epsilon Eri may have a less massive planet of m~0.2 m_J at a similar distance of 55-65 AU. This conclusion is testable: Each resonant feature is stationary in the reference frame co-rotating with the planet, but it is not so for the observer at Earth. Therefore, if our interpretation of asymmetric clumps in circumstellar disks as dynamical resonant structures is correct, the above pattern revolves around the star with an angular velocity of (1.2-1.6) deg/yr (Vega) and (0.6-0.8) deg/yr (Epsilon Eri) -- a prediction that can be tested on a timescale of several years.

Using numerical simulations, we examine the structure of a cometary population near a massive planet, such as a giant planet of the Solar system, starting with one-planet approximation (the Sun plus one planet). By studying the distributions of comets in semimajor axis, eccentricity, pericenter, and apocenter distances, we have revealed several interesting features in these distributions. The most remarkable ones include (i) spatial accumulation of comets near the planetary orbit (which we call the `cometary belt') and (ii) avoidance of resonant orbits by comets. Then we abandon one-planet approximation and examine as to how a cometary belt is modified when the influence of all four giant planets is taken into consideration. To this end, we simulate a stationary distribution of comets, which results from the gravitational scattering of the Kuiper belt objects on the four giant planets and accounts for the effects of mean motion resonances. Accounting for the influence of four giant planets makes the cometary belts overlapping, but nevertheless keeping almost all their basic features found in one-planet approximation. In particular, the belts maintain the gaps in the (a,e)- and (a,i)-space similar to the Kirkwood gaps in the main asteroid belt. We conclude that the large-scale structure of the Solar system is featured by the four cometary belts expected to contain 20-30 millions of scattered comets, and only a tiny fraction of them is currently visible as Jupiter-, Saturn-, etc. family comets.

Non-steady and eruptive phenomena in quasars are thought to be associated with the Eddington or super-Eddington luminous stage. Although there is no lack in hypotheses about the total duration of such a stage, the latter remains essentially unknown. We calculate the duration of quasar luminous phase in dependence upon the initial mass of a newborn massive black hole (MBH) by comparing the observed luminosity- and redshift distributions of quasars with mass distribution of the central MBHs in normal galactic nuclei. It is assumed that, at the quasar stage, each MBH goes through a single (or recurrent) phase(s) of accretion with, or close to, the Eddington luminosity. The mass distributions of quasars is found to be connected with that of MBHs residing in normal galaxies by a one-to-one corrrespondence through the entire mass range of the inferred MBHs if the accretion efficiency of mass-to-energy transformation $\eta \sim 0.1$.

We develop a physical model of the zodiacal cloud incorporating the real dust sources of asteroidal, cometary, and kuiperoidal origin. Using the inferred distribution of the zodiacal dust, we compute its thermal emission and scattering at several wavelengths (1.25, 5, and 20 $\mu$m) as a function of NGST location assumed to be at 1 AU or 3 AU. Areas on the sky with a minimum of zodiacal light are determined.

The determination of velocities of stars from precise Doppler measurements is described here using relativistic theory of astronomical reference frames so as to determine the Keplerian and post-Keplerian parameters of binary systems. We apply successive Lorentz transformations and the relativistic equation of light propagation to establish the exact treatment of Doppler effect in binary systems both in special and general relativity theories. As a result, the Doppler shift is a sum of (1) linear in $c^{-1}$ terms, which include the ordinary Doppler effect and its variation due to the secular radial acceleration of the binary with respect to observer; (2) terms proportional to $c^{-2}$, which include the contributions from the quadratic Doppler effect caused by the relative motion of binary star with respect to the Solar system, motion of the particle emitting light and diurnal rotational motion of observer, orbital motion of the star around the binary's barycenter, and orbital motion of the Earth; and (3) terms proportional to $c^{-2}$, which include the contributions from redshifts due to gravitational fields of the star, star's companion, Galaxy, Solar system, and the Earth. After parameterization of the binary's orbit we find that the presence of periodically changing terms in the Doppler schift enables us disentangling different terms and measuring, along with the well known Keplerian parameters of the binary, four additional post-Keplerian parameters, including the inclination angle of the binary's orbit, $i$. We briefly discuss feasibility of practical implementation of these theoretical results, which crucially depends on further progress in the technique of precision Doppler measurements.

In part I, using an effective computational approach, we have reconstructed the population of dust sources between Jupiter and Neptune. Here, in part II, we present the results on distribution of dust produced by 157 real sources (100 Jupiter-family comets with semi-major axes $a<20$ AU, 51 Kuiper belt, and 6 Centaur objects) as well as 211 fictitious sources taken from our computed sample. The following processes that influence the dust particle dynamics are taken into account: 1) gravitational scattering on four giant planets; 2) planetary resonances; and 3) the Poynting-Robertson (P-R) and solar wind drags. A file consisting of $0.9\times 10^6$ particle positions has been computed to simulate the dust distribution in the outer parts of the Solar system. We find that this distribution is highly non-uniform, with most of the dust concentrating into four belts associated with the orbits of four giant planets, with a sharp rise (depending on the size of particles) at the innermost part of the ring. As in I, we reveal a rich and sophisticated resonant structure of these belts containing families of resonances and gaps. A dissipative nature of the P-R drag results in specific features of particle's capture into, and evolution in, the resonances. Based on our simulations, we expect a new, quasi-stationary dust population to exist in the belts near Jupiter and Saturn, which is highly inclined and possesses large eccentricities. This population is basically non-resonant and is an important addition to otherwise resonant dust belts. The simulated dust is likely the main source of the zodiacal light in the outer Solar system, which will be analyzed in our further work.

We employ an efficient numerical approach to simulate a stationary distribution of test objects, which results from their gravitational scattering on the four giant planets, with accounting for effects of mean motion resonances. Using the observed distribution of the Kuiper belt objects, we reconstruct, in the space of orbital coordinates, the distribution function $n(a,e,i)$ for the population of minor bodies beyond Jupiter. We confirm that thousands of large yet cold comets and Centaurs might be located between the orbits of Jupiter and Neptune. Moreover, we find as an important result that they are concentrated into four circumsolar belts, with a highly non-uniform and well structured distribution of the objects. This huge yet unrevealed population, with only a few of its representatives presently known, is expected to have, like our simulations demonstrate, a rich resonant structure containing both density maxima and gaps. The resonant structure is formed due to gravitational perturbations, i.e. in a non-dissipative way. If plotted in the ($a,e,i$)-space of orbital coordinates, the belts contain gaps (including those between resonant groups), quite similar to the Kirkwood gaps in the main asteroid belt. An appreciable fraction of the test bodies reveals, for some time, an accumulation near (rather than in) the resonances, both interior and exterior, with the giant planets. An accompanying paper considers the population simulated in this work as the major source of dust in the outer Solar system. The simple but fast and efficient numerical approach employed in this work would allow the reader for applying it to many other problems of his/her interest.

Recently we have elaborated a new cosmological model of gamma-ray burst (GRB) origin (1998, ApJ 502, 192), which employs the dynamical evolution of central dense stellar clusters in the galactic nuclei. Those clusters inevitably contain a large fraction of compact stellar remnants (CSRs), such as neutron stars (NSs) and stellar mass black holes (BHs), and close encounters between them result in radiative captures into short-living binaries, with subsequent merging of the components, thereby producing GRBs (typically at large distances from the nucleus). In the present paper, we calculate the redshift distribution of the rate of GRBs produced by close encounters of NSs in distant galactic nuclei. To this end, the following steps are undertaken: (i) we establish a connection between the parameters of the fast evolving central stellar clusters (i.e. those for which the time of dynamical evolution exceeds the age of the Universe) with masses of the forming central supermassive black holes (SMBHs) using a dynamical evolution model; (ii) we connect these masses with the inferred mass distributions of SMBHs in the galactic nuclei and the redshift distribution of quasars by assuming a certain `Eddington luminosity phase' in their activity; (iii) we incorporate available observational data on the redshift distribution of quasars as well as a recently found correlation between the masses of galaxies and their central SMBHs. The resulting redshift distribution of the GRB rate, which accounts for both fast and slowly evolving galactic nuclei is consistent with that inferred from the BATSE data if the fraction of fast evolving galactic nuclei is in the range $0.016-0.16$.

We discuss mechanisms for angular momentum transport in the clumpy medium of the circumnuclear disk at the Galactic center. The viscosity due to clump-clump collisions is found to be less than some critical viscosity; this meets the conditions at which a collective mode of nonaxisymmetric shear perturbations in the disk is able to grow until going into the saturation regime where fully developed turbulence is established. We find that the angular momentum transfer due to this turbulent viscosity turns out to comparable to the transport due to magnetic torques. Taken together, the turbulent and magnetic transfer of angular momentum are able to provide the inflow of mass into the central parsec with a rate of about 10^{-2} M_{sun}/yr, consistent with the available data.

A new cosmological scenario for the origin of gamma ray bursts (GRBs) is proposed. In our scenario, a highly evolved central core in the dense galactic nucleus is formed containing a subsystem of compact stellar remnants (CSRs), such as neutron stars and black holes. Those subsystems result from the dynamical evolution of dense central stellar clusters in the galactic nuclei through merging of stars, thereby forming (as has been realized by many authors) the short-living massive stars and then CSRs. We estimate the rate of random CSR collisions in the evolved galactic nuclei by taking into account, similar to Quinlan & Shapiro (1987), the dissipative encounters of CSRs, mainly due to radiative losses of gravitational waves, which results in the formation of intermediate short-living binaries, with further coalescence of the companions to produce GRBs. We also consider how the possible presence of a central supermassive black hole, formed in a highly evolved galactic nucleus, influences the CSR binary formation. This scenario does not postulate ad hoc a required number of tight binary neutron stars in the galaxies. Instead, it gives, for the most realistic parameters of the evolved nuclei, the expected rate of GRBs consistent with the observed one, thereby explaining the GRB appearance in a natural way of the dynamical evolution of galactic nuclei. In addition, this scenario provides an opportunity for a cosmological GRB recurrence, previously considered to be a distinctive feature of GRBs of a local origin only. We also discuss some other observational tests of the proposed scenario.

Just as the ordinary Doppler effect serves as a tool to measure radial velocities of celestial objects, so can the relativistic Doppler effect be implemented to measure a combination of radial and transverse velocities by using recent improvements in observing techniques. A key element that makes a further use of this combination feasible is the periodicity in changes of the orbital velocity direction for the source. Two cases are considered: (i) a binary star; and (ii) a solitary star with the planetary companion. It is shown that, in case (i), several precision Doppler measurements employing the gas absorption cell technique would determine both the total orbital velocity and the inclination angle of the binary orbit disentangled from the peculiar velocity of the system. The necessary condition for that is the measured, at least with a modest precision, proper motion and distance to the system.

We point out that a high number density of stars in the core of a dense star cluster such as the central stellar cluster at the Galactic center, where many stars possess strong stellar winds, should result in collisions of those winds. The wind collisions in the dense stellar core at the Galactic center would result in production of strong X-ray flares with the rate of $\sim 10^{-4}~ (N_w/10^3)^2 / yr$ and duration of $\sim 1$ week, where $N_w$ is the number of the wind producing stars in the core. Presence of a massive black hole would enhance the stellar density around it and would make collisions of the winds in the core substantially more frequent. Collisions of the stellar winds in the cluster have also a number of interesting observable implications, including generation of $\gamma$-rays by particles accelerated by the shocks from the colliding winds. These processes are also expected to be relevant to compact regions of intense star formation elsewhere.

The effect of solar or stellar radiation on dust particles' trajectories (the Poynting-Robertson drag) has been studied by a number of authors and applied to interplanetary dust dynamics in numerical computations. Meanwhile some important features of dust flows can be studied analytically by implementing our novel hydrodynamical approach to use the continuity equation written in the particle's orbital elements as coordinates (Gor'kavyi, Ozernoy, & Mather 1997). By employing this approach and integrating the continuity equation, we are able to find two integrals of motion when the Poynting-Robertson drag dominates the dissipative forces in the dust flow. These integrals of motion enable us to explore basic characteristics of dust flows from any sources in the Solar system (such as asteroids, comets, Kuiper belt, etc.) or in another planetary system. In particular, we have reproduced the classical solution $n(r)\propto r^{-1}$ that approximately represents the overall distribution of dust in the Solar system. We have also investigated various factors that could be responsible for the deviations of the power law index in $n(r)\propto r^{\delta}$ from $\delta=-1$, including the influences of the orbital characteristics of dust sources, the evaporation of dust particles, as well as mixtures of dust particles of both asteroidal and cometary origin. We have calculated the masses and number densities of asteroidal and cometary components of the zodiacal cloud at different distances from the Sun.

We model the evolution of the selected types of X-ray sources (both with BHs and NSs) for the star formation burst at the Galactic center. For the currently assumed starburst age of 4-7 millions years, our results are consistent with the recent {\it GRANAT} X-ray observations of the Galactic center and predict a substantial number of X-ray binaries with black holes. In addition we made a very simple model of the spatial distribution of the X-ray sources, which is also consistent with the {\it GRANAT} observations. The results can be applied to other starformation regions.

X-ray binary population at the Galactic center is proposed here to form as a result of a starburst, which is synthesized by using the ``Scenario Machine'' code of binary star evolution. For the currently assumed starburst age of 4-7 millions years, our results are consistent with the recent {\it GRANAT} X-ray observations of the Galactic center and predict a substantial number of X-ray binaries with black holes.