Source author record

Michael Pürrer

Michael Pürrer 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

20works
4topics
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

20 published item(s)

preprint2021arXiv

An aligned-spin neutron-star--black-hole waveform model based on the effective-one-body approach and numerical-relativity simulations

After the discovery of gravitational waves from binary black holes (BBHs) and binary neutron stars (BNSs) with the LIGO and Virgo detectors, neutron-star--black-holes (NSBHs) are the natural next class of binary systems to be observed. In this work, we develop a waveform model for aligned-spin neutron-star--black-holes (NSBHs) combining a BBH baseline waveform (available in the effective-one-body approach) with a phenomenological description of tidal effects (extracted from numerical-relativity simulations), and correcting the amplitude during the late inspiral, merger and ringdown to account for the NS tidal disruption. We calibrate the amplitude corrections using NSBH waveforms obtained with the SpEC and the SACRA codes. The model was calibrated using simulations with NS masses in the range $1.2-1.4 M_\odot$, tidal deformabilities up to $4200$ (for a 1.2 $M_\odot$ NS), and dimensionless BH spin magnitude up to 0.9. Based on the simulations used, and on checking that sensible waveforms are produced, we recommend our model to be employed with NS mass in the range $1\mbox{--}3 M_\odot$, tidal deformability $0\mbox{--}5000$, and BH spin magnitude up to $0.9$. We also validate our model against two new, highly accurate NSBH waveforms with BH spin 0.9 and mass ratios 3 and 4, characterized by tidal disruption, produced with SpEC, and find very good agreement. We find that it will be challenging for the advanced LIGO-Virgo--detector network at design sensitivity to distinguish different source classes. We perform parameter-estimation on a synthetic numerical-relativity signal in zero noise to study parameter biases. Finally, we reanalyze GW170817, with the hypothesis that it is a NSBH. We do not find evidence to distinguish the BNS and NSBH hypotheses, however the posterior for the mass ratio is shifted to less equal masses under the NSBH hypothesis. [Abstract abridged for arxiv].

preprint2020arXiv

Frequency domain reduced order model of aligned-spin effective-one-body waveforms with higher-order modes

We present a frequency domain reduced order model (ROM) for the aligned-spin effective-one-body (EOB) model for binary black holes (BBHs) SEOBNRv4HM that includes the spherical harmonics modes $(\ell, |m|) = (2,1),(3,3),(4,4),(5,5)$ beyond the dominant $(\ell, |m|) = (2,2)$ mode. These higher modes are crucial to accurately represent the waveform emitted from asymmetric BBHs. We discuss a decomposition of the waveform, extending other methods in the literature, that allows us to accurately and efficiently capture the morphology of higher mode waveforms. We show that the ROM is very accurate with median (maximum) values of the unfaithfulness against SEOBNRv4HM lower than $0.001\% (0.03\%)$ for total masses in $[2.8,100] M_\odot$. For a total mass of $M = 300 M_\odot$ the median (maximum) value of the unfaithfulness increases up to $0.004\% (0.17\%)$. This is still at least an order of magnitude lower than the estimated accuracy of SEOBNRv4HM compared to numerical relativity simulations. The ROM is two orders of magnitude faster in generating a waveform compared to SEOBNRv4HM. Data analysis applications typically require $\mathcal{O}(10^6-10^8)$ waveform evaluations for which SEOBNRv4HM is in general too slow. The ROM is therefore crucial to allow the SEOBNRv4HM waveform to be used in searches and Bayesian parameter inference. We present a targeted parameter estimation study that shows the improvements in measuring binary parameters when using waveforms that includes higher modes and compare against three other waveform models.

preprint2020arXiv

Multipolar Effective-One-Body Waveforms for Precessing Binary Black Holes: Construction and Validation

