BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Gaurav Khanna

Gaurav Khanna appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
gr-qcphysics.comp-phastro-ph.HEhep-thastro-phastro-ph.GAhep-phDistributed, Parallel, and Cluster Computingastro-ph.COEmerging TechnologiesPerformance

Catalog footprint

What is connected

38works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

38 published item(s)

preprint2026arXiv

Phenomenology and origin of late-time tails in eccentric binary black hole mergers

We investigate the late-time tail behavior in gravitational waves from merging eccentric binary black holes (BBH) using black hole perturbation theory. For simplicity, we focus only on the dominant quadrupolar mode of the radiation. We demonstrate that such tails become more prominent as eccentricity increases. Exploring the phenomenology of the tails in both spinning and non-spinning eccentric binaries, with the spin magnitude varying from $χ=-0.6$ to $χ=+0.6$ and eccentricity as high as $e=0.98$, we find that these tails can be well approximated by a slowly decaying power law. We study the power law for varying systems and find that the power law exponent lies close to the theoretically expected value $-4$. Finally, using both plunge geodesic and radiation-reaction-driven orbits, we perform a series of numerical experiments to understand the origin of the tails in BBH simulations. Our results suggest that the late-time tails are strongly excited in eccentric BBH systems when the smaller black hole is in the neighborhood of the apocenter, as opposed to any structure in the strong field of the larger black hole. Our analysis framework is publicly available through the \texttt{gwtails} Python package.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Source-Driven Tails in Kerr Spacetime: Nonlinear effects in Late-Time Behavior

We present the long-duration time-domain simulations of scalar-field tails in Kerr spacetimes driven by \emph{outgoing} multipolar sources. Extending the recent work in the literature from Schwarzschild to rotating black holes, we evolve sources with $\ell'=\{0,1,2,3,4\}$ on backgrounds with dimensionless spin $a/M=\{0.0, 0.8, 1.0\}$ and extract the late-time decay rates of measured modes $\ell\le4$ for a nonlinearity-inspired outgoing source with a $1/r^2$ fall-off. In all cases we find the inverse power-law index $p_{\ell\ell'}$ to be larger than the source-free Price law values by one unit, i.e. $p^{\text{sourced}}_{\ell\ell'} = p^{\text{Price}}_{\ell\ell'} + 1$. We also include a power-law index value computation for a similar source-driven gravitational wave case $(\ell,m)=(4,4)$ and confirm closely related results in the recent literature.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Measuring quasi-normal mode amplitudes with misaligned binary black hole ringdowns

In recent work, we examined how different modes in the ringdown phase of a binary coalescence are excited as a function of the final plunge geometry. At least in the large mass ratio limit, we found a clean mapping between angles describing the plunge and the amplitude of different quasi-normal modes (QNMs) which constitute the ringdown. In this study, we use that mapping to construct a waveform model expressed as a sum of QNMs where the mode amplitudes and phases are determined by the source plunge parameters. We first generate a large number of calibration waveforms and interpolate between fits of each mode amplitude and phase up to $\ell \leq 8$ and $\ell - |m| \leq 4$. The density of our calibration data allows us to resolve important features such as phase transition discontinuities at large misalignments. Using our ringdown waveform model, we then perform Bayesian parameter estimation with added white Gaussian noise to demonstrate that, in principle, the mode amplitudes can be measured and used to constrain the plunge geometry. We find that inferences are substantially improved by incorporating prior information constraining mode excitation, which motivates work to understand and characterize how the QNM excitation depends on the coalescence geometry. These results are part of a broader effort to map the mode excitation from arbitrary masses and spins, which will be useful for characterizing ringdown waves in upcoming gravitational-wave measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Scalar and Gravitational Transient "Hair" for Near-Extremal Black Holes

