Related papers: Improved fast-rotating black hole evolution simula…
Solving the 4-d Einstein equations as evolution in time requires solving equations of two types: the four elliptic initial data (constraint) equations, followed by the six second order evolution equations. Analytically the constraint…
Active galactic nucleus (AGN) feedback from accreting supermassive black holes (SMBHs) is an essential ingredient of galaxy formation simulations. The orbital evolution of SMBHs is affected by dynamical friction that cannot be predicted…
Binary black holes with spins that are aligned with the orbital angular momentum do not precess. However, post-Newtonian calculations predict that "up-down" binaries, in which the spin of the heavier (lighter) black hole is aligned…
We present techniques for successfully performing numerical relativity simulations of binary black holes with fourth-order accuracy. Our simulations are based on a new coding framework which currently supports higher order finite…
In this study, we examine the impacts of black hole spin, Lorentz-violating parameter, and the scalar field's mass on quasinormal modes (QNMs) of rotating Einstein-Bumblebee black holes, including computations up to the second-order…
We investigate the future evolution of Nariai black hole which is extremal limit of Schwarzschild-de Sitter one in modified gravity. The perturbations equations around Nariai black hole are derived in static and cosmological patches for…
Approximate solutions to the Einstein field equations are a valuable tool to investigate gravitational phenomena. An important aspect of any approximation is to investigate and quantify its regime of validity. We present a study that…
We present a full 3D numerical evolution code to study neutron stars in massive-scalar-tensor theories. The code is embedded in the Einstein Toolkit framework and its implementation constitutes a modified version of the…
This paper investigates the dynamic behavior of static, spherically symmetric black holes within the Einstein-Bumblebee gravity model with a cosmological constant, focusing on scalar field perturbations. Through separation of the angular…
We consider unconstrained evolution schemes for the hyperboloidal initial value problem in numerical relativity as a promising candidate for the optimally efficient numerical treatment of radiating compact objects. Here, spherical symmetry…
We present techniques for long-term, stable, and accurate evolutions of multiple-black-hole spacetimes using the `moving puncture' approach with fourth- and eighth-order finite difference stencils. We use these techniques to explore…
We investigate the evolution of black holes on orbits with small inclinations ($i < 2^\circ$) to the gaseous discs of active galactic nuclei. We perform 3D adiabatic hydrodynamic simulations within a shearing frame, studying the damping of…
We present long term evolutions of a single black hole of mass $M$ with the BSSN system using pseudospectral methods. For our simulations we use the SGRID code where the BSSN system is implemented in its standard second order in space form.…
We calculate the first dynamical evolutions of merging black hole-neutron star binaries that construct the combined black hole-neutron star spacetime in a general relativistic framework. We treat the metric in the conformal flatness…
We present cosmological-scale numerical simulations of an evolving universe in full general relativity (GR) and introduce a new numerical tool, {\sc CosmoGRaPH}, which employs the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formalism on a…
In this work we present 3D numerical relativity simulations of thick accretion disks around tilted Kerr BH. We investigate the evolution of three different initial disk models with a range of initial black hole spin magnitudes and tilt…
Three dimensional (3D) numerical evolutions of static black holes with excision are presented. These evolutions extend to about 8000M, where M is the mass of the black hole. This degree of stability is achieved by using growth-rate…
We present a new many-parameter family of hyperbolic representations of Einstein's equations, which we obtain by a straightforward generalization of previously known systems. We solve the resulting evolution equations numerically for a…
We construct a numerical relativity code based on the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formulation for the gravitational quadratic $f(R)$ Starobinsky model. By removing the assumption that the determinant of the conformal 3-metric…
The standard approach to the numerical evolution of black hole data using the ADM formulation with maximal slicing and vanishing shift is extended to non-symmetric black hole data containing black holes with linear momentum and spin by…