Numerical Relativity Waveforms Research Articles

Hybrid waveforms for generic precessing binaries for gravitational-wave data analysis

We construct hybrid binary black holes merger waveforms using analytical model waveforms for the early inspiral phase and numerical relativity waveforms for late inspiral to merger and post merger phases. To hybridize analytical and numerical waveforms, we first perform a 3-dimensional rotation to align the instantaneous orbital planes associated with the two waveforms at some fiducial frequency, we then find appropriate phase and time translations that maximize the overlap of the two waveforms in a hybridization interval. We discuss the accuracy and limitations for hybrids constructed by this procedure in the context of LIGO-Virgo-KAGRA observations. Our goal is to hybridize waveforms for more generic precessing binaries and construct longer waveforms that are sufficiently accurate for the parameter estimation techniques for upcoming LIGO observations.

Black hole spectroscopy for KAGRA future prospect in O5

Ringdown gravitational waves of compact binary mergers are an important target to test general relativity. The main components of the ringdown waveform after merger are black hole quasinormal modes. In general relativity, all multipolar quasinormal modes of a black hole should give the same values of black hole parameters. Although the observed binary black hole events so far are not significant enough to perform the test with ringdown gravitational waves, it is expected that the test will be achieved in third generation detectors. The Japanese gravitational wave detector KAGRA, called bKAGRA for the current configuration, has started observation, and discussions for the future upgrade plans have also started. In this study, we consider which KAGRA upgrade plan is the best to detect the subdominant quasinormal modes of black holes in the aim of testing general relativity. We use a numerical relativity waveform as injected signals that contains two multipolar modes and analyze each mode by matched filtering. Our results suggest that the plan FDSQZ, which improves the sensitivity of KAGRA for broad frequency range, is the most suitable configuration for black hole spectroscopy.

Precessing numerical relativity waveform surrogate model for binary black holes: A Gaussian process regression approach

Gravitational wave astrophysics relies heavily on the use of matched filtering both to detect signals in noisy data from detectors, and to perform parameter estimation on those signals. Matched filtering relies upon prior knowledge of the signals expected to be produced by a range of astrophysical systems, such as binary black holes. These waveform signals can be computed using numerical relativity techniques, where the Einstein field equations are solved numerically, and the signal is extracted from the simulation. Numerical relativity simulations are, however, computationally expensive, leading to the need for a surrogate model which can predict waveform signals in regions of the physical parameter space which have not been probed directly by simulation. We present a method for producing such a surrogate using Gaussian process regression which is trained directly on waveforms generated by numerical relativity. This model returns not just a single interpolated value for the waveform at a new point, but a full posterior probability distribution on the predicted value. This model is therefore an ideal component in a Bayesian analysis framework, through which the uncertainty in the interpolation can be taken into account when performing parameter estimation of signals.

Validating the effective-one-body numerical-relativity waveform models for spin-aligned binary black holes along eccentric orbits

Effective-one-body (EOB) numerical-relativity (NR) waveform models for spin-aligned binary black holes (BBHs), known as the SEOBNR waveform models, are based on the EOB theoretical framework and NR simulations. SEOBNR models have played an important role in the LIGO scientific collaboration (LSC) gravitational wave (GW) data analysis for both signal search and parameter estimation. SEOBNR models for quasi-circular orbits have evolved through version 1 to version 4 by extending their validity domain and including more NR results. Along another direction, we recently extended SEOBNRv1 model to SEOBNRE model which is valid for spin-aligned BBH coalescence along eccentric orbits. In this paper we validate this theoretical waveform model by comparing them against the numerical relativity simulation bank, Simulating eXtreme Spacetimes (SXS) catalog. In total, 278 NR waveforms are investigated which include binaries with large eccentricity; large spin and large mass ratio. Our SEOBNRE can model the NR waveforms quite well. The fitting factor for most of the 278 waveforms is larger than 99\%. It indicates that the SEOBNRE model could be used as template waveforms for eccentric spin-aligned BBH coalescence.

Search for intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network

