Python 2D Jean Inestabilidad en simulación de galaxia espiral

Tengo ~ 1000 instantáneas de una simulación de galaxia espiral (desde perturbaciones de nubes moleculares t = 0 a t> 1Gyr. Lo que debo hacer es determinar si alguna posición en cualquiera de los tiempos de la instantánea cumple con los criterios de inestabilidad de Jeans .

La simulación es 2D, y el script de instantánea puede leer la posición, la velocidad y la densidad de la superficie ∑, pero no la temperatura. La simulación tampoco tiene gas, en cambio, la simulación se modeló dando a cada estrella individual una masa adicional que corresponde al gas. Miré en línea y varios sitios indican la masa de Jeans en 2d igual a c$_{s}$$^4$ / G * ∑. Pero no tengo idea de cómo determinar la velocidad del sonido c . Además, no estoy realmente interesado en la masa de Jeans, más aún si cualquier posición puede colapsar para formar estrellas, es decir, satisfacer la inestabilidad de Jeans.$^2$$_{s}$

Primero una advertencia, luego un truco:

1) Si no ejecutó esta simulación, si no es un astrónomo teórico capacitado de nivel universitario o superior, y si esta simulación no tiene el tipo de datos que necesita para responder la pregunta que está haciendo, no debe confiar los resultados en absoluto .

Se necesitan años de estudio para comprender los efectos numéricos extraños (cosas que suceden en la simulación debido a la forma en que está configurada que nunca sucedería en la realidad), la física de fondo y un montón de suposiciones y aproximaciones que se incluyen en este tipo de cosas.

A menudo decimos 'basura adentro, basura afuera'. Las respuestas que obtenga de esto probablemente no sean muy esclarecedoras porque dicen más sobre la simulación que el espacio real.

2) Dicho esto, si solo estás haciendo esto por diversión o para aprender las habilidades, ¡GENIAL! ¡Sigue jugando! Y una suposición muy básica que puede hacer para obtener algunos números (eso será bastante incorrecto) es que todo el gas en el espacio es de ~ 10K.