Inside-out Collapse Solution

In 1977 Shu found another self-similar solution which is realized after a central protostar with infinitesimal mass is formed in the singular isothermal sphere solution. The gas begins to accrete to the protostar. Outside the region where the accretion occurs, the initial SIS is kept unchanged, since the SIS is a hydrostatic solution. And the front of accretion expands radially outward in time. Since the inflow region expand outwardly, he called it the inside-out collapse solution. In Figure 4.10, the evolution is shown for density and inflow velocity. This solution gives the evolution of a protostellar core formed in a subcritical cloud/cloud core.

Using equations (4.79) and (4.80) and assuming the inflow velocity should reduce at large radius, we obtain

$$V(\xi)=-\frac{A-2}{\xi}-\frac{(A-2)(A-6)}{6\xi^3}+\cdots,$$ (4.83)

$$\Omega(\xi)=\frac{A}{\xi^2}-\frac{A(A-2)}{2\xi^4}+\cdots.$$ (4.84)

Since $\xi\rightarrow \infty$ means $t\rightarrow t_0$ (if $r$ is finite), $\Omega(\xi) \rightarrow {A}/{\xi^2}$ means that

$$\rho(r,t_0)=\frac{Ac_{\rm is}^2}{4\pi G r^2}.$$ (4.85)

Comparing with the SIS, when $A=2$ this gives the SIS and when $A>2$ this gives a density distribution in which the pressure is inefficient and the cloud is contracting. The solution with $A>2$ is obtained by a procedure as (1) at a sufficiently large radius $\xi_1$, calculate $V(\xi_1)$ and $\Omega(\xi_1)$. (2) from these values, integrate equations (4.73) and (4.74) inwardly. Figure 4.10 show the solution of this type. The solution with $A>2$ inflow speed is accelerated towards the center. Decreasing $A$ ( $A\rightarrow 2$), it is shown that an outer part $\xi\mbox{\raisebox{0.3ex}{$>$}\hspace{-1.1em} \raisebox{-0.7ex}{$\sim$}} 1$ reaches $V\rightarrow 0$. For $A=2^+$, the solution reaches the singular line $V=-\xi+1$ at $\xi=1$ ($V=0$) ^4.2. Since $V=0$ and $\Omega=2$ at $\xi=1$, this solution with $A=2^+$ converges to the SIS at $\xi=1$. This means that if there is an infinitesimally small amount of excess mass at the center of SIS, the accretion begins from the center while outside a radius the cloud is left static. The inner part of the solution $\xi \mbox{\raisebox{0.3ex}{$<$}\hspace{-1.1em} \raisebox{-0.7ex}{$\sim$}} 1$, $V$ and $\Omega$ are well expressed as $V \propto \xi^{1/2}$ and $\Omega\propto \xi^{3/2}$.

Figure 4.10: Self-similar solution which shows the inside-out collapse. (Left:) Infall velocity ( $-V=-v/c_{\rm is}$) is plotted against the similarity variable $\xi\equiv r/c_{\rm is}(t-t_0)$. Three curves correspond to models $A=2.4$, 2.1, and 2.01. (Right:) Density ( $\Omega =4\pi G \rho (t-t_0)^2$) is plotted against the similarity variable $\xi\equiv r/c_{\rm is}(t-t_0)$.

The power-law distributions of $V \propto \xi^{1/2}$ and $\Omega\propto \xi^{3/2}$ are explained as follows: Conservation of the total energy of an inflowing gas shell is expressed as

$$\frac{v^2}{2}-\frac{GM_{r}}{r}=-\frac{GM_{r_0}}{r_0},$$ (4.86)

where $r_0$ denotes the initial radius of a gas shell and $M_{r_0}=M_r$ represents the mass inside the gas shell. Neglecting $GM_{r_0}/r_0$ compared with the term $GM_{r}/r$, this gives

$$v\simeq \left( \frac{2GM_r}{r}\right)^{1/2}\propto r^{-1/2},$$ (4.87)

where we assumed a major part of $M_r$ comes from the mass of a protostar $M_*$, that is, $M_r=M_*+\int_0^r \rho 4\pi r^2 dr \simeq M_*$ Since the average density inside the radius $r$ $\bar{\rho}(<r)=\int_0^r \rho 4\pi r^2 dr /\int_0^r 4\pi r^2 dr =3 \rho(r)$ for SIS distribution, the time necessary for a gas shell to reach the center is proportional to $\propto \bar{\rho}(<r)^{-1/2}\propto \rho(r)^{-1/2}\propto r^1\propto (t-t_0)^1$, where we used a fact that the front of accretion expands with a constant speed $c_s{\rm is}$. This means that the time necessary for the gas shell to travel from the radius of the accretion wave front to the center is proportional to $t-t_0$. Since the mass of the shell which begins accretion in a unit time is equal to $\rho(r) 4\pi r^2 c_{\rm is}$ and is constant irrespective of $t-t_0$. These two facts indicate that the mass accretion rate is constant in time. That is,

$$\dot{M}=4\pi r^2 \rho v={\rm const}.$$ (4.88)

Since the mass swept by the sound wave per unit time is equal to $\dot{M}_{sw}= 4\pi r^2 \rho_{\rm SIS} c_s = 2c_s^3/G$, the accretion rate is proportional to $c_s^3/G$. Shu (1977) obtained the accretion rate

$$\dot{M}=0.975 \frac{c_s^3}{G},$$ (4.89)

for SIS $A=2^+$. Using this equation, equations (4.87) and (4.88) indicate that the spatial density distribution is expressed by a power-law as

$$\rho(r)\propto r^{-3/2},$$ (4.90)

which is valid for the region except for the vicinity of the front of the accretion.