Gravitational wave astronomy has been firmly established with the detection of gravitational waves from the merger of ten stellar mass binary black holes and a neutron star binary. This paper reports on the all-sky search for gravitational waves from intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network. The search uses three independent algorithms: two based on matched filtering of the data with waveform templates of gravitational wave signals from compact binaries, and a third, model-independent algorithm that employs no signal model for the incoming signal. No intermediate mass black hole binary event was detected in this search. Consequently, we place upper limits on the merger rate density for a family of intermediate mass black hole binaries. In particular, we choose sources with total masses $M=m_1+m_2\in[120,800]$M$_\odot$ and mass ratios $q = m_2/m_1 \in[0.1,1.0]$. For the first time, this calculation is done using numerical relativity waveforms (which include higher modes) as models of the real emitted signal. We place a most stringent upper limit of $0.20$~Gpc$^{-3}$yr$^{-1}$ (in co-moving units at the 90% confidence level) for equal-mass binaries with individual masses $m_{1,2}=100$M$_\odot$ and dimensionless spins $\chi_{1,2}= 0.8$ aligned with the orbital angular momentum of the binary. This improves by a factor of $\sim 5$ that reported after Advanced LIGO's first observing run.

Physics of eccentric binary black hole mergers: A numerical relativity perspective

Gravitational wave observations of eccentric binary black hole mergers will provide unequivocal evidence for the formation of these systems through dynamical assembly in dense stellar environments. The study of these astrophysically motivated sources is timely in view of electromagnetic observations, consistent with the existence of stellar mass black holes in the globular cluster M22 and in the Galactic center, and the proven detection capabilities of ground-based gravitational wave detectors. In order to get insights into the physics of these objects in the dynamical, strong-field gravity regime, we present a catalog of 89 numerical relativity waveforms that describe binary systems of non-spinning black holes with mass-ratios $1\leq q \leq 10$, and initial eccentricities as high as $e_0=0.18$ fifteen cycles before merger. We use this catalog to quantify the loss of energy and angular momentum through gravitational radiation, and the astrophysical properties of the black hole remnant, including its final mass and spin, and recoil velocity. We discuss the implications of these results for gravitational wave source modeling, and the design of algorithms to search for and identify eccentric binary black hole mergers in realistic detection scenarios.

Fusing numerical relativity and deep learning to detect higher-order multipole waveforms from eccentric binary black hole mergers

We determine the mass-ratio, eccentricity and binary inclination angles that maximize the contribution of the higher-order waveform multipoles $(\ell, \, |m|)= \{(2,\,2),\, (2,\,1),\, (3,\,3),\, (3,\,2), \, (3,\,1),\, (4,\,4),\, (4,\,3),\, (4,\,2),\,(4,\,1)\}$ for the gravitational wave detection of eccentric binary black hole mergers. We carry out this study using numerical relativity waveforms that describe non-spinning black hole binaries with mass-ratios $1\leq q \leq 10$, and orbital eccentricities as high as $e_0=0.18$ fifteen cycles before merger. For stellar-mass, asymmetric mass-ratio, binary black hole mergers, and assuming LIGO's Zero Detuned High Power configuration, we find that in regions of parameter space where black hole mergers modeled with $\ell=|m|=2$ waveforms have vanishing signal-to-noise ratios, the inclusion of $(\ell, \, |m|)$ modes enables the observation of these sources with signal-to-noise ratios that range between 30\% to 45\% the signal-to-noise ratio of optimally oriented binary black hole mergers modeled with $\ell=|m|=2$ numerical relativity waveforms. Having determined the parameter space where $(\ell, \, |m|)$ modes are important for gravitational wave detection, we construct waveform signals that describe these astrophysically motivate scenarios, and demonstrate that these topologically complex signals can be detected and characterized in real LIGO noise with deep learning algorithms.

Characterization of numerical relativity waveforms of eccentric binary black hole mergers