As gravitational-wave detectors become more sensitive, we will access a greater variety of signals emitted by compact binary systems, shedding light on their astrophysical origin and environment. A key physical effect that can distinguish among formation scenarios is the misalignment of the spins with the orbital angular momentum, causing the spins and the binary's orbital plane to precess. To accurately model such systems, it is crucial to include multipoles beyond the dominant quadrupole. Here, we develop the first multipolar precessing waveform model in the effective-one-body (EOB) formalism for the inspiral, merger and ringdown (IMR) of binary black holes: SEOBNRv4PHM. In the nonprecessing limit, the model reduces to SEOBNRv4HM, which was calibrated to numerical-relativity (NR) simulations, and waveforms from perturbation theory. We validate SEOBNRv4PHM by comparing it to the public catalog of 1405 precessing NR waveforms of the Simulating eXtreme Spacetimes (SXS) collaboration, and also to new 118 precessing NR waveforms, which span mass ratios 1-4 and spins up to 0.9. We stress that SEOBNRv4PHM is not calibrated to NR simulations in the precessing sector. We compute the unfaithfulness against the 1523 SXS precessing NR waveforms, and find that, for $94\%$ ($57\%$) of the cases, the maximum value, in the total mass range $20-200 M_\odot$, is below $3\%$ ($1\%$). Those numbers become $83\%$ ($20\%$) when using the IMR, multipolar, precessing phenomenological model IMRPhenomPv3HM. We investigate the impact of such unfaithfulness values with two parameter-estimation studies on synthetic signals. We also compute the unfaithfulness between those waveform models and identify in which part of the parameter space they differ the most. We validate them also against the multipolar, precessing NR surrogate model NRSur7dq4, and find that the SEOBNRv4PHM model outperforms IMRPhenomPv3HM.

preprint2019arXiv

Regression methods in waveform modeling: a comparative study

Gravitational-wave astronomy of compact binaries relies on theoretical models of the gravitational-wave signal that is emitted as binaries coalesce. These models do not only need to be accurate, they also have to be fast to evaluate in order to be able to compare millions of signals in near real time with the data of gravitational-wave instruments. A variety of regression and interpolation techniques have been employed to build efficient waveform models, but no study has systematically compared the performance of these regression methods yet. Here we provide such a comparison of various techniques, including polynomial fits, radial basis functions, Gaussian process regression and artificial neural networks, specifically for the case of gravitational waveform modeling. We use all these techniques to regress analytical models of non-precessing and precessing binary black hole waveforms, and compare the accuracy as well as computational speed. We find that most regression methods are reasonably accurate, but efficiency considerations favour in many cases the most simple approach. We conclude that sophisticated regression methods are not necessarily needed in standard gravitational-wave modeling applications, although problems with higher complexity than what is tested here might be more suitable for machine-learning techniques and more sophisticated methods may have side benefits.

preprint2016arXiv

Fast and Accurate Inference on Gravitational Waves from Precessing Compact Binaries

Inferring astrophysical information from gravitational waves emitted by compact binaries is one of the key science goals of gravitational-wave astronomy. In order to reach the full scientific potential of gravitational-wave experiments we require techniques to mitigate the cost of Bayesian inference, especially as gravitational-wave signal models and analyses become increasingly sophisticated and detailed. Reduced order models (ROMs) of gravitational waveforms can significantly reduce the computational cost of inference by removing redundant computations. In this paper we construct the first reduced order models of gravitational-wave signals that include the effects of spin-precession, inspiral, merger, and ringdown in compact object binaries, and which are valid for component masses describing binary neutron star, binary black hole and mixed binary systems. This work utilizes the waveform model known as "IMRPhenomPv2". Our ROM enables the use of a fast reduced order quadrature (ROQ) integration rule which allows us to approximate Bayesian probability density functions at a greatly reduced computational cost. We find that the ROQ rule can be used to speed up inference by factors as high as 300 without introducing systematic bias. This corresponds to a reduction in computational time from around half a year to a half a day, for the longest duration/lowest mass signals. The ROM and ROQ rule are available with the main inference library of the LIGO Scientific Collaboration, LALInference.

preprint2015arXiv

Can we measure individual black-hole spins from gravitational-wave observations?

