Source author record

Chad R. Galley

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
19topics
4close collaborators

Actions

Connect this record

Log in to claim

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

25 published item(s)

preprint2020arXiv

Constraining the parameters of GW150914 & GW170104 with numerical relativity surrogates

Gravitational-wave detectors have begun to observe coalescences of heavy black holes at a consistent pace for the past few years. Accurate models of gravitational waveforms are essential for unbiased and precise estimation of source parameters. Recently developed surrogate models based on high-accuracy numerical relativity (NR) simulations are ideal for constraining physical parameters of heavy black hole merger events. In this paper, we first demonstrate the viability of these multi-modal surrogates as reliable parameter estimation tools. We show that NR surrogates can extract additional information from GW data that is inaccessible to traditional models, by analyzing a set of synthetic signals with the NR surrogate and other approximants. We also consider the case of two of the earliest binary black holes detected by the LIGO observatories: GW150914 and GW170104. We reanalyze their data with fully-precessing NR-surrogate templates and freely provide the resulting posterior samples as supplemental material. We find that our refined analysis is able to extract information from sub-dominant GW harmonics in data, and therefore better resolve the degeneracy in measuring source luminosity distance and orbital inclination for both events. We estimate the sources of both events to be 20-25% further away than was previously estimated. Our analyses also constrain their orbital orientation more tightly around face-on or face-off configurations than before. Additionally, for GW150914 we constrain the effective inspiral spin more tightly around zero. This work is one of the first to unambiguously extract sub-dominant GW mode information from real events. It is also a first step toward eliminating the approximations used in semi-analytic waveform models from GW parameter estimation. It also motivates that NR surrogates be extended to cover more of the binary black hole parameter space.

preprint2016arXiv

Fast and efficient evaluation of gravitational waveforms via reduced-order spline interpolation

Numerical simulations of merging black hole binaries produce the most accurate gravitational waveforms. The availability of hundreds of these numerical relativity (NR) waveforms, often containing many higher spherical harmonic modes, allows one to study many aspects of gravitational waves. Amongst these are the response of data analysis pipelines, the calibration of semi-analytical models, the building of reduced-order surrogates, the estimation of the parameters of detected gravitational waves, and the composition of public catalogs of NR waveform data. The large number of generated NR waveforms consequently requires efficient data storage and handling, especially since many more waveforms will be generated at an increased rate in the forthcoming years. In addition, gravitational wave data analyses often require the NR waveforms to be interpolated and uniformly resampled at high sampling rates. Previously, this resulted in very large data files (up to $\sim$ several GB) in memory-intensive operations, which is unfeasible when confronted with hundreds of multi-modal NR waveforms. To handle these challenges, we present a simple and efficient method to significantly \emph{compress} the original waveform data sets while accurately reproducing the original data via spline interpolation. The method is generically applicable to relatively smooth, one-dimensional datasets and uses a greedy algorithm to determine the most relevant subset of the full data such that a spline interpolant of a specified polynomial degree will represent the original data to within a requested point-wise tolerance. We find significant compression of the original NR data sets presented here. These compressed data sets can then be evaluated fast and efficiently and resampled as desired.

preprint2016arXiv

The tail effect in gravitational radiation-reaction: time non-locality and renormalization group evolution

