Dynamical Collapse

In 1969, Larson (1969) and Penston (1969) found a self-similar solution which is suited for the dynamical contraction. Figure 4.8 is a radial density distribution for a spherical collapse of an isothermal cloud, where the cloud has a four-times larger mass than that of the Bonnor-Ebert mass. Although the figure is taken from a recent numerical study by Ogino et al (1999), a similar solution was obtained in Larson (1969). We can see that the solution has several characteristic points as follows:

- The cloud evolves in a self-similar way. That is, the spatial distribution of the density (left) at is well fitted by that at after shifting in the and the directions. As for the infall velocity spatial distribution, only a shift in the direction is needed.
- The density distribution in the envelope, which is fitted by , is almost unchanged. Only the central part of the cloud (high-density region) contracts.
- The time before the core formation epoch (the core formation time is defined as the time at which the central density increases greatly) is a good indicator to know how high the central density is. That is, reading from the figure, at () the central density reaches and at () the density is equal to . This shows the maximum (central) density is approximately proportional to , which is reasonable from the description of the free-fall time .

(4.63) | |||

(4.64) | |||

(4.65) |

where represents the mass included in the radius and denotes the mass of the protostar. A self-similar solution which has a form

should be found, where and are functions only on . For example, equation (4.66) asks the shape of the density distribution is the same after resizing of equation (4.69) and re-normalizing as . Since

(4.70) |

(4.71) |

Equations (4.73) and (4.74) have a singular point at which or . Since the point of const moves with , the flow velocity relative to this const is equal to . Thus the singular point at which corresponds to a sonic point. Therefore, since the flow has to pass the sonic point smoothly, the rhs of equations (4.73) and (4.74) have to be equal to zero at the singular point . This gives at the sonic pont ,

(4.75) |

(4.76) |

(4.77) |

These equations (4.72), (4.73) and (4.74) have an analytic solution

(4.78) |

This shows that for sufficiently large radius the gas flows with a constant inflow velocity .

This has a solution in which the density and the infall velocity should be regular with
reaching the center ().
Such kind of solution is plotted in Figure 4.9 (left) and
Figure 4.9 (right) with .
This time evolution is expected from the self-similar solution.
This shows that

(4.81) | |||

(4.82) |

Reaching the outer boundary the numerical solution (Fig.4.8) differs from the self-similar solution (Fig.4.9). For example, is reduced to zero in the numerical simulations, while it reaches a finite value 3.28 in the self-similar solution. And as for the density distribution, drops near the outer boundary in the numerical simulations while it decreases proportional to . However, in the region except for the vicinity of the outer boundary the self-similar solution expresses well the dynamical collapse of the spherical isothermal cloud. This solution gives the evolution of a pre-protostellar core formed in a supercritical cloud/cloud core.