Search

Charge collection efficiency (CE) η 0 under small signal conditions, corresponding to a uniform field in the detector medium, has been widely used in evaluating the performance of photoconductive detectors. The present paper answers the question, ‘What is the variance of the collected charge in an integrating detector as a function of photoinjection level and what are the errors if we continue to use the small signal equations?’ The variance σ 0 2 in η 0 under small signals has been theoretically derived in the literature and has been a key factor in the detective quantum efficiency modeling of integrating detectors based on various semiconductors. σ 0 2 is a noise source and can degrade the detector performance under incomplete charge collection. The statistical variance σ 0 2 in the CE η 0, under small signals and the variance σ r 2 in the CE ηr under an arbitrary injection level r (injected charge divided by charge on the electrodes) have been studied using the Monte Carlo simulation model developed in this work to evaluate the difference between σ r 2 and σ 0 2 from small to large signals. Initial injection of electron and hole pairs and their subsequent transport and trapping in the presence of an electric field, which is calculated from the Poisson equation, is used to calculate the photocurrent. Each injected carrier is tracked as it moves in the semiconductor until it is either trapped or reaches the collection electrode. Trapped carriers do not contribute to the photocurrent but continue to contribute to the field through the Poisson equation. The instantaneous photocurrent i ph(t) is calculated from the drift of the free carriers through the Shockley–Ramo theorem. i ph(t) is integrated over the duration of the photocurrent to calculate the total collected charge and hence the CE ηr . The variance σ r 2 in ηr is found from multiple simulations of ηr . The ηr and σ r 2 have been generated over varying charge injection ratios r, the electron and hole ranges μτ, mean photoinjection depths δ and drift mobility ratios b. At full injection, the deviation Δ σ r 2 of the CE variance σ r 2 from the uniform field case σ 0 2 ( i . e . Δ σ r 2 = σ r 2 − σ 0 2 ) may be as much as 40% larger or 20% lower than the small signal model prediction. This study provides the extent of errors involved in the variance of the CE in non-uniform fields and quantifies the increase in errors that can occur in high injection cases. In practice, typical injection ratios are less than 0.2, which means that the magnitude of percentage error Δ σ r 2 / σ 0 2 is less than 5%.

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., 24 pages, 18 figures. Abstract abridged

We present a detailed investigation into the properties of GW170729, the gravitational wave with the most massive and distant source confirmed to date. We employ an extensive set of waveform models, including new improved models that incorporate the effect of higher-order waveform modes which are particularly important for massive systems. We find no indication of spin-precession, but the inclusion of higher-order modes in the models results in an improved estimate for the mass ratio of $(0.3-0.8)$ at the 90\% credible level. Our updated measurement excludes equal masses at that level. We also find that models with higher-order modes lead to the data being more consistent with a smaller effective spin, with the probability that the effective spin is greater than zero being reduced from $99\%$ to $94\%$. The 90\% credible interval for the effective spin parameter is now $(-0.01-0.50)$. Additionally, the recovered signal-to-noise ratio increases by $\sim0.3$ units compared to analyses without higher-order modes. We study the effect of common spin priors on the derived spin and mass measurements, and observe small shifts in the spins, while the masses remain unaffected. We argue that our conclusions are robust against systematic errors in the waveform models. We also compare the above waveform-based analysis which employs compact-binary waveform models to a more flexible wavelet- and chirplet-based analysis. We find consistency between the two, with overlaps of $\sim 0.9$, typical of what is expected from simulations of signals similar to GW170729, confirming that the data are well-described by the existing waveform models. Finally, we study the possibility that the primary component of GW170729 was the remnant of a past merger of two black holes and find this scenario to be indistinguishable from the standard formation scenario., Comment: 14 pages, 8 figures, final published version, samples available at https://git.ligo.org/katerina.chatziioannou/gw170729hm_datarelease