We use the effective field theory (EFT) framework to calculate the tail effect in gravitational radiation reaction, which enters at 4PN order in the dynamics of a binary system. The computation entails a subtle interplay between the near (or potential) and far (or radiation) zones. In particular, we find that the tail contribution to the effective action is non-local in time, and features both a dissipative and a `conservative' term. The latter includes a logarithmic ultraviolet (UV) divergence, which we show cancels against an infrared (IR) singularity found in the (conservative) near zone. The origin of this behavior in the long-distance EFT is due to the point-particle limit -shrinking the binary to a point- which transforms a would-be infrared singularity into an ultraviolet divergence. This is a common occurrence in an EFT approach, which furthermore allows us to use renormalization group (RG) techniques to resum the resulting logarithmic contributions. We then derive the RG evolution for the binding potential and total mass/energy, and find agreement with the results obtained imposing the conservation of the (pseudo) stress-energy tensor in the radiation theory. While the calculation of the leading tail contribution to the effective action involves only one diagram, five are needed for the one-point function. This suggests logarithmic corrections may be easier to incorporate in this fashion. We conclude with a few remarks on the nature of these IR/UV singularities, the (lack of) ambiguities recently discussed in the literature, and the completeness of the analytic Post-Newtonian framework.

preprint2015arXiv

Fast and accurate prediction of numerical relativity waveforms from binary black hole coalescences using surrogate models

Simulating a binary black hole (BBH) coalescence by solving Einstein's equations is computationally expensive, requiring days to months of supercomputing time. Using reduced order modeling techniques, we construct an accurate surrogate model, which is evaluated in a millisecond to a second, for numerical relativity (NR) waveforms from non-spinning BBH coalescences with mass ratios in $[1, 10]$ and durations corresponding to about $15$ orbits before merger. We assess the model's uncertainty and show that our modeling strategy predicts NR waveforms {\em not} used for the surrogate's training with errors nearly as small as the numerical error of the NR code. Our model includes all spherical-harmonic ${}_{-2}Y_{\ell m}$ waveform modes resolved by the NR code up to $\ell=8.$ We compare our surrogate model to Effective One Body waveforms from $50$-$300 M_\odot$ for advanced LIGO detectors and find that the surrogate is always more faithful (by at least an order of magnitude in most cases).

preprint2015arXiv

Slimplectic Integrators: Variational Integrators for General Nonconservative Systems

Symplectic integrators are widely used for long-term integration of conservative astrophysical problems due to their ability to preserve the constants of motion; however, they cannot in general be applied in the presence of nonconservative interactions. In this Letter, we develop the "slimplectic" integrator, a new type of numerical integrator that shares many of the benefits of traditional symplectic integrators yet is applicable to general nonconservative systems. We utilize a fixed time-step variational integrator formalism applied to the principle of stationary nonconservative action developed in Galley, 2013; Galley, Tsang & Stein, 2014. As a result, the generalized momenta and energy (Noether current) evolutions are well-tracked. We discuss several example systems, including damped harmonic oscillators, Poynting-Robertson drag, and gravitational radiation reaction, by utilizing our new publicly available code to demonstrate the slimplectic integrator algorithm. Slimplectic integrators are well-suited for integrations of systems where nonconservative effects play an important role in the long-term dynamical evolution. As such they are particularly appropriate for cosmological or celestial N-body dynamics problems where nonconservative interactions, e.g. gas interactions or dissipative tides, can play an important role.

preprint2014arXiv

A sparse representation of gravitational waves from precessing compact binaries

Many relevant applications in gravitational wave physics share a significant common problem: the seven-dimensional parameter space of gravitational waveforms from precessing compact binary inspirals and coalescences is large enough to prohibit covering the space of waveforms with sufficient density. We find that by using the reduced basis method together with a parametrization of waveforms based on their phase and precession, we can construct ultra-compact yet high-accuracy representations of this large space. As a demonstration, we show that less than $100$ judiciously chosen precessing inspiral waveforms are needed for $200$ cycles, mass ratios from $1$ to $10$ and spin magnitudes $\le 0.9$. In fact, using only the first $10$ reduced basis waveforms yields a maximum mismatch of $0.016$ over the whole range of considered parameters. We test whether the parameters selected from the inspiral regime result in an accurate reduced basis when including merger and ringdown; we find that this is indeed the case in the context of a non-precessing effective-one-body model. This evidence suggests that as few as $\sim 100$ numerical simulations of binary black hole coalescences may accurately represent the seven-dimensional parameter space of precession waveforms for the considered ranges.

preprint2014arXiv

Fast prediction and evaluation of gravitational waveforms using surrogate models

[Abridged] We propose a solution to the problem of quickly and accurately predicting gravitational waveforms within any given physical model. The method is relevant for both real-time applications and in more traditional scenarios where the generation of waveforms using standard methods can be prohibitively expensive. Our approach is based on three offline steps resulting in an accurate reduced-order model that can be used as a surrogate for the true/fiducial waveform family. First, a set of m parameter values is determined using a greedy algorithm from which a reduced basis representation is constructed. Second, these m parameters induce the selection of m time values for interpolating a waveform time series using an empirical interpolant. Third, a fit in the parameter dimension is performed for the waveform's value at each of these m times. The cost of predicting L waveform time samples for a generic parameter choice is of order m L + m c_f online operations where c_f denotes the fitting function operation count and, typically, m << L. We generate accurate surrogate models for Effective One Body (EOB) waveforms of non-spinning binary black hole coalescences with durations as long as 10^5 M, mass ratios from 1 to 10, and for multiple harmonic modes. We find that these surrogates are three orders of magnitude faster to evaluate as compared to the cost of generating EOB waveforms in standard ways. Surrogate model building for other waveform models follow the same steps and have the same low online scaling cost. For expensive numerical simulations of binary black hole coalescences we thus anticipate large speedups in generating new waveforms with a surrogate. As waveform generation is one of the dominant costs in parameter estimation algorithms and parameter space exploration, surrogate models offer a new and practical way to dramatically accelerate such studies without impacting accuracy.

preprint2014arXiv

Self-force via Green functions and worldline integration

A compact object moving in curved spacetime interacts with its own gravitational field. This leads to both dissipative and conservative corrections to the motion, which can be interpreted as a self-force acting on the object. The original formalism describing this self-force relied heavily on the Green function of the linear differential operator that governs gravitational perturbations. However, because the global calculation of Green functions in non-trivial black hole spacetimes has been an open problem until recently, alternative methods were established to calculate self-force effects using sophisticated regularization techniques that avoid the computation of the global Green function. We present a method for calculating the self-force that employs the global Green function and is therefore closely modeled after the original self-force expressions. Our quantitative method involves two stages: (i) numerical approximation of the retarded Green function in the background spacetime; (ii) evaluation of convolution integrals along the worldline of the object. This novel approach can be used along arbitrary worldlines, including those currently inaccessible to more established computational techniques. Furthermore, it yields geometrical insight into the contributions to self-interaction from curved geometry (back-scattering) and trapping of null geodesics. We demonstrate the method on the motion of a scalar charge in Schwarzschild spacetime. This toy model retains the physical history-dependence of the self-force but avoids gauge issues and allows us to focus on basic principles. We compute the self-field and self-force for many worldlines including accelerated circular orbits, eccentric orbits at the separatrix, and radial infall. This method, closely modeled after the original formalism, provides a promising complementary approach to the self-force problem.

preprint2014arXiv

The principle of stationary nonconservative action for classical mechanics and field theories

We further develop a recently introduced variational principle of stationary action for problems in nonconservative classical mechanics and extend it to classical field theories. The variational calculus used is consistent with an initial value formulation of physical problems and allows for time-irreversible processes, such as dissipation, to be included at the level of the action. In this formalism, the equations of motion are generated by extremizing a nonconservative action $\mathcal{S}$, which is a functional of a doubled set of degrees of freedom. The corresponding nonconservative Lagrangian contains a potential $K$ which generates nonconservative forces and interactions. Such a nonconservative potential can arise in several ways, including from an open system interacting with inaccessible degrees of freedom or from integrating out or coarse-graining a subset of variables in closed systems. We generalize Noether's theorem to show how Noether currents are modified and no longer conserved when $K$ is non-vanishing. Consequently, the nonconservative aspects of a physical system are derived solely from $K$. We show how to use the formalism with examples of nonconservative actions for discrete systems including forced damped harmonic oscillators, radiation reaction on an accelerated charge, and RLC circuits. We present examples for nonconservative classical field theories. Our approach naturally allows for irreversible thermodynamic processes to be included in an unconstrained variational principle. We present the nonconservative action for a Navier-Stokes fluid including the effects of viscous dissipation and heat diffusion, as well as an action that generates the Maxwell model for viscoelastic materials, which can be easily generalized to more realistic rheological models. We show that the nonconservative action can be derived as the classical limit of a more complete quantum theory.

preprint2013arXiv

Gravitational self-force in the ultra-relativistic limit: The 'large-N' expansion

We study the gravitational self-force using the effective field theory formalism. We show that in the ultra-relativistic limit γ\to \infty, with γthe boost factor, many simplifications arise. Drawing parallels with the large N limit in quantum field theory, we introduce the parameter 1/N = 1/γ^2 and show that the effective action admits a well defined expansion in powers of λ= Nε, at each order in 1/N, where ε= E_m/M and E_m=γm is the (kinetic) energy of the small mass. Moreover, we show that diagrams with nonlinear bulk interactions first enter at O(λ^2/N^2) and only diagrams with nonlinearities in the worldline couplings, which are significantly easier to compute, survive in the large N/ultra-relativistic limit. Finally, we derive the self-force to O(λ^4/N) and provide expressions for some conservative quantities for circular orbits.

preprint2013arXiv

The classical mechanics of non-conservative systems

Hamilton's principle of stationary action lies at the foundation of theoretical physics and is applied in many other disciplines from pure mathematics to economics. Despite its utility, Hamilton's principle has a subtle pitfall that often goes unnoticed in physics: it is formulated as a boundary value problem in time but is used to derive equations of motion that are solved with initial data. This subtlety can have undesirable effects. I present a formulation of Hamilton's principle that is compatible with initial value problems. Remarkably, this leads to a natural formulation for the Lagrangian and Hamiltonian dynamics of generic non-conservative systems, thereby filling a long-standing gap in classical mechanics. Thus dissipative effects, for example, can be studied with new tools that may have application in a variety of disciplines. The new formalism is demonstrated by two examples of non-conservative systems: an object moving in a fluid with viscous drag forces and a harmonic oscillator coupled to a dissipative environment.

preprint2012arXiv

Caustic echoes from a Schwarzschild black hole

We present the first numerical construction of the scalar Schwarzschild Green function in the time-domain, which reveals several universal features of wave propagation in black hole spacetimes. We demonstrate the trapping of energy near the photon sphere and confirm its exponential decay. The trapped wavefront propagates through caustics resulting in echoes that propagate to infinity. The arrival times and the decay rate of these caustic echoes are consistent with propagation along null geodesics and the large l-limit of quasinormal modes. We show that the four-fold singularity structure of the retarded Green function is due to the well-known action of a Hilbert transform on the trapped wavefront at caustics. A two-fold cycle is obtained for degenerate source-observer configurations along the caustic line, where the energy amplification increases with an inverse power of the scale of the source. Finally, we discuss the tail piece of the solution due to propagation within the light cone, up to and including null infinity, and argue that, even with ideal instruments, only a finite number of echoes can be observed. Putting these pieces together, we provide a heuristic expression that approximates the Green function with a few free parameters. Accurate calculations and approximations of the Green function are the most general way of solving for wave propagation in curved spacetimes and should be useful in a variety of studies such as the computation of the self-force on a particle.

preprint2012arXiv

Non-sky-averaged sensitivity curves for space-based gravitational-wave observatories

(abridged) The signal-to-noise ratio (SNR) is used in gravitational-wave observations as the basic figure of merit for detection confidence and, together with the Fisher matrix, for the amount of physical information that can be extracted from a detected signal. SNRs are usually computed from a sensitivity curve, which describes the gravitational-wave amplitude needed by a monochromatic source of given frequency to achieve a threshold SNR. For interferometric space-based detectors similar to LISA, which are sensitive to long-lived signals and have constantly changing position and orientation, exact SNRs need to be computed on a source-by-source basis. For convenience, most authors prefer to work with sky-averaged sensitivities, accepting inaccurate SNRs for individual sources and giving up control over the statistical distribution of SNRs for source populations. In this paper, we describe a straightforward end-to-end recipe to compute the non-sky-averaged sensitivity of interferometric space-based detectors of any geometry: in effect, we derive error bars for the sky-averaged sensitivity curve, which provide a stringent statistical interpretation for previously unqualified statements about sky-averaged SNRs. As a worked-out example, we consider isotropic and Galactic-disk populations of monochromatic sources, as observed with the "classic LISA" configuration. We confirm that the (standard) inverse-rms average sensitivity for the isotropic population remains the same whether or not the LISA orbits are included in the computation. However, detector motion tightens the distribution of sensitivities, so for 50% of sources the sensitivity is within 30% of its average. For the Galactic-disk population, the average and the distribution of the sensitivity for a moving detector turn out to be similar to the isotropic case.

preprint2012arXiv

Radiation reaction at 3.5 post-Newtonian order in effective field theory

We derive the radiation reaction forces on a compact binary inspiral through 3.5 order in the post-Newtonian expansion using the effective field theory approach. We utilize a recent formulation of Hamilton's variational principle that rigorously extends the usual Lagrangian and Hamiltonian formalisms to dissipative systems, including the inspiral of a compact binary from the emission of gravitational waves. We find agreement with previous results, which thus provides a non-trivial confirmation of the extended variational principle. The results from this work nearly complete the equations of motion for the generic inspiral of a compact binary with spinning constituents through 3.5 post-Newtonian order, as derived entirely with effective field theory, with only the spin-orbit corrections to the potential at 3.5 post-Newtonian remaining.

preprint2012arXiv

Reduced Basis representations of multi-mode black hole ringdown gravitational waves

We construct compact and high accuracy Reduced Basis (RB) representations of single and multiple quasinormal modes (QNMs). The RB method determines a hierarchical and relatively small set of the most relevant waveforms. We find that the exponential convergence of the method allows for a dramatic compression of template banks used for ringdown searches. Compressing a catalog with a minimal match $\MMm=0.99$, we find that the selected RB waveforms are able to represent {\em any} QNM, including those not in the original bank, with extremely high accuracy, typically less than $10^{-13}$. We then extend our studies to two-mode QNMs. Inclusion of a second mode is expected to help with detection, and might make it possible to infer details of the progenitor of the final black hole. We find that the number of RB waveforms needed to represent any two-mode ringdown waveform with the above high accuracy is {\em smaller} than the number of metric-based, one-mode templates with $\MMm=0.99$. For unconstrained two-modes, which would allow for consistency tests of General Relativity, our high accuracy RB has around $10^4$ {\em fewer} waveforms than the number of metric-based templates for $\MMm=0.99$. The number of RB elements grows only linearly with the number of multipole modes versus exponentially with the standard approach, resulting in very compact representations even for many multiple modes. The results of this paper open the possibility of searches of multi-mode ringdown gravitational waves.

preprint2012arXiv

Theory of optomechanics: Oscillator-field model of moving mirrors

In this paper we present a model for the kinematics and dynamics of optomechanics which describe the coupling between an optical field, here modeled by a massless scalar field, and the internal (e.g., determining its reflectivity) and mechanical (e.g., displacement) degrees of freedom of a moveable mirror. As opposed to implementing boundary conditions on the field we highlight the internal dynamics of the mirror which provides added flexibility to describe a variety of setups relevant to current experiments. The inclusion of the internal degrees of freedom in this model allows for a variety of optical activities of mirrors from those exhibiting broadband reflective properties to the cases where reflection is suppressed except for a narrow band centered around the characteristic frequency associated with the mirror's internal dynamics. After establishing the model and the reflective properties of the mirror we show how appropriate parameter choices lead to useful optomechanical models such as the well known Barton-Calogeracos model [G. Barton and A. Calogeracos, Ann. Phys. 238, 227 (1995)] and the important yet lesser explored nonlinear models (e.g., $Nx$ coupling) for small photon numbers $N$, which present models based on side-band approximations [H. Kimble et al., Phys. Rev. D 65, 022002 (2001)] cannot cope with. As a simple illustrative application we consider classical radiation pressure cooling with this model. To expound its theoretical structure and physical meanings we connect our model to field-theoretical models using auxiliary fields and the ubiquitous Brownian motion model of quantum open systems. Finally we describe the range of applications of this model, from a full quantum mechanical treatment of radiation pressure cooling, quantum entanglement between macroscopic mirrors, to the backreaction of Hawking radiation on black hole evaporation in a moving mirror analog.

preprint2012arXiv

Towards beating the curse of dimensionality for gravitational waves using Reduced Basis

Using the Reduced Basis approach, we efficiently compress and accurately represent the space of waveforms for non-precessing binary black hole inspirals, which constitutes a four dimensional parameter space (two masses, two spin magnitudes). Compared to the non-spinning case, we find that only a {\it marginal} increase in the (already relatively small) number of reduced basis elements is required to represent any non-precessing waveform to nearly numerical round-off precision. Most parameters selected by the algorithm are near the boundary of the parameter space, leaving the bulk of its volume sparse. Our results suggest that the full eight dimensional space (two masses, two spin magnitudes, four spin orientation angles on the unit sphere) may be highly compressible and represented with very high accuracy by a remarkably small number of waveforms, thus providing some hope that the number of numerical relativity simulations of binary black hole coalescences needed to represent the entire space of configurations is not intractable. Finally, we find that the {\it distribution} of selected parameters is robust to different choices of seed values starting the algorithm, a property which should be useful for indicating parameters for numerical relativity simulations of binary black holes. In particular, we find that the mass ratios $m_1/m_2$ of non-spinning binaries selected by the algorithm are mostly in the interval $[1,3]$ and that the median of the distribution follows a power-law behavior $\sim (m_1/m_2)^{-5.25}$.

preprint2011arXiv

A nonlinear scalar model of extreme mass ratio inspirals in effective field theory II. Scalar perturbations and a master source

The motion of a small compact object (SCO) in a background spacetime is investigated further in a class of model nonlinear scalar field theories having a perturbative structure analogous to the General Relativistic description of extreme mass ratio inspirals (EMRIs). We derive regular expressions for the scalar perturbations generated by the SCO's motion valid through third order in $ε$, the size of the SCO to the background curvature length scale. Our expressions are compared to those calculated through second order in $ε$ by Rosenthal in [E. Rosenthal, CQG 22, S859 (2005)] and found to agree but our procedure for regularizing the scalar perturbations is considerably simpler. Following the Detweiler-Whiting (DW) scheme, we use our regular expressions for the field and derive the regular self-force corrections through third order. We find agreement with our previous derivation based on a variational principle of an effective action for the worldline associated with the SCO thus demonstrating the internal consistency of our formalism. This also explicitly demonstrates that the DW decomposition of Green's functions is a valid and practical method of self force computation at higher orders in perturbation theory and, as we show in an appendix, at all orders in perturbation theory. Finally, we identify a master source from which all other physically relevant quantities are derivable. Knowing the master source perturbatively allows one to construct the waveform measured by an observer, the regular part of the field on the worldline, the regular part of the self force, and orbital quantities such as shifts of the innermost stable circular orbit, etc. The existence of a master source together with the regularization methods implemented in this series should be indispensable for derivations of higher-order gravitational self force corrections.

preprint2011arXiv

Reduced basis catalogs for gravitational wave templates

We introduce a reduced basis approach as a new paradigm for modeling, representing and searching for gravitational waves. We construct waveform catalogs for non-spinning compact binary coalescences, and we find that for accuracies of 99% and 99.999% the method generates a factor of about $10-10^5$ fewer templates than standard placement methods. The continuum of gravitational waves can be represented by a finite and comparatively compact basis. The method is robust under variations in the noise of detectors, implying that only a single catalog needs to be generated.

preprint2010arXiv

A nonlinear scalar model of extreme mass ratio inspirals in effective field theory I. Self force through third order

The motion of a small compact object in a background spacetime is investigated in the context of a model nonlinear scalar field theory. This model is constructed to have a perturbative structure analogous to the General Relativistic description of extreme mass ratio inspirals (EMRIs). We apply the effective field theory approach to this model and calculate the finite part of the self force on the small compact object through third order in the ratio of the size of the compact object to the curvature scale of the background (e.g., black hole) spacetime. We use well-known renormalization methods and demonstrate the consistency of the formalism in rendering the self force finite at higher orders within a point particle prescription for the small compact object. This nonlinear scalar model should be useful for studying various aspects of higher-order self force effects in EMRIs but within a comparatively simpler context than the full gravitational case. These aspects include developing practical schemes for higher order self force numerical computations, quantifying the effects of transient resonances on EMRI waveforms and accurately modeling the small compact object's motion for precise determinations of the parameters of detected EMRI sources.

preprint2010arXiv

Finite size corrections to the radiation reaction force in classical electrodynamics

We introduce an effective field theory approach that describes the motion of finite size objects under the influence of electromagnetic fields. We prove that leading order effects due to the finite radius $R$ of a spherically symmetric charge is order $R^2$ rather than order $R$ in any physical model, as widely claimed in the literature. This scaling arises as a consequence of Poincaré and gauge symmetries, which can be shown to exclude linear corrections. We use the formalism to calculate the leading order finite size correction to the Abraham-Lorentz-Dirac force.

preprint2010arXiv

Statistical constraints on binary black hole inspiral dynamics

We perform a statistical analysis of the binary black hole problem in the post-Newtonian approximation by systematically sampling and evolving the parameter space of initial configurations for quasi-circular inspirals. Through a principal component analysis of spin and orbital angular momentum variables we systematically look for uncorrelated quantities and find three of them which are highly conserved in a statistical sense, both as functions of time and with respect to variations in initial spin orientations. We also look for and find the variables that account for the largest variations in the problem. We present binary black hole simulations of the full Einstein equations analyzing to what extent these results might carry over to the full theory in the inspiral and merger regimes. Among other applications these results should be useful both in semi-analytical and numerical building of templates of gravitational waves for gravitational wave detectors.

preprint2009arXiv

Radiation reaction and gravitational waves in the effective field theory approach

We compute the contribution to the Lagrangian from the leading order (2.5 post-Newtonian) radiation reaction and the quadrupolar gravitational waves emitted from a binary system using the effective field theory (EFT) approach of Goldberger and Rothstein. We use an initial value formulation of the underlying (quantum) framework to implement retarded boundary conditions and describe these real-time dissipative processes. We also demonstrate why the usual scattering formalism of quantum field theory inadequately accounts for these. The methods discussed here should be useful for deriving real-time quantities (including radiation reaction forces and gravitational wave emission) and hereditary terms in the post-Newtonian approximation (including memory, tail and other causal, history-dependent integrals) within the EFT approach. We also provide a consistent formulation of the radiation sector in the equivalent effective field theory approach of Kol and Smolkin.

preprint2005arXiv

Self-Force with a Stochastic Component from Radiation Reaction of a Scalar Charge Moving in Curved Spacetime

We give a quantum field theoretical derivation of the scalar Abraham-Lorentz-Dirac (ALD) equation and the self-force for a scalar charged particle interacting with a quantum scalar field in curved spacetime. We regularize the causal Green's function using a quasi-local expansion in the spirit of effective field theory and obtain a regular expression for the self-force. The scalar ALD equation obtained in this way for the classical motion of the particle checks with the equation obtained by Quinn earlier \cite{Quinn}. We further derive a scalar ALD-Langevin equation with a classical stochastic force accounting for the effect of quantum fluctuations in the field, which causes small fluctuations on the particle trajectory. This equation will be useful for the study of stochastic motion of charges under the influence of both quantum and classical noise sources, derived either self-consistently (as done here) or put in by hand (with warnings). We show the possibility of secular effects from such stochastic influences on the trajectory that may impact on the present calculations of gravitational waveform templates.