Miren Radia

We present MHDuet , an open source evolution code for general relativistic magnetohydrodynamics with neutrino transport. The code solves the full set of Einstein equations coupled to a relativistic, magnetized fluid with an M1 neutrino radiation scheme using advanced techniques, including adaptive mesh and large eddy simulation techniques, to achieve high accuracy. The Simflowny platform generates the code from a high-level specification of the computational system, producing code that runs with either the structured adaptive mesh refinement application infrastructure (SAMRAI) or adaptive mesh refinement framework for exascale computing (AMReX) infrastructure. The choice of AMReX enables compilation and execution on graphical processing units, running an order of magnitude faster than on central processing units at the node level. We validate the code against benchmark tests, reproducing previous results obtained with the SAMRAI infrastructure, and demonstrate its capabilities with simulations of neutron stars employing realistic tabulated equations of state. Resolution studies clearly demonstrate convergence faster than second order in the grid spacing. Scaling tests reveal excellent strong and weak scaling performance when running on graphical processing units. The goal of the code is to provide a powerful tool for studying the dynamics of compact objects within multi-messenger astrophysics.

We present a generalisation of the curative initial data construction derived for equal-mass compact binaries in Helferet al(2019Phys. Rev.D99044046; 2022Class. Quantum Grav.39074001) to arbitrary mass ratios. We demonstrate how these improved initial data avoid substantial spurious artifacts in the collision dynamics of unequal-mass boson-star binaries in the same way as has previously been achieved with the simpler method restricted to the equal-mass case. We employ the improved initial data to explore in detail the impact of phase offsets in the coalescence of equal- and unequal-mass boson star binaries.

In this work we study the long-lived post-merger gravitational wave signature of a boson-star binary coalescence. We use full numerical relativity to simulate the post-merger and track the gravitational afterglow over an extended period of time. We implement recent innovations for the binary initial data, which significantly reduce spurious initial excitations of the scalar field profiles, as well as a measure for the angular momentum that allows us to track the total momentum of the spatial volume, including the curvature contribution. Crucially, we find the afterglow to last much longer than the spin-down timescale. This prolonged gravitational wave afterglow provides a characteristic signal that may distinguish it from other astrophysical sources.

We demonstrate the flexibility and utility of the Berger-Rigoutsos Adaptive Mesh Refinement (AMR) algorithm used in the open-source numerical relativity code GRChombo for generating gravitational waveforms from binary black-hole inspirals, and for studying other problems involving non-trivial matter configurations. We show that GRChombo can produce high quality binary black-hole waveforms through a code comparison with the established numerical relativity code Lean. We also discuss some of the technical challenges involved in making use of full AMR (as opposed to, e.g. moving box mesh refinement), including the numerical effects caused by using various refinement criteria when regridding. We suggest several "rules of thumb" for when to use different tagging criteria for simulating a variety of physical phenomena. We demonstrate the use of these different criteria through example evolutions of a scalar field theory. Finally, we also review the current status and general capabilities of GRChombo.

We demonstrate the flexibility and utility of the Berger–Rigoutsos adaptive mesh refinement (AMR) algorithm used in the open-source numerical relativity (NR) code GRChombo for generating gravitational waveforms from binary black-hole (BH) inspirals, and for studying other problems involving non-trivial matter configurations. We show that GRChombo can produce high quality binary BH waveforms through a code comparison with the established NR code Lean. We also discuss some of the technical challenges involved in making use of full AMR (as opposed to, e.g. moving box mesh refinement), including the numerical effects caused by using various refinement criteria when regridding. We suggest several ‘rules of thumb’ for when to use different tagging criteria for simulating a variety of physical phenomena. We demonstrate the use of these different criteria through example evolutions of a scalar field theory. Finally, we also review the current status and general capabilities of GRChombo.

Through numerical simulations of boson-star head-on collisions, we explore the quality of binary initial data obtained from the superposition of single-star spacetimes. Our results demonstrate that evolutions starting from a plain superposition of individual boosted boson-star spacetimes are vulnerable to significant unphysical artefacts. For equal-mass binaries, these difficulties can be overcome with a simple modification of the initial data suggested in Helferet al(2019Phys. Rev. D99044046) for collisions of oscillations. While we specifically consider massive complex scalar field boson star models of very high and low compactness, we conjecture that this vulnerability be also present in other kinds of exotic compact systems and hence needs to be addressed.