Source author record

Stephen R. Green

Stephen R. Green 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

26works
12topics
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

26 published item(s)

preprint2025arXiv

GW200208_222617 as an eccentric black-hole binary merger: properties and astrophysical implications

Detecting orbital eccentricity in a stellar-mass black-hole merger would point to a non-isolated formation channel. Eccentric binaries can form in dense stellar environments such as globular clusters or active galactic nuclei, or from triple stellar systems in the Galactic field. However, confidently measuring eccentricity is challenging -- short signals from high-mass eccentric mergers can mimic spin-induced precession, making the two effects hard to disentangle. This degeneracy weakens considerably for longer-duration signals. Here, GW200208_222617 provides a rare opportunity. Originating from a relatively low-mass binary with source-frame chirp mass $\sim20$ M$_\odot$, its gravitational-wave signal spanned $\sim14$ orbital cycles in band, with no indication of data quality issues. Previous analyses for quasi-circular binaries found no evidence for spin precession, and multiple subsequent studies found the data to favour an eccentric merger despite notable technical differences. All in all, we believe GW200208_222617 is the black-hole merger event from GWTC-3 with the least ambiguous detection of eccentricity. We present a critical discussion of properties and astrophysical interpretation of GW200208_222617 as an eccentric black-hole merger using models of field triples, globular clusters, and active galactic nuclei. We find that if GW200208_222617 was indeed eccentric, its origin is consistent with a field triple or globular cluster. Formation in the inner regions of an active galactic nucleus is disfavoured. The outer regions of such a disk remain a viable origin for GW200208_222617; we demonstrate how future detections of eccentric mergers formed in such environments could be powerful tools for constraining the disk geometry.

preprint2022arXiv

Gamma-radiation sky maps from compact binaries

We study sky maps and light curves of gamma-ray emission from neutron stars in compact binaries, and in isolation. We briefly review some gamma-ray emission models, and reproduce sky maps from a standard isolated pulsar in the Separatrix Layer model. We consider isolated pulsars with several variations of a dipole magnetic field, including superpositions, and predict their gamma-ray emission. Our results provide new heuristics on what can and cannot be inferred about the magnetic field configuration of pulsars from high-energy observations. We find that typical double-peak light curves can be produced by pulsars with significant multipole structure beyond a single dipole. For binary systems, we also present a simple approximation that is useful for rapid explorations of binary magnetic field structure. Finally, we predict the gamma-ray emission pattern from a compact black hole-neutron star binary moments before merger by applying the Separatrix Layer model to data simulated in full general relativity; we find that face-on observers receive little emission, equatorial observers see one broad peak, and more generic observers typically see two peaks.

preprint2022arXiv

Using machine learning to parametrize postmerger signals from binary neutron stars

There is growing interest in the detection and characterization of gravitational waves from postmerger oscillations of binary neutron stars. These signals contain information about the nature of the remnant and the high-density and out-of-equilibrium physics of the postmerger processes, which would complement any electromagnetic signal. However, the construction of binary neutron star postmerger waveforms is much more complicated than for binary black holes: (i) there are theoretical uncertainties in the neutron-star equation of state and other aspects of the high-density physics, (ii) numerical simulations are expensive and available ones only cover a small fraction of the parameter space with limited numerical accuracy, and (iii) it is unclear how to parametrize the theoretical uncertainties and interpolate across parameter space. In this work, we describe the use of a machine-learning method called a conditional variational autoencoder (CVAE) to construct postmerger models for hyper/massive neutron star remnant signals based on numerical-relativity simulations. The CVAE provides a probabilistic model, which encodes uncertainties in the training data within a set of latent parameters. We estimate that training such a model will ultimately require $\sim 10^4$ waveforms. However, using synthetic training waveforms as a proof-of-principle, we show that the CVAE can be used as an accurate generative model and that it encodes the equation of state in a useful latent representation.

preprint2021arXiv

New metric reconstruction scheme for gravitational self-force calculations

