Related papers: $\texttt{Jipole}$: A Differentiable $\texttt{ipole…
With Advanced LIGO detecting the gravitational waves emitted from a pair of merging black holes in late 2015, we have a new perspective into the strong field regime of binary black hole systems. Event horizons are the defining features of…
Worldwide very long baseline radio interferometry arrays are expected to obtain horizon-scale images of supermassive black hole candidates as well as of relativistic jets in several nearby active galactic nuclei. This motivates the…
We compare the fully nonlinear and perturbative evolution of nonrotating black holes with odd-parity distortions utilizing the perturbative results to interpret the nonlinear results. This introduction of the second polarization…
This article outlines our derivation of the second order perturbations to a Schwarzschild black hole, highlighting our use of, and necessary reliance on, computer algebra. The particular perturbation scenario that is presented here is the…
Accretion disks, swirling structures of matter spiraling into black holes, play a pivotal role in our understanding of binary star systems and their intricate evolutionary processes. While current models often simplify these complex…
Ray tracing radiative transfer is a powerful method for comparing theoretical models of black hole accretion flows and jets with observations. We present a public code, grtrans, for carrying out such calculations in the Kerr metric,…
Detection and characterization of extended structures is a crucial goal in high contrast imaging. However, these structures face challenges in data reduction, leading to over-subtraction from speckles and self-subtraction with most existing…
We introduce the radiative transfer postprocessing code Subsweep. The code is based on the method of transport sweeps, in which the exact solution to the scattering-less radiative transfer equation is computed in a single pass through the…
The present article is concerned scattered data approximation for higher dimensional data sets which exhibit an anisotropic behavior in the different dimensions. Tailoring sparse polynomial interpolation to this specific situation, we…
The EHT collaboration released in 2019 the first horizon-scale images of a black hole accretion flow, opening a novel route for plasma physics comprehension and gravitational tests. Although the present unresolved images deeply depend on…
When numerically solving Einstein's equations for the evolution of binary black holes, physical imperfections in the initial data manifest as a transient, high-frequency pulse of ''junk radiation.'' This unphysical signal must be removed…
Differentiable rendering is a technique used in an important emerging class of visual computing applications that involves representing a 3D scene as a model that is trained from 2D images using gradient descent. Recent works (e.g. 3D…
We present results from a new code for binary black hole evolutions using the moving-puncture approach, implementing finite differences in generalised coordinates, and allowing the spacetime to be covered with multiple communicating…
Spherically symmetric (1D) black-hole spacetimes are considered as a test for numerical relativity. A finite difference code, based in the hyperbolic structure of Einstein's equations with the harmonic slicing condition is presented.…
We derive algorithms for higher order derivative computation of the rectangular $QR$ and eigenvalue decomposition of symmetric matrices with distinct eigenvalues in the forward and reverse mode of algorithmic differentiation (AD) using…
We present a new class of 3D black hole initial data sets for numerical relativity. These data sets go beyond the axisymmetric, ``gravity wave plus rotating black hole'' single black hole data sets by creating a dynamic, distorted hole with…
Black hole ringdowns in extensions of General Relativity (GR) generically exhibit two distinct signatures: (1) theory-dependent shifts in the standard black-hole quasinormal modes, and (2) additional modes arising from extra fundamental…
Polarized synchrotron emission is a fundamental process in high-energy astrophysics, particularly in the environments around black holes and pulsars. Accurate modeling of this emission requires precise computation of the emission,…
Ptychography in a reflection geometry shows great promise for non-destructive imaging of 3-dimensional nanostructures at the surface of a thick substrate. A major challenge to obtain high quality reflection-ptychographic images under…
We show how forward-mode automatic differentiation (AD) can be employed within larger reverse-mode computations to dynamically differentiate broadcast operations in a GPU-friendly manner. Our technique fully exploits the broadcast…