Related papers: Hyperboloidal discontinuous time-symmetric numeric…
In this paper we study the finite element approximation of systems of second-order nonlinear hyperbolic equations. The proposed numerical method combines a $hp$-version discontinuous Galerkin finite element approximation in the time…
We study the prospects of Gaussian processes (GP), a machine learning (ML) algorithm, as a tool to reconstruct the Hubble parameter $H(z)$ with two upcoming gravitational wave missions, namely the evolved Laser Interferometer Space Antenna…
We present a new method designed to avoid numerical challenges that have impeded calculation of the Lorenz gauge self-force acting on a compact object inspiraling into a Kerr black hole. This type of calculation is valuable in creating…
Local spectral similarity (LSS) algorithm has been developed for detecting homogeneous areas and edges in hyperspectral images (HSIs). The proposed algorithm transforms the 3-D data cube (within a spatial window) into a spectral similarity…
We put forward a new Bayesian modeling strategy for spatiotemporal count data that enables efficient posterior sampling. Most previous models for such data decompose logarithms of the response Poisson rates into fixed effects and spatial…
High order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution $\Delta x$) with only a moderate increase to computational expense. Significant…
We present an acceleration method for sequences of large-scale linear systems, such as the ones arising from the numerical solution of time-dependent partial differential equations coupled with algebraic constraints. We discuss different…
The kinetic energy of a local system of objects placed in a curved spacetime is gained by the subsequent acceleration of the object following the more contracted region of spacetime. Normally this happens near massive gravitating stars.…
The foundations are laid for the numerical computation of the actual worldline for a particle orbiting a black hole and emitting gravitational waves. The essential practicalities of this computation are here illustrated for a scalar…
We study linear gravitational perturbations of Schwarzschild spacetime by solving numerically Regge-Wheeler-Zerilli equations in time domain using hyperboloidal surfaces and a compactifying radial coordinate. We stress the importance of…
Supermassive black hole binaries (SMBHBs) are expected to emit continuous gravitational waves in the pulsar timing array (PTA) frequency band ($10^{-9}$--$10^{-7}$ Hz). The development of data analysis techniques aimed at efficient…
Self-force theory is the leading method of modeling extreme-mass-ratio inspirals (EMRIs), key sources for the gravitational-wave detector LISA. It is well known that for an accurate EMRI model, second-order self-force effects are critical,…
Spacetime discontinuous Galerkin (SDG) finite element methods are used to solve such PDEs involving space and time variables arising from wave propagation phenomena in important applications in science and engineering. To support an…
We present the algebraic foundations of the symmetric Zassenhaus algorithm and some of its variants. These algorithms have proven effective in devising higher-order methods for solving the time-dependent Schr\"{o}dinger equation in the…
The future space based gravitational wave detector LISA (Laser Interferometer Space Antenna) will observe millions of Galactic binaries constantly present in the data stream. A small fraction of this population (of the order of several…
Time-Delay Interferometry (TDI) is the data processing technique needed for generating interferometric combinations of data measured by the multiple Doppler readouts available onboard the three LISA spacecraft. Within the space of all…
We anticipate that the data acquired by the Laser Interferometer Space Antenna (LISA) will be dominated by the gravitational wave signals from several astrophysical populations. The analysis of these data is a new challenge and is the main…
High-accuracy gravitational-wave modeling demands going beyond linear, first-order perturbation theory. Particularly motivated by the need for second-order perturbative models of extreme-mass-ratio inspirals and black hole ringdowns, we…
The present study seeks to investigate mathematical structures of a multi-frequency subspace migration weighted by the natural logarithmic function for imaging of thin electromagnetic inhomogeneities from measured far-field pattern. To this…
Equal-arm interferometric detectors of gravitational radiation allow phase measurements many orders of magnitude below the intrinsic phase stability of the laser injecting light into their arms. This is because the noise in the laser light…