Inspirals of stellar-mass objects into massive black holes will be important sources for the space-based gravitational-wave detector LISA. Modelling these systems requires calculating the metric perturbation due to a point particle orbiting a Kerr black hole. Currently, the linear perturbation is obtained with a metric reconstruction procedure that puts it in a "no-string" radiation gauge which is singular on a surface surrounding the central black hole. Calculating dynamical quantities in this gauge involves a subtle procedure of "gauge completion" as well as cancellations of very large numbers. The singularities in the gauge also lead to pathological field equations at second perturbative order. In this paper we re-analyze the point-particle problem in Kerr using the corrector-field reconstruction formalism of Green, Hollands, and Zimmerman (GHZ). We clarify the relationship between the GHZ formalism and previous reconstruction methods, showing that it provides a simple formula for the "gauge completion". We then use it to develop a new method of computing the metric in a more regular gauge: a Teukolsky puncture scheme. This scheme should ameliorate the problem of large cancellations, and by constructing the linear metric perturbation in a sufficiently regular gauge, it should provide a first step toward second-order self-force calculations in Kerr. Our methods are developed in generality in Kerr, but we illustrate some key ideas and demonstrate our puncture scheme in the simple setting of a static particle in Minkowski spacetime.

preprint2021arXiv

Nonlinear effects in the black hole ringdown: absorption-induced mode excitation

Gravitational-wave observations of black hole ringdowns are commonly used to characterize binary merger remnants and to test general relativity. These analyses assume linear black hole perturbation theory, in particular that the ringdown can be described in terms of quasinormal modes even for times approaching the merger. Here we investigate a nonlinear effect during the ringdown, namely how a mode excited at early times can excite additional modes as it is absorbed by the black hole. This is a third-order secular effect: the change in the black-hole mass causes a shift in the mode spectrum, so that the original mode is projected onto the new ones. Using nonlinear simulations, we study the ringdown of a spherically-symmetric scalar field around an asymptotically anti-de Sitter black hole, and we find that this "absorption-induced mode excitation" (AIME) is the dominant nonlinear effect. We show that this effect takes place well within the nonadiabatic regime, so we can analytically estimate it using a sudden mass-change approximation. Adapting our estimation technique to asymptotically-flat Schwarzschild black holes, we expect AIME to play a role in the analysis and interpretation of current and future gravitational wave observations.

preprint2020arXiv

Complete parameter inference for GW150914 using deep learning

The LIGO and Virgo gravitational-wave observatories have detected many exciting events over the past five years. As the rate of detections grows with detector sensitivity, this poses a growing computational challenge for data analysis. With this in mind, in this work we apply deep learning techniques to perform fast likelihood-free Bayesian inference for gravitational waves. We train a neural-network conditional density estimator to model posterior probability distributions over the full 15-dimensional space of binary black hole system parameters, given detector strain data from multiple detectors. We use the method of normalizing flows---specifically, a neural spline normalizing flow---which allows for rapid sampling and density estimation. Training the network is likelihood-free, requiring samples from the data generative process, but no likelihood evaluations. Through training, the network learns a global set of posteriors: it can generate thousands of independent posterior samples per second for any strain data consistent with the prior and detector noise characteristics used for training. By training with the detector noise power spectral density estimated at the time of GW150914, and conditioning on the event strain data, we use the neural network to generate accurate posterior samples consistent with analyses using conventional sampling techniques.

preprint2020arXiv

Teukolsky formalism for nonlinear Kerr perturbations

