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.