We introduce a method to quantify the initial eccentricity, gravitational wave frequency, and mean anomaly of numerical relativity simulations that describe non-spinning black holes on moderately eccentric orbits. We demonstrate that this method provides a robust characterization of eccentric binary black hole mergers with mass-ratios $q\leq10$ and eccentricities $e_0\leq0.2$ fifteen cycles before merger. We quantify the circularization rate of a variety of eccentric numerical relativity waveforms introduced in [1] by computing overlaps with their quasi-circular counterparts, finding that $50M$ before merger they attain overlaps ${\cal{O}}\geq0.99$, furnishing evidence for the circularization of moderately eccentric binary black hole mergers with mass-ratios $q\leq10$. We also quantify the importance of including higher-order waveform modes for the characterization of eccentric binary black hole mergers. Using two types of numerical waveforms, one that includes $(\ell, \, |m|)= \{(2,\,2),\, (2,\,1),\, (3,\,3),\, (3,\,2), \, (3,\,1),\, (4,\,4),\, (4,\,3),\, (4,\,2),\,(4,\,1)\}$ and one that only includes the $\ell=|m|=2$ mode, we find that the overlap between these two classes of waveforms is as low as ${\cal{O}}=0.89$ for $q=10$ eccentric binary black hole mergers, underscoring the need to include higher-order waveform modes for the description of these gravitational wave sources. We discuss the implications of these findings for future source modeling and gravitational wave detection efforts.

Surrogate model of hybridized numerical relativity binary black hole waveforms

Numerical relativity (NR) simulations provide the most accurate binary black hole gravitational waveforms, but are prohibitively expensive for applications such as parameter estimation. Surrogate models of NR waveforms have been shown to be both fast and accurate. However, NR-based surrogate models are limited by the training waveforms' length, which is typically about 20 orbits before merger. We remedy this by hybridizing the NR waveforms using both post-Newtonian and effective one body waveforms for the early inspiral. We present NRHybSur3dq8, a surrogate model for hybridized nonprecessing numerical relativity waveforms, that is valid for the entire LIGO band (starting at $20~\text{Hz}$) for stellar mass binaries with total masses as low as $2.25\,M_{\odot}$. We include the $\ell \leq 4$ and $(5,5)$ spin-weighted spherical harmonic modes but not the $(4,1)$ or $(4,0)$ modes. This model has been trained against hybridized waveforms based on 104 NR waveforms with mass ratios $q\leq8$, and $|\chi_{1z}|,|\chi_{2z}| \leq 0.8$, where $\chi_{1z}$ ($\chi_{2z}$) is the spin of the heavier (lighter) BH in the direction of orbital angular momentum. The surrogate reproduces the hybrid waveforms accurately, with mismatches $\lesssim 3\times10^{-4}$ over the mass range $2.25M_{\odot} \leq M \leq 300 M_{\odot}$. At high masses ($M\gtrsim40M_{\odot}$), where the merger and ringdown are more prominent, we show roughly two orders of magnitude improvement over existing waveform models. We also show that the surrogate works well even when extrapolated outside its training parameter space range, including at spins as large as 0.998. Finally, we show that this model accurately reproduces the spheroidal-spherical mode mixing present in the NR ringdown signal.

Gravitational waveforms for high spin and high mass-ratio binary black holes: A synergistic use of numerical-relativity codes

Observation and characterisation of gravitational waves from binary black holes requires accurate knowledge of the expected waveforms. The late inspiral and merger phase of the waveform is obtained through direct numerical integration of the full 3-dimensional Einstein equations. The Spectral Einstein Code (SpEC) utilizes a multi-domain pseudo-spectral method tightly adapted to the geometry of the black holes; it is computationally efficient and accurate, but--for high mass-ratios and large spins--sometimes requires manual fine-tuning for the merger-phase of binaries. The Einstein Toolkit (ET) employs finite difference methods and the moving puncture technique; it is less computationally efficient, but highly robust. For some mergers with high mass ratio and large spins, the efficient numerical algorithms used in SpEC have failed, whereas the simpler algorithms used in the ET were successful. Given the urgent need of testing the accuracy of waveform models currently used in LIGO and Virgo inference analyses for high mass ratios and spins, we present here a synergistic approach to numerical-relativity: We combine SpEC and ET waveforms into complete inspiral-merger-ringdown waveforms, taking advantage of the computational efficiency of the pseudo-spectral code during the inspiral, and the robustness of the finite-difference code at the merger. We validate our method against a case where complete waveforms from both codes are available, compute three new hybrid numerical-relativity waveforms, and compare them with analytical waveform models currently used in LIGO and Virgo science. All the waveforms and the hybridization code are publicly available.

Systematic effects from black hole-neutron star waveform model uncertainties on the neutron star equation of state

