The Evolution of Mass Loaded Supernova Remnants. II. Temperature Dependent Mass Injection Rates

J.M. Pittard,S. J. Arthur,J.E. Dyson,S.A.E.G. Falle,T.W. Hartquist,M.I. Knight,M. Pexton

DOI: https://doi.org/10.1051/0004-6361%3A20030157

2003-03-12

Abstract:We investigate the evolution of spherically symmetric supernova remnants in which mass loading takes place due to conductively driven evaporation of embedded clouds. Numerical simulations reveal significant differences between the evolution of conductively mass loaded and the ablatively mass loaded remnants studied in Paper I. A main difference is the way in which conductive mass loading is extinguished at fairly early times, once the interior temperature of the remnant falls below ~10 million K. Thus, at late times remnants that ablatively mass load are dominated by loaded mass and thermal energy, while those that conductively mass load are dominated by swept-up mass and kinetic energy. Simple approximations to the remnant evolution, complementary to those in Paper I, are given.

Astrophysics

What problem does this paper attempt to address?

### Problems the paper attempts to solve

This paper aims to study the evolution of supernova remnants (SNR) under the assumption of spherical symmetry, especially when mass loading is taken into account. Specifically, the author focuses on the mass loading caused by conductive evaporation driven by heat conduction and compares it with the mass loading caused by hydrodynamic ablation studied in the previous paper.

#### Main problems
1. **Influence of mass loading mechanisms**: The paper explores the different influences of different mass loading mechanisms (such as conductive evaporation and hydrodynamic ablation) on the evolution of supernova remnants.
2. **Temperature - dependent mass injection rate**: It studies how temperature affects the mass injection rate, especially when the internal temperature of the remnant drops below about \(10^7\) K and conductive evaporation stops.
3. **Differences in remnant evolution**: It compares the different effects of conductive evaporation and hydrodynamic ablation on the structure and evolution path of the remnant, especially in the late stage. The former is mainly dominated by swept - up mass and kinetic energy, while the latter is dominated by loaded mass and thermal energy.

### Specific content
- **Physical background of conductive evaporation**: It is mentioned in the paper that in many astrophysical systems, the hot dilute phase and the cold dense phase coexist, and heat - conduction evaporation is an important process. However, this process may be suppressed by the magnetic field.
- **Numerical simulation results**: Through numerical simulation, the author finds that the mass loading of conductive evaporation is significant in the early stage, but as the internal temperature of the remnant drops below \(10^7\) K, this mass loading will rapidly weaken and finally stop. This is in contrast to the situation of hydrodynamic ablation, which continues throughout the evolution process.
- **Simplified approximation**: In order to better understand these complex behaviors, the author provides simplified approximate formulas to describe the evolution of the remnant, especially suitable for the situation in multi - phase media.

### Conclusion
The main contribution of this paper lies in revealing the importance and uniqueness of conductive evaporation as a mass loading mechanism in the evolution of supernova remnants. The research shows that conductive evaporation has a significant impact on the structure and evolution of the remnant in the early stage, but its role gradually weakens in the late stage, and swept - up mass and kinetic energy become the dominant factors. This forms a clear contrast with the remnants dominated by hydrodynamic ablation.

### Related formulas
Some of the key formulas involved in the paper are as follows:
- Evaporation rate formula:
\[
\dot{m} = 4\pi a^2 n T^{1/2} (k \mu m_H)^{1/2} \phi F(\sigma_0)
\]
where \(a\) is the radius of the cloud, \(n\) and \(T\) are the density and temperature of the surrounding hot phase respectively, \(\mu m_H\) is the average mass per particle, \(\phi\) is the efficiency factor, and \(F(\sigma_0)\) is a function defined in Cowie & McKee (1977).

- Saturation parameter:
\[
\sigma_0 = \left( \frac{T}{1.54\times 10^7} \right)^2 \frac{1}{n a_{pc} \phi}
\]
where \(a_{pc} = a/(1 \text{ pc})\) is the radius of the cloud in parsecs.

- Mass loading rate of conductive evaporation:
\[
q =
\begin{cases}
316 f' q_0 & \text{if } T > 10^7 \text{K} \\
f' q_0 \left( \frac{T}{T_6} \right)^{5/2} & \text{if } 10^5 \text{K} < T < 10^7 \text{K} \\
0 & \text{if } T < 10^5 \text{K}
\end{cases}
\]
where \(T_6 = 10^6 \text{K}\), \(q_0\)