BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Stephen R. Lau

Stephen R. Lau contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
gr-qc

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
1topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 published item(s)

preprint2012arXiv

Fast evaluation of asymptotic waveforms from gravitational perturbations

In the context of blackhole perturbation theory, we describe both exact evaluation of an asymptotic waveform from a time series recorded at a finite radial location and its numerical approximation. From the user's standpoint our technique is easy to implement, affords high accuracy, and works for both axial (Regge-Wheeler) and polar (Zerilli) sectors. Our focus is on the ease of implementation with publicly available numerical tables, either as part of an existing evolution code or a post-processing step. Nevertheless, we also present a thorough theoretical discussion of asymptotic waveform evaluation and radiation boundary conditions, which need not be understood by a user of our methods. In particular, we identify (both in the time and frequency domains) analytical asymptotic waveform evaluation kernels, and describe their approximation by techniques developed by Alpert, Greengard, and Hagstrom. This paper also presents new results on the evaluation of far-field signals for the ordinary (acoustic) wave equation. We apply our method to study late-time decay tails at null-infinity, "teleportation" of a signal between two finite radial values, and luminosities from extreme-mass-ratio binaries. Through numerical simulations with the outer boundary as close in as r = 30M, we compute asymptotic waveforms with late-time t^{-4} decay (l = 2 perturbations), and also luminosities from circular and eccentric particle-orbits that respectively match frequency domain results to relative errors of better than 10^{-12} and 10^{-9}. Furthermore, we find that asymptotic waveforms are especially prone to contamination by spurious junk radiation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Implicit-explicit (IMEX) evolution of single black holes

Numerical simulations of binary black holes---an important predictive tool for the detection of gravitational waves---are computationally expensive, especially for binaries with high mass ratios or with rapidly spinning constituent holes. Existing codes for evolving binary black holes rely on explicit timestepping methods, for which the timestep size is limited by the smallest spatial scale through the Courant-Friedrichs-Lewy condition. Binary inspiral typically involves spatial scales (the spatial resolution required by a small or rapidly spinning hole) which are orders of magnitude smaller than the relevant (orbital, precession, and radiation-reaction) timescales characterizing the inspiral. Therefore, in explicit evolutions of binary black holes, the timestep size is typically orders of magnitude smaller than the relevant physical timescales. Implicit timestepping methods allow for larger timesteps, and they often reduce the total computational cost (without significant loss of accuracy) for problems dominated by spatial rather than temporal error, such as for binary-black-hole inspiral in corotating coordinates. However, fully implicit methods can be difficult to implement for nonlinear evolution systems like the Einstein equations. Therefore, in this paper we explore implicit-explicit (IMEX) methods and use them for the first time to evolve black-hole spacetimes. Specifically, as a first step toward IMEX evolution of a full binary-black-hole spacetime, we develop an IMEX algorithm for the generalized harmonic formulation of the Einstein equations and use this algorithm to evolve stationary and perturbed single-black-hole spacetimes. Numerical experiments explore the stability and computational efficiency of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Sparse spectral-tau method for the three-dimensional helically reduced wave equation on two-center domains

We describe a multidomain spectral-tau method for solving the three-dimensional helically reduced wave equation on the type of two-center domain that arises when modeling compact binary objects in astrophysical applications. A global two-center domain may arise as the union of Cartesian blocks, cylindrical shells, and inner and outer spherical shells. For each such subdomain, our key objective is to realize certain (differential and multiplication) physical-space operators as matrices acting on the corresponding set of modal coefficients. We achieve sparse banded realizations through the integration "preconditioning" of Coutsias, Hagstrom, Hesthaven, and Torres. Since ours is the first three-dimensional multidomain implementation of the technique, we focus on the issue of convergence for the global solver, here the alternating Schwarz method accelerated by GMRES. Our methods may prove relevant for numerical solution of other mixed-type or elliptic problems, and in particular for the generation of initial data in general relativity.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Discontinuous Galerkin method for the spherically reduced BSSN system with second-order operators

We present a high-order accurate discontinuous Galerkin method for evolving the spherically-reduced Baumgarte-Shapiro-Shibata-Nakamura (BSSN) system expressed in terms of second-order spatial operators. Our multi-domain method achieves global spectral accuracy and long-time stability on short computational domains. We discuss in detail both our scheme for the BSSN system and its implementation. After a theoretical and computational verification of the proposed scheme, we conclude with a brief discussion of issues likely to arise when one considers the full BSSN system.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Persistent junk solutions in time-domain modeling of extreme mass ratio binaries

