Related papers: YNOGK: A new public code for calculating null geod…
We consider the motion of test particles and light rays in a static cylindrically symmetric conformal spacetime given by Said et al [1]. We derive the equations of motion and present their analytical solutions in terms of the Weierstrass…
We derive alternate and new closed-form analytic solutions for the non-equatorial eccentric bound trajectories, $\{ \phi \left( r, \theta \right)$, $\ t \left( r, \theta \right),\ r \left( \theta \right) \}$, around a Kerr black hole by…
In the next few years, the near-infrared interferometer GRAVITY will be able to observe the Galactic center. Astrometric data will be obtained with an anticipated accuracy of 10 $\mu$as. To analyze these future data, we have developed a…
We present a novel numerical method for solving the elliptic partial differential equation problem for the electrostatic potential with piecewise constant conductivity. We employ an integral equation approach for which we derive a system of…
We revisit a little known theorem due to Beltrami, through which the integration of the geodesic equations of a curved manifold is accomplished by a method which, even if inspired by the Hamilton-Jacobi method, is purely geometric. The…
We present a first numerical implementation of a new scheme by Pound et al. that enables the calculation of the gravitational self-force in Kerr spacetime from a reconstructed metric-perturbation in a radiation gauge. The numerical task of…
This paper describes the Fortran 77 code SIMU, version 1.1, designed for numerical simulations of observational relations along the past null geodesic in the Lemaitre-Tolman-Bondi (LTB) spacetime. SIMU aims at finding scale invariant…
We derive the approximate analytical solutions of the bound timelike geodesic orbits in the effective-one-body (EOB) frame with extreme-mass ratio limit. The analytical solutions are expressed in terms of the elliptic integrals using Mino…
We calculate the transformation and inverse transformation, in the form of Taylor expansions, from arbitrary coordinates to Fermi-Walker coordinates in tubular neighborhoods of arbitrary timelike paths for general spacetimes. Explicit…
The null geodesics that describe photon orbits in the spacetime of a rotating electrically charged black hole (Kerr-Newman) are solved exactly including the contribution from the cosmological constant. We derive elegant closed form…
We derive energy-norm a posteriori error bounds for an Euler time-stepping method combined with various spatial discontinuous Galerkin schemes for linear parabolic problems. For accessibility, we address first the spatially semidiscrete…
Fully implicit Runge-Kutta (IRK) methods have many desirable accuracy and stability properties as time integration schemes, but high-order IRK methods are not commonly used in practice with large-scale numerical PDEs because of the…
The numerical integration of particle trajectories in curved spacetimes is fundamental for obtaining realistic models of the particle dynamics around massive compact objects such as black holes and neutron stars. Generalized algorithms…
Dirac's constraint analysis and the symplectic structure of geodesic equations are obtained for the general cylindrically symmetric stationary spacetime. For this metric, using the obtained first order Lagrangian, the geodesic equations of…
We explore the properties of test-particle orbits in "bumpy" spacetimes - stationary, reflection-symmetric, asymptotically flat solutions of Einstein equations that have a non-Kerr (anomalous) higher-order multipole-moment structure but can…
We study strong gravitational lensing in rotating space-times which can be thought of as realistic galactic models in General Relativity. To this end, using the Newman-Janis algorithm, we first obtain a rotating version of a static galactic…
The purpose of the present article is to show that the Newman-Janis and Newman et al algorithm used to derive the Kerr and Kerr-Newman metrics respectively, automatically leads to the extension of the initial non negative polar radial…
The overlap Dirac operator at nonzero quark chemical potential involves the computation of the sign function of a non-Hermitian matrix. In this talk we present an iterative method, first proposed by us in Ref. [1], which allows for an…
We consider a pointed curve $(X,P)$ which is given by the Weierstrass normal form, $y^r + A_{1}(x) y^{r-1} + A_{2}(x) y^{r-2} +\cdots + A_{r-1}(x) y + A_{r}(x)$ where $x$ is an affine coordinate on $\mathbb{P}^1$, the point $\infty$ on $X$…
Bearing in mind the application to the collapsar models of gamma-ray bursts (GRBs), we develop a numerical scheme and code for estimating the deposition of energy and momentum due to the neutrino pair annihilation ($\nu + {\bar \nu}…