We identify various contributors of systematic effects in the measurement of the neutron star (NS) tidal deformability and quantify their magnitude for several types of neutron star - black hole (NSBH) binaries. Gravitational waves from NSBH mergers contain information about the components' masses and spins as well as the NS equation of state. Extracting this information requires comparison of the signal in noisy detector data with theoretical templates derived from some combination of post-Newtonian (PN) approximants, effective one-body (EOB) models and %analytic fits to numerical relativity (NR) simulations. The accuracy of these templates is limited by errors in the NR simulations, by the approximate nature of the PN/EOB waveforms, and by the hybridization procedure used to combine them. In this paper, we estimate the impact of these errors by constructing and comparing a set of PN-NR hybrid waveforms, for the first time with NR waveforms from two different codes, namely, SpEC and SACRA, for such systems. We then attempt to recover the parameters of the binary using two non-precessing template approximants. We find that systematic errors are too large for tidal effects to be accurately characterized for any realistic NS equation of state model. We conclude that NSBH waveform models must be significantly improved if they are to be useful for the extraction of NS equation of state information or even for distinguishing NSBH systems from binary black holes.

Matter imprints in waveform models for neutron star binaries: Tidal and self-spin effects

The combined observation of gravitational and electromagnetic waves from the coalescence of two neutron stars marks the beginning of multi-messenger astronomy with gravitational waves (GWs). The development of accurate gravitational waveform models is a crucial prerequisite to extract information about the properties of the binary system that generated a detected GW signal. In binary neutron star systems (BNS), tidal effects also need to be incorporated in the modeling for an accurate waveform representation. Building on previous work [Phys.Rev.D96 121501], we explore the performance of inspiral-merger waveform models that are obtained by adding a numerical relativity (NR) based approximant for the tidal part of the phasing (NRTidal) to existing models for nonprecessing and precessing binary black hole systems (SEOBNRv4, PhenomD and PhenomPv2), as implemented in the LSC Algorithm Library Suite. The resulting BNS waveforms are compared and contrasted to target waveforms hybridizing NR waveforms, covering the last approx. 10 orbits up to merger and extending through the postmerger phase, with inspiral waveforms calculated from 30Hz obtained with TEOBResumS. The latter is a state-of-the-art effective-one-body waveform model that blends together tidal and spin effects. We probe that the combination of the PN-based self-spin terms and of the NRTidal description is necessary to obtain minimal mismatches (< 0.01) and phase differences (< 1 rad) with respect to the target waveforms. However, we also discuss possible improvements and drawbacks of the NRTidal approximant in its current form, since we find that it tends to overestimate the tidal interaction with respect to the TEOBResumS model during the inspiral.

Time-domain effective-one-body gravitational waveforms for coalescing compact binaries with nonprecessing spins, tides, and self-spin effects

We present TEOBResumS, a new effective-one-body (EOB) waveform model for nonprecessing (spin-aligned) and tidally interacting compact binaries.Spin-orbit and spin-spin effects are blended together by making use of the concept of centrifugal EOB radius. The point-mass sector through merger and ringdown is informed by numerical relativity (NR) simulations of binary black holes (BBH) computed with the SpEC and BAM codes. An improved, NR-based phenomenological description of the postmerger waveform is developed.The tidal sector of TEOBResumS describes the dynamics of neutron star binaries up to merger and incorporates a resummed attractive potential motivated by recent advances in the post-Newtonian and gravitational self-force description of relativistic tidal interactions. Equation-of-state dependent self-spin interactions (monopole-quadrupole effects) are incorporated in the model using leading-order post-Newtonian results in a new expression of the centrifugal radius. TEOBResumS is compared to 135 SpEC and 19 BAM BBH waveforms. The maximum unfaithfulness to SpEC data $\bar{F}$ -- at design Advanced-LIGO sensitivity and evaluated with total mass $M$ varying between $10M_\odot \leq M \leq 200 M_\odot$ --is always below $2.5 \times 10^{-3}$ except for a single outlier that grazes the $7.1 \times 10^{-3}$ level. When compared to BAM data, $\bar{F}$ is smaller than $0.01$ except for a single outlier in one of the corners of the NR-covered parameter space, that reaches the $0.052$ level.TEOBResumS is also compatible, up to merger, to high end NR waveforms from binary neutron stars with spin effects and reduced initial eccentricity computed with the BAM and THC codes. The model is designed to generate accurate templates for the analysis of LIGO-Virgo data through merger and ringdown. We demonstrate its use by analyzing the publicly available data for GW150914.

