Source author record

Richard A. Matzner

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

ResearcherUnclaimed source record

gr-qc astro-ph hep-ph astro-ph.CO astro-ph.SR hep-th math-ph math.MP

Catalog footprint

What is connected

12works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2020arXiv

Gauge Fixing and Constrained Dynamics

We review the Dirac formalism for dealing with constraints in a canonical Hamiltonian formulation and discuss gauge freedom and display constraints for gauge theories in a general context. We introduce the Dirac bracket and show that it provides a consistent method to remove any gauge freedom present. We discuss stability in evolution of gauge theories and show that fixing all gauge freedom is sufficient to ensure well-posedness for a large class of gauge theories. Electrodynamics provides examples of the methods outlined for general gauge theories. Future work will apply the formalism, and results derived here, to General Relativity.

preprint2013arXiv

Error-analysis and comparison to analytical models of numerical waveforms produced by the NRAR Collaboration

The Numerical-Relativity-Analytical-Relativity (NRAR) collaboration is a joint effort between members of the numerical relativity, analytical relativity and gravitational-wave data analysis communities. The goal of the NRAR collaboration is to produce numerical-relativity simulations of compact binaries and use them to develop accurate analytical templates for the LIGO/Virgo Collaboration to use in detecting gravitational-wave signals and extracting astrophysical information from them. We describe the results of the first stage of the NRAR project, which focused on producing an initial set of numerical waveforms from binary black holes with moderate mass ratios and spins, as well as one non-spinning binary configuration which has a mass ratio of 10. All of the numerical waveforms are analysed in a uniform and consistent manner, with numerical errors evaluated using an analysis code created by members of the NRAR collaboration. We compare previously-calibrated, non-precessing analytical waveforms, notably the effective-one-body (EOB) and phenomenological template families, to the newly-produced numerical waveforms. We find that when the binary's total mass is ~100-200 solar masses, current EOB and phenomenological models of spinning, non-precessing binary waveforms have overlaps above 99% (for advanced LIGO) with all of the non-precessing-binary numerical waveforms with mass ratios <= 4, when maximizing over binary parameters. This implies that the loss of event rate due to modelling error is below 3%. Moreover, the non-spinning EOB waveforms previously calibrated to five non-spinning waveforms with mass ratio smaller than 6 have overlaps above 99.7% with the numerical waveform with a mass ratio of 10, without even maximizing on the binary parameters.

preprint2011arXiv

Measuring emission coordinates in a pulsar-based relativistic positioning system

A relativistic deep space positioning system has been proposed using four or more pulsars with stable repetition rates. (Each pulsar emits pulses at a fixed repetition period in its rest frame.) The positioning system uses the fact that an event in spacetime can be fully described by emission coordinates: the proper emission time of each pulse measured at the event. The proper emission time of each pulse from four different pulsars---interpolated as necessary---provides the four spacetime coordinates of the reception event in the emission coordinate system. If more than four pulsars are available, the redundancy can improve the accuracy of the determination and/or resolve degeneracies resulting from special geometrical arrangements of the sources and the event. We introduce a robust numerical approach to measure the emission coordinates of an event in any arbitrary spacetime geometry. Our approach uses a continuous solution of the eikonal equation describing the backward null cone from the event. The pulsar proper time at the instant the null cone intersects the pulsar world line is one of the four required coordinates. The process is complete (modulo degeneracies) when four pulsar world lines have been crossed by the light cone. The numerical method is applied in two different examples: measuring emission coordinates of an event in Minkowski spacetime using pulses from four pulsars stationary in the spacetime; and measuring emission coordinates of an event in Schwarzschild spacetime using pulses from four pulsars freely falling toward a static black hole. These numerical simulations are merely exploratory, but with improved resolution and computational resources the method can be applied to more pertinent problems. For instance one could measure the emission coordinates, and therefore the trajectory, of the Earth.

preprint2011arXiv

Super-Extremal Spinning Black Holes via Accretion

A Kerr black hole with mass $M$ and angular momentum $J$ satisfies the extremality inequality $|J| \le M^2$. In the presence of matter and/or gravitational radiation, this bound needs to be reformulated in terms of local measurements of the mass and the angular momentum directly associated with the black hole. The isolated and dynamical horizon framework provides such quasi-local characterization of black hole mass and angular momentum. With this framework, it is possible in axisymmetry to reformulate the extremality limit as $|J| \le 2\,M_H^2$, with $M_H$ the irreducible mass of the black hole computed from its apparent horizon area and $J$ obtained using approximate rotational Killing vectors on the apparent horizon. The $|J| \le 2\,M_H^2$ condition is also equivalent to requiring a non-negative black hole surface gravity. We present numerical experiments of an accreting black hole that temporarily violates this extremality inequality. The initial configuration consists of a single, rotating black hole surrounded by a thick, shell cloud of negative energy density. For these numerical experiments, we introduce a new matter-without-matter evolution method.