We study the existence and nature of Aretakis "hair" and its potentially observable imprint at a finite distance from the horizon (Ori-coefficient) in near-extremal black hole backgrounds. Specifically, we consider the time evolution of horizon penetrating scalar and gravitational perturbations with compact support on near-extremal Reissner-Nordstrom (NERN) and Kerr (NEK). We do this by numerically solving the Teukolsky equation and determining the Aretakis charge values on the horizon and at a finite distance from the black hole. We demonstrate that these values are no longer strictly conserved in the non-extremal case; however, their decay rate can be arbitrarily slow as the black hole approaches extremality allowing for the possibility of their observation as a transient hair.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Surrogate Model for Gravitational Wave Signals from Comparable- to Large- Mass-Ratio Black Hole Binaries

Gravitational wave signals from compact astrophysical sources such as those observed by LIGO and Virgo require a high-accuracy, theory-based waveform model for the analysis of the recorded signal. Current inspiral-merger-ringdown models are calibrated only up to moderate mass ratios, thereby limiting their applicability to signals from high-mass ratio binary systems. We present EMRISur1dq1e4, a reduced-order surrogate model for gravitational waveforms of 13,500M in duration and including several harmonic modes for non-spinning black hole binary systems with mass-ratios varying from 3 to 10,000 thus vastly expanding the parameter range beyond the current models. This surrogate model is trained on waveform data generated by point-particle black hole perturbation theory (ppBHPT) both for large mass-ratio and comparable mass-ratio binaries. We observe that the gravitational waveforms generated through a simple application of ppBHPT to the comparable mass-ratio cases agree remarkably (and surprisingly) well with those from full numerical relativity after a rescaling of the ppBHPT's total mass parameter. This observation and the EMRISur1dq1e4 surrogate model will enable data analysis studies in the high-mass ratio regime, including potential intermediate mass-ratio signals from LIGO/Virgo and extreme-mass ratio events of interest to the future space-based observatory LISA.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Basis Function Method for Numerical Loop Quantum Cosmology: The Schwarzschild Black Hole Interior

Loop quantum cosmology is a symmetry-reduced application of loop quantum gravity that has led to the resolution of classical singularities such as the big bang, and those at the center of black holes. This can be seen through numerical simulations involving the quantum Hamiltonian constraint that is a partial difference equation. The equation allows one to study the evolution of sharply-peaked Gaussian wave packets that generically exhibit a quantum "bounce" or a non-singular passage through the classical singularity, thus offering complete singularity resolution. In addition, von-Neumann stability analysis of the difference equation -- treated as a stencil for a numerical solution that steps through the triad variables -- yields useful constraints on the model and the allowed space of states. In this paper, we develop a new method for the numerical solution of loop quantum cosmology models using a set of basis functions that offer a number of advantages over computing a solution by stepping through the triad variables. We use the Corichi and Singh model for the Schwarzschild interior as the main case study in this effort. The main advantage of this new method is computational efficiency and the ease of parallelization. In addition, we also discuss how the stability analysis appears in the context of this new approach.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Arrival Times of Gravitational Radiation Peaks for Binary Inspiral

Modeling of gravitational waves (GWs) from binary black hole inspiral brings together early post-Newtonian waveforms and late quasinormal ringing waveforms. Attempts to bridge the two limits without recourse to numerical relativity involve predicting the time of the peak GW amplitude. This prediction will require solving the question of why the peak of the "source," i.e., the peak of the binary angular velocity, does not correspond to the peak of the GW amplitude. We show here that this offset can be understood as due to the existence two distinct components of the radiation: the "direct" radiation analogous to that in flat spacetime, and "scattered" radiation associated with curved spacetime. The time dependence of these two components, and of their relative phases determines the location of the peak amplitude. We use a highly simplified model to clarify the twocomponent nature of the source, then demonstrate that the explanation is valid also for an extreme mass ratio binary inspiral.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Black hole binary inspiral: Analysis of the plunge