We develop a formalism to treat higher order (nonlinear) metric perturbations of the Kerr spacetime in a Teukolsky framework. We first show that solutions to the linearized Einstein equation with nonvanishing stress tensor can be decomposed into a pure gauge part plus a zero mode (infinitesimal perturbation of the mass and spin) plus a perturbation arising from a certain scalar ("Debye-Hertz") potential, plus a so-called "corrector tensor." The scalar potential is a solution to the spin $-2$ Teukolsky equation with a source. This source, as well as the tetrad components of the corrector tensor, are obtained by solving certain decoupled ordinary differential equations involving the stress tensor. As we show, solving these ordinary differential equations reduces simply to integrations in the coordinate $r$ in outgoing Kerr-Newman coordinates, so in this sense, the problem is reduced to the Teukolsky equation with source, which can be treated by a separation of variables ansatz. Since higher order perturbations are subject to a linearized Einstein equation with a stress tensor obtained from the lower order perturbations, our method also applies iteratively to the higher order metric perturbations, and could thus be used to analyze the nonlinear coupling of perturbations in the near-extremal Kerr spacetime, where weakly turbulent behavior has been conjectured to occur. Our method could also be applied to the study of perturbations generated by a pointlike body traveling on a timelike geodesic in Kerr, which is relevant to the extreme mass ratio inspiral problem.

preprint2019arXiv

Excited hairy black holes: dynamical construction and level transitions

We study the dynamics of unstable Reissner-Nordström anti-de Sitter black holes under charged scalar field perturbations in spherical symmetry. We unravel their general behavior and approach to the final equlibrium state. In the first part of this work, we present a numerical analysis of massive charged scalar field quasinormal modes. We identify the known mode families - superradiant modes, zero-damped modes, AdS modes, and the near-horizon mode - and we track their migration under variation of the black hole and field parameters. We show that the zero-damped modes become superradiantly unstable for large RNAdS with large gauge coupling; the leading unstable mode is identified with the near-horizon condensation instability. In the second part, we present results of numerical simulations of perturbed large RNAdS, showing the nonlinear development of these unstable modes. For generic initial conditions, charge and mass are ransferred from the black hole to the scalar field, until an equilibrium solution with a scalar condensate is reached. We use results from the linear analysis, however, to select special initial data corresponding to an unstable overtone mode. We find that these data evolve to produce a new equilibrium state - an excited hairy black hole with the scalar condensate in an overtone configuration. This state is, however, unstable, and the black hole eventually decays to the generic end state. Nevertheless, this demonstrates the potential relevance of overtone modes as transients in black hole dynamics.

preprint2016arXiv

A Simple, Heuristic Derivation of our "No Backreaction" Results

We provide a simple discussion of our results on the backreaction effects of density inhomogeneities in cosmology, without mentioning one-parameter families or weak limits. Emphasis is placed on the manner in which "averaging" is done and the fact that one is solving Einstein's equation. The key assumptions and results that we rigorously derived within our original mathematical framework are thereby explained in a heuristic way.

preprint2016arXiv

Nonlinear Evolution and Final Fate of Charged Anti-de Sitter Black Hole Superradiant Instability

We describe the full nonlinear development of the superradiant instability for a charged massless scalar field, coupled to general relativity and electromagnetism, in the vicinity of a Reissner-Nordstrom-AdS black hole. The presence of the negative cosmological constant provides a natural context for considering perfectly reflecting boundary conditions and studying the dynamics as the scalar field interacts repeatedly with the black hole. At early times, small superradiant perturbations grow as expected from linearized studies. Backreaction then causes the black hole to lose charge and mass until the perturbation becomes nonsuperradiant, with the final state described by a stable hairy black hole. For large gauge coupling, the instability extracts a large amount of charge per unit mass, resulting in greater entropy increase. We discuss the implications of the observed behavior for the general problem of superradiance in black hole spacetimes.

preprint2015arXiv

Conserved quantities and dual turbulent cascades in Anti-de Sitter spacetime

We consider the dynamics of a spherically symmetric massless scalar field coupled to general relativity in Anti--de Sitter spacetime in the small-amplitude limit. Within the context of our previously developed two time framework (TTF) to study the leading self-gravitating effects, we demonstrate the existence of two new conserved quantities in addition to the known total energy $E$ of the modes: The particle number $N$ and Hamiltonian $H$ of our TTF system. Simultaneous conservation of $E$ and $N$ implies that weak turbulent processes undergo dual cascades (direct cascade of $E$ and inverse cascade of $N$ or vice versa). This partially explains the observed dynamics of 2-mode initial data. In addition, conservation of $E$ and $N$ limits the region of phase space that can be explored within the TTF approximation and in particular rules out equipartion of energy among the modes for general initial data. Finally, we discuss possible effects of conservation of $N$ and $E$ on late time dynamics.