We present a prescription for computing gravitational waveforms for the inspiral, merger and ringdown of non-spinning eccentric binary black hole systems. The inspiral waveform is computed using the post-Newtonian expansion and the merger waveform is computed by interpolating a small number of quasi-circular NR waveforms. The use of circular merger waveforms is possible because eccentric binaries circularize in the last few cycles before the merger, which we demonstrate up to mass ratio $q = m_1/m_2 = 3$. The complete model is calibrated to 23 numerical relativity (NR) simulations starting ~20 cycles before the merger with eccentricities $e_\text{ref} \le 0.08$ and mass ratios $q \le 3$, where $e_\text{ref}$ is the eccentricity ~7 cycles before the merger. The NR waveforms are long enough that they start above 30 Hz (10 Hz) for BBH systems with total mass $M \ge 80 M_\odot$ ($230 M_\odot$). We find that, for the sensitivity of advanced LIGO at the time of its first observing run, the eccentric model has a faithfulness with NR of over 97% for systems with total mass $M \ge 85 M_\odot$ across the parameter space ($e_\text{ref} \le 0.08, q \le 3$). For systems with total mass $M \ge 70 M_\odot$, the faithfulness is over 97% for $e_\text{ref} \lesssim 0.05$ and $q \le 3$. The NR waveforms and the Mathematica code for the model are publicly available.

We improve the accuracy of the effective-one-body (EOB) waveforms that were employed during the first observing run of Advanced LIGO for binaries of spinning, nonprecessing black holes by calibrating them to a set of 141 numerical-relativity (NR) waveforms. The NR simulations expand the domain of calibration towards larger mass ratios and spins, as compared to the previous EOBNR model. Merger-ringdown waveforms computed in black-hole perturbation theory for Kerr spins close to extremal provide additional inputs to the calibration. For the inspiral-plunge phase, we use a Markov-chain Monte Carlo algorithm to efficiently explore the calibration space. For the merger-ringdown phase, we fit the NR signals with phenomenological formulae. After extrapolation of the calibrated model to arbitrary mass ratios and spins, the (dominant-mode) EOBNR waveforms have faithfulness --- at design Advanced-LIGO sensitivity --- above $99\%$ against all the NR waveforms, including 16 additional waveforms used for validation, when maximizing only on initial phase and time. This implies a negligible loss in event rate due to modeling for these binary configurations. We find that future NR simulations at mass ratios $\gtrsim 4$ and double spin $\gtrsim 0.8$ will be crucial to resolve discrepancies between different ways of extrapolating waveform models. We also find that some of the NR simulations that already exist in such region of parameter space are too short to constrain the low-frequency portion of the models. Finally, we build a reduced-order version of the EOBNR model to speed up waveform generation by orders of magnitude, thus enabling intensive data-analysis applications during the upcoming observation runs of Advanced LIGO., 27 pages, 15 figures

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, 51 pages, 10 figures; published version

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., Presented at Amaldi 9

We describe the Einstein Toolkit, a community-driven, freely accessible computational infrastructure intended for use in numerical relativity, relativistic astrophysics, and other applications. The Toolkit, developed by a collaboration involving researchers from multiple institutions around the world, combines a core set of components needed to simulate astrophysical objects such as black holes, compact objects, and collapsing stars, as well as a full suite of analysis tools. The Einstein Toolkit is currently based on the Cactus Framework for high-performance computing and the Carpet adaptive mesh refinement driver. It implements spacetime evolution via the BSSN evolution system and general-relativistic hydrodynamics in a finite-volume discretization. The toolkit is under continuous development and contains many new code components that have been publicly released for the first time and are described in this article. We discuss the motivation behind the release of the toolkit, the philosophy underlying its development, and the goals of the project. A summary of the implemented numerical techniques is included, as are results of numerical test covering a variety of sample astrophysical problems., 62 pages, 20 figures