preprint2010arXiv

Final Mass and Spin of Merged Black Holes and the Golden Black Hole

We present results on the mass and spin of the final black hole from mergers of equal mass, spinning black holes. The study extends over a broad range of initial orbital configurations, from direct plunges to quasi-circular inspirals to more energetic orbits (generalizations of Newtonian elliptical orbits). It provides a comprehensive search of those configurations that maximize the final spin of the remnant black hole. We estimate that the final spin can reach a maximum spin $a/M_h \approx 0.99\pm 0.01$ for extremal black hole mergers. In addition, we find that, as one increases the orbital angular momentum from small values, the mergers produce black holes with mass and spin parameters $\lbrace M_h/M, a/M_h \rbrace$ ~spiraling around the values $\lbrace \hat M_h/M, \hat a/M_h \rbrace$ of a {\it golden} black hole. Specifically, $(M_h-\hat M_h)/M \propto e^{\pm B\,ϕ}\cosϕ$ and $(a-\hat a)/M_h \propto e^{\pm C\,ϕ}\sinϕ$, with $ϕ$ a monotonically growing function of the initial orbital angular momentum. We find that the values of the parameters for the \emph{golden} black hole are those of the final black hole obtained from the merger of a binary with the corresponding spinning black holes in a quasi-circular inspiral.

preprint2009arXiv

Area Invariance of Apparent Horizons under Arbitrary Boosts

It is a well known analytic result in general relativity that the 2-dimensional area of the apparent horizon of a black hole remains invariant regardless of the motion of the observer, and in fact is independent of the $ t=constant $ slice, which can be quite arbitrary in general relativity. Nonetheless the explicit computation of horizon area is often substantially more difficult in some frames (complicated by the coordinate form of the metric), than in other frames. Here we give an explicit demonstration for very restricted metric forms of (Schwarzschild and Kerr) vacuum black holes. In the Kerr-Schild coordinate expression for these spacetimes they have an explicit Lorentz-invariant form. We consider {\it boosted} versions with the black hole moving through the coordinate system. Since these are stationary black hole spacetimes, the apparent horizons are two dimensional cross sections of their event horizons, so we compute the areas of apparent horizons in the boosted space with (boosted) $ t = constant $, and obtain the same result as in the unboosted case. Note that while the invariance of area is generic, we deal only with black holes in the Kerr-Schild form, and consider only one particularly simple change of slicing which amounts to a boost. Even with these restrictions we find that the results illuminate the physics of the horizon as a null surface and provide a useful pedagogical tool. As far as we can determine, this is the first explicit calculation of this type demonstrating the area invariance of horizons. Further, these calculations are directly relevant to transformations that arise in computational representation of moving black holes. We present an application of this result to initial data for boosted black holes.

preprint2008arXiv

The Volume Inside a Black Hole

The horizon (the surface) of a black hole is a null surface, defined by those hypothetical "outgoing" light rays that just hover under the influence of the strong gravity at the surface. Because the light rays are orthogonal to the spatial 2-dimensional surface at one instant of time, the surface of the black hole is the same for all observers (i.e. the same for all coordinate definitions of "instant of time"). This value is 4*(pi)* (2Gm/c^2)^2 for nonspinning black holes, with G= Newton's constant, c= speed of light, and m= mass of the black hole. The 3-dimensional spatial volume inside a black hole, in contrast, depends explicitly on the definition of time, and can even be time dependent, or zero. We give examples of the volume found inside a standard, nonspinning spherical black hole, for several different standard time-coordinate definitions. Elucidating these results for the volume provides a new pedagogical resource of facts already known in principle to the relativity community, but rarely worked out.

preprint2007arXiv

Gravitational recoil from spinning binary black hole mergers

The inspiral and merger of binary black holes will likely involve black holes with both unequal masses and arbitrary spins. The gravitational radiation emitted by these binaries will carry angular as well as linear momentum. A net flux of emitted linear momentum implies that the black hole produced by the merger will experience a recoil or kick. Previous studies have focused on the recoil velocity from unequal mass, non-spinning binaries. We present results from simulations of equal mass but spinning black hole binaries and show how a significant gravitational recoil can also be obtained in these situations. We consider the case of black holes with opposite spins of magnitude $a$ aligned/anti-aligned with the orbital angular momentum, with $a$ the dimensionless spin parameters of the individual holes. For the initial setups under consideration, we find a recoil velocity of $V = 475 \KMS a$. Supermassive black hole mergers producing kicks of this magnitude could result in the ejection from the cores of dwarf galaxies of the final hole produced by the collision.