preprint2015arXiv

Coupled Oscillator Model for Nonlinear Gravitational Perturbations

Motivated by the gravity/fluid correspondence, we introduce a new method for characterizing nonlinear gravitational interactions. Namely we map the nonlinear perturbative form of the Einstein equation to the equations of motion of a collection of nonlinearly-coupled harmonic oscillators. These oscillators correspond to the quasinormal or normal modes of the background spacetime. We demonstrate the mechanics and the utility of this formalism within the context of perturbed asymptotically anti-de Sitter black brane spacetimes. We confirm in this case that the boundary fluid dynamics are equivalent to those of the hydrodynamic quasinormal modes of the bulk spacetime. We expect this formalism to remain valid in more general spacetimes, including those without a fluid dual. In other words, although borne out of the gravity/fluid correspondence, the formalism is fully independent and it has a much wider range of applicability. In particular, as this formalism inspires an especially transparent physical intuition, we expect its introduction to simplify the often highly technical analytical exploration of nonlinear gravitational dynamics.

preprint2015arXiv

Islands of stability and recurrence times in AdS

We study the stability of anti-de Sitter (AdS) spacetime to spherically symmetric perturbations of a real scalar field in general relativity. Further, we work within the context of the "two time framework" (TTF) approximation, which describes the leading nonlinear effects for small amplitude perturbations, and is therefore suitable for studying the weakly turbulent instability of AdS---including both collapsing and non-collapsing solutions. We have previously identified a class of quasi-periodic (QP) solutions to the TTF equations, and in this work we analyze their stability. We show that there exist several families of QP solutions that are stable to linear order, and we argue that these solutions represent islands of stability in TTF. We extract the eigenmodes of small oscillations about QP solutions, and we use them to predict approximate recurrence times for generic non-collapsing initial data in the full (non-TTF) system. Alternatively, when sufficient energy is driven to high-frequency modes, as occurs for initial data far from a QP solution, the TTF description breaks down as an approximation to the full system. Depending on the higher order dynamics of the full system, this often signals an imminent collapse to a black hole.

preprint2015arXiv

Superradiant instabilities of asymptotically anti-de Sitter black holes

We study the linear stability of asymptotically anti-de Sitter black holes in general relativity in spacetime dimension $d\ge4$. Our approach is an adaptation of the general framework of Hollands and Wald, which gives a stability criterion in terms of the sign of the canonical energy, $\mathcal{E}$. The general framework was originally formulated for static or stationary and axisymmetric black holes in the asymptotically flat case, and the stability analysis for that case applies only to axisymmetric perturbations. However, in the asymptotically anti-de Sitter case, the stability analysis requires only that the black hole have a single Killing field normal to the horizon and there are no restrictions on the perturbations (apart from smoothness and appropriate behavior at infinity). For an asymptotically anti-de Sitter black hole, we define an ergoregion to be a region where the horizon Killing field is spacelike; such a region, if present, would normally occur near infinity. We show that for black holes with ergoregions, initial data can be constructed such that $\mathcal{E}<0$, so all such black holes are unstable. To obtain such initial data, we first construct an approximate solution to the constraint equations using the WKB method, and then we use the Corvino-Schoen technique to obtain an exact solution. We also discuss the case of charged asymptotically anti-de Sitter black holes with generalized ergoregions.

preprint2015arXiv

Universality of non-equilibrium dynamics of CFTs from holography