Binary black hole coalescence has its peak of gravitational wave generation during the "plunge," the transition from quasicircular early motion to late quasinormal ringing. Although advances in numerical relativity have provided plunge waveforms, there is still no intuitive or phenomenological understanding of plungecomparable to that of the early and late stages. Here we make progress in developing such understanding by focusing on the excitation of quasinormal ringing (QNR) during the plunge. We rely on insights of the linear mathematics of the particle perturbation model for the extreme mass limit. Our analysis, based on the Fourier domain Green function, and a simple initial model, point to the crucial role played by the kinematics near the "light ring" (the circular photon orbit) in determining the excitation of QNR. That insight is then shown to successfully explain Schwarzschild QNR found with evolution codes. Lastly, a phenomenological explanation is given for the underlying importance of the light ring.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gravitational waves from a plunge into a nearly extremal Kerr black hole

We study numerically in the time domain the linearized gravitational waves emitted from a plunge into a nearly extremal Kerr black hole by solving the inhomogeneous Teukolsky equation. We consider spinning black holes for which the specific spin angular momentum $a/M=1-ε$, and we consider values of $ε\geq 10^{-6}$. We find an effective transient behavior for the quasi-normal ringdown: the early phase of the quasi-normal ringdown is governed by a decay according to inverse time, with frequency equaling twice the black hole's horizon frequency. The smaller $ε$ the later the transition from this transient inverse time decay to exponential decay. Such sources, if exist, may be interesting potential sources for terrestrial or space borne gravitational wave observatories.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Scientific Computing Using Consumer Video-Gaming Hardware Devices

Commodity video-gaming hardware (consoles, graphics cards, tablets, etc.) performance has been advancing at a rapid pace owing to strong consumer demand and stiff market competition. Gaming hardware devices are currently amongst the most powerful and cost-effective computational technologies available in quantity. In this article, we evaluate a sample of current generation video-gaming hardware devices for scientific computing and compare their performance with specialized supercomputing general purpose graphics processing units (GPGPUs). We use the OpenCL SHOC benchmark suite, which is a measure of the performance of compute hardware on various different scientific application kernels, and also a popular public distributed computing application, Einstein@Home in the field of gravitational physics for the purposes of this evaluation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Black holes in anti-de Sitter: quasinormal modes, tails and tales of flat spacetime

Black holes in asymptotically anti-de Sitter (AdS) spacetimes have been the subject of intense scrutiny, including detailed frequency-domain analysis and full nonlinear evolutions. Remarkably, studies of linearized perturbations in the time-domain are scarce or non-existing. We close this gap by evolving linearized scalar wavepackets in the background of rotating BHs in AdS spacetimes. Our results show a number of interesting features. Small BHs in AdS behave as asymptotically flat BHs for early/intermediate times, displaying the same ringdown modes and power-law tails. As the field bounces back and forth between the horizon and the timelike boundary it "thermalizes" and the modes of AdS settle in. Finally, we have indications that wavepackets in the vicinity of fastly spinning BHs grow exponentially in time, signalling a superradiant instability of the geometry previously reported through a frequency-domain analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Brief note on high-multipole Kerr tails

In this note we reconsider the late-time, power-law decay rate of scalar fields in a Kerr space-time background. We implement a number of mathematical and computational enhancements to our time-domain, (2+1)D Teukolsky evolution code and are able to obtain reliable decay rates for multipoles as high as $\ell=16$. Our numerical results suggest full agreement with the proposed decay expressions in recent work for both the finite distance, and null infinity cases. We also study the same in the context of extremal Kerr space-time and find that the same results hold, except that the horizon tails follow the null infinity expressions instead.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Self-force gravitational waveforms for extreme and intermediate mass ratio inspirals. III: Spin-orbit coupling revisited