preprint2004arXiv

Extended Lifetime in Computational Evolution of Isolated Black Holes

Solving the 4-d Einstein equations as evolution in time requires solving equations of two types: the four elliptic initial data (constraint) equations, followed by the six second order evolution equations. Analytically the constraint equations remain solved under the action of the evolution, and one approach is to simply monitor them ({\it unconstrained} evolution). The problem of the 3-d computational simulation of even a single isolated vacuum black hole has proven to be remarkably difficult. Recently, we have become aware of two publications that describe very long term evolution, at least for single isolated black holes. An essential feature in each of these results is {\it constraint subtraction}. Additionally, each of these approaches is based on what we call "modern," hyperbolic formulations of the Einstein equations. It is generally assumed, based on computational experience, that the use of such modern formulations is essential for long-term black hole stability. We report here on comparable lifetime results based on the much simpler ("traditional") $\dot g$ - $\dot K$ formulation. We have also carried out a series of {\it constrained} 3-d evolutions of single isolated black holes. We find that constraint solution can produce substantially stabilized long-term single hole evolutions. However, we have found that for large domains, neither constraint-subtracted nor constrained $\dot g$ - $\dot K$ evolutions carried out in Cartesian coordinates admit arbitrarily long-lived simulations. The failure appears to arise from features at the inner excision boundary; the behavior does generally improve with resolution.

preprint1998arXiv

Light Propagation in inhomogeneous Universes

Using a multi-plane lensing method that we have developed, we follow the evolution of light beams as they propagate through inhomogeneous universes. We use a P3M code to simulate the formation and evolution of large-scale structure. The resolution of the simulations is increased to sub-Megaparsec scales by using a Monte Carlo method to locate galaxies inside the computational volume according to the underlying particle distribution. The galaxies are approximated by isothermal spheres, with each morphological type having its own distribution of masses and core radii. The morphological types are chosen in order to reproduce the observed morphology-density relation. This algorithm has an effective resolution of 9 orders of magnitudes in length, from the size of superclusters down to the core radii of the smallest galaxies. We consider cold dark matter models normalized to COBE, and perform a large parameter survey by varying the cosmological parameters Omega_0, lambda_0, H_0, and n (the tilt of the primordial power spectrum). The values of n are chosen by imposing particular values or sigma_8, the rms mass fluctuation at a scale of 8/h Mpc. We use the power spectrum given by Bunn & White. This is the largest parameter survey ever done is this field.

preprint1996arXiv

Analysis of ``Gauge Modes'' in Linearized Relativity

By writing the complete set of $3 + 1$ (ADM) equations for linearized waves, we are able to demonstrate the properties of the initial data and of the evolution of a wave problem set by Alcubierre and Schutz. We show that the gauge modes and constraint error modes arise in a straightforward way in the analysis, and are of a form which will be controlled in any well specified convergent computational discretization of the differential equations.

preprint1994arXiv

Spin Dynamics of the LAGEOS Satellite in Support of a Measurement of the Earth's Gravitomagnetism

LAGEOS is an accurately-tracked, dense spherical satellite covered with 426 retroreflectors. The tracking accuracy is such as to yield a medium term (years to decades) inertial reference frame determined via relatively inexpensive observations. This frame is used as an adjunct to the more difficult and data intensive VLBI absolute frame measurements. There is a substantial secular precession of the satellite's line of nodes consistent with the classical, Newtonian precession due to the non-sphericity of the earth. Ciufolini has suggested the launch of an identical satellite (LAGEOS-3) into an orbit supplementary to that of LAGEOS-1: LAGEOS-3 would then experience an equal and opposite classical precession to that of LAGEOS-1. Besides providing a more accurate real-time measurement of the earth's length of day and polar wobble, this paired-satellite experiment would provide the first direct measurement of the general relativistic frame-dragging effect. Of the five dominant error sources in this experiment, the largest one involves surface forces on the satellite, and their consequent impact on the orbital nodal precession. The surface forces are a function of the spin dynamics of the satellite. Consequently, we undertake here a theoretical effort to model the spin ndynamics of LAGEOS. In this paper we present our preliminary results.