Motivated by a low-energy effective description of gauge theory/string theory duality, we conjecture that the dynamics of $SO(4)$-invariant states in a large class of four-dimensional conformal gauge theories on $S^3$ with non-equal central charges $c\ne a$ are universal on time scales $t_{\rm universal}\propto ({\cal E}-{\cal E}_{\rm vacuum})^{-1}$, in the limit where the energy ${\cal E}\to {\cal E}_{\rm vacuum}$. We show that low-energy excitations in $c\ne a$ CFTs do not thermalize in this limit. The holographic universality conjecture then implies that within the Einstein-scalar field system (dual to theories with $c=a$), $AdS_5$ is stable to spherically symmetric perturbations against formation of trapped surfaces within time scales $t_{\rm universal}$.

preprint2014arXiv

Holographic Thermalization, stability of AdS, and the Fermi-Pasta-Ulam-Tsingou paradox

For a real massless scalar field in general relativity with a negative cosmological constant, we uncover a large class of spherically symmetric initial conditions that are close to AdS, but whose numerical evolution does not result in black hole formation. According to the AdS/CFT dictionary, these bulk solutions are dual to states of a strongly interacting boundary CFT that fail to thermalize at late times. Furthermore, as these states are not stationary, they define dynamical CFT configurations that do not equilibrate. We develop a two-timescale perturbative formalism that captures both direct and inverse cascades of energy and agrees with our fully nonlinear evolutions in the appropriate regime. We also show that this formalism admits a large class of quasi-periodic solutions. Finally, we demonstrate a striking parallel between the dynamics of AdS and the classic Fermi-Pasta-Ulam-Tsingou problem.

preprint2014arXiv

How well is our universe described by an FLRW model?

Extremely well! In the $Λ$CDM model, the spacetime metric, $g_{ab}$, of our universe is approximated by an FLRW metric, $g_{ab}^{(0)}$, to about 1 part in $10^4$ or better on both large and small scales, except in the immediate vicinity of very strong field objects, such as black holes. However, derivatives of $g_{ab}$ are not close to derivatives of $g_{ab}^{(0)}$, so there can be significant differences in the behavior of geodesics and huge differences in curvature. Consequently, observable quantities in the actual universe may differ significantly from the corresponding observables in the FLRW model. Nevertheless, as we shall review here, we have proven general results showing that---within the framework of our approach to treating backreaction---the large matter inhomogeneities that occur on small scales cannot produce significant effects on large scales, so $g_{ab}^{(0)}$ satisfies Einstein's equation with the averaged stress-energy tensor of matter as its source. We discuss the flaws in some other approaches that have suggested that large backreaction effects may occur. As we also will review here, with a suitable "dictionary," Newtonian cosmologies provide excellent approximations to cosmological solutions to Einstein's equation (with dust and a cosmological constant) on all scales. Our results thereby provide strong justification for the mathematical consistency and validity of the $Λ$CDM model within the context of general relativistic cosmology.

preprint2013arXiv

A Holographic Path to the Turbulent Side of Gravity

We study the dynamics of a 2+1 dimensional relativistic viscous conformal fluid in Minkowski spacetime. Such fluid solutions arise as duals, under the "gravity/fluid correspondence", to 3+1 dimensional asymptotically anti-de Sitter (AAdS) black brane solutions to the Einstein equation. We examine stability properties of shear flows, which correspond to hydrodynamic quasinormal modes of the black brane. We find that, for sufficiently high Reynolds number, the solution undergoes an inverse turbulent cascade to long wavelength modes. We then map this fluid solution, via the gravity/fluid duality, into a bulk metric. This suggests a new and interesting feature of the behavior of perturbed AAdS black holes and black branes, which is not readily captured by a standard quasinormal mode analysis. Namely, for sufficiently large perturbed black objects (with long-lived quasinormal modes), nonlinear effects transfer energy from short to long wavelength modes via a turbulent cascade within the metric perturbation. As long wavelength modes have slower decay, this lengthens the overall lifetime of the perturbation. We also discuss various implications of this behavior, including expectations for higher dimensions, and the possibility of predicting turbulence in more general gravitational scenarios.

preprint2013arXiv

Dynamic and Thermodynamic Stability of Relativistic, Perfect Fluid Stars

