Source author record

David Garfinkle

David Garfinkle 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

66works
14topics
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

66 published item(s)

preprint2026arXiv

Asymptotically Euclidean Solutions of the Constraint Equations with Prescribed Asymptotics

We demonstrate that in constructing asymptotically flat vacuum initial data sets in General Relativity via the conformal method, certain asymptotic structures may be prescribed a priori through the specified seed data, including the ADM momentum components, the leading- and next-to-leading-order decay rates, and the anisotropy in the metric's mass term, yielding a recipe to construct initial data sets with desired asymptotics. We numerically construct a simple explicit example of an initial data set, with stronger asymptotics than have been obtained in previous work, such that the evolution of this initial data set does not exhibit the conjectured antipodal symmetry between future and past null infinity.

preprint2023arXiv

An experiment to measure electromagnetic memory

We describe an experiment to measure the electromagnetic analog of gravitational wave memory, the so-called electromagnetic memory. Whereas gravitational wave memory is a residual displacement of test masses, electromagnetic memory is a residual velocity (i.e. kick) of test charges. The source of gravitational wave memory is energy that is not confined to any bounded spatial region: in the case of binary black hole mergers the emitted energy of gravitational radiation as well as the recoil energy of the final black hole. Similarly, electromagnetic memory requires a source whose charges are not confined to any bounded spatial region. While particle beams can provide unbounded charges, their currents are too small to be practical for such an experiment. Instead we propose a short microwave pulse applied to the center of a long dipole antenna. In this way the measurement of the kick can be done quickly enough that the finite size of the antenna does not come into play and it acts for our purposes the same as if it were an infinite antenna.

preprint2022arXiv

A non-trivial PT-symmetric continuum Hamiltonian and its Eigenstates and Eigenvalues

In this paper, a non-trivial system governed by a continuum PT-symmetric Hamiltonian is discussed. We show that this Hamiltonian is iso-spectral to the simple harmonic oscillator. We find its eigenfunctions and the path in the complex plane along which these functions form an orthonormal set. We also find the hidden symmetry operator, ${\cal C}$, for this system. All calculations are performed analytically and without approximation.

preprint2022arXiv

A numerical stability analysis of mean curvature flow of noncompact hypersurfaces with Type-II curvature blowup: II

In previous work [GIKW21], we have presented evidence from numerical simulations that the Type-II singularities of mean curvature flow (MCF) of rotationally-symmetric, complete, noncompact embedded hypersurfaces constructed in [IW19, IWZ21] are stable. More precisely, it is shown in that paper that for small rotationally-symmetric perturbations of initial embeddings near the "tip", numerical simulations of MCF of such initial embeddings develop the same Type-II singularities with the same "bowl soliton" blowup behaviors in a neighborhood of the singularity. It is also shown in that work that for small rotationally-symmetric perturbations of the initial embeddings that are sufficiently far away from the tip, MCF develops Type-I "neckpinch" singularities. In this work, we again use numerical simulations to show that MCF subject to initial perturbations that are not rotationally symmetric behaves asymptotically like it does for rotationally-symmetric perturbations. In particular, if we impose sinusoidal angular dependence on the initial embeddings, we find that for perturbations near the tip, evolutions by MCF asymptotically lose their angular dependence -- becoming round -- and develop Type-II bowl soliton singularities. As well, if we impose sinusoidal angular dependence on the initial embeddings for perturbations sufficiently far from the tip, the angular dependence again disappears as Type-I neckpinch singularities develop. The numerical analysis carried out in this work is an adaptation of the "overlap" method introduced in [GIKW21] and permits angular dependence.

preprint2022arXiv

Gravitational wave memory and the wave equation

Gravitational wave memory and its electromagnetic analog are shown to be straightforward consequences of the wave equation. From Maxwell's equations one can derive a wave equation for the electric field, while from the Bianchi identity one can derive a wave equation for the Riemann tensor in linearized gravity. Memory in both cases is derived from the structure of the source of those wave equations.

preprint2021arXiv

Can Thorne-$\dot{\textrm{Z}}$ytkow Objects source GW190814-type events?