Tracking Black Hole Kicks from Gravitational-Wave Observations.

Coalescing binary black holes emit anisotropic gravitational radiation. This causes a net emission of linear momentum that produces a gradual acceleration of the source. As a result, the final remnant black hole acquires a characteristic velocity known as recoil velocity or gravitational kick. The symmetries of gravitational wave emission are reflected in the interactions of the gravitational wave modes emitted by the binary. In this Letter, we make use of the rich information encoded in the higher-order modes of the gravitational wave emission to infer the component of the kick along the line of sight (or radial kick). We do this by performing parameter inference on simulated signals given by numerical relativity waveforms for nonspinning binaries using numerical relativity templates of aligned-spin (nonprecessing) binary black holes. We find that for suitable sources, namely those with mass ratio q≥2 and total mass M∼100 M_{⊙}, and for modest radial kicks of 120 km/s, the 90% credible intervals of our posterior probability distributions can exclude a zero kick at a signal-to-noise ratio of 15, using a single Advanced LIGO detector working at its early sensitivity. The measurement of a nonzero radial kick component would provide the first observational signature of net transport of linear momentum by gravitational waves away from their source.

Enriching the symphony of gravitational waves from binary black holes by tuning higher harmonics

For the first time, we construct an inspiral-merger-ringdown waveform model within the effective-one-body formalism for spinning, nonprecessing binary black holes that includes gravitational modes beyond the dominant $(\ell,|m|) = (2,2)$ mode, specifically $(\ell,|m|)=(2,1),(3,3),(4,4),(5,5)$. Our multipolar waveform model incorporates recent (resummed) post-Newtonian results for the inspiral and information from 157 numerical-relativity simulations, and 13 waveforms from black-hole perturbation theory for the (plunge-)merger and ringdown. We quantify the improved accuracy including higher-order modes by computing the faithfulness of the waveform model against the numerical-relativity waveforms used to construct the model. We define the faithfulness as the match maximized over time, phase of arrival, gravitational-wave polarization and sky position of the waveform model, and averaged over binary orientation, gravitational-wave polarization and sky position of the numerical-relativity waveform. When the waveform model contains only the $(2,2)$ mode, we find that the averaged faithfulness to numerical-relativity waveforms containing all modes with $\ell \leq$ 5 ranges from $90\%$ to $99.9\%$ for binaries with total mass $20-200 M_\odot$ (using the Advanced LIGO's design noise curve). By contrast, when the $(2,1),(3,3),(4,4),(5,5)$ modes are also included in the model, the faithfulness improves to $99\%$ for all but four configurations in the numerical-relativity catalog, for which the faithfulness is greater than $98.5\%$. Using our results, we also develop also a (stand-alone) waveform model for the merger-ringdown signal, calibrated to numerical-relativity waveforms, which can be used to measure multiple quasi-normal modes. The multipolar waveform model can be extended to include spin-precession, and will be employed in upcoming observing runs of Advanced LIGO and Virgo.

Detection and characterization of spin-orbit resonances in the advanced gravitational wave detectors era

In this paper, we test the performance of templates in detection and characterization of Spin-orbit resonant (SOR) binaries. We use precessing SEOBNRv3 waveforms as well as {\it four} numerical relativity (NR) waveforms to model GWs from SOR binaries and filter them through IMRPhenomD, SEOBNRv4 (non-precessing) and IMRPhenomPv2 (precessing) approximants. We find that IMRPhenomD and SEOBNRv4 recover only $\sim70\%$ of injections with fitting factor (FF) higher than 0.97 (or 90\% of injections with ${\rm FF} >0.9$).However, using the sky-maxed statistic, IMRPhenomPv2 performs magnificently better than their non-precessing counterparts with recovering $99\%$ of the injections with FFs higher than 0.97. Interestingly, injections with $\Delta \phi = 180^{\circ}$ have higher FFs ($\Delta \phi$ is the angle between the components of the black hole spins in the plane orthogonal to the orbital angular momentum) as compared to their $\Delta \phi =0^{\circ}$ and generic counterparts. This implies that we will have a slight observation bias towards $\Delta \phi=180^{\circ}$ SORs while using non-precessing templates for searches. All template approximants are able to recover most of the injected NR waveforms with FFs $>0.95$. For all the injections including NR, the error in estimating chirp mass remains below $<10\%$ with minimum error for $\Delta \phi = 180^{\circ}$ resonant binaries. The symmetric mass ratio can be estimated with errors below $15\%$. The effective spin parameter $\chi_{\rm eff}$ is measured with maximum absolute error of 0.13. The in-plane spin parameter $\chi_p$ is mostly underestimated indicating that a precessing signal will be recovered as a relatively less precessing signal. Based on our findings, we conclude that we not only need improvements in waveform models towards precession and non-quadrupole modes but also better search strategies for precessing GW signals.

Inferring black-hole orbital dynamics from numerical-relativity gravitational waveforms

Binary-black-hole dynamics cannot be related to the resulting gravitational-wave signal by a constant retarded time. This is due to the nontrivial dynamical spacetime curvature between the source and the signal. In a numerical-relativity simulation there is also some ambiguity in the black-hole dynamics, which depend on the gauge (coordinate) choices used in the numerical solution of Einstein's equations. It has been shown previously that a good approximation to the direction of the binary's time-dependent orbital angular momentum $\stackrel{^}{\mathbf{L}}(t)$ can be calculated from the gravitational-wave signal. This is done by calculating the direction that maximizes the quadrupolar ($\ensuremath{\ell}=2$, $|m|=2$) emission. The direction depends on whether we use the Weyl scalar ${\ensuremath{\psi}}_{4}$ or the gravitational-wave strain $h$, but these directions are nonetheless invariant for a given binary configuration. We treat the ${\ensuremath{\psi}}_{4}$-based direction as a proxy to $\stackrel{^}{\mathbf{L}}(t)$. We investigate how well the binary's orbital phase, ${\ensuremath{\phi}}_{\mathrm{orb}}(t)$, can also be estimated from the signal. For this purpose we define a quantity $\mathrm{\ensuremath{\Phi}}(t)$ that agrees well with ${\ensuremath{\phi}}_{\mathrm{orb}}(t)$. One application is to studies that involve injections of numerical-relativity waveforms into gravitational-wave detector data.

Total-variation methods for gravitational-wave denoising: Performance tests on Advanced LIGO data

We assess total-variation methods to denoise gravitational-wave signals in real noise conditions, by injecting numerical-relativity waveforms from core-collapse supernovae and binary black hole mergers in data from the first observing run of Advanced LIGO. This work is an extension of our previous investigation where only Gaussian noise was used. Since the quality of the results depends on the regularization parameter of the model, we perform an heuristic search for the value that produces the best results. We discuss various approaches for the selection of this parameter, either based on the optimal, mean, or multiple values, and compare the results of the denoising upon these choices. Moreover, we also present a machine-learning-informed approach to obtain the Lagrange multiplier of the method through an automatic search. Our results provide further evidence that total-variation methods can be useful in the field of Gravitational-Wave Astronomy as a tool to remove noise.

Measuring the properties of nearly extremal black holes with gravitational waves

Characterizing the properties of black holes is one of the most important science objectives for gravitational-wave observations. Astrophysical evidence suggests that black holes that are nearly extremal (i.e. spins near the theoretical upper limit) might exist and thus might be among the merging black holes observed with gravitational waves. In this paper, we explore how well current gravitational wave parameter estimation methods can measure the spins of rapidly spinning black holes in binaries. We simulate gravitational-wave signals using numerical-relativity waveforms for nearly-extremal, merging black holes. For simplicity, we confine our attention to binaries with spins parallel or antiparallel with the orbital angular momentum. We find that recovering the holes' nearly extremal spins is challenging. When the spins are nearly extremal and parallel to each other, the resulting parameter estimates do recover spins that are large, though the recovered spin magnitudes are still significantly smaller than the true spin magnitudes. When the spins are nearly extremal and antiparallel to each other, the resulting parameter estimates recover the small effective spin but incorrectly estimate the individual spins as nearly zero. We study the effect of spin priors and argue that a commonly used prior (uniform in spin magnitude and direction) hinders unbiased recovery of large black-hole spins.

Eccentric binary black hole inspiral-merger-ringdown gravitational waveform model from numerical relativity and post-Newtonian theory

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.