In the context of metric perturbation theory for non-spinning black holes, extreme mass ratio binary (EMRB) systems are described by distributionally forced master wave equations. Numerical solution of a master wave equation as an initial boundary value problem requires initial data. However, because the correct initial data for generic-orbit systems is unknown, specification of trivial initial data is a common choice, despite being inconsistent and resulting in a solution which is initially discontinuous in time. As is well known, this choice leads to a "burst" of junk radiation which eventually propagates off the computational domain. We observe another unintended consequence of trivial initial data: development of a persistent spurious solution, here referred to as the Jost junk solution, which contaminates the physical solution for long times. This work studies the influence of both types of junk on metric perturbations, waveforms, and self-force measurements, and it demonstrates that smooth modified source terms mollify the Jost solution and reduce junk radiation. Our concluding section discusses the applicability of these observations to other numerical schemes and techniques used to solve distributionally forced master wave equations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

IMEX evolution of scalar fields on curved backgrounds

Inspiral of binary black holes occurs over a time-scale of many orbits, far longer than the dynamical time-scale of the individual black holes. Explicit evolutions of a binary system therefore require excessively many time steps to capture interesting dynamics. We present a strategy to overcome the Courant-Friedrichs-Lewy condition in such evolutions, one relying on modern implicit-explicit ODE solvers and multidomain spectral methods for elliptic equations. Our analysis considers the model problem of a forced scalar field propagating on a generic curved background. Nevertheless, we encounter and address a number of issues pertinent to the binary black hole problem in full general relativity. Specializing to the Schwarzschild geometry in Kerr-Schild coordinates, we document the results of several numerical experiments testing our strategy.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Differential Forms and Wave Equations for General Relativity

Recently, Choquet-Bruhat and York and Abrahams, Anderson, Choquet-Bruhat, and York (AACY) have cast the 3+1 evolution equations of general relativity in gauge-covariant and causal ``first-order symmetric hyperbolic form,'' thereby cleanly separating physical from gauge degrees of freedom in the Cauchy problem for general relativity. A key ingredient in their construction is a certain wave equation which governs the light-speed propagation of the extrinsic curvature tensor. Along a similar line, we construct a related wave equation which, as the key equation in a system, describes vacuum general relativity. Whereas the approach of AACY is based on tensor-index methods, the present formulation is written solely in the language of differential forms. Our approach starts with Sparling's tetrad-dependent differential forms, and our wave equation governs the propagation of Sparling's 2-form, which in the ``time-gauge'' is built linearly from the ``extrinsic curvature 1-form.'' The tensor-index version of our wave equation describes the propagation of (what is essentially) the Arnowitt-Deser-Misner gravitational momentum.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

On the Canonical Reduction of Spherically Symmetric Gravity

In a thorough paper Kuchar has examined the canonical reduction of the most general action functional describing the geometrodynamics of the maximally extended Schwarzschild geometry. This reduction yields the true degrees of freedom for (vacuum) spherically symmetric general relativity. The essential technical ingredient in Kuchar's analysis is a canonical transformation to a certain chart on the gravitational phase space which features the Schwarzschild mass parameter $M_{S}$, expressed in terms of what are essentially Arnowitt-Deser-Misner variables, as a canonical coordinate. In this paper we discuss the geometric interpretation of Kuchar's canonical transformation in terms of the theory of quasilocal energy-momentum in general relativity given by Brown and York. We find Kuchar's transformation to be a ``sphere-dependent boost to the rest frame," where the ``rest frame'' is defined by vanishing quasilocal momentum. Furthermore, our formalism is general enough to cover the case of (vacuum) two-dimensional dilaton gravity. Therefore, besides reviewing Kuchař's original work for Schwarzschild black holes from the framework of hyperbolic geometry, we present new results concerning the canonical reduction of Witten-black-hole geometrodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Spinors and the reference point of quasilocal energy

This paper investigates the relationship between the quasilocal energy of Brown and York and certain spinorial expressions for gravitational energy constructed from the Witten-Nester integral. A key feature of the Brown-York method for defining quasilocal energy is that it allows for the freedom to assign the reference point of the energy. When possible, it is perhaps most natural to reference the energy against flat space, i.e. assign flat-space the zero value of energy. It is demonstrated that the Witten-Nester integral when evaluated on solution spinors to the Sen-Witten equation (obeying appropriate boundary conditions) is essentially the Brown-York quasilocal energy with a reference point determined by the Sen-Witten spinors. For the case of round spheres in the Schwarzschild geometry, these spinors determine the flat-space reference point. A similar viewpoint is proposed for the Schwarzschild-case quasilocal energy of Dougan and Mason.

0 citations0 reviewsNo code linkedOpen access