Though the result obtained in the previous section is simple and in good agreement with the experimental result, the mechanism which determines the value of is not at all clear. In this section, we try to determine with more approximate approaches and to obtain qualitative understanding. First, we obtain the growth rate and growth mode of a typical planetesimal, using the collision probability obtained in section 2, and then derive the stationary distribution from them.
The growth rate of a planetesimal is obtained by integrating equation (5) over the range as Safronov1969,BargePellat1991 e.g.
By substituting equations (5) through (8) into equation (21) and neglecting the dependence on of and r and that on m of , we have
For the power-law mass distribution of equation (1), equation (22) is reduced to
where is the cut-off mass of the distribution function at the low-mass end. These two cases correspond to two different modes of accretion. If , the main mode of the growth is to eat small planetesimals, since the growth rate is determined by the low end of the integration over mass. If , on the other hand, the growth is driven by the collision with particles with similar masses. In other words, the growth is hierarchical.
Both the numerical results and the result obtained in the previous section suggest that . So we consider the case . The stationary state is realized if the mass flux does not depend on either the time or the mass, which is expressed as
Note that here the condition of the stationary distribution is the constant flux of the mass nm, and not the constant flux of the number of planetesimals n. This is because the dominant mode of the growth of the planetesimals is the hierarchical growth. The number of planetesimals would become smaller as they grow, since growth is driven by collisions of planetesimals of similar sizes.
If collisions are dominantly hierarchical, strictly speaking it does not make sense to use the growth equation expressed as the differential equation (21), since the change of the mass of planetesimals is not continuous. However, in this case we can use the mass doubling timescale
and consider the discrete mass bins with the ranges , , , ... , ... In this picture, collisions take place only between particles in the same mass bins. The product of the collision moves to the next mass bin.
In this case, the stationary state is realized when the incoming and outgoing flux for each mass bin is [is are] equal. For mass bin k, the total mass of planetesimals lost per unit time is
where the subscript k denotes the mass bin k and is the total mass of planetesimals in the mass bin k. Note that , since the width of the mass bin is proportional to m itself.
The incoming mass to bin k is what is lost from bin k-1. Thus, the condition for the stationary state is that the outgoing mass is the same for all mass bins,
Equations (25) through (27) leads to equation (24) as the condition for the stationary state. Thus, no matter [no matter whether] we assume the [delete ``the''] continuous growth of the mass or discontinuous jump[s] in the mass, the condition for the stationary state is the same.
From equation (1), (23), and (24), we obtain
This value is the same as what [delete ``what''] we obtained in section 2.
If we assume the [the a] constant flux of the number of planetesimals, we obtain , which is significantly smaller.
The reason why has a large negative value is now clear. As can be seen in equation (23), the growth rate of planetesimals has [a] strong positive dependence on their mass. This is, of course, why the runaway growth occurs (). [delete extra set of parentheses around Wetherill and Stewart (1989)] This strong dependence requires that the number of planetesimal goes down quickly as its mass increases, if a stationary distribution is to be realized. The strong dependence of the growth rate on the mass comes from the size, scale h[e]ight and the effect of gravitational focusing, all of them [them which] enhance the growth rate for larger mass.
From equations (23) and (29), we obtain
which implies the mass of a planetesimal can reach infinity in a finite time. Of course, this cannot happen in reality, since the number of planetesimals with this infinite mass would go below unity. In addition, at [a] certain point the most massive particle grows so massive that its effect dominates the velocity dispersion of less massive particles, resulting in the slowing down of the growth . This slowing down causes the massive planetesimals of similar mass to be formed in roughly equal orbital separation. Thus, these massive particles cannot collide with each other . Nevertheless, equation (30) implies that the stationary distribution of the mass for a wide range of mass[es] can be realized in a short time.
In the case of the two-dimensional space, we would obtain
However, with this formula, the total mass diverges if the high-mass end is taken infinite. Therefore, for the total mass to conserve, there must be a cut off [cut off cutoff] at the high-mass end, and no stationary solution is possible.
Note that the existence of the stationary solution is directly related to the occurrence of the runaway growth. The runaway implies that the heavy particles can become heavier by themselves, without being affected by lighter particles. Thus, the power-law behavior is observed because there is no characteristic scale. In the case of the two-dimensional simulation, the evolution of the heaviest particle is slower than that of lighter particles, and therefore all particles tend to have the [the a] mass similar to that of the heaviest particle. The power-law behavior is not observed in this case.