The peeling theorem of general relativity predicts that the Weyl curvature scalars Psi_n (n=0...4), when constructed from a suitable null tetrad in an asymptotically flat spacetime, fall off asymptotically as r^(n-5) along outgoing radial null geodesics. This leads to the interpretation of Psi_4 as outgoing gravitational radiation at large distances from the source. We have performed numerical simulations in full general relativity of a binary black hole inspiral and merger, and have computed the Weyl scalars in the standard tetrad used in numerical relativity. In contrast with previous results, we observe that all the Weyl scalars fall off according to the predictions of the theorem., 7 pages, 3 figures, published version

We quantify the consistency of numerical-relativity black-hole-binary waveforms for use in gravitational-wave (GW) searches with current and planned ground-based detectors. We compare previously published results for the $(\ell=2,| m | =2)$ mode of the gravitational waves from an equal-mass nonspinning binary, calculated by five numerical codes. We focus on the 1000M (about six orbits, or 12 GW cycles) before the peak of the GW amplitude and the subsequent ringdown. We find that the phase and amplitude agree within each code's uncertainty estimates. The mismatch between the $(\ell=2,| m| =2)$ modes is better than $10^{-3}$ for binary masses above $60 M_{\odot}$ with respect to the Enhanced LIGO detector noise curve, and for masses above $180 M_{\odot}$ with respect to Advanced LIGO, Virgo and Advanced Virgo. Between the waveforms with the best agreement, the mismatch is below $2 \times 10^{-4}$. We find that the waveforms would be indistinguishable in all ground-based detectors (and for the masses we consider) if detected with a signal-to-noise ratio of less than $\approx14$, or less than $\approx25$ in the best cases., Comment: 17 pages, 9 figures. Version accepted by PRD

The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise., Comment: Proceedings of Numerical Relativity and Data Analysis NRDA 2008 (Syracuse NY, August 2008) - 14 pages, 1 figure

We discuss results that have been obtained from the implementation of the initial round of testbeds for numerical relativity which was proposed in the first paper of the Apples with Apples Alliance. We present benchmark results for various codes which provide templates for analyzing the testbeds and to draw conclusions about various features of the codes. This allows us to sharpen the initial test specifications, design a new test and add theoretical insight., Corrected version

Tremendous progress has been made towards the solution of the binary-black-hole problem in numerical relativity. The waveforms produced by numerical relativity will play a role in gravitational wave detection as either test beds for analytic template banks or as template banks themselves. As the parameter space explored by numerical relativity expands, the importance of quantifying the effect that each parameter has on first the detection of gravitational waves and then the parameter estimation of their sources increases. In light of this, we present a study of equal-mass, spinning binary-black-hole evolutions through matched filtering techniques commonly used in data analysis. We study how the match between two numerical waveforms varies with numerical resolution, initial angular momentum of the black holes, and the inclination angle between the source and the detector. This study is limited by the fact that the spinning black-hole binaries are oriented axially and the waveforms only contain approximately two and a half orbits before merger. We find that for detection purposes, spinning black holes require the inclusion of the higher harmonics in addition to the dominant mode, a condition that becomes more important as the black-hole spins increase. In addition, we conduct a preliminary investigation of how well a template of fixed spin and inclination angle can detect target templates of arbitrary but nonprecessing spin and inclination for the axial case considered here.

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

We propose the creation of an on-line journal that integrates with a component framework. By means of the framework, simulations referenced in the paper can be run (or re-run) on journal supercomputers at will, allowing verification of the results or deeper analysis of the paper's problem space. The on-line journal would provide standardized job launch and control mechanisms, data formats, and address code portability issues. Through the use of the framework and its components authors can make use of built-in productivity and readability enhancing features (such as automatic parallelism), leverage code published in previous works, and gain stature as their published modules are used by others. We feel that Cactus is an ideal candidate for the on-line journal, therefore we propose the name “The Prickly Pear Archive” for the on-line journal, taking the name from a species of cactus that can be used to make paper.