We consider perfect fluid bodies (stars) in general relativity, characterized by particle number density and entropy per particle. A star is said to be in dynamic equilibrium if it is a stationary, axisymmetric solution to the Einstein-fluid equations with circular flow. We prove that for a star in dynamic equilibrium, the necessary and sufficient condition for thermodynamic equilibrium (extremum of total entropy S) is constancy of angular velocity (Ω), redshifted temperature, and redshifted chemical potential. Friedman previously identified positivity of canonical energy, E, as a criterion for dynamic stability and argued that all rotating stars are dynamically unstable to sufficiently nonaxisymmetric perturbations (the CFS instability), so our main focus is on axisymmetric stability. We show that for a star in dynamic equilibrium, mode stability holds with respect to all axisymmetric perturbations if E is positive on a certain subspace, V, of axisymmetric Lagrangian perturbations that, in particular, have no Lagrangian change in angular momentum density. Conversely, if E fails to be positive on V, then there exist perturbations that can't become asymptotically stationary at late times. We further show that for a star in thermodynamic equilibrium, the canonical energy in the rotating frame, E_r, is related to second order changes in ADM mass, M, and angular momentum, J, by E_r = δ^2 M - Ωδ^2 J. Thus, positivity of E_r for perturbations with δJ = 0 is a necessary condition for thermodynamic stability (local maximum of S). For axisymmetric perturbations we have E = E_r, so a necessary condition for thermodynamic stability with respect to axisymmetric perturbations is positivity of E on all perturbations with δJ = 0, not merely on the perturbations in V. Many of our results are in close parallel to results of Hollands and Wald for the theory of black holes.

preprint2013arXiv

Examples of backreaction of small scale inhomogeneities in cosmology

In previous work, we introduced a new framework to treat large scale backreaction effects due to small scale inhomogeneities in general relativity. We considered one-parameter families of spacetimes for which such backreaction effects can occur, and we proved that, provided the weak energy condition on matter is satisfied, the leading effect of small scale inhomogeneities on large scale dynamics is to produce a traceless effective stress-energy tensor that itself satisfies the weak energy condition. In this work, we illustrate the nature of our framework by providing two explicit examples of one-parameter families with backreaction. The first, based on previous work of Berger, is a family of polarized vacuum Gowdy spacetimes on a torus, which satisfies all of the assumptions of our framework. As the parameter approaches its limiting value, the metric uniformly approaches a smooth background metric, but spacetime derivatives of the deviation of the metric from the background metric do not converge uniformly to zero. The limiting metric has nontrivial backreaction from the small scale inhomogeneities, with an effective stress-energy that is traceless and satisfies the weak energy condition, in accord with our theorems. Our second one-parameter family consists of metrics which have a uniform Friedmann-Lemaitre-Robertson-Walker limit. This family satisfies all of our assumptions with the exception of the weak energy condition for matter. In this case, the limiting metric has an effective stress-energy tensor which is not traceless. We emphasize the importance of imposing energy conditions on matter in studies of backreaction.

preprint2012arXiv

Constraints on String Cosmology

String theory contains sources like orientifold planes that support higher derivative interactions. These interactions make possible static flux compactifications which are forbidden in supergravity. They can also lead to violations of the strong energy condition (SEC) which is needed for an accelerating universe. We examine how large a violation is possible in the context of the heterotic string compactified to four dimensions. We find that de Sitter solutions are still not possible but that classically forbidden anti-de Sitter solutions are possible.

preprint2012arXiv

Newtonian and Relativistic Cosmologies