The LIGO-Virgo collaboration reported in their third run the coalescence event GW190814 involving a 2.6 $M_{\odot}$ object with a 23 $M_{\odot}$ black hole. In this letter we study the conditions under which Thorne-$\dot{\textrm{Z}}$ytkow objects (T$\dot{\textrm{Z}}$Os) can be connected to that type of events. We evaluate first the rate of appearance of T$\dot{\textrm{Z}}$Os in the local Universe. Under the assumption that T$\dot{\textrm{Z}}$Os eventually become low mass gap black holes we evaluate how those black holes end up in binaries with other stellar mass black holes and compare to the reported rate for GW190814-type of events (1-23 $\textrm{Gpc}^{-3} \textrm{yr}^{-1}$). We find that T$\dot{\textrm{Z}}$Os in dense stellar clusters can not explain the LIGO-Virgo rate without a T$\dot{\textrm{Z}}$Os population in the field providing a dominant contribution. We also find that T$\dot{\textrm{Z}}$Os formed within hierarchical triple systems in the field with the third more distant star being the progenitor of a stellar mass black hole, may be able to give a rate comparable to that of GW190814-type events. In that case, future observations should discover mergers between stellar mass and low mass gap black holes, with the lower mass spanning the entire low mass gap range.

preprint2020arXiv

Consistent Blandford-Znajek Expansion

The Blandford-Znajek mechanism is the continuous extraction of energy from a rotating black hole via plasma currents flowing on magnetic field lines threading the horizon. In the discovery paper, Blandford and Znajek demonstrated the mechanism by solving the equations of force-free electrodynamics in a perturbative expansion valid at small black hole spin. Attempts to extend this perturbation analysis to higher order have encountered inconsistencies.We overcome this problem using the method of matched asymptotic expansions, taking care to resolve all of the singular surfaces (light surfaces) in the problem. Working with the monopole field configuration, we show explicitly how the inconsistencies are resolved in this framework and calculate the field configuration to one order higher than previously known. However, there is no correction to the energy extraction rate at this order. These results confirm the basic consistency of the split monopole at small spin and lay a foundation for further perturbative studies of the Blandford-Znajek mechanism.

preprint2020arXiv

Cosmological initial data for numerical relativity

We find initial data for numerical relativity simulations of inhomogeneous cosmologies. This involves treating an exceptional case of the general relativity constraint equations. We devise analytic and numerical methods to treat this exceptional case. We apply the analytic method to the standard case of cosmology with a single scalar field. The numerical method is applied to the two-field ekpyrotic cosmology.

preprint2016arXiv

A selection rule for transitions in PT-symmetric quantum theory

Carl Bender and collaborators have developed a quantum theory governed by Hamiltonians that are PT-symmetric rather than Hermitian. To implement this theory, the inner product was redefined to guarantee positive norms of eigenstates of the Hamiltonian. In the general case, which includes arbitrary time-dependence in the Hamiltonian, a modification of the Schrodinger equation is necessary as shown by Gong and Wang to conserve probability. In this paper, we derive the following selection rule: transitions induced by time dependence in a PT-symmetric Hamiltonian cannot occur between normalized states of differing PT-norm. We show three examples of this selection rule in action: two matrix models and one in the continuum.

preprint2016arXiv

Numerical Relativity Beyond Astrophysics

Though the main applications of computer simulations in relativity are to astrophysical systems such as black holes and neutron stars, nonetheless there are important applications of numerical methods to the investigation of general relativity as a fundamental theory of the nature of space and time. This paper gives an overview of some of these applications. In particular we cover (i) investigations of the properties of spacetime singularities such as those that occur in the interior of black holes and in big bang cosmology. (ii) investigations of critical behavior at the threshold of black hole formation in gravitational collapse. (iii) investigations inspired by string theory, in particular analogs of black holes in more than 4 spacetime dimensions and gravitational collapse in spacetimes with a negative cosmological constant.

preprint2015arXiv

Gravitational wave memory in de Sitter spacetime

We examine gravitational wave memory in the case where sources and detector are in an expanding cosmology. For simplicity, we treat the case where the cosmology is de Sitter spacetime, and discuss the possibility of generalizing our results to the case of a more realistic cosmology. We find results very similar to those of gravitational wave memory in an asymptotically flat spacetime, but with the magnitude of the effect multiplied by a redshift factor.

preprint2015arXiv

Gravitational Waves and Their Memory in General Relativity