Measurements of black-hole spins from gravitational-wave observations of black-hole binaries with ground-based detectors are expected to be hampered by partial degeneracies in the gravitational-wave phasing: between the two component spins, and between the spins and the binary's mass ratio, at least for signals that are dominated by the binary's inspiral. Through the merger and ringdown, however, a different set of degeneracies apply. This suggests the possibility that, if the inspiral, merger and ringdown are all within the sensitive frequency band of a detector, we may be able to break these degeneracies and more accurately measure both spins. In this work we investigate our ability to measure individual spins for non-precessing binaries, for a range of configurations and signal strengths, and conclude that in general the spin of the larger black hole will be measurable (at best) with observations from Advanced LIGO and Virgo. This implies that in many applications waveform models parameterized by only one \emph{effective spin} will be sufficient. Our work does not consider precessing binaries or sub-dominant harmonics, although we provide some arguments why we expect that these will not qualitatively change our conclusions.

preprint2015arXiv

Comparison of subdominant gravitational wave harmonics between post-Newtonian and numerical relativity calculations and construction of multi-mode hybrids

Gravitational waveforms which describe the inspiral, merger and ringdown of coalescing binaries are usually constructed by synthesising information from perturbative descriptions, in particular post-Newtonian theory and black-hole perturbation theory, with numerical solutions of the full Einstein equations. In this paper we discuss the "glueing" of numerical and post-Newtonian waveforms to produce hybrid waveforms which include subdominant spherical harmonics ("higher order modes"), and focus in particular on the process of consistently aligning the waveforms, which requires a comparison of both descriptions and a discussion of their imprecisions. We restrict to the non-precessing case, and illustrate the process using numerical waveforms of up to mass ratio $q=18$ produced with the BAM code, and publicly available waveforms from the SXS catalogue. The results also suggest new ways of analysing finite radius errors in numerical simulations.

preprint2015arXiv

Frequency domain reduced order model of aligned-spin effective-one-body waveforms with generic mass-ratios and spins

I provide a frequency domain reduced order model (ROM) for the aligned-spin effective-one-body (EOB) model "SEOBNRv2" for data analysis with second and third generation ground based gravitational wave (GW) detectors. SEOBNRv2 models the dominant mode of the GWs emitted by the coalescence of black hole (BH) binaries. The large physical parameter space (dimensionless spins $-1 \leq χ_i \leq 0.99$ and symmetric mass-ratios $0.01 \leq η\leq 0.25$) requires sophisticated reduced order modeling techniques, including patching in the parameter space and in frequency. I find that the time window over which the inspiral-plunge and the merger-ringdown waveform in SEOBNRv2 are connected is discontinuous when the spin of the deformed Kerr BH $χ=0.8$ or the symmetric mass-ratio $η\sim 0.083$. This discontinuity increases resolution requirements for the ROM. The ROM can be used for compact binary systems with total masses of $2 M_\odot$ or higher for the advanced LIGO (aLIGO) design sensitivity and a $10$ Hz lower cutoff frequency. The ROM has a worst mismatch against SEOBNRv2 of $\sim 1\%$, but in general mismatches are better than $\sim 0.1\%$. The ROM is crucial for key data analysis applications for compact binaries, such as GW searches and parameter estimation carried out within the LIGO Scientific Collaboration (LSC).

preprint2015arXiv

Frequency-domain gravitational waves from non-precessing black-hole binaries. I. New numerical waveforms and anatomy of the signal

In this paper we discuss the anatomy of frequency-domain gravitational-wave signals from non-precessing black-hole coalescences with the goal of constructing accurate phenomenological waveform models. We first present new numerical-relativity simulations for mass ratios up to 18 including spins. From a comparison of different post-Newtonian approximants with numerical-relativity data we select the uncalibrated SEOBNRv2 model as the most appropriate for the purpose of constructing hybrid post-Newtonian/numerical-relativity waveforms, and we discuss how we prepare time-domain and frequency-domain hybrid data sets. We then use our data together with results in the literature to calibrate simple explicit expressions for the final spin and radiated energy. Equipped with our prediction for the final state we then develop a simple and accurate merger-ringdown-model based on modified Lorentzians in the gravitational wave amplitude and phase, and we discuss a simple method to represent the low frequency signal augmenting the TaylorF2 post-Newtonian approximant with terms corresponding to higher orders in the post-Newtonian expansion. We finally discuss different options for modelling the small intermediate frequency regime between inspiral and merger-ringdown. A complete phenomenological model based on the present work is presented in a companion paper.

