  Nearly all of the previous gravitational wave (GW) searches in the LIGO-Virgo
data included GW waveforms with only the dominant quadrupole mode, i.e.,
omitting higher-order harmonics which are predicted by general relativity.
Based on the techniques developed in Wadekar et al. [1,2], we improve the IAS
pipeline by ($i$) introducing higher harmonics in the GW templates, ($ii$)
downweighting noise transients ('glitches') to improve the search sensitivity
to high-mass and high-redshift binary black hole (BBH) mergers. We find 14 new
BBH mergers with $0.53\leq p_{\rm astro}\leq 0.88$ on running our pipeline over
the public LIGO-Virgo data from the O3 run (we use the detection threshold as
$p_{\rm astro}>0.5$ following the approach of other pipelines). We also broadly
recover the high-significance events from earlier catalogs, except some which
were either vetoed or fell below our SNR threshold for trigger collection.
  A few notable properties of our new candidate events are as follows. At
$>95$\% credibility, 4 candidates have total masses in the IMBH range (i.e.,
above 100 $M_\odot$), and 9 candidates have $z>0.5$. 9 candidates have median
mass of the primary BH falling roughly within the pair instability mass gap,
with the highest primary mass being $300_{+60}^{-120} M_\odot$. 5 candidates
have median mass ratio $q < 0.5$. Under a prior uniform in effective spin
$\chi_{\rm eff}$, 6 candidates have $\chi_{\rm eff} > 0$ at $>95\%$
credibility. We also find that including higher harmonics in our search raises
the significance of a few previously reported marginal events (e.g.,
GW190711_030756). While our new candidate events have modest false alarm rates
($\gtrsim 1.6 $/yr), a population inference study including these can better
inform the parameter space of BHs corresponding to the pair instability mass
gap, high redshifts, positive effective spins and asymmetric mass ratios.
