Brilliantov et al. (PNAS, 2015) propose a model for the size distribution ~R^-3 for small particles with radius R and ~exp(-(R/Rc)^3) for large particles, where Rc=5.5 m. In 1989 Longaretti found analytically ~R^-3 for small particles and R^-6 for large ones (2). The law R^-6 also describes moonlets with size ~ 0.1-1 km. Cut-off law from Brilliantov et al. model does not describe moonlets and requires new mechanism for the origin of ~ 1 km size bodies. This model does not take into account the key effects of self-gravitation of large particles and differential rotation of the rings. Longaretti used a more accurate model of destruction of particles: "A relative increase of the erosion/destruction rate of the large particles must take place, because these particles have relative velocities of collision larger than the dispersion velocity, due to the differential Keplerian motion". The new model does not explain the difference between the rings and the satellites and the authors suggest that their calculations are applicable to the rings of the Chariklo and Chiron. In fact, the Chariklo rings are not examples of planetary rings, but proto-satellite disks. Simplified physical model that does not describe moonlets and the outer boundary of the rings makes the new model not a step forward but a step backwards in comparison with the papers of 20-30 years old.

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.

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.

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.