The first- and second-order dissipative self force and the first order conservative self force are applied together with spin-orbit coupling to the quasi-circular motion of a test mass in the spacetime of a Schwarzschild black hole, for extreme or intermediate mass ratios. The partial dephasing of the gravitational waveform (at the order that is independent of the system's mass ratio) due to the self force is compared with that of spin-orbit coupling. We find that accurate waveforms for parameter estimation need to include both effects. Specifically, we find a particular value for the spin parameter such that the spin-orbit effect cancels out the self-force effect on the waveform. Exclusion of dephasing effects that are independent of the mass ratio therefore might lead to a non-perturbative error in the estimation of the system's parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Intermediate behavior of Kerr tails

The numerical investigation of wave propagation in the asymptotic domain of Kerr spacetime has only recently been possible thanks to the construction of suitable hyperboloidal coordinates. The asymptotics revealed an apparent puzzle in the decay rates of scalar fields: the late-time rates seemed to depend on whether finite distance observers are in the strong field domain or far away from the rotating black hole, an apparent phenomenon dubbed "splitting". We discuss far-field "splitting" in the full field and near-horizon "splitting" in certain projected modes using horizon-penetrating, hyperboloidal coordinates. For either case we propose an explanation to the cause of the "splitting" behavior, and we determine uniquely decay rates that previous studies found to be ambiguous or immeasurable. The far-field "splitting" is explained by competition between projected modes. The near-horizon "splitting" is due to excitation of lower multipole modes that back excite the multipole mode for which "splitting" is observed. In both cases "splitting" is an intermediate effect, such that asymptotically in time strong field rates are valid at all finite distances. At any finite time, however, there are three domains with different decay rates whose boundaries move outwards during evolution. We then propose a formula for the decay rate of tails that takes into account the inter--mode excitation effect that we study.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Small mass plunging into a Kerr black hole: Anatomy of the inspiral-merger-ringdown waveforms

We numerically solve the Teukolsky equation in the time domain to obtain the gravitational-wave emission of a small mass inspiraling and plunging into the equatorial plane of a Kerr black hole. We account for the dissipation of orbital energy using the Teukolsky frequency-domain gravitational-wave fluxes for circular, equatorial orbits, down to the light-ring. We consider Kerr spins $-0.99 \leq q \leq 0.99$, and compute the inspiral-merger-ringdown (2,2), (2,1), (3,3), (3,2), (4,4), and (5,5) modes. We study the large-spin regime, and find a great simplicity in the merger waveforms, thanks to the extremely circular character of the plunging orbits. We also quantitatively examine the mixing of quasinormal modes during the ringdown, which induces complicated amplitude and frequency modulations in the waveforms. Finally, we explain how the study of small mass-ratio black-hole binaries helps extending effective-one-body models for comparable-mass, spinning black-hole binaries to any mass ratio and spin magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Black hole binary inspiral and trajectory dominance

Gravitational waves emitted during the inspiral, plunge and merger of a black hole binary carry linear momentum. This results in an astrophysically important recoil to the final merged black hole, a ``kick'' that can eject it from the nucleus of a galaxy. In a previous paper we showed that the puzzling partial cancellation of an early kick by a late antikick, and the dependence of the cancellation on black hole spin, can be understood from the phenomenology of the linear momentum waveforms. Here we connect that phenomenology to its underlying cause, the spin-dependence of the inspiral trajectories. This insight suggests that the details of plunge can be understood more broadly with a focus on inspiral trajectories.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Mode coupling mechanism for late-time Kerr tails

We consider the decay rate for scalar fields in Kerr spacetime. We consider pure initial (azimuthal) multipoles $\ell'$ with respect to the class which includes Boyer-Lindquist coordinates, and focus attention on the decay rate of the multipole $\ell$. We use an iterative method proposed by Gleiser, Price, and Pullin, and identify the mode coupling mechanism through the iterations in powers of the square of the Kerr black hole's specific angular momentum that gives rise to a decay rate formula recently proposed by Zenginoğlu, Khanna, and Burko. Modes $\ell$ may be excited through different channels, each leading to its own decay rate. The asymptotic decay rate of the mode $\ell$ is the slowest of the decay rate of the various channels. In some cases, more than one channel leads to the same decay rate, and then the amplitude of the mode is the sum of the amplitudes of the partial fields generated by the individual channels. We also show that one may identify the asymptotically-dominant channel of mode excitations, and obtain approximate results for the mode of interest by studying the dominant channel. The results of the dominant channel approximation approach the full-mode results at late times, and their difference approaches zero quadratically in inverse time.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modeling the horizon-absorbed gravitational flux for equatorial-circular orbits in Kerr spacetime

We propose an improved analytical model for the horizon-absorbed gravitational-wave energy flux of a small body in circular orbit in the equatorial plane of a Kerr black hole. Post-Newtonian (PN) theory provides an analytical description of the multipolar components of the absorption flux through Taylor expansions in the orbital frequency. Building on previous work, we construct a mode-by-mode factorization of the absorbed flux whose Taylor expansion agrees with current PN results. This factorized form significantly improves the agreement with numerical results obtained with a frequency-domain Teukolsky code, which evolves through a sequence of circular orbits up to the photon orbit. We perform the comparison between model and numerical data for dimensionless Kerr spins $-0.99 \leq q \leq 0.99$ and for frequencies up to the light ring of the Kerr black hole. Our proposed model enforces the presence of a zero in the flux at an orbital frequency equal to the frequency of the horizon, as predicted by perturbation theory. It also reproduces the expected divergence of the flux close to the light ring. Neither of these features are captured by the Taylor-expanded PN flux. Our proposed absorption flux can also help improve models for the inspiral, merger, ringdown of small mass-ratio binary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Self-force gravitational waveforms for extreme and intermediate mass ratio inspirals. II: Importance of the second-order dissipative effect

We consider the importance of the second-order dissipative self force for gravitational wave dephasing for an extreme or intermediate mass ratio system moving along a quasi-circular Schwarzschild orbit. For the first-order self force we use the fully relativistic force in the Lorenz gauge for eternally circular geodesics. The second-order self force is modeled by its 3.5 post Newtonian counterpart. We evolve the system using the osculating orbits method, and obtain the gravitational waveforms, whose phase includes all the terms - within our approximation (and using the self force along circular geodesics) - that are independent of the system's mass ratio. The partial dephasing due to the second-order dissipative self force is substantially smaller than that of the first-order conservative self force, although they are both at the same order in the mass ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

High Accuracy Gravitational Waveforms from Black Hole Binary Inspirals Using OpenCL

There is a strong need for high-accuracy and efficient modeling of extreme-mass-ratio binary black hole systems because these are strong sources of gravitational waves that would be detected by future observatories. In this article, we present sample results from our Teukolsky EMRI code: a time-domain Teukolsky equation solver (a linear, hyperbolic, partial differential equation solver using finite-differencing), that takes advantage of several mathematical and computational enhancements to efficiently generate long-duration and high-accuracy EMRI waveforms. We emphasize here the computational advances made in the context of this code. Currently there is considerable interest in making use of many-core processor architectures, such as Nvidia and AMD graphics processing units (GPUs) for scientific computing. Our code uses the Open Computing Language (OpenCL) for taking advantage of the massive parallelism offered by modern GPU architectures. We present the performance of our Teukolsky EMRI code on multiple modern processor architectures and demonstrate the high level of accuracy and performance it is able to achieve. We also present the code's scaling performance on a large supercomputer i.e. NSF's XSEDE resource: Keeneland (a 201 TeraFLOP/s, 120-node HP SL390 system with 240 Intel Xeon 5660 CPUs and 360 NVIDIA Fermi M2070 graphics processors, with the nodes connected by an InfiniBand QDR network).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Modeling multipolar gravitational-wave emission from small mass-ratio mergers

Using the effective-one-body (EOB) formalism and a time-domain Teukolsky code, we generate inspiral, merger, and ringdown waveforms in the small mass-ratio limit. We use EOB inspiral and plunge trajectories to build the Teukolsky equation source term, and compute full coalescence waveforms for a range of black hole spins. By comparing EOB waveforms that were recently developed for comparable mass binary black holes to these Teukolsky waveforms, we improve the EOB model for the (2,2), (2,1), (3,3), and (4,4) modes. Our results can be used to quickly and accurately extract useful information about merger waves for binaries with spin, and should be useful for improving analytic models of such binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical Techniques in Loop Quantum Cosmology

In this article, we review the use of numerical techniques to obtain solutions for the quantum Hamiltonian constraint in loop quantum cosmology (LQC). First, we summarize the basic features of LQC, and describe features of the constraint equations to solve - generically, these are difference (rather than differential) equations. Important issues such as differing quantization methods, stability of the solutions, the semi-classical limit, and the relevance of lattice refinement in the difference equations are discussed. Finally, the cosmological models already considered in the literature are listed, along with typical features in these models and open issues.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A multi-domain hybrid method for head-on collision of black holes in particle limit

A hybrid method is developed based on the spectral and finite-difference methods for solving the inhomogeneous Zerilli equation in time-domain. The developed hybrid method decomposes the domain into the spectral and finite-difference domains. The singular source term is located in the spectral domain while the solution in the region without the singular term is approximated by the higher-order finite-difference method. The spectral domain is also split into multi-domains and the finite-difference domain is placed as the boundary domain. Due to the global nature of the spectral method, a multi-domain method composed of the spectral domains only does not yield the proper power-law decay unless the range of the computational domain is large. The finite-difference domain helps reduce boundary effects due to the truncation of the computational domain. The multi-domain approach with the finite-difference boundary domain method reduces the computational costs significantly and also yields the proper power-law decay. Stable and accurate interface conditions between the finite-difference and spectral domains and the spectral and spectral domains are derived. For the singular source term, we use both the Gaussian model with various values of full width at half maximum and a localized discrete $δ$-function. The discrete $δ$-function was generalized to adopt the Gauss-Lobatto collocation points of the spectral domain. The gravitational waveforms are measured. Numerical results show that the developed hybrid method accurately yields the quasi-normal modes and the power-law decay profile. The numerical results also show that the power-law decay profile is less sensitive to the shape of the regularized $δ$-function for the Gaussian model than expected. The Gaussian model also yields better results than the localized discrete $δ$-function.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An Exploration of OpenCL for a Numerical Relativity Application

Currently there is considerable interest in making use of many-core processor architectures, such as Nvidia and AMD graphics processing units (GPUs) for scientific computing. In this work we explore the use of the Open Computing Language (OpenCL) for a typical Numerical Relativity application: a time-domain Teukolsky equation solver (a linear, hyperbolic, partial differential equation solver using finite-differencing). OpenCL is the only vendor-agnostic and multi-platform parallel computing framework that has been adopted by all major processor vendors. Therefore, it allows us to write portable source-code and run it on a wide variety of compute hardware and perform meaningful comparisons. The outcome of our experimentation suggests that it is relatively straightforward to obtain order-of-magnitude gains in overall application performance by making use of many-core GPUs over multi-core CPUs and this fact is largely independent of the specific hardware architecture and vendor. We also observe that a single high-end GPU can match the performance of a small-sized, message-passing based CPU cluster.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gravitational Wave Emission from the Single-Degenerate Channel of Type Ia Supernovae

The thermonuclear explosion of a C/O white dwarf as a Type Ia supernova (SN Ia) generates a kinetic energy comparable to that released by a massive star during a SN II event. Current observations and theoretical models have established that SNe Ia are asymmetric, and therefore--like SNe II--potential sources of gravitational wave (GW) radiation. We perform the first detailed calculations of the GW emission for a SN Ia of any type within the single-degenerate channel. The gravitationally-confined detonation (GCD) mechanism predicts a strongly-polarized GW burst in the frequency band around 1 Hz. Third-generation spaceborne GW observatories currently in planning may be able to detect this predicted signal from SNe Ia at distances up to 1 Mpc. If observable, GWs may offer a direct probe into the first few seconds of the SNe Ia detonation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Null infinity waveforms from extreme-mass-ratio inspirals in Kerr spacetime

We describe the hyperboloidal compactification for Teukolsky equations in Kerr spacetime. We include null infinity on the numerical grid by attaching a hyperboloidal layer to a compact domain surrounding the rotating black hole and the orbit of an inspiralling point particle. This technique allows us to study, for the first time, gravitational waveforms from large- and extreme-mass-ratio inspirals in Kerr spacetime extracted at null infinity. Tests and comparisons of our results with previous calculations establish the accuracy and efficiency of the hyperboloidal layer method.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical modeling of gravitational wave sources accelerated by OpenCL

In this work, we make use of the OpenCL framework to accelerate an EMRI modeling application using the hardware accelerators -- Cell BE and Tesla CUDA GPU. We describe these compute technologies and our parallelization approach in detail, present our performance results, and then compare them with those from our previous implementations based on the native CUDA and Cell SDKs. The OpenCL framework allows us to execute identical source-code on both architectures and yet obtain strong performance gains that are comparable to what can be derived from the native SDKs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Systematics of black hole binary inspiral kicks and the slowness approximation

During the inspiral and merger of black holes, the interaction of gravitational wave multipoles carries linear momentum away, thereby providing an astrophysically important recoil, or "kick" to the system and to the final black hole remnant. It has been found that linear momentum during the last stage (quasinormal ringing) of the collapse tends to provide an "antikick" that in some cases cancels almost all the kick from the earlier (quasicircular inspiral) emission. We show here that this cancellation is not due to peculiarities of gravitational waves, black holes, or interacting multipoles, but simply to the fact that the rotating flux of momentum changes its intensity slowly. We show furthermore that an understanding of the systematics of the emission allows good estimates of the net kick for numerical simulations started at fairly late times, and is useful for understanding qualitatively what kinds of systems provide large and small net kicks.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Testing the Cosmic Censorship Conjecture with point particles: the effect of radiation reaction and the self-force

A classical thought-experiment to destroy black holes was envisaged by Wald in 1974: it consists of throwing particles with large angular momentum into an extremal black hole, checking whether their capture can over-spin the black hole past the extremal limit and create a naked singularity. Wald showed that in the test-particle limit, particles that would be otherwise capable of producing naked singularities are simply scattered. Recently Jacobson and Sotiriou showed that if one considers instead a black hole that is almost, but not exactly extremal, then in the absence of backreaction effects particle capture could indeed over-spin the spacetime above the Kerr limit. Here we analyze back-reaction effects and show that for some of the trajectories giving rise to naked singularities, radiative effects can be neglected. However, for these orbits the conservative self-force is important, and seems to have the right sign to prevent the formation of naked singularities.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Binary black hole merger gravitational waves and recoil in the large mass ratio limit

Spectacular breakthroughs in numerical relativity now make it possible to compute spacetime dynamics in almost complete generality, allowing us to model the coalescence and merger of binary black holes with essentially no approximations. The primary limitation of these calculations is now computational. In particular, it is difficult to model systems with large mass ratio and large spins, since one must accurately resolve the multiple lengthscales which play a role in such systems. Perturbation theory can play an important role in extending the reach of computational modeling for binary systems. In this paper, we present first results of a code which allows us to model the gravitational waves generated by the inspiral, merger, and ringdown of a binary system in which one member of the binary is much more massive than the other. This allows us to accurately calibrate binary dynamics in the large mass ratio regime. We focus in this analysis on the recoil imparted to the merged remnant by these waves. We closely examine the "antikick", an anti-phase cancellation of the recoil arising from the plunge and ringdown waves, described in detail by Schnittman et al. We find that, for orbits aligned with the black hole spin, the antikick grows as a function of spin. The total recoil is smallest for prograde coalescence into a rapidly rotating black hole, and largest for retrograde coalescence. Amusingly, this completely reverses the predicted trend for kick versus spin from analyses that only include inspiral information.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Late-time Kerr tails: generic and non-generic initial data sets, "up" modes, and superposition

Three interrelated questions concerning Kerr spacetime late-time scalar-field tails are considered numerically, specifically the evolutions of generic and non-generic initial data sets, the excitation of "up" modes, and the resolution of an apparent paradox related to the superposition principle. We propose to generalize the Barack-Ori formula for the decay rate of any tail multipole given a generic initial data set, to the contribution of any initial multipole mode. Our proposal leads to a much simpler expression for the late-time power law index. Specifically, we propose that the late-time decay rate of the $Y_{\ell m}$ spherical harmonic multipole moment because of an initial $Y_{\ell' m}$ multipole is independent of the azimuthal number $m$, and is given by $t^{-n}$, where $n=\ell'+\ell+1$ for $\ell<\ell'$ and $n=\ell'+\ell+3$ for $\ell\ge\ell'$. We also show explicitly that the angular symmetry group of a multipole does not determine its late-time decay rate.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Test bodies and naked singularities: is the self-force the cosmic censor?

Jacobson and Sotiriou showed that rotating black holes could be spun-up past the extremal limit by the capture of non-spinning test bodies, if one neglects radiative and self-force effects. This would represent a violation of the Cosmic Censorship Conjecture in four-dimensional, asymptotically flat spacetimes. We show that for some of the trajectories giving rise to naked singularities, radiative effects can be neglected. However, for these orbits the conservative self-force is important, and seems to have the right sign to prevent the formation of naked singularities.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

An exploration of CUDA and CBEA for a gravitational wave data-analysis application (Einstein@Home)

We present a detailed approach for making use of two new computer hardware architectures -- CBEA and CUDA -- for accelerating a scientific data-analysis application (Einstein@Home). Our results suggest that both the architectures suit the application quite well and the achievable performance in the same software developmental time-frame, is nearly identical.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The CBE Hardware Accelerator for Numerical Relativity: A Simple Approach

Hardware accelerators (such as the Cell Broadband Engine) have recently received a significant amount of attention from the computational science community because they can provide significant gains in the overall performance of many numerical simulations at a low cost. However, such accelerators usually employ a rather unfamiliar and specialized programming model that often requires advanced knowledge of their hardware design. In this article, we demonstrate an alternate and simpler approach towards managing the main complexities in the programming of the Cell processor, called software caching. We apply this technique to a numerical relativity application: a time-domain, finite-difference Kerr black hole perturbation evolver, and present the performance results. We obtain gains in the overall performance of generic simulations that are close to the theoretical maximum that can be obtained through our parallelization approach.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Close limit of grazing black hole collisions: non-spinning holes

Using approximate techniques we study the final moments of the collision of two (individually non-spinnning) black holes which inspiral into each other. The approximation is based on treating the whole space-time as a single distorted black hole. We obtain estimates for the radiated energy, angular momentum and waveforms for the gravitational waves produced in such a collision. The results can be of interest for analyzing the data that will be forthcoming from gravitational wave interferometric detectors, like the LIGO, GEO, LISA, VIRGO and TAMA projects.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

The imposition of Cauchy data to the Teukolsky equation III: The rotating case

We solve the problem of expressing the Weyl scalars $ψ$ that describe gravitational perturbations of a Kerr black hole in terms of Cauchy data. To do so we use geometrical identities (like the Gauss-Codazzi relations) as well as Einstein equations. We are able to explicitly express $ψ$ and $\partial _tψ$ as functions only of the extrinsic curvature and the three-metric (and geometrical objects built out of it) of a generic spacelike slice of the spacetime. These results provide the link between initial data and $ψ$ to be evolved by the Teukolsky equation, and can be used to compute the gravitational radiation generated by two orbiting black holes in the close limit approximation. They can also be used to extract waveforms from spacetimes completely generated by numerical methods.

0 citations0 reviewsNo code linkedOpen access