Velocidad de caída en supernovas de colapso del núcleo

En este artículo Transporte de neutrinos en supernovas con colapso del núcleo , en la descripción del mecanismo de supernova con colapso del núcleo, se afirma que

La velocidad de la materia que cae en el núcleo aumenta a medida que nos alejamos del centro del núcleo.

Inicialmente pensé que esta velocidad es la velocidad de caída libre. Dado que el tiempo de caída libre varía en 1/sqrt (densidad) y la densidad disminuye a medida que nos movemos hacia afuera desde el centro del núcleo, el tiempo de caída libre será mayor y, por lo tanto, la velocidad de caída libre será menor.

Pero la velocidad de caída muestra una tendencia opuesta. ¿Qué establece esta 'velocidad de caída' y por qué aumenta a medida que nos alejamos del centro central?

¡Interesante pregunta! ¿Será que la materia que cae está luchando contra la presión de radiación de los neutrinos, etc...?
@ProfRob en algún otro texto encontré que esto sucede al comienzo del colapso del núcleo.
Creo que tienes razón y he estado bajo un malentendido.
por lo tanto, el tiempo de caída libre será mayor y, por lo tanto, la velocidad de caída libre será menor. ¿ Por qué? ¿La velocidad de caída libre es constante?

Respuestas (1)

Es posible que desee consultar "Estructura estelar y evolución" de Kippenhahn, Weigert y Weiss. El párrafo 36.3.1 "Soluciones simples de colapso" contiene justo lo que está buscando. En particular, la figura 36.5 muestra el perfil de velocidad en función de la coordenada de masa radial. He reproducido la figura aquí:

perfil de velocidad en función de la coordenada de masa radial

El razonamiento es más o menos el siguiente: puede modelar el núcleo colapsado como un gas degenerado relativista con ecuación de estado

PAG = k ρ 4 / 3

También necesitas la ecuación de continuidad de masa.

d metro d r = 4 π r 2 ρ
y la ecuación de Euler
d v r d t + 1 ρ PAG r + Φ r = 0

Si define un factor de escala a y una coordenada radial z tal que r = a ( t ) z , puede ver que el sistema de ecuaciones se desacopla en una ecuación para a ( t ) que se pueden integrar fácilmente, y una ecuación para ρ ( z ) :

1 z 2 d d z ( z 2 d w d z ) + w 3 = λ

Que es muy similar a la ecuación de Lane-Emden. Aquí w ( z ) se define tal que ρ ( z ) = ρ C w 3 ( z ) , como en Lane-Emden. El parámetro λ es una constante que mide de alguna manera la desviación del equilibrio hidrostático. De hecho, si λ = 0 , se obtiene la ecuación de equilibrio de Lane-Emden.

Al igual que la ecuación de Lane-Emden, esta ecuación admite soluciones sólo hasta una masa total METRO < METRO C , dónde METRO C no es mucho más alta que la masa de Chandrasekhar.

El resultado es que si la masa del núcleo que colapsa es mayor que METRO C , todo el núcleo no puede ser descrito por la solución de la ecuación. La parte interna del núcleo, hasta METRO C seguirá la ecuación esbozada y por lo tanto caerá homólogamente con v r = a ˙ a r : la velocidad radial es mayor cuanto más lejos del centro. La parte externa del núcleo, no descrita por la ecuación, estará casi en caída libre.

Además, la masa de Chandrasekhar es proporcional a m mi 2 . Pero los electrones son capturados durante el colapso, por lo tanto m mi aumenta y METRO C disminuye Esto significa que durante el colapso se encogerá la parte del núcleo que colapsa homólogamente.

Esto fue solo un resumen del proceso, si desea obtener más detalles, le recomiendo encarecidamente que eche un vistazo a Kippenhahn.

tl; dr: La parte interna del núcleo colapsa homólogamente con v r , mientras que la parte exterior está casi en caída libre con v r 1 / 2


Infall supersónico y algunas referencias más

Impulsado por el comentario de Lekha, he encontrado más referencias sobre el tema. Kippenhahn(1) basa su argumento en Goldreich(3) y Yahil(4)(5). La figura proviene en cambio de Müller(2).

Siguiendo a Yahil(4), dado que la velocidad de caída no puede exceder la velocidad de caída libre v F F 2 GRAMO metro r , la masa de la parte interna del núcleo se puede determinar encontrando dónde la velocidad de colapso homólogo v r es igual v F F . Cálculos numéricos muestran que en la parte exterior ρ r 3 = C o norte s t . Razonando a grandes rasgos, esto significa que metro ( r ) yo o gramo ( r ) y v F F ρ 1 / 6 . Pero también la velocidad del sonido va como C s = PAG ρ ρ 1 / 6 , por lo tanto v C s C o norte s t . De hecho, los resultados de la simulación muestran que la velocidad de caída del núcleo externo está constantemente alrededor de Mach 2. Tenga en cuenta que el modelo de Yahil no asume un politropo con γ = 4 3 y por lo tanto es más general que la aquí expuesta.

Referencias:

(1) Kippenhahn, R., Weigert, A. y Weiss, A., Stellar Structure and Evolution . 2012. págs. 462-465. doi:10.1007/978-3-642-30304-3.

(2) Müller, E. "Métodos computacionales para el flujo de fluidos astrofísicos, Curso avanzado 27 de Saas Fee". 1997. ed. por LeVeque, RJ, Mihalas, D., Dorfi, EA, Müller, E. (Springer, Berlin Heidelberg), pág. 343

(3) Goldreich, P. y Weber, SV, “Núcleos estelares que colapsan homólogamente”, The Astrophysical Journal , vol. 238, págs. 991–997, 1980. doi:10.1086/158065.

(4) Yahil, A. y Lattimer, JM, “Supernovas para peatones”, en Supernovas: A Survey of Current Research , 1982, vol. 90, págs. 53–70. ( No pude encontrar este documento, pero se hace referencia en el siguiente )

(5) Yahil, A., “Colapso estelar autosimilar”, The Astrophysical Journal , vol. 265, págs. 1047–1055, 1983. doi:10.1086/160746.

Entonces, 'La velocidad de la materia que cae en el núcleo aumenta a medida que nos alejamos del centro del núcleo'. debe ser cierto para el núcleo interno, ¿verdad? Pero el autor escribe: "Eventualmente, la velocidad de caída excede la velocidad del sonido local, es decir, la caída se vuelve supersónica". Entonces, ¿cómo se vuelve supersónica la velocidad de la parte exterior del núcleo, que está en caída libre?
@Lekha sí, la parte exterior del núcleo es supersónica, porque la densidad de la parte exterior es significativamente menor hasta el punto de que la velocidad del sonido es menor que la velocidad de caída libre. He buscado el tema y parece que estos resultados se derivan principalmente del cálculo numérico. Agregaré un párrafo y alguna referencia adicional a la respuesta.