Cosmological N-body simulations are now being performed using Newtonian gravity on scales larger than the Hubble radius. It is well known that a uniformly expanding, homogeneous ball of dust in Newtonian gravity satisfies the same equations as arise in relativistic FLRW cosmology, and it also is known that a correspondence between Newtonian and relativistic dust cosmologies continues to hold in linearized perturbation theory in the marginally bound/spatially flat case. Nevertheless, it is far from obvious that Newtonian gravity can provide a good global description of an inhomogeneous cosmology when there is significant nonlinear dynamical behavior at small scales. We investigate this issue in the light of a perturbative framework that we have recently developed, which allows for such nonlinearity at small scales. We propose a relatively straightforward "dictionary"---which is exact at the linearized level---that maps Newtonian dust cosmologies into general relativistic dust cosmologies, and we use our "ordering scheme" to determine the degree to which the resulting metric and matter distribution solve Einstein's equation. We find that Einstein's equation fails to hold at "order 1" at small scales and at "order $ε$" at large scales. We then find the additional corrections to the metric and matter distribution needed to satisfy Einstein's equation to these orders. While these corrections are of some interest in their own right, our main purpose in calculating them is that their smallness should provide a criterion for the validity of the original dictionary (as well as simplified versions of this dictionary). We expect that, in realistic Newtonian cosmologies, these additional corrections will be very small; if so, this should provide strong justification for the use of Newtonian simulations to describe relativistic cosmologies, even on scales larger than the Hubble radius.

preprint2011arXiv

A new framework for analyzing the effects of small scale inhomogeneities in cosmology

We develop a new, mathematically precise framework for treating the effects of nonlinear phenomena occurring on small scales in general relativity. Our approach is an adaptation of Burnett's formulation of the "shortwave approximation", which we generalize to analyze the effects of matter inhomogeneities as well as gravitational radiation. Our framework requires the metric to be close to a "background metric", but allows arbitrarily large stress-energy fluctuations on small scales. We prove that, within our framework, if the matter stress-energy tensor satisfies the weak energy condition (i.e., positivity of energy density in all frames), then the only effect that small scale inhomogeneities can have on the dynamics of the background metric is to provide an "effective stress-energy tensor" that is traceless and has positive energy density---corresponding to the presence of gravitational radiation. In particular, nonlinear effects produced by small scale inhomogeneities cannot mimic the effects of dark energy. We also develop "perturbation theory" off of the background metric. We derive an equation for the "long-wavelength part" of the leading order deviation of the metric from the background metric, which contains the usual terms occurring in linearized perturbation theory plus additional contributions from the small-scale inhomogeneities. Under various assumptions concerning the absence of gravitational radiation and the non-relativistic behavior of the matter, we argue that the "short wavelength" deviations of the metric from the background metric near a point $x$ should be accurately described by Newtonian gravity, taking into account only the matter lying within a "homogeneity lengthscale" of $x$. Finally, we argue that our framework should provide an accurate description of the actual universe.

preprint2010arXiv

Numerical Parameter Survey of Nonradiative Black Hole Accretion -- Flow Structure and Variability of the Rotation Measure

We conduct a survey of numerical simulations to probe the structure and appearance of non-radiative black hole accretion flows like the supermassive black hole at the Galactic centre. We find a generic set of solutions, and make specific predictions for currently feasible rotation measure (RM) observations, which are accessible to current instruments including the EVLA, GMRT and ALMA. The slow time variability of the RM is a key quantitative signature of this accretion flow. The time variability of RM can be used to quantitatively measure the nature of the accretion flow, and to differentiate models. Sensitive measurements of RM can be achieved using RM synthesis or using pulsars. Our energy conserving ideal magneto-hydrodynamical simulations, which achieve high dynamical range by means of a deformed-mesh algorithm, stretch from several Bondi radii to about one thousandth of that radius, and continue for tens of Bondi times. Magnetized flows which lack outward convection possess density slopes around -1, almost independent of physical parameters, and are more consistent with observational constraints than are strongly convective flows We observe no tendency for the flows to become rotationally supported in their centres, or to develop steady outflow. We support these conclusions with formulae which encapsulate our findings in terms of physical and numerical parameters. We discuss the relation of these solutions to other approaches. The main potential uncertainties are the validity of ideal MHD and the absence of a fully relativistic inner boundary condition. The RM variability predictions are testable with current and future telescopes.