General relativity explains gravitational radiation from binary black hole or neutron star mergers, from core-collapse supernovae and even from the inflation period in cosmology. These waves exhibit a unique effect called memory or Christodoulou effect, which in a detector like LIGO or LISA shows as a permanent displacement of test masses and in radio telescopes like NANOGrav as a change in the frequency of pulsars' pulses. It was shown that electromagnetic fields and neutrino radiation enlarge the memory. Recently it has been understood that the two types of memory addressed in the literature as `linear' and `nonlinear' are in fact two different phenomena. The former is due to fields that do not and the latter is due to fields that do reach null infinity.

preprint2014arXiv

A perturbative and gauge invariant treatment of gravitational wave memory

We present a perturbative treatment of gravitational wave memory. The coordinate invariance of Einstein's equations leads to a type of gauge invariance in perturbation theory. As with any gauge invariant theory, results are more clear when expressed in terms of manifestly gauge invariant quantities. Therefore we derive all our results from the perturbed Weyl tensor rather than the perturbed metric. We derive gravitational wave memory for the Einstein equations coupled to a general energy-momentum tensor that reaches null infinity.

preprint2014arXiv

Examination of a simple example of gravitational wave memory

We examine a simple example of gravitational wave memory due to the decay of a point particle into two point particles. In the case where one of the decay products is null, there are two types of memory: a null memory due to the null particle and an ordinary memory due to the recoiling timelike particle. In the case where both decay products are timelike, there is only ordinary memory. However, this ordinary memory can mimic the null memory in the limit where one of the decay products has a large velocity.

preprint2013arXiv

An electromagnetic analog of gravitational wave memory

We present an electromagnetic analog of gravitational wave memory. That is, we consider what change has occurred to a detector of electromagnetic radiation after the wave has passed. Rather than a distortion in the detector, as occurs in the gravitational wave case, we find a residual velocity (a "kick") to the charges in the detector. In analogy with the two types of gravitational wave memory ("ordinary" and "nonlinear") we find two types of electromagnetic kick.

preprint2013arXiv

Killing-Yano tensors in spaces admitting a hypersurface orthogonal Killing vector

Methods are presented for finding Killing-Yano tensors, conformal Killing-Yano tensors, and conformal Killing vectors in spacetimes with a hypersurface orthogonal Killing vector. These methods are similar to a method developed by the authors for finding Killing tensors. In all cases one decomposes both the tensor and the equation it satisfies into pieces along the Killing vector and pieces orthogonal to the Killing vector. Solving the separate equations that result from this decomposition requires less computing than integrating the original equation. In each case, examples are given to illustrate the method.

preprint2013arXiv

Linear Stability Analysis and the Speed of Gravitational Waves in Dynamical Chern-Simons Modified Gravity

We perform a linear stability analysis of dynamical Chern-Simons modified gravity in the geometric optics approximation and find that it is linearly stable on the backgrounds considered. Our analysis also reveals that gravitational waves in the modified theory travel at the speed of light in Minkowski spacetime. However, on a Schwarzschild background the characteristic speed of propagation along a given direction splits into two modes, one subluminal and one superluminal. The width of the splitting depends on the azimuthal components of the propagation vector, is linearly proportional to the mass of the black hole, and decreases with the third inverse power of the distance from the black hole. Radial propagation is unaffected, implying that as probed by gravitational waves the location of the event horizon of the spacetime is unaltered. The analysis further reveals that when a high frequency, pure gravitational wave is scattered from a black hole, a scalar wave of comparable amplitude is excited, and vice-versa.

preprint2013arXiv

Matters of Gravity, The Newsletter of the Topical Group on Gravitation of the American Physical Society, Volume 42, Fall 2013

GGR News: we hear that..., by David Garfinkle What's new in LIGO, by David Shoemaker News from NSF, by Pedro Marronetti Citation counts and indicies: Beware of bad data, by Clifford Will Research Briefs: Results from Planck, by William Jones Conference Reports: Gravity- Strings and Higher Dimensions, by Simon Ross Capra 16, by Leor Barack Reduced Order Modeling in General Relativity, by Michele Vallisneri

preprint2013arXiv

Neutrino Radiation Showing a Christodoulou Memory Effect in General Relativity

We describe neutrino radiation in general relativity by introducing the energy-momentum tensor of a null fluid into the Einstein equations. Investigating the geometry and analysis at null infinity, we prove that a component of the null fluid enlarges the Christodoulou memory effect of gravitational waves. The description of neutrinos in general relativity as a null fluid can be regarded as a limiting case of a more general description using the massless limit of the Einstein-Vlasov system. The present authors with co-authors have work in progress to generalize the results of this paper using this more general description. Gigantic neutrino bursts occur in our universe in core-collapse supernovae and in the mergers of neutron star binaries.

preprint2013arXiv

Nonperturbative analysis of the evolution of cosmological perturbations through a nonsingular bounce

In bouncing cosmology, the primordial fluctuations are generated in a cosmic contraction phase before the bounce into the current expansion phase. For a nonsingular bounce, curvature and anisotropy grow rapidly during the bouncing phase, raising questions about the reliability of perturbative analysis. In this paper, we study the evolution of adiabatic perturbations in a nonsingular bounce by nonperturbative methods including numerical simulations of the nonsingular bounce and the covariant formalism for calculating nonlinear perturbations. We show that the bounce is disrupted in regions of the universe with significant inhomogeneity and anisotropy over the background energy density, but is achieved in regions that are relatively homogeneous and isotropic. Sufficiently small perturbations, consistent with observational constraints, can pass through the nonsingular bounce with negligible alteration from nonlinearity. We follow scale invariant perturbations generated in a matter-like contraction phase through the bounce. Their amplitude in the expansion phase is determined by the growing mode in the contraction phase, and the scale invariance is well preserved across the bounce.

preprint2012arXiv

Matters of Gravity, The Newsletter of the Topical Group on Gravitation of the American Physical Society, Volume 39, Winter 2012

GGR News: News From NSF, by David Garfinkle We hear that..., by David Garfinkle GGR Program at the APS meeting in Atlanta, GA, by David Garfinkle Research Briefs: AdS instability, by David Garfinkle Conference reports: Isenberg Fest, by Robert M. Wald COSMO 11, by Carlos Martins Astro-GR 6 2011, by Sascha Husa et al

preprint2012arXiv

Matters of Gravity, The Newsletter of the Topical Group on Gravitation of the American Physical Society, Volume 40, Fall 2012

GGR News: we hear that..., by David Garfinkle 100 years ago, by David Garfinkle New publisher and new book, by Vesselin Petkov Research Briefs: Dark matter news, by Katie Freese LARES satellite, by Richard Matzner Conference reports: Workshop on gravitational wave bursts, by Pablo Laguna JoshFest, by Ed Glass Electromagnetic and gravitational wave astronomy, by Sean McWilliams Bits, branes, and black holes, by Ted Jacobson and Don Marolf

preprint2012arXiv

On Field Theory Thermalization from Gravitational Collapse

Motivated by its field theory interpretation, we study gravitational collapse of a minimally coupled massless scalar field in Einstein gravity with a negative cosmological constant. After demonstrating the accuracy of the numerical algorithm for the questions we are interested in, we investigate various aspects of the apparent horizon formation. In particular, we study the time and radius of the apparent horizon formed as functions of the initial Gaussian profile for the scalar field. We comment on several aspects of the dual field theory picture.

preprint2012arXiv

Rapid Thermalization in Field Theory from Gravitational Collapse

Motivated by the duality with thermalization in field theory, we study gravitational collapse of a minimally coupled massless scalar field in Einstein gravity with a negative cosmological constant. We investigate the system numerically and establish that for small values of the initial amplitude of the scalar field there is no black hole formation, rather, the scalar field performs an oscillatory motion typical of geodesics in AdS. For large enough values of the amplitude of the scalar field we find black hole formation which we detect numerically as the emergence of an apparent horizon. Using the time of formation as an estimate for thermalization in the field theory we conclude that thermalization occurs very rapidly, close to the causal bound for a very wide range of black hole masses. We further study the thermalization time in more detail as a function of the amplitude and the width of the initial Gaussian scalar field profile and detect a rather mild structure.

preprint2011arXiv

A positive energy theorem for Einstein-aether and Hořava gravity

Energy positivity is established for a class of solutions to Einstein-aether theory and the IR limit of Hořava gravity within a certain range of coupling parameters. The class consists of solutions where the aether 4-vector is divergence free on a spacelike surface to which it is orthogonal (which implies that the surface is maximal). In particular, this result holds for spherically symmetric solutions at a moment of time symmetry.

preprint2011arXiv

Gravitational collapse of k-essence

We perform numerical simulations of the gravitational collapse of a k-essence scalar field. When the field is sufficiently strongly gravitating, a black hole forms. However, the black hole has two horizons: a light horizon (the ordinary black hole horizon) and a sound horizon that traps k-essence. In certain cases the k-essence signals can travel faster than light and the sound horizon is inside the light horizon. Under those circumstances, k-essence signals can escape from the black hole. Eventually, the two horizons merge and the k-essence signals can no longer escape.

preprint2011arXiv

Matters of Gravity, The Newsletter of the Topical Group on Gravitation of the American Physical Society, Volume 38, Fall 2011

GGR News: New version of LISA, by Karsten Danzmann We hear that, by David Garfinkle Research Briefs: Finally, results from Gravity Probe B, by Clifford M. Will Discovery of the highest mass neutron star, by Wynn C.G. Ho Conference reports: Numerical relativity beyond astrophysics, by Carsten Gundlac Cold materials, hot nuclei and black holes, by Vijay Balasubramanian Benasque workshop on gravity, by Oscar J.C. Dias

preprint2011arXiv

Non-Stationary Dark Energy Around a Black Hole

Numerical simulations of the accretion of test scalar fields with non-standard kinetic terms (of the k-essence type) onto a Schwarzschild black hole are performed. We find a full dynamical solution for the spherical accretion of a Dirac-Born-Infeld type scalar field. The simulations show that the accretion eventually settles down to a well known stationary solution. This particular analytical steady state solution maintains two separate horizons. The standard horizon is for the usual particles propagating with the limiting speed of light, while the other sonic horizon is for the k-essence perturbations propagating with the speed of sound around this accreting background. For the case where the k-essence perturbations propagate superluminally, we show that one can send signals from within a black hole during the approach to the stationary solution. We also find that a ghost condensate model settles down to a stationary solution during the accretion process.

preprint2010arXiv

Matters of Gravity, The Newsletter of the Topical Group on Gravitation of the American Physical Society, Volume 36, Fall 2010

GGR News: -we hear that..., by David Garfinkle -Network of gravitational-wave detectors, by Stan Whitcomb Research Briefs: -New tests of general relativity, by Quentin Bailey Conference Reports: -Theory meets data analysis, by Steve Detweiler -Ascona conference, by Simon Ross -Condensed Matter at AdS/CFT/Strings 2010, by Christopher Herzog

preprint2009arXiv

Matters of Gravity, The Newsletter of the Topical Group on Gravitation of the American Physical Society, Volume 34, Fall 2009

GGR News: New Matters of Gravity Correspondents, by David Garfinkle we hear that..., by David Garfinkle Research Briefs: Numerical Relativity Studies of Black Hole Kicks, by Pablo Laguna Conference Reports: AbhayFest, by Chris Beetle Eastern Gravity Meeting, by Jeff Winicour Benasque Workshop on Gravity, by Roberto Emparan and Veronika Hubeny Workshop on the Fluid-Gravity Correspondence, by Mukund Rangamani

preprint2009arXiv

The motion of galaxy clusters in inhomogeneous cosmologies

Lemaitre-Tolman-Bondi inhomogeneous spacetimes can be used as a cosmological model to account for the type Ia supernova data. However, such models also give rise to large velocities of galaxy clusters with respect to the cosmic microwave background. Those velocities can be measured using the kinematic Sunyaev-Zeldovich effect. This paper presents a calculation of galaxy cluster velocities as a function of redshift for such a model.

preprint2009arXiv

The Parallelometer: a mechanical device to study curvature

A simple mechanical device is introduced, the parallelometer, that can be used to measure curvatures of surfaces. The device can be used as a practical illustration of parallel transport of a vector and to study Berry phase shift when it is carried along a loop on the surface. Its connection to the Foucault pendulum is discussed. The experimental results can be successfully compared with the theoretical expectations. The experiment is inexpensive and conceptually easy to perform and understand for a beginner.

preprint2006arXiv

Matters of Gravity, The Newsletter of the Topical Group in Gravitation of the American Physical Society, Volume 28, Fall 2006

Research Briefs: Singularity Avoidance in Canonical Quantum Gravity, by Viqar Husain What's New in LIGO, by David Shoemaker Conference reports: Scanning New Horizons: GR Beyond 4 dimensions, by Donald Marolf Quantum Gravity in the Americas III, by Jorge Pullin New Frontiers in Numerical Relativity, by Luciano Rezzolla Teaching General Relativity to Undergraduates, by Greg Comer Ninth Capra Meeting on Radiation Reaction, by Lior Burko

preprint2002arXiv

Harmonic coordinate method for simulating generic singularities

This paper presents both a numerical method for general relativity and an application of that method. The method involves the use of harmonic coordinates in a 3+1 code to evolve the Einstein equations with scalar field matter. In such coordinates, the terms in Einstein's equations with the highest number of derivatives take a form similar to that of the wave equation. The application is an exploration of the generic approach to the singularity for this type of matter. The preliminary results indicate that the dynamics as one approaches the singularity is locally the dynamics of the Kasner spacetimes.

preprint2002arXiv

Perturbations of an exact solution for 2+1 dimensional critical collapse

We find the perturbation spectrum of a family of spherically symmetric and continuously self-similar (CSS) exact solutions that appear to be relevant for the critical collapse of scalar field matter in 2+1 spacetime dimensions. The rate of exponential growth of the unstable perturbation yields the critical exponent. Our results are compared to the numerical simulations of Pretorius and Choptuik and are inconclusive: We find a CSS solution with exactly one unstable mode, which suggests that it may be the critical solution, but another CSS solution which has three unstable modes fits the numerically found critical solution better.

preprint1998arXiv

Scale invariance and critical gravitational collapse

We examine ways to write the Choptuik critical solution as the evolution of scale invariant variables. It is shown that a system of scale invariant variables proposed by one of the authors does not evolve periodically in the Choptuik critical solution. We find a different system, based on maximal slicing. This system does evolve periodically, and may generalize to the case of axisymmetry or of no symmetry at all.

preprint1998arXiv

Scaling of curvature in sub-critical gravitational collapse

We perform numerical simulations of the gravitational collapse of a spherically symmetric scalar field. For those data that just barely do not form black holes we find the maximum curvature at the position of the central observer. We find a scaling relation between this maximum curvature and distance from the critical solution. The scaling relation is analogous to that found by Choptuik for black hole mass for those data that do collapse to form black holes. We also find a periodic wiggle in the scaling exponent.

preprint1997arXiv

Choptuik scaling and the scale invariance of Einstein's equation

The relationship of Choptuik scaling to the scale invariance of Einstein's equation is explored. Ordinary dynamical systems often have limit cycles: periodic orbits that are the asymptotic limit of generic solutions. We show how to separate Einstein's equation into the dynamics of the overall scale and the dynamics of the ``scale invariant'' part of the metric. Periodicity of the scale invariant part implies periodic self-similarity of the spacetime. We also analyze a toy model that exhibits many of the features of Choptuik scaling.

preprint1996arXiv

New Algorithm for Mixmaster Dynamics

We present a new numerical algorithm for evolving the Mixmaster spacetimes. By using symplectic integration techniques to take advantage of the exact Taub solution for the scattering between asymptotic Kasner regimes, we evolve these spacetimes with higher accuracy using much larger time steps than previously possible. The longer Mixmaster evolution thus allowed enables detailed comparison with the Belinskii, Khalatnikov, Lifshitz (BKL) approximate Mixmaster dynamics. In particular, we show that errors between the BKL prediction and the measured parameters early in the simulation can be eliminated by relaxing the BKL assumptions to yield an improved map. The improved map has different predictions for vacuum Bianchi Type IX and magnetic Bianchi Type VI$_0$ Mixmaster models which are clearly matched in the simulation.

preprint1996arXiv

Solution generating with perfect fluids

We apply a technique, due to Stephani, for generating solutions of the Einstein-perfect fluid equations. This technique is similar to the vacuum solution generating techniques of Ehlers, Harrison, Geroch and others. We start with a ``seed'' solution of the Einstein-perfect fluid equations with a Killing vector. The seed solution must either have (i) a spacelike Killing vector and equation of state P=rho or (ii) a timelike Killing vector and equation of state rho+3P=0. The new solution generated by this technique then has the same Killing vector and the same equation of state. We choose several simple seed solutions with these equations of state and where the Killing vector has no twist. The new solutions are twisting versions of the seed solutions.

preprint1994arXiv

Choptuik scaling in null coordinates

A numerical simulation is performed of the gravitational collapse of a spherically symmetric scalar field. The algorithm uses the null initial value formulation of the Einstein-scalar equations, but does {\it not} use adaptive mesh refinement. A study is made of the critical phenomena found by Choptuik in this system. In particular it is verified that the critical solution exhibits periodic self-similarity. This work thus provides a simple algorithm that gives verification of the Choptuik results.