Neste trabalho implementamos melhorias no código BAM para simulações de sistemas astrofísicos no contexto da Relatividade Numérica. Tais melhorias consistem na introdução da infraestrutura computacional necessária ao emprego de equações de estado (EoSs) oriundas de teorias microscópicas para descrever de maneira realística a matéria estelar em regimes de altas densidades e temperatura finita. Além disso, mecanismos de interação fraca dominados por neutrinos são considerados através do chamado Neutrinos Leakage Scheme (NLS) e do formalismo M1 para as equações de transporte de radiação.
Com as novas implementações, realizamos simulações-teste com estrelas TOV empregando o NLS e reproduzimos com sucesso o colapso gravitacional induzido por neutrinos.
Também realizamos pela primeira vez em BAM simulações das etapas finais da colisão de sistemas binários de estrelas de nêutrons (BNS) utilizando o novo ferramental numérico implementado neste trabalho a fim de validar nossos métodos pela comparação com resultados reportados. Nossos resultados encontram-se publicados na Ref. [1].
A seguir, realizamos simulações explorando o espaço de parâmetros dos sistemas binários com diferentes razões de massa e spins iniciais, obtendo bom acordo com os resultados reportados na literatura em termos dos estados finais do sistema e propriedades da matéria ejetada.
Por fim, realizamos simulações-teste de BNSs com o esquema M1 para o tratamento dos neutrinos a fim de demonstrar a robustez de nossos métodos. Obtivemos evoluções numéricas estáveis, a despeito do substancial incremento de complexidade das simulações, e cujas propriedades preditas nas fases de pós-colisão são bastante convincentes. Num futuro breve realizaremos simulações de maior resolução com o esquema M1 com o intuito de estudar em detalhes colisões genéricas de BNSs.