preprint2015arXiv

Frequency-domain gravitational waves from non-precessing black-hole binaries. II. A phenomenological model for the advanced detector era

We present a new frequency-domain phenomenological model of the gravitational-wave signal from the inspiral, merger and ringdown of non-precessing (aligned-spin) black-hole binaries. The model is calibrated to 19 hybrid effective-one-body--numerical-relativity waveforms up to mass ratios of 1:18 and black-hole spins of $|a/m| \sim 0.85$ ($0.98$ for equal-mass systems). The inspiral part of the model consists of an extension of frequency-domain post-Newtonian expressions, using higher-order terms fit to the hybrids. The merger-ringdown is based on a phenomenological ansatz that has been significantly improved over previous models. The model exhibits mismatches of typically less than 1\% against all 19 calibration hybrids, and an additional 29 verification hybrids, which provide strong evidence that, over the calibration region, the model is sufficiently accurate for all relevant gravitational-wave astronomy applications with the Advanced LIGO and Virgo detectors. Beyond the calibration region the model produces physically reasonable results, although we recommend caution in assuming that \emph{any} merger-ringdown waveform model is accurate outside its calibration region. As an example, we note that an alternative non-precessing model, SEOBNRv2 (calibrated up to spins of only 0.5 for unequal-mass systems), exhibits mismatch errors of up to 10\% for high spins outside its calibration region. We conclude that waveform models would benefit most from a larger number of numerical-relativity simulations of high-aligned-spin unequal-mass binaries.

preprint2015arXiv

Impact of gravitational radiation higher order modes on single aligned-spin gravitational wave searches for binary black holes

Current template-based gravitational wave searches for compact binary coalescences (CBC) use waveform models that neglect the higher order modes content of the gravitational radiation emitted, considering only the quadrupolar $(\ell,|m|)=(2,2)$ modes. We study the effect of such a neglection for the case of aligned-spin CBC searches for equal-spin (and non-spinning) binary black holes in the context of two versions of Advanced LIGO: the upcoming 2015 version, known as early Advanced LIGO (eaLIGO) and its Zero-Detuned High Energy Power version, that we will refer to as Advanced LIGO (AdvLIGO). In addition, we study the case of a non-spinning search for initial LIGO (iLIGO). We do this via computing the effectualness of the aligned-spin SEOBNRv1 ROM waveform family, which only considers quadrupolar modes, towards hybrid post-Newtonian/Numerical Relativity waveforms which contain higher order modes. We find that for all LIGO versions, losses of more than $10\%$ of events occur for mass ratio $q\geq6$ and $M \geq 100M_\odot$ due to the neglection of higher modes. Moreover, for iLIGO and eaLIGO, losses notably increase up to $(39,23)\%$ respectively for the highest mass $(220M_\odot)$ and mass ratio ($q=8$) studied. For the case of early AdvLIGO, losses of $10\%$ occur for $M>50M_\odot$ and $q\geq6$. Neglection of higher modes leads to observation-averaged systematic parameter biases towards lower spin, total mass and chirp mass. For completeness, we perform a preliminar, non-exhaustive comparison of systematic biases to statistical errors. We find that, for a given SNR, systematic biases dominate over statistical errors at much lower total mass for eaLIGO than for AdvLIGO.

preprint2015arXiv

Measuring intermediate mass black hole binaries with advanced gravitational wave detectors

We perform a systematic study to explore the accuracy with which the parameters of intermediate-mass black-hole binary systems can be measured from their gravitational wave (GW) signatures using second-generation GW detectors. We make use of the most recent reduced-order models containing inspiral, merger and ringdown signals of aligned-spin effective-one-body waveforms (SEOBNR) to significantly speed up the calculations. We explore the phenomenology of the measurement accuracies for binaries with total masses between 50 and 500 $M_\odot$ and mass ratios between 0.1 and 1. We find that (i) at total masses below ~200 $M_\odot$, where the signal-to-noise-ratio is dominated by the inspiral portion of the signal, the chirp mass parameter can be accurately measured; (ii) at higher masses, the information content is dominated by the ringdown, and total mass is measured more accurately; (iii) the mass of the lower-mass companion is poorly estimated, especially at high total mass and more extreme mass ratios; (iv) spin cannot be accurately measured for our injection set with non-spinning components. Most importantly, we find that for binaries with non-spinning components at all values of the mass ratio in the considered range and at network signal-to-noise ratio of 15, analyzed with spin-aligned templates, the presence of an intermediate-mass black hole with mass >100 $M_\odot$ can be confirmed with 95% confidence in any binary that includes a component with a mass of 130 $M_\odot$ or greater.

