Entonces, me estoy interesando en el colapso gravitatorio de grandes nubes moleculares que se fragmentan y colapsan en múltiples estrellas. ¿Existen paquetes de simulación académica conocidos y populares que los simulen con fines de investigación? ¿O hay documentos o referencias que pueda usar para escribir los míos?
¡Esto suena muy ambicioso! Creo que su forma de abordar este tema podría ser echar un vistazo a las obras de dos autores: Matthew Bate y Mark Krumholz .
La razón por la que sugiero estos dos nombres es que son representativos de los dos enfoques básicos para simular el colapso de las nubes moleculares (que tiene poco que ver con las simulaciones de n-cuerpos). Estos enfoques se conocen como hidrodinámica de partículas suavizadas (SPH) y refinamiento de malla adaptable (AMR).
En pocas palabras, SPH trata el fluido como partículas discretas cuyas propiedades se suavizan sobre el tamaño de un núcleo; la resolución se controla mediante el número de partículas en la simulación y la longitud del suavizado. AMR es un enfoque más "tradicional" basado en cuadrículas en el que el tamaño de la cuadrícula varía de forma adaptativa con el tiempo y el espacio para obtener la resolución requerida.
Un rápido google revela una serie de códigos públicos, por ejemplo, para SPH hay SPHYSICS o GRLAB ; para AMR hay ENZO o PARAMESH .
Una de las principales actividades en las simulaciones de modelos en estos días es tratar de incorporar "retroalimentación" en la formación de estrellas. es decir, los vientos y la radiación de las estrellas de gran masa recién nacidas y los flujos de salida de las protoestrellas. Tales cosas no forman parte de los códigos estándar.