Parameter Estimation with Targeted Eccentric Numerical-Relativity Simulations for GW200208_22 and GW190620

We have analyzed LVK gravitational wave events that show some evidence of eccentricity from TEOBResumS modeling parameter estimations and have confronted them independently with full numerical generated waveforms from our bank of nearly 2000 simulations of binary black holes. We have used RIFT for Bayesian parameter estimation and found that GW200208_22 KDE estimates favor eccentricities e₂₀= 0.191^+0.123_-0.174 upon entering the LVK band at 20 Hz within a 90% confidence limit. We employed 39 new targeted NR simulations and found a top improved likelihood lnL matching waveform, compared to model-based analysis, with an eccentricity at 20Hz, e₂₀ = 0.200, thus reinforcing the eccentric hypothesis of the binary. We have also used our full bank of numerical waveforms on GW190620 finding that it favors eccentricities at 10 Hz in the range of 0 ≤ e₁₀ ≤ 0.3. New specifically targeted simulations will be required to narrow this eccentricity range.