preprint2014arXiv

A simple model of complete precessing black-hole-binary gravitational waveforms

The construction of a model of the gravitational-wave (GW) signal from generic configurations of spinning-black-hole binaries, through inspiral, merger and ringdown, is one of the most pressing theoretical problems in the build-up to the era of GW astronomy. We present the first such model in the frequency domain, "PhenomP", which captures the basic phenomenology of the seven-dimensional parameter space of binary configurations with only three key physical parameters. Two of these (the binary's mass ratio and an effective total spin parallel to the orbital angular momentum, which determines the inspiral rate) define an underlying non-precessing-binary model. The non-precessing-binary waveforms are then "twisted up" with approximate expressions for the precessional motion, which require only one additional physical parameter, an effective precession spin, $χ_p$. All other parameters (total mass, sky location, orientation and polarisation, and initial phase) can be specified trivially. The model is constructed in the frequency domain, which will be essential for efficient GW searches and source measurements. We have tested the model's fidelity for GW applications by comparison against hybrid post-Newtonian-numerical-relativity waveforms at a variety of configurations --although we did not use these numerical simulations in the construction of the model. Our model can be used to develop GW searches, to study the implications for astrophysical measurements, and as a simple conceptual framework to form the basis of generic-binary waveform modelling in the advanced-detector era.

preprint2014arXiv

Frequency domain reduced order models for gravitational waves from aligned-spin compact binaries

Black-hole binary coalescences are one of the most promising sources for the first detection of gravitational waves. Fast and accurate theoretical models of the gravitational radiation emitted from these coalescences are highly important for the detection and extraction of physical parameters. Spinning effective-one-body (EOB) models for binaries with aligned spins have been shown to be highly faithful, but are slow to generate and thus have not yet been used for parameter estimation studies. I provide a frequency-domain singular value decomposition (SVD)-based surrogate reduced order model that is thousands of times faster for typical system masses and has a faithfulness mismatch of better than $\sim 0.1\%$ with the original SEOBNRv1 model for advanced LIGO detectors. This model enables parameter estimation studies up to signal-to-noise ratios (SNRs) of 20 and even up to SNR 50 for total masses below $50 M_\odot$. This article discusses various choices for approximations and interpolation over the parameter space that can be made for reduced order models of spinning compact binaries, provides a detailed discussion of errors arising in the construction and assesses the fidelity of such models.

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.

preprint2013arXiv

Testing the validity of the single-spin approximation in inspiral-merger-ringdown waveforms

Gravitational-wave signals from black-hole binaries with non-precessing spins are described by four parameters -- each black hole's mass and spin. It has been shown that the dominant spin effects can be modeled by a \emph{single} spin parameter, leading to the development of several \emph{three-parameter} waveform models. Previous studies indicate that these models should be adequate for gravitational-wave detection. In this paper we focus on the systematic biases that would result from using them to estimate binary parameters, and consider a one-parameter family of configurations at one choice of mass ratio and effective single spin. We find that for low-mass binaries within that family of configurations, where the observable waveform is dominated by the inspiral, the systematic bias in all physical parameters is smaller than the parameter uncertainty due to degeneracies between the mass ratio and the spins, at least up to signal-to-noise ratios (SNRs) of 50. For higher-mass binaries, where the merger and ringdown make a greater contribution to the observed signal, the bias in the mass ratio is comparable to its uncertainty at SNRs of only $\sim$30, and the bias in the measurement of the total spin is \emph{larger} than the uncertainty defined by the 90% confidence region even at an SNR of only 10. Although this bias may be mitigated in future models by a better choice of single-effective-spin parameter, these results suggest that it may be possible to accurately measure \emph{both} black-hole spins in intermediate-mass binaries.

preprint2012arXiv

An efficient iterative method to reduce eccentricity in numerical-relativity simulations of compact binary inspiral

We present a new iterative method to reduce eccentricity in black-hole-binary simulations. Given a good first estimate of low-eccentricity starting momenta, we evolve puncture initial data for ~4 orbits and construct improved initial parameters by comparing the inspiral with post-Newtonian calculations. Our method is the first to be applied directly to the gravitational-wave (GW) signal, rather than the orbital motion. The GW signal is in general less contaminated by gauge effects, which, in moving-puncture simulations, limit orbital-motion-based measurements of the eccentricity to an uncertainty of $Δe \sim 0.002$, making it difficult to reduce the eccentricity below this value. Our new method can reach eccentricities below $10^{-3}$ in one or two iteration steps; we find that this is well below the requirements for GW astronomy in the advanced detector era. Our method can be readily adapted to any compact-binary simulation with GW emission, including black-hole-binary simulations that use alternative approaches, and neutron-star-binary simulations. We also comment on the differences in eccentricity estimates based on the strain $h$, and the Newman-Penrose scalar $Ψ_4$.

preprint2012arXiv

The NINJA-2 catalog of hybrid post-Newtonian/numerical-relativity waveforms for non-precessing black-hole binaries

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search and parameter-estimation algorithms using numerically generated waveforms, and to foster closer collaboration between the numerical relativity and data analysis communities. The first NINJA project used only a small number of injections of short numerical-relativity waveforms, which limited its ability to draw quantitative conclusions. The goal of the NINJA-2 project is to overcome these limitations with long post-Newtonian - numerical relativity hybrid waveforms, large numbers of injections, and the use of real detector data. We report on the submission requirements for the NINJA-2 project and the construction of the waveform catalog. Eight numerical relativity groups have contributed 63 hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. We summarize the techniques used by each group in constructing their submissions. We also report on the procedures used to validate these submissions, including examination in the time and frequency domains and comparisons of waveforms from different groups against each other. These procedures have so far considered only the $(\ell,m)=(2,2)$ mode. Based on these studies we judge that the hybrid waveforms are suitable for NINJA-2 studies. We note some of the plans for these investigations.

preprint2005arXiv

News from Critical Collapse: Bondi Mass, Tails and Quasinormal Modes

We discuss critical gravitational collapse on the threshold of apparent horizon formation as a model both for the discussion of global aspects of critical collapse and for numerical studies in a compactified context. For our matter model we choose a self-gravitating massless scalar field in spherical symmetry, which has been studied extensively in the critical collapse literature. Our evolution system is based on Bondi coordinates, the mass function is used as an evolution variable to ensure regularity at null infinity. We compute radiation quantities like the Bondi mass and news function and find that they reflect the DSS behavior. Surprisingly, the period of radiation at null infinity is related to the formal result for the leading quasi-normal mode of a black hole with rapidly decreasing mass. Furthermore, our investigations shed some light on global versus local issues in critical collapse, and the validity and usefulness of the concept of null infinity when predicting detector signals.

preprint2000arXiv

Type II Critical Collapse of a Self-Gravitating Nonlinear $σ$-Model

We report on the existence and phenomenology of type II critical collapse within the one-parameter family of SU(2) $σ$-models coupled to gravity. Numerical investigations in spherical symmetry show discretely self-similar (DSS) behavior at the threshold of black hole formation for values of the dimensionless coupling constant $\ccbeta$ ranging from 0.2 to 100; at 0.18 we see small deviations from DSS. While the echoing period $Δ$ of the critical solution rises sharply towards the lower limit of this range, the characteristic mass scaling has a critical exponent $γ$ which is almost independent of $\ccbeta$, asymptoting to $0.1185 \pm 0.0005$ at large $\ccbeta$. We also find critical scaling of the scalar curvature for near-critical initial data. Our numerical results are based on an outgoing-null-cone formulation of the Einstein-matter equations, specialized to spherical symmetry. Our numerically computed initial-data critical parameters $p^*$ show 2nd order convergence with the grid resolution, and after compensating for this variation in $p^*$, our individual evolutions are uniformly 2nd order convergent even very close to criticality.