We applied our earlier numerical method to compute the thermal emission from a dusty resonant structure around Vega and compared the new observational evidence for a ring arc at 95 AU near Vega with our modeling. Our high resolution simulations indicate that Vega may have an outerm ost planet at a distance of 90-100 AU with coordinates near 18 35'16.25" (right ascension, B1950) and 38 44'15" (declination). Another symmetrical position, 18 35'14.75" and 38 44'32.5", of the planet is also possible. Radii of the observed clumps (and planet's orbital radius) may be smaller, e.g. ~60-75 AU. In this case the planet stays approximately on the line between the two possible positions above. From our modeling we make following conclusions: 1.Dust clumps rotate around Vega exactly wi th the planet orbital speed, 0.6-1.2 deg/year, as function of orbital radius of parent planet. 2. Orbital radius of the planet is approximately the same as the orbital radius of clumps (90-100 AU from observation by Koerner, Sargent & Ostroff, 2001 or 60-75 AU from observation by Wilner, Holman, Ho & Kuchner, 2002). 3. The most dense clump is ahead of the planet in the direction of the orbital motion. 4. The orbit of the planet may be circular or with a small eccentricity. 5. Absence of dust i n the inner part of Vega planetary system is a signature of inner planet(s) with orbital radii <50-60 AU.

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.

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.

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.

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.

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.

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.

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.