\section{Introduction}

Solving the two-body problem in General Relativity (GR) remains a challenge of both theoretical interest and astrophysical relevance.
Albeit an analytical solution is lacking, advances in numerical relativity
(NR) in the past decades provided the first numerical evolutions of
merging compact objects \cite{Pretorius:2005gq,Campanelli:2005dd,Baker:2005vv}, as
well as catalogs of waveforms \cite{Mroue:2013xna,Jani:2016wkt,Healy:2017psd,Dietrich:2018phi, Boyle:2019kee}. On the
analytical side of the problem, approximations to the binary motion and gravitational
radiation, via expansions in one or more small parameters, have been applied to different
domains of validity~\cite{Blanchet:2013haa,Barack:2018yvs,Tiec:2014lba}, providing us
with a variety of waveform models.

The effective-one-body (EOB) framework is a synergistic approach that allows us to resum information from several analytical approximations.
NR-calibrated inspiral-merger-ringdown models based on EOB theory \cite{Pan:2011gk,Taracchini:2012ig,Taracchini:2013rva,Bohe:2016gbl,Cotesta:2018fcv} were employed by LIGO-Virgo experiments to detect gravitational waves (GWs) and infer astrophysical and cosmological information from them~\cite{Abbott:2016blz,Abbott:2016nmj,Abbott:2017vtc,Abbott:2017oio,Abbott:2017gyy,TheLIGOScientific:2017qsa,LIGOScientific:2018mvr,TheLIGOScientific:2016wfe,LIGOScientific:2018mvr,TheLIGOScientific:2016src}. In view of the expected increase in the signal-to-noise ratio of signals detected with upcoming LIGO-Virgo runs, and next generation detectors in space (LISA~\cite{Audley:2017drz}) and on Earth (Einstein Telescope~\cite{Punturo:2010zz} and Cosmic Explorer~\cite{Evans:2016mbw}),
it is important and timely to include more physics and build more accurate waveforms in the EOB approach.

Historically, physical information from the two-body problem has
mostly entered EOB theory via the post-Newtonian (PN) expansion
\cite{Buonanno:1998gg, Buonanno:2000ef, Damour:2000we}, valid for
bound orbits at large distances and for velocities smaller than the
speed of light $v^2/c^2\sim\text{G}M/r c^2 \ll 1$ (here $M=m_1+m_2$
is the total mass, with $m_1$ the mass of the primary and $m_2$ the
mass of the secondary body).  PN conservative-dynamics information has
so far been calculated up to fourth order, in the nonspinning
  case, using the Arnowitt-Deser-Misner (ADM)
\cite{Jaranowski:2015lha,Damour:2014jta,Damour:2016abl}, Fokker
\cite{Bernard:2015njp,Bernard:2016wrg,Bernard:2017bvn} and
effective-field-theory approaches \cite{Foffa:2019rdf,Foffa:2019yfl} (which were also employed to determine the 5PN gravitational interaction in the static limit~\cite{Foffa:2019hrb,Blumlein:2019zku}). In the quasi-circular-orbit limit,
4PN information has been successfully included in the EOB dynamics in
the form of an expansion in the inverse radius $u \equiv \text{G}M/r
c^2 \ll 1$ and in the momenta $\boldsymbol{p}^2$, with exact
dependence on the symmetric mass ratio $\nu=m_{1}m_{2}/M^{2}$
\cite{Damour:2015isa}.
Further resummations of this PN expansion form the core of the EOB waveform models \cite{Damour:2008gu,Damour:2009kr,Barausse:2009xi,Pan:2011gk,Hinderer:2016eia}.
Post-Minkowskian (PM) information, valid in the weak field $\text{G}M/r c^2 \ll 1$, but for all velocities $v^2/c^2\leq 1$, has also provided valuable insight in the structure of EOB Hamiltonians, for both spinning and non-spinning bound systems \cite{Damour:2017zjx, Damour:2016pm,Vines:2017pm,Vines:2018gqi,Antonelli:2019ytb}.

The self-force (SF) program, initiated in Refs.~\cite{Mino:1996nk,Quinn:1996am}  and based on an expansion of Einstein's equations in the small mass ratio (SMR) $q=m_2/m_1$,
has been successful in the calculation of the gravitational SF of a small body around Schwarzschild \cite{Barack:2007tm,Barack:2010tm},  and recently Kerr black-holes \cite{Shah:2012gu,vandeMeent:2015lxa,vandeMeent:2016pee,vandeMeent:2017bcc}, to first order in the mass ratio and for generic bound orbits. The results, corroborated by the use of several gauges and numerical techniques (see, e.g., Ref.~\cite{Barack:2018yvs} and references therein), have been already used to evolve extreme-mass-ratio-inspirals (EMRIs) around a Schwarzschild black-hole \cite{Warburton:2011fk,Osburn:2015duj} and they represent a key input for EMRI waveform modeling schemes recently developed \cite{vandeMeent:2018rms} and under development \cite{Hinderer:2008dm}.

As the SF program employs different gauge-dependent schemes to obtain
its results \cite{Barack:2018yvs}, it is paramount to be able to check
results via gauge-invariant quantities, such as the innermost-stable--circular-orbit (ISCO)-shift
\cite{Barack:2009ey}, periastron advance
\cite{Barack:2010ny,LeTiec:2011bk,vandeMeent:2016hel}, spin-precession
\cite{Dolan:2013roa,Kavanagh:2017wot,Bini:2014ica,Bini:2018ylh},
tidal invariants \cite{Dolan:2014pja,Bini:2014zxa} and the Detweiler redshift
\cite{Detweiler:2008ft,Barack:2011ed,Akcay:2015pza,Kavanagh:2015lva,Shah:2013uya,Johnson-McDaniel:2015vva}.
For a particle with four-velocity $\tilde{u}^\alpha$ normalized in an
effective metric $\tilde{g}_{\alpha\beta}=g^\text{(0)}_{\alpha\beta}+h^{\text{R}}_{\alpha\beta}$
[i.e., moving around a Schwarzschild background
$g^\text{(0)}_{\alpha\beta}$ perturbed by a regularized metric
$h_{\alpha\beta}^{\text{R}}$ and such that
$\tilde{g}_{\alpha\beta}\tilde{u}^\alpha\tilde{u}^\beta=-1+\mathcal{O}(\nu)$],
the Detweiler redshift is defined as the ratio between proper time
measured in an orbit around the effective metric
$\tilde{g}_{\alpha\beta}$, $d\tilde{\tau}$, and coordinate time,
$d\tau$ \footnote{As pointed out in Ref.~\cite{Barack:2018yvs}, $z$
  does not correspond to the gravitational redshift due to the use of
  the regularized perturbation $h^{\text{R}}_{\alpha\beta}$ in its
  definition. It does only in the full geometry, e.g., including a
  singular metric $h^{\text{S}}_{\alpha\beta}$ at the location of the
  particle such that the body perturbation is $h_{\alpha\beta}\equiv
  h^{\text{R}}_{\alpha\beta}+h^{\text{S}}_{\alpha\beta}$. A sounder
  physical description can be obtained if the small companion is a
  black hole, since the Detweiler redshift can then be linked to the
  surface gravity $\kappa$ of the small
  body~\cite{Zimmerman:2016ajr}. }:
$z \equiv (\tilde{u}^t)^{-1}=d\tilde{\tau}/d\tau$
\cite{Detweiler:2008ft,Barack:2018yvs}.  Recently, the
Detweiler redshift has been used for cross-cultural studies between
approximations to the two-body problem in GR
\cite{LeTiec:2011ab,Tiec:2014lba}, and it has provided an important
benchmark to check PN and SMR results in the small-mass-ratio and
weak-field domain, in which both PN and SMR frameworks are expected to
be valid. This synergistic program has been extended
to NR simulations of equal--mass-ratio binaries with the computation of
the Detweiler redshift in Ref.~\cite{Zimmerman:2016ajr}.

As pointed out in Ref.~\cite{Damour:2009sf}, gauge-invariant SMR
quantities such as the Detweiler redshift can be also used to inform
the conservative sector of EOB Hamiltonians
\cite{LeTiec:2011dp,Damour:2009sf,Barausse:2011dq,LeTiec:2011bk}.
There are two ways in which this valuable information
 could be incorporated into the EOB approach: it can be either used to partially determine high-order PN coefficients of EOB Hamiltonians \cite{kavanaghetal:2015,Hopper:2015icj,Kavanagh:2016idg,Bini:2013rfa,Bini:2014nfa,Bini:2015bla, Bini:2015bfb,Bini:2015xua,Bini:2016dvs,Bini:2016qtx} or it can be used to resum PN dynamics around the test-body limit \cite{LeTiec:2011dp,Barausse:2011dq,Akcay:2012ea}. Here, we focus on the latter approach.

Currently available EOB Hamiltonians informed with the Detweiler redshift cannot be reliably evolved near the Schwarzschild light-ring (LR) radius, i.e., $r=3\text{G}M/c^2$. Such an issue, hereafter called the \textit{LR-divergence} problem,  appears as a coordinate singularity of the effective Hamiltonian at the Schwarzschild LR \cite{LeTiec:2011dp,Akcay:2012ea}. In this paper we address the problem and, adopting a different EOB gauge, we obtain a Hamiltonian with SMR information that exhibits no divergence at the LR radius. This result allows us to use the precious near-LR, strong-field information from SF calculations in the evolutions of EOB Hamiltonians.

The organization of the paper is as follows.
In Sec.~\ref{sec:prelim} we review the LR-divergence arising from informing the conservative sector of standard EOB Hamiltonians with the Detweiler redshift and we discuss how a different EOB gauge (introduced in Ref.~\cite{Damour:2017zjx} in the context of PM calculations) helps to solve the issue.
In Sec.~\ref{sec:EOBSMR}, we inform the conservative sector of EOB Hamiltonians in the alternative gauge with circular-orbit information from the Detweiler redshift, and with both non-circular-orbit and higher-order-in-mass-ratio information from the PN approximation.
In Sec.~\ref{sec:evolutions}, we evolve quasi-circular inspirals from this LR-divergence-free Hamiltonian and show that the evolution of the orbital separation crosses the LR radius without encountering singularities. Moreover, we perform systematic comparisons against NR predictions of phase and binding energy for non-spinning systems with mass ratios $1/10 \leq q \leq 1$. We conclude in Sec.~\ref{sec:conclusions}. In Appendix~\ref{sec:appred} we present high-precision fits to the Detweiler redshift with improved data in the strong field.
We use geometric units G=$c$=1 throughout the paper.

\section{On gauges and the light-ring divergence}
\label{sec:prelim}

We begin by noting some conventions to be used in the following sections. In the present paper, we do not consider spinning systems; we denote the reduced mass by $\mu=(m_1 m_2)/M$ and the total mass by $M=m_1+m_2$.  We work with generalized (polar) coordinates ${q}_a\equiv (r,\phi)$ in the orbital plane, with canonically conjugate momenta $p_a\equiv(p_r,p_\phi)$, and we often employ the mass-reduced inverse orbital separation $u \equiv M/r$ and the mass-reduced momenta $\hat{p}_r \equiv p_r/\mu$ and $\hat{p}_\phi\equiv p_\phi/(M\mu)$.

\subsection{The light-ring divergence}

In the EOB approach, the real two-body motion is mapped to the
effective motion of a test body in an effective \textit{deformed}
Schwarzschild spacetime with coordinates $(t,r,\theta,\phi)$, with the deformation parameter being the symmetric mass ratio $\nu$.
The mapping can be obtained via a dictionary between the action integrals $I_a=(2\pi)^{-1}\oint
{p}_a\, d{q}_a$ of a two-body system in the
center-of-mass frame and those of a test-body moving in the effective
metric $g_{\mu\nu}^{\text{eff}}$. Considering orbits in the equatorial plane $\theta=\pi/2$, identifying the radial and angular
action integrals of real and effective systems, i.e., setting
$I_\text{r}^{\text{real}}=I_\text{r}^{\text{eff}}$ and
$I_\phi^{\text{real}}=I_\phi^{\text{eff}}$, the EOB approach allows a
simple relation between the real
$H_{\text{EOB}}(r,p_r,p_\phi,\nu)$ and effective
$H_{\text{eff}}(r,p_r,p_\phi,\nu)$ Hamiltonians
\cite{Buonanno:1998gg}:
\begin{equation}\label{enmap}
H_{\text{EOB}} \equiv M \hat H_{\text{EOB}} =M\sqrt{1+2\nu\bigg(\frac{H_{\text{eff}}}{\mu}-1\bigg)}\,.
\end{equation}
$H_{\text{eff}}$ describes the motion of a test body with mass $\mu$ and is determined by a mass-shell constraint of the form~\cite{Damour:2000we}
\begin{equation}\label{HamJac}
g^{\mu\nu}_\text{eff}p_\mup_\nu+\mu^2+Q(r,p_r,p_\phi,\nu)
=0,
\end{equation}
where the effective metric is given by
\begin{align}\label{effmetric}
ds^{2}=-A(r,\nu)dt^{2}+[A(r,\nu)\bar D(r,\nu)]^{-1}dr^{2}+r^{2} d\Omega^{2}\,,
\end{align}
with the potentials $A(r,\nu)$ and $\bar D(r,\nu)$ depending on the orbital separation $r$ and the symmetric mass ratio $\nu$. In terms of the inverse radius $u=M/r$, they reduce to $A_0(u)=1-2u$ and $\bar D_0=1$ in the test particle limit ($\nu\rightarrow0$). Inserting the inverse of the metric~\eqref{effmetric} into Eq.~\eqref{HamJac}, and using $p_\mu=(-H_\text{eff},p_r,p_\theta=0,p_\phi)$, the mass-reduced effective Hamiltonian $\hat{H}_{\text{eff}}\equiv H_{\text{eff}}/\mu$ is found to be \cite{Damour:2000we}
\begin{align}\label{HeffDJS}
\hat{H}_{\text{eff}}^2&=A(u,\nu)\big[1+\hat{p}_{\phi}^{2}u^{2}+A(u,\nu)\bar D(u,\nu)\hat{p}_{r}^{2}
\nonumber\\
&\quad+\hat{Q}(u,\hat{p}_r,\hat{p}_\phi,\nu)\big],
\end{align}
with $\hat Q\equiv Q/\mu^2$.  The non-geodesic function $ Q$ in
Eq.~\eqref{HamJac} has been introduced to extend the EOB Hamiltonian
through 3PN order without changing the mapping (\ref{enmap}) (for a geodesic
one-body motion at 3PN order with an energy map different from
(\ref{enmap}) see Appendix A in Ref.~\cite{Damour:2000we}). Its
mass-reduced form $\hat{Q}(u,\hat{p}_r,\hat{p}_\phi,\nu)$ in Eq.~\eqref{HeffDJS} generically depends on both the mass-recuded
radial momentum $\hat{p}_{r}$ and the mass-recuded angular momentum $\hat{p}_{\phi}$. Reference \cite{Damour:2000we}
  showed that at 3PN order $\hat{Q}$ must be
fourth order in the momenta, and that the non-geodesic term
is not uniquely fixed. By setting some of the free parameters to zero, it is possible
to make the function  $\hat{Q}(u,\hat{p}_r,\hat{p}_\phi,\nu)$ depend only on the radial momentum
[i.e., $\hat Q(u,\hat{p}_r,\hat{p}_\phi,\nu)\rightarrow\hat
Q(u,\hat{p}_r,\nu)$]. Since 2000, this choice of $\hat{Q}$
has been adopted in several EOB papers [although see Refs.~\cite{Damour:2002vi,Buonanno:2007pf} for alternative choices
of $\hat{Q}$]. Henceforth, we shall
denote the $\hat{Q}$ function that only depends on the radial momentum
as $\hat{Q}^{\rm DJS}({u},\hat{p}_r,\nu)$, after the
initials of the three authors of
Ref.~\cite{Damour:2000we}. We refer to the DJS EOB Hamiltonian as
the Hamiltonian that uses the $\hat{Q}^{\rm DJS}({u},\hat{p}_r,\nu)$ function.
Note that in this gauge, the angular momentum $\hat{p}_{\phi}$ only appears in the second term in brackets in
Eq.~\eqref{HeffDJS}. Moreover, in the circular orbit limit ($\hat{p}_{r}=0$) the
conservative dynamics information is fully described by the $A(u,\nu)$
potential in this gauge, as found at 2PN order~\cite{Buonanno:1998gg}.
The 4PN expressions for $A(u,\nu)$, $\bar D(u,\nu)$
and $\hat Q^{\rm DJS}({u},\hat{p}_{r},\nu)$ in the DJS gauge, for quasi-circular orbits,
are obtained mapping Eq.~\eqref{enmap} to the 4PN-expanded Hamiltonian and can be
found in Ref.~\cite{Damour:2015isa}.

The first efforts to incorporate SMR quantities in EOB Hamiltonians sought to do so
using the gauge of Eq.~\eqref{HeffDJS} with $\hat Q(u,\hat p_{r},\hat p_{\phi},\nu)\rightarrow \hat Q^{\rm DJS}(u,\hat p_{r},\nu)$
\cite{Barausse:2011dq,Akcay:2015pjz,Barack:2010ny,Akcay:2012ea}. In this gauge, the function $A(u, \nu)$, having the complete dynamical information for circular orbits, allows a linear-in-$\nu$ expansion about the Schwarzschild limit:
\begin{align}\label{ADnuexp}
&A(u, \nu)=1-2u+\nu a (u)+\mathcal{O}(\nu^2)\,.
\end{align}
The $a(u)$ function resums the complete circular-orbit PN dynamics in linear order in $\nu$.
References~\cite{LeTiec:2011dp,Barausse:2011dq} obtained an expression for $a(u)$ employing the linear-in-$\nu$ correction to the Detweiler redshift. Notably, the Detweiler redshift is expanded around the Schwarzschild background, $z(x)=\sqrt{1-3x}+\nu \Delta z(x)+\mathcal{O}(\nu^2)$ [where $x\equiv (M\Omega)^{2/3}$ is the gauge-independent inverse radius], and the $\Delta z$ correction is linked to $a(u)$ via the first law of binary black-hole mechanics \cite{LeTiec:2011ab}.
The resulting expression reads:

\begin{equation}\label{Asf}
a(u)= \Delta z(u)\sqrt{1-3u}-u\bigg(1+\frac{1-4u}{\sqrt{1-3u}}\bigg)\,.
\end{equation}
In Eq.~\eqref{Asf}, $\Delta z$ depends on the gauge-dependent inverse radius $u$, rather than its gauge-independent counterpart $x$. This is only correct if we restrict to first order in $\nu$, since $x=u+\mathcal{O}(\nu)$.
The quantity $\Delta z(x)$, has been fitted with data extending to the LR \cite{Akcay:2012ea}, allowing precious strong-field information to enter the EOB dynamics.

The form of $a(u)$ is suggestive of trouble arising at the Schwarzschild light ring, i.e., at $u_{\text{LR}}=1/3$, where the second term in Eq.~\eqref{Asf} diverges. In principle, this divergence might be tamed by the behaviour of the redshift $\Delta z(u)$ appearing in the first term in brackets, but data for the redshift up to the LR show that this is not the case and that $a(u)$ indeed diverges there~\cite{Akcay:2012ea}. This is worrisome, as $a(u)$ directly enters the effective Hamiltonian and, via the energy map, the EOB-resummed dynamics. The EOB dynamics thus contains a divergence \textit{for generic orbits} (e.g., for any value of $\hat{p}_{\phi}$ and $\hat{p}_{r}$).
It was pointed out in Ref.~\cite{Akcay:2012ea} that the LR-divergence is a phase-space coordinate singularity that arises due to the use of the DJS gauge, and that can be solved adopting a different gauge in which the function $\hat Q$ grows as $\hat Q\propto \hat{p}_\phi^3$ when $\hat{p}_\phi\rightarrow\infty$ and $\hat{p}_r\rightarrow0$.

It is worth mentioning that the argument in Ref.~\cite{Akcay:2012ea} stems from a similar LR divergence that has appeared when including tidal effects
in the EOB approach \cite{Bini:2012gu}. Tidal effects enter the potential $A(u)$ via a correction in a tidal expansion akin to Eq.~\eqref{ADnuexp}: $A(u)=A_{\text{2pp}}+\mu_\text{T} a_\text{T}(u,\nu)+\mathcal{O}(\mu_\text{T}^2)$, where $A_{\text{2pp}}$ is the two point-particle (pp) EOB potential \cite{Bini:2012gu} and $\mu_\text{T}$ the small tidal parameter. It has been found in Ref.~\cite{Bini:2012gu} that, in the extreme-mass-ratio limit and for circular orbits, the first-order correction scales as $a_\text{T}(u,\nu)\propto (1-3u)^{-1}$ when $u\rightarrow u_\text{LR}$. An alternative EOB Hamiltonian that includes dynamical tides without introducing poles at the LR has been introduced in Ref.~\cite{Steinhoff:2016rfi}; this has been achieved by abandoning the DJS gauge (see, e.g., their Appendix D).

\subsection{The post-Schwarzschild effective-one-body gauge}

Reference~\cite{Damour:2017zjx} has shown that it is possible to obtain a different EOB gauge, hereafter the \textit{post-Schwarzschild} (PS) gauge, solving Eq.~\eqref{HamJac} with the Schwarzschild limit of the metric~\eqref{effmetric}.
The mass-reduced effective Hamiltonian thus obtained has the following form:
\begin{equation}\label{HeffEG}
\hat{H}_{\text{eff}}^{\text{PS}}=\sqrt{\hat{H}_{\text{S}}^{2}+(1-2u)\hat{Q}^{\text{PS}}(u,\nu,\hat{H}_{\text{S}})}\,,
\end{equation}
where $\hat{H}_{\text{S}}$ is the Schwarzschild Hamiltonian:
\begin{equation}\label{Hs}
\hat{H}_{\text{S}}(u,\hat{p}_{r},\hat{p}_{\phi})=\sqrt{(1-2u)\left [1+\hat{p}_{\phi}^{2}u^{2}+(1-2u)\hat{p}_{r}^{2}\right ]}\,.
\end{equation}

In Ref.~\cite{Damour:2017zjx}, the PS function $\hat{Q}^{\text{PS}}$ has been derived to 2PM order via a scattering-angle calculation and to 3PN order via a canonical transformation from the DJS Hamiltonian at 3PN. In Ref.~\cite{Antonelli:2019ytb}, these calculations have been extended to 3PM and 4PN orders, respectively (the latter only in the near-circular orbit limit).

It is noticed that, in PS EOB Hamiltonians, all the information on the two-body problem with $\nu \neq 0$ is contained in $\hat{Q}^{\text{PS}}(u,\nu,\hat{H}_{\text{S}})$. This feature and the fact that circular-orbit dynamics is contained also in the $\hat{Q}$ function, significantly differentiate PS Hamiltonians from DJS ones. The PS gauge is uniquely fixed resumming the angular and radial momenta into the Schwarzschild Hamiltonian \eqref{Hs}. The powers of such momenta are furthermore not bound in any way, due to the generic functional dependence of $\hat{Q}^{\text{PS}}(u,\nu,\hat{H}_{\text{S}})$ on $\hat{H}_{\text{S}}$. In principle, then, arbitrary powers of $\hat{p}_{\phi}$ are contained in $\hat{Q}^{\text{PS}}(u,\nu,\hat{H}_{\text{S}})$ via $\hat{H}_{\text{S}}$. In particular, differently from $\hat{Q}^{\text{DJS}}(u,\nu,\hat{p}_r)$,  powers of momentum enter at second order
in $\hat{Q}^{\text{PS}}(u,\nu,\hat{H}_{\text{S}})$ instead of fourth order.

The unconstrained dependence of $\hat{Q}^{\text{PS}}$ on $\hat{H}_{\text{S}}$ makes the use of PS Hamiltonians very appealing in the context of our work. It was shown in Ref.~\cite{Damour:2017zjx} that, in the high energy limit for which $\hat{p}_{\phi}\rightarrow\infty$,
the LR-divergence can be captured by the coefficient of a term proportional to $\hat{H}_{\text{S}}^3$. This result is in agreement with a point made in the conclusions of Ref.~\cite{Akcay:2012ea}. As it approaches the LR radius, the effective mass moving in a deformed-Schwarzschild background described by Eqs.~\eqref{ADnuexp} and~\eqref{Asf} has a divergent-energy behaviour that must be removed with an appropriate energy-corrected mass-ratio parameter $\tilde\nu=\nu \hat{H}_{\text{S}}$.
In the next section, building from this knowledge and making use of a simple ansatz for $\hat{Q}^{\text{PS}}(u,\nu,\hat{H}_{\text{S}})$, we construct a Hamiltonian in the PS gauge that contains information from $\Delta z$, while remaining analytic at the LR.

\section{Conservative dynamics of post-Schwarzschild Hamiltonians}
\label{sec:EOBSMR}

\subsection{Information from circular orbits}\label{sec:preamble}
In this section, we link the conservative sector of the PS EOB Hamiltonian to the SMR contribution to $\Delta z$.
Following Ref.~\cite{Barausse:2011dq}, we do so matching, at fixed frequency, the circular orbit binding energy at linear order in $\nu$ from the EOB Hamiltonian with the binding energy in the same limit from SF results. The latter is obtained in Ref.~\cite{LeTiec:2011dp} and is a consequence of the first law of binary-black-hole mechanics. As a function of
$\Delta z$ and the gauge-invariant inverse radius $x$, it reads \cite{LeTiec:2011dp}:

\begin{equation}\label{EsEsf}
\hat{E}^{\text{SF}}_{\text{bind}}=\frac{1-2x}{\sqrt{1-3x}}-1+\nu \hat{E}_{\text{SMR}}(x,\Delta z,\Delta z')+\mathcal{O}(\nu^2)\,,
\end{equation}
\begin{align}\label{Ezold}
\hat{E}_{\text{SMR}}(x,\Delta z,\Delta z')=&-1+\sqrt{1-3 x}-\frac{x}{3}
\Delta z'(x)\nonumber\\
&+\frac{\Delta z(x)}{2}+\frac{(7-24 x) x}{6
	(1-3 x)^{3/2}}\,.
\end{align}
The prime denotes differentiation with respect to $x$.
We find it useful to rewrite the redshift as:

\begin{equation}\label{zsmr}
\Delta z(x)=\frac{\zzz^{(0)}(x)}{1-3x}+\frac{\zzz^{(1)}(x)}{\sqrt{1-3x}}+\frac
\zzz^{(2)}(x)}{1-3x}\ln E_\text{S}^{-2}(x)
\,,
\end{equation}
 In the above expression, we have defined $E_\text{S}(x)\equiv (1-2x)/\sqrt{1-3x}$.
In Appendix~\ref{sec:appred}, $\zzz^{(0)}(x)$, $\zzz^{(1)}(x)$ and $\zzz^{(2)}(x)$ are fitted to high-precision SF data and such to be analytic at the LR.
Equation~\eqref{Ezold} then reads:
\begin{widetext}
\begin{align}\label{Esf}
\hat{E}_{\text{SMR}}=&\sqrt{1-3 x}-1+\frac{(7-24 x) x}{6 (1-3 x)^{3/2}}
+\frac{1}{2(1-3x)}
\bigg[\zzz^{(0)}(x)+\zzz^{(1)}(x)\sqrt{1-3 x}
+\zzz^{(2)}(x)\ln E_\text{S}^{-2}(x)\bigg]
\nonumber\\
&
-\frac{x}{3(1-3x)}\bigg\{\frac{3\zzz^{(0)}(x)}{1-3 x}+\frac{3 \zzz^{(1)}(x)}{2 \sqrt{1-3 x}}+
\bigg[\frac{1-6 x}{(1-2x)(1-3x)}+\frac{3\ln E_\text{S}^{-2}(x)}{(1-3x)}\bigg]\zzz^{(2)}(x)
\nonumber\\
&
+(\zzz^{(0)})'(x)+ \sqrt{1-3 x}\,(\zzz^{(1)})'(x)+ (\zzz^{(2)})'(x)\ln E_\text{S}^{-2}(x)\bigg\}\,.
\end{align}
\end{widetext}
We next consider the PS EOB Hamiltonian $H_{\text{EOB}}$, with an ansatz for $\hat Q^{\rm PS}$ reading:

\begin{equation}\label{HSMR}
\hat{Q}^{\text{PS}}_{\text{SMR}}(u,\nu,\hat{H}_{\text{S}})=\nu\big[f_{0}(u)\hat{H}_{\text{S}}^{5}+f_{1}(u)\hat{H}_{\text{S}}^{2}+f_{2}(u)\hat{H}_{\text{S}}^{3}\ln\hat{H}^{-2}_{\text{S}}\big]\,.
\end{equation}
In the rest of this section, when matching to the SMR results, we limit to circular orbits; thus we use $\hat{H}_{\text{S}}(u,\hat{p}_r=0,\hat{p}_\phi)$ in Eq.~(\ref{HSMR}). The role of the $\hat{H}_{\text{S}}^5$ term is to capture the global divergence $(1-3x)^{-3/2}$ of Eq.~\eqref{Esf}\footnote{In principle, a $\hat{H}_{\text{S}}^3$ term will suffice to capture the divergence. However, we find that this minimal choice leads to evolutions that are not well behaved for systems with comparable masses.}, while the second term $\hat{H}_{\text{S}}^2$ is devised to incorporate the $\sqrt{1-3x}$ terms appearing in the numerator of the same equation, which would make the Hamiltonian imaginary after the light ring. The term proportional to $\ln\hat{H}_{\text{S}}^{-2}$ incorporates the logs in the fit that would make the Hamiltonian non-smooth at the light ring.
Setting $p_r=0$ and using:
\begin{equation}
\dot{p}_r=-\frac{\partial  H_{\text{EOB}}}{dr}(r,p_r=0,p_\phi^{\text{circ}},\nu)=0\,,
\end{equation}
the (mass-reduced) circular-orbit momentum $\hat{p}_\phi^{\text{circ}}$ as a function of the inverse radius $u$ is determined at linear order in $\nu$ (with $f_i'(u)=df_i/du$):
\begin{align}\label{pu}
\hat{p}_{\phi}^{\text{circ}}(u,\nu)=&\frac{1}{\sqrt{u(1-3u)}}+\nu\,
\frac{(1-2 u)^2 }{4 (1-3
	u)^3 \sqrt{u}}\times\nonumber\\
&\bigg[2(1-2u)^3f_{0}(u)+2(1-3u)^{3/2}f_1 (u)
\nonumber\\
&+2(1-2u)(1-3u)f_2(u)\ln E_\text{S}^{-2} (u)\nonumber\\
&-(1-2u)^4f'_{0}(u)-(1-2u)(1-3u)^{3/2}f'_{1}(u)\nonumber\\
&-(1-2u)^2(1-3u)\ln E_\text{S}^{-2}(u)f'_2(u)\bigg]+\mathcal{O}(\nu^2)\,.
\end{align}
We further use the relation:
\begin{equation}\label{omega}
\Omega=\frac{\partial H_{\text{EOB}}}{dp_\phi}(r,p_r=0,p_{\phi}^{\text{circ}},\nu)\,,
\end{equation}
and exploit its link to the gauge-independent inverse radius $x$ given by $x=(M\Omega)^{2/3}$.
Inserting Eq.~\eqref{pu} in Eq.~\eqref{omega} and inverting the obtained expression at linear order in $\nu$,
we establish a link between the gauge-dependent $u$ and the gauge-independent $x$ inverse radii:
\begin{widetext}
\begin{align}\label{ux}
u^{\text{circ}}(x,\nu) = &x+
\frac{x\,\nu}{6 (1-3
	x)^{3/2}}
\bigg\{
4-20x+24 x^2-(4-12x) \sqrt{1-3 x}-10 (1-2 x)^4 f_0(x)\nonumber\\
&-4 \sqrt{1-3 x} \left(1-5x+6 x^2\right) f_1(x)+\big[4 - 28 x + 64 x^2 - 48 x^3 -(6-42x+96x^2-72x^3)\ln E_\text{S}^{-2}(x)\big] f_2(x) \nonumber\\
&+\left(1-10x+40x^2-80x^3+80x^4-32x^5\right) f_0'(x)+\sqrt{1-3 x}\left(1-7x+16x^2-12x^3\right)
f_1'(x)\nonumber\\
&+\left(1-9x+30x^2-44x^3+24x^4\right)\ln E_\text{S}^{-2}(x) f_2'(x)
\bigg\}+\mathcal{O}(\nu^2)\,.
\end{align}
\end{widetext}
To calculate the (mass-reduced) gauge-invariant, circular-orbit binding energy at linear order in $\nu$ from $H_{\text{EOB}}$, we employ the definition :
\begin{equation}\label{Ebind}
\hat{E}^{\text{EOB}}_{\text{bind}} \equiv  (H_{\text{EOB}}-M)/\mu\,,
\end{equation}
Inserting Eqs.~\eqref{pu} and \eqref{ux} in $H_{\text{EOB}}$ and retaining only terms up to first order in the mass ratio, we get:
\begin{widetext}
	\begin{align}\label{Eeob}
\hat{E}^{\text{EOB}}_{\text{bind}}(x,\nu)=&\frac{1-2x}{\sqrt{1-3 x}}-1-\frac{\nu}{6(1-3x)^3}
\bigg\{
6-55x+170 x^2-189x^3+36 x^4-\sqrt{1-3 x} \left(6-46 x+108 x^2-72 x^3\right)	\nonumber\\
&-\left(3-7x-18 x^2\right) (1-2 x)^4 f_0(x)-(1-3x)^{3/2} \left(3-16 x+20 x^2\right)
f_1(x)\nonumber\\
&+(1-3x)(1-2 x)^2  \big[2x (1-6 x)-\left(3-9x-6 x^2\right) \ln E_\text{S}^{-2}(x)\big]f_2(x)+2 x (1-2 x)^5 (1-3 x) f_0'(x) 	\nonumber\\
&+2x (1-3 x)^{5/2} (1-2 x)^2 f_1'(x)+2x (1-3 x)^2 (1-2 x)^3 \ln E_\text{S}^{-2}(x) f_2'(x)
\bigg\}+\mathcal{O}(\nu^2)\,.
	\end{align}
\end{widetext}
Matching Eq.~\eqref{EsEsf} [with correction given by Eq.~\eqref{Esf}] and Eq.~\eqref{Eeob}, we obtain differential equations to be solved for $f_0(x)$, $f_1(x)$ and $f_2(x)$.
Further splitting the $f_i$ coefficients as follows:
\begin{align}
&f_0(x)=\tilde{f}_0(x)+\sum_{i=0}^{i=2}f_0^{(i)}(x)\Delta z^{(i)}(x)\label{X0ans1}
\\
&f_1(x)=\tilde{f}_1(x)+\sum_{i=0}^{i=2}f_1^{(i)}(x)\Delta z^{(i)}(x)\label{X1ans1}\,,\\
&f_2(x)=\tilde{f}_2(x)+\sum_{i=0}^{i=2}f_2^{(i)}(x)\Delta z^{(i)}(x)\label{X2ans1}\,,
\end{align}
and imposing that the Hamiltonian coefficients be analytic at the LR radius (i.e., that they do not contain  $\sqrt{1-3x}$ or $\ln E_\text{S}^{-2}(x)$ terms),
we obtain the following non-zero solutions\footnote{Similarly to what is done in Eq.~\cite{Barausse:2011dq}, we impose that the PN expansion cannot admit half-integer powers of $x$. This allows us to set all constants of integration to zero.}:
\begin{subequations}
\begin{align}
&\tilde{f}_0(x)=-\frac{x (1-3x)\left(1-4x
\right)}{(1-2x)^5}\,,\\
&\tilde{f}_1(x)=-\frac{x
}{(1-2 x)^2}\,,\\
&f_0^{(0)}(x)=\frac{1-3x}{(1-2x)^5}\,,\\
&f_1^{(1)}(x)=\frac{1}{(1-2x)^2}\,,\\
&f_2^{(2)}(x)=\frac{1}{(1-2x)^3}\,.
\end{align}
\end{subequations}
The $f_{i}(x)$ coefficients are readily found via Eqs.~\eqref{X0ans1}, ~\eqref{X1ans1} and ~\eqref{X2ans1} and then inserted in the non-geodesic term in the effective Hamiltonian \eqref{HSMR} to obtain:
\begin{align}\label{QSMR}
\frac{\hat{Q}^{\text{PS}}_{\text{SMR}}}{\nu}(u,\nu,\hat{H}_{\text{S}})=&(1-3u)\bigg[\frac{\zzz^{(0)}(u)}{(1-2u)^5}-\frac{(1-4 u)\text{ }u}{(1-2 u)^5}\bigg]\hat{H}_{\text{S}}^5\nonumber\\
&+\bigg[\frac{\zzz^{(1)}(u)}{(1-2 u)^2}-\frac{u}{(1-2
	u)^2}\bigg]\hat{H}_{\text{S}}^2 \nonumber\\
&+\frac{\zzz^{(2)}(u)}{(1-2u)^3}\hat{H}_{\text{S}}^{3}\ln\hat{H}_{\text{S}}^{-2}\,.
\end{align}
We see that the resulting Hamiltonian concisely resums the complete circular-orbit PN dynamics at linear order in $\nu$. The non-geodesic function $\hat{Q}^{\text{PS}}_{\text{SMR}}$ does not contain any term divergent at the LR, as $\zzz^{(0)}(u)$, $\zzz^{(1)}(u)$ and $\zzz^{(2)}(u)$ are constructed to be analytic there.

\subsection{Information from non-circular orbits and from higher orders in the mass ratio}\label{sec:noncirc}

The calculation in Sec.~\ref{sec:preamble} is carried out in the circular-orbit limit at linear order in the mass ratio. However, it is  possible to include more physical information to the Hamiltonian, coming both from non-circular-orbit terms and from terms at higher orders in the mass ratio. For instance, self-force information for mildly eccentric orbits can be obtained via the SMR correction to the periastron advance $\rho_\text{SF}$ \cite{Barack:2010ny}, which can then be linked to the EOB potentials. This was the strategy used in Refs.~\cite{Damour:2009sf,Barausse:2011dq} to obtain an expression for the potential $\bar D(r)$ in terms of $\Delta z(u)$ and $\rho_\text{SF}(u)$ and introduce non-circular SF data into the EOB Hamilonian up to the Schwarzschild ISCO (i.e., $u_\text{ISCO}=1/6$). Alternatively, one can exploit the generalized redshift~\cite{Barack:2011ed} and link it to $\bar D(r)$, as done in Refs.~\cite{Tiec:2015cxa,Akcay:2015pjz}.
Here, we insert generic-orbit PN information in our Hamiltonian and leave the inclusion of non-circular SMR information in $\hat Q^{\rm PS}$  to future work.

Post-Schwarzschild EOB Hamiltonians with PN information from generic-orbits have been already considered in the literature. For example, the PS Hamiltonians at 3PN order has been investigated in Ref.~\cite{Damour:2017zjx}.
Using the PN parameters $Y \equiv (\hat H^{2}_{\text{S}}-1)\sim \mathcal{O}(1/c^2)$ and $u$, its expression is given by:
\begin{align}
\label{QPScoeff}
\hat Q^{\text{PS}}_{\text{3PN}}=&3\nu u^{2}Y+5\nu u^{3}
\nonumber\\
&+\bigg(3\nu-\frac{9}{4}\nu^{2}\bigg)u^{2}Y^{2}+\bigg(27\nu-\frac{23}{4}\nu^{2}\bigg)u^{3}Y\nonumber\\
&+\bigg(\frac{175}{3}\nu-\frac{41\pi^{2}}{32}\nu-\frac{7}{2}\nu^{2}\bigg)u^{4}\,.
\end{align}
As discussed, the above Hamiltonian contains two-body information that is not captured by the calculation leading to $\hat{Q}^{\text{PS}}_{\text{SMR}}$ and that we wish to add to it.

To this end, we consider a mixed SMR-3PN non-geodesic function of the following form:

\begin{equation}\label{Qtot}
\hat Q^{\text{PS}}_{\text{SMR-3PN}}=\hat Q^{\text{PS}}_{\text{SMR}}+\Delta\hat Q^{\text{PS}}\,,
\end{equation}
where $\hat Q^{\text{PS}}_{\text{SMR}}$ is given by Eq.~\eqref{QSMR} and contains all the circular-orbit terms at linear order in $\nu$, while $\Delta\hat Q^{\text{PS}}$ is fixed demanding that it contains all the additional PN information from Eq.~\eqref{QPScoeff}, in such a way not to contribute to the linear-in-$\nu$ binding energy in the circular-orbit limit.

We opt to further split $\Delta\hat Q^{\text{PS}}$ into two contributions: $\Delta\hat Q^{\text{PS}}_{\text{extra}}$ collects the extra terms up to 3PN order (including both non-circular 3PN terms at linear order in $\nu$ and $\nu^2$ terms), while $\Delta\hat Q^{\text{PS}}_{\text{count.}}$ is a counterterm whose functionality is explained below. We then have:
\begin{equation}\label{Q4pnsplit}
\Delta\hat Q^{\text{PS}} \equiv \Delta\hat Q^{\text{PS}}_{\text{extra}}-\Delta\hat Q^{\text{PS}}_{\text{count.}}\,.
\end{equation}
The former contribution is readily obtained calculating the difference between Eq.~\eqref{QPScoeff} and the 3PN expansion of Eq.~\eqref{QSMR}\footnote{That is, Eq.~\eqref{QSMR} is expanded in the PN parameters $u$ and $Y=\hat H_{\text{S}}^2-1$. The redshift functions $\zzz^{(0)}(u)$, $\zzz^{(1)}(u)$ and $\zzz^{(2)}(u)$ also need to be PN expanded: their expressions are obtained matching the 3PN expansion of the redshift from Ref.~\cite{Kavanagh:2016idg} and Eq.~\eqref{zsmr}.}. The result reads:
\begin{align}\label{Qpn1}
\Delta\hat Q^{\text{PS}}_{\text{extra}}=& 3\nu u^2 Y+\bigg(3\nu-\frac{9}{4} \nu^2 \bigg) u^2Y^2+3 \nu u^3 \nonumber\\
&+\bigg(22\nu-\frac{23}{4} \nu^2 \bigg) u^3Y+\bigg(16 \nu-\frac{7}{2} \nu^2 \bigg)  u^4\,.
\end{align}
In the PS gauge $\hat Q^{\text{PS}}$ depends on momenta via $\hat H_\text{S}(u,\hat{p}_{r},\hat{p}_{\phi})$, which cannot be separated into circular and non-circular orbit contributions.
Because of that, the linear-in-$\nu$ portion of Eq.~\eqref{Qpn1} contributes to the linear-in-$\nu$ binding energy for circular orbits.
Therefore, the addition of $\Delta\hat Q^{\text{PS}}_{\text{extra}}$ to $\Delta\hat Q^{\text{PS}}_{\text{SMR}}$ spoils the matching between EOB and SF binding energies for circular orbits at linear order in the mass ratio guaranteed by the sole presence of $\Delta\hat Q^{\text{PS}}_{\text{SMR}}$.

The matching between the two binding energies can be maintained with a particular choice of the second contribution to Eq.~\eqref{Q4pnsplit}, i.e., $\Delta\hat Q^{\text{PS}}_{\text{count.}}$. We choose a counterterm that starts at 4PN, in order not to spoil the agreement at 3PN for generic orbits guaranteed by Eq.~\eqref{Qpn1}:
\begin{equation}\label{Qpn2}
	\Delta\hat Q^{\text{PS}}_{\text{count.}}=\nu \big[q_{(3,2)}u^3Y^2+q_{(4,1)}u^4Y+q_{(5,0)}u^5\big]\,.
\end{equation}
We impose that the linear-in-$\nu$ binding energy from $\Delta\hat Q^{\text{PS}}$ from Eq.~\eqref{Q4pnsplit} [calculated as done for Eq.\eqref{Eeob} in Sec.~\ref{sec:preamble}] vanishes and we obtain:
\begin{equation}\label{qcoeff}
q_{(3,2)}=9\,;\, q_{(4,1)}=96\,;\, q_{(5,0)}=112\,.
\end{equation}

The final PN correction $\Delta\hat Q^{\text{PS}}$ thus contains all the extra information from generic orbits at 3PN that is not captured by $\hat Q^{\text{PS}}_{\text{SMR}}$, without contributing to the linear in mass ratio binding energy for circular orbits.
The exercise above can be repeated at one PN order higher to obtain $\Delta\hat Q^{\text{PS}}$ at 4PN starting from the 4PN EOB Hamiltonian in the PS gauge.~\cite{Antonelli:2019ytb}. Such a computation does not present major differences from the calculation above: the only feature changing is the counterterm, which needs to start at 5PN and include logarithmic terms. We have decided not to include $\Delta\hat Q^{\text{PS}}$ at 4PN in this paper, as the 4PN Hamiltonian from which it is constructed is only valid for near-circular orbits. The $\Delta\hat Q^{\text{PS}}$ at 3PN that we obtain here is instead valid for generic orbits.

\section{Inspirals in effective-one-body theory}
\label{sec:evolutions}

\subsection{Plunging through the light ring with small mass-ratio Hamiltonians}

\begin{figure}[t]
	\includegraphics[width=\columnwidth]{figuresSMR/Revolq10.pdf}
	\caption{\textbf{Plunges through the light-ring radius:} the evolved orbital separation for the SMR Hamiltonians is presented. The effective masses of models $H_{\text{SMR}}^{\rm EOB,PS}$ and $H_{\text{SMR-3PN}}^{\rm EOB,PS}$
	plunge through the LR radius $r_S=3M$. Conversely, the plunge of the effective mass of $H_{\text{SMR}}^{\rm EOB}$ presents unphysical features associated to the LR-divergence.
	}
	\label{fig:orbitsep}
\end{figure}

In this section, we evolve the EOB Hamiltonians constructed in Secs.~\ref{sec:preamble} and ~\ref{sec:noncirc} [i.e., Eq.~\eqref{HeffEG} with non-geodesic functions~\eqref{QSMR} and~\eqref{Qtot}], and the EOB Hamiltonian with SMR information in the DJS gauge. We refer to them as $H_{\text{SMR}}^{\rm EOB,PS}$, $H_{\text{SMR-3PN}}^{\rm EOB,PS}$ and $H_{\text{SMR}}^{\rm EOB}$, see Table~\ref{table:models}.

The EOB approach comprises of a conservative sector, discussed in detail in Sec.~\ref{sec:prelim}, and a dissipative sector, responsible for the slow GW-driven inspiral of the compact bodies towards merger. The basic set of equations for inspiraling orbits in the EOB framework are the Hamilton equations augmented with a radiation-reaction force $\mathcal{F}_{\text{RR}}$.
In terms of a generic mass-reduced EOB Hamiltonian $\hat{H}^{\rm EOB}(\hat r,\hat{p}_{r_{*}},\hat{p}_\phi)$, the equations read \cite{Buonanno:2000ef,Buonanno:2007pf,Damour:2007xr,Pan:2011gk}:
\begin{subequations}
\begin{align}
& \frac{d\hat r}{d\hat t}=\frac{A(\hat r)}{\sqrt{D(\hat r)}}\frac{\partial \hat{H}^{\rm EOB} }{\partial \hat{p}_{r_{*}}}\,, \label{HH1}\\
&\frac{d\phi}{d\hat t}=\frac{\partial \hat{H}^{\rm EOB}}{\partial \hat{p}_{\phi}}\,, \label{HH2}\\
&\frac{d\hat{p}_{ r_{*}}}{d\hat t}=-\frac{A(\hat r)}{\sqrt{D(\hat r)}}\frac{\partial \hat{H}^{\rm EOB}}{\partial \hat r}+\mathcal{F}_{\text{RR}}\frac{\hat{p}_{r_{*}}}{\hat{p}_{\phi}}\,, \label{HH3}\\
&\frac{d\hat{p}_{\phi}}{d\hat t}=\mathcal{F}_{\text{RR}}\,, \label{HH4}
\end{align}
\end{subequations}
where we have introduced the mass-reduced radius $\hat r \equiv r/M$ and coordinate time $\hat t \equiv t/M$ and used the mass-reduced radial momentum $\hat{p}_{r_{*}}$ conjugate to the radius $r_{*}$ in tortoise coordinates,
defined for generic potentials $A(\hat r)$ and $D(\hat r)$\footnote{Here $D(\hat r)$ is the inverse of $\bar D(\hat r)$ mentioned in Sec.~\ref{sec:prelim}.} by:
\begin{equation}\label{prprstar}
\frac{d\hat r_{*}}{d\hat r} \equiv \frac{\sqrt{D(\hat r)}}{A(\hat r)}
=\frac{\hat{p}_{r}}{\hat{p}_{r_{*}}}
\,.
\end{equation}
In the evolution of the EOB Hamiltonian in the DJS gauge we use the PN-expanded expressions for $A(\hat r)$, $D(\hat r)$ and $\hat Q^{\text{DJS}}$ at the required PN order~\cite{Buonanno:1998gg,Damour:2000we,Damour:2015isa}, whereas we use their test-body limits in the evolutions of Hamiltonians in the PS gauge\footnote{The effective Hamiltonian in the PS gauge \eqref{HeffEG} is obtained solving the Hamilton-Jacobi equations with the Schwarzschild metric. The $A(\hat r)$ and $D(\hat r)$ are therefore fixed by their Schwarzschild limits.}. The Hamiltonians in both gauges depend on $\hat{p}_{ r_{*}}$, rather than $\hat{p}_{ r}$.

The radiation reaction force $\mathcal{F}_{\text{RR}}$ drives the inspiral of the system and it contains semi-analytical two-body information \cite{Damour:2007xr,Damour:2008gu,Pan:2010hz}. In this paper, we employ its non-Keplerian form (with $\hat \Omega \equiv d\phi/d\hat t=M\Omega$):
\begin{equation}
^{\text{nK}}\mathcal{F}_{\text{RR}}=-\frac{1}{\nu\hat \Omega}\frac{dE}{dt}\,,
\end{equation}
where $dE/dt$ is the GW flux for quasi-circular orbits~\cite{Damour:2008gu}:
\begin{equation}
\frac{dE}{dt}=\frac{\hat\Omega^2}{8\pi}\sum_{l=2}^{l_\text{max}=8}\sum_{m=l-2}^{l} m^2 \big|\hat rh_{lm}\big|^2\,.
\end{equation}
The modes $h_{lm}$ are built from PN theory, but resummed multiplicatively (see e.g., Ref.~\cite{Damour:2008gu}). Here, we use the resummation of the (non-spinning) modes and flux presented in Ref.~\cite{Pan:2011gk} (which coincides with the state-of-the-art modes and flux used in the EOB waveform model for
LIGO/Virgo data-analsyis~\cite{Bohe:2016gbl}, when spins are set to zero). We do not include
the ``next-to-quasi-circular'' (NQC) coefficients~\cite{Bohe:2016gbl}, or any calibration parameter obtained
imposing better agreement with numerical-relativity waveforms. Our main motivation here is to compare how well the conservative
EOB-dynamics of SMR models compare to PN ones and with NR.

The result of the evolved orbital separations $\hat r$ of both DJS and PS  Hamiltonians  for $q=1/10$ are reported in Fig.~\ref{fig:orbitsep}.
Focusing on the evolution in the DJS case, it is seen that the pole in the conservative part of the DJS Hamiltonian affects the motion of the effective body close to the LR radius. That is, $H_{\text{SMR}}^{\rm EOB}$ diverges at $\hat r_{\text{S}}^{\text{LR}}=3$, at which point it acts as an infinite potential barrier that the effective mass cannot cross. Conversely, the effective mass plunges through the Schwarzschild LR radius in the cases of $H_{\text{SMR}}^{\rm EOB,PS}$ and $H_{\text{SMR-3PN}}^{\rm EOB,PS}$. This finding confirms that there is no unphysical behaviour at the LR radius for SMR Hamiltonians in the PS gauge.
To conclude, we also notice that the horizons of the $H_{\text{SMR}}^{\rm EOB,PS}$ and $H_{\text{SMR-3PN}}^{\rm EOB,PS}$ models (red and blue dots) are quite close to the LR. Such a large deformation from the Schwarzschild background, while not presenting an issue by itself, could pose problems in the evolution of the EOB dynamics after the LR radius and, thus, in the modelling of EOB waveforms and frequencies during the transition between plunge and merger-ringdown phases.

\subsection{Comparisons against numerical relativity}
\label{sec:energetics}

\begin{table}
	\caption{{\bf Two-body EOB Hamiltonians evolved.} We summarize the EOB Hamiltonians evolved.
		\label{table:models}}
	\begin{ruledtabular}
		\begin{tabular}{lp{5.2cm}p{2.5cm}}
			\addlinespace[.1cm]
			$H_{\text{SMR}}^{\rm EOB,PS}$
			& SMR Hamiltonian in PS gauge	& This paper \\
			\midrule
			\addlinespace[.1cm]
			$H_{\text{SMR-3PN}}^{\rm EOB,PS}$
			& SMR-3PN Hamiltonian in PS gauge & This paper \\
			\midrule
			\addlinespace[.1cm]
			$H_{\text{SMR}}^{\rm EOB}$
			& SMR Hamiltonian in the DJS gauge (with LR divergence)	& \cite{Barausse:2011dq} \\
			\midrule
			\addlinespace[.1cm]
			$H_{n{\rm PN}}^{\rm EOB,PS}$
			&   $n$PN  Hamiltonian in PS gauge & \cite{Damour:2017zjx}	\\
			\midrule
			\addlinespace[.1cm]
			$H_{n\rm PN}^{\rm EOB}$
			&	$n$PN  Hamiltonian in DJS gauge& \cite{Buonanno:1998gg,Damour:2000we} 	\\
			\addlinespace[.1cm]
		\end{tabular}
	\end{ruledtabular}
\end{table}

\begin{table}
	\caption{{\bf Set of non-spinning NR simulations and alignment time-windows.}  We list the SXS IDs, the mass ratios $q$ and the number of orbital cycles $N_{\rm orb}^{\rm merg}$ from the beginning of the simulation up to the binary black-hole merger (peak of $h^{\rm NR}_{22}$), as reported in the SXS catalog. We further include the time $t^\text{alig}_\text{in}$ at which the alignment procedure starts, the time $t^\text{alig}_\text{fin}$ at which it ends (in units of $M$) and the estimated NR error at merger $\Delta \phi^\text{merg}_\text{NR}$ (in radians).
		\label{table:setofNR}}
	\begin{ruledtabular}
		\begin{tabular}{ccc  ccc}
			\toprule
			SXS ID:						&	q$^{-1}$	& $N_{\rm orb}^{\rm merg}$	&    $t^\text{alig}_\text{in}$					&	$t^\text{alig}_\text{fin}$		&$\Delta \phi^\text{merg}_\text{NR}$ 	\\
			\midrule
			\addlinespace[.02cm]
			\hline
			\addlinespace[.1cm]
			0180	&	1	& 28.18	& 820	&	2250   & $\pm$0.25 \\
			1222	&	2	& 28.76	& 1000	&	2555  & $\pm$1.26 \\
			1221	&	3	& 27.18	& 1800	&	3000  & $\pm$0.21 \\
			1220	&	4	& 26.26	& 1800	&	3000   & $\pm$1.82 \\
			0056	&	5	& 28.81	& 1500	&	3000   & $\pm$0.39 \\
			0181	&	6	& 26.47	& 1000	&	2500   & $\pm$0.01 \\
			0298	&	7	& 19.68	& 780	&	2180   & $\pm$0.10 \\
			0063	&	8	& 25.83	& 1140	&	2540   & $\pm$0.85 \\
			0301	&	9	& 18.93	& 780	&	2180   & $\pm$0.13 \\
			0303	&	10	& 19.27	& 700	&	1900   & $\pm$0.49 \\

		\end{tabular}
	\end{ruledtabular}
\end{table}

\begin{table*}
	\caption{{\bf Details of the dephasing comparison.} We report the dephasing (in radians) of the SMR and 3PN models in both gauges at 8 and 4 GW cycles before NR merger, as found using the time-windows of Table~\ref{table:setofNR}. We also report the corresponding estimated NR error, which we denote by $\Delta \phi_\text{NR}$. The error for each NR simulation is estimated taking the phase differences between the highest two resolutions of the NR simulation (at fixed extrapolation order) and between two successive extrapolation orders (at fixed resolution), and adding them in quadrature.
		\label{table:simulations0}}
	\begin{ruledtabular}
		\begin{tabular}{cc | cccc c | cccc c  cccc c| ccc ccc }
			\toprule
			&\multirow{2}{*}{$q^{-1}$}
			&\multicolumn{4}{c}{8 GW cycles before merger
			& \multirow{3}{*}{$\Delta\phi_{\text{NR}}$ }
			&\multicolumn{4}{c}{4 GW cycles before merger }
			& \multirow{3}{*}{$\Delta\phi_{\text{NR}}$ }
		    \\

		    &

		    & $\Delta \phi_{\text{SMR}}^{\rm EOB,PS}$
		    & $\Delta \phi_{\text{SMR-3PN}}^{\rm EOB,PS}$
		    & $\Delta \phi_{\text{3PN}}^{\rm EOB,PS}$
		    & $\Delta \phi_{\text{3PN}}^{\rm EOB}$
		    &

		    & $\Delta \phi_{\text{SMR}}^{\rm EOB,PS}$
		    & $\Delta \phi_{\text{SMR-3PN}}^{\rm EOB,PS}$
		    & $\Delta \phi_{\text{3PN}}^{\rm EOB,PS}$
		    & $\Delta \phi_{\text{3PN}}^{\rm EOB}$
		    &

			\\
			\midrule
			\addlinespace[.02cm]
			\hline
			\addlinespace[.1cm]
			&	1	& ~0.111	& -0.033	&	-0.971   & ~0.032 & $\pm$0.032 & ~0.352 & -0.012 & -2.630 & ~0.084 & $\pm$0.056\\
			&	2	& ~0.112	& -0.061	&	-1.342   & -0.023 & $\pm$0.105 & ~0.512 & -0.021 & -5.586 & -0.043 & $\pm$0.224\\
		    & 	3	& ~0.050	& -0.021	&	-0.617   & -0.023 & $\pm$0.093 & ~0.111 & -0.026 & -1.209 & -0.048 & $\pm$0.144\\
			&	4	& ~0.046	& -0.038	&	-0.859   & -0.078 & $\pm$0.203 & ~0.187 & -0.041 & -2.540 & -0.212 & $\pm$0.372\\
			&	5	& ~0.037	& -0.034	&	-0.846   & -0.086 & $\pm$0.023 & ~0.125 & -0.044 & -2.077 & -0.211 & $\pm$0.064\\
			&	6	& -0.035	& -0.064	&	-0.433   & -0.093 & $\pm$0.006 & -0.041 & -0.082 & -0.599 & -0.126 & $\pm$0.007\\
			&	7	& ~0.024	& -0.009	&	-0.462   & -0.070 & $\pm$0.001 & ~0.092 & -0.003 & -1.403 & -0.211 & $\pm$0.009\\
			&	8	& ~0.021	& -0.021	&	-0.676   & -0.107 & $\pm$0.057 & ~0.076 & -0.025 & -1.660 & -0.260 & $\pm$0.155\\
			&	9	& ~0.017	& -0.005	&	-0.368   & -0.068 & $\pm$0.002 & ~0.063 & -0.005 & -1.185 & -0.220 & $\pm$0.012\\
			&	10	& ~0.022	& -0.001	&	-0.413   & -0.076 & $\pm$0.033 & ~0.070 & -0.004 & -1.245 & -0.233 & $\pm$0.083\\

		\end{tabular}
	\end{ruledtabular}
\end{table*}

\begin{table}
	\caption{{\bf Alternative alignment time-windows.}Time-windows (in units of $M$) employed for Fig.~\ref{fig:dphivsq}: here, $t^\text{alig}_\text{in}$ is the time corresponding to 34 GW cycles before merger for each NR simulation, whereas $t^\text{alig}_\text{fin}$ is chosen to encompass 10 GW cycles. The time at merger is given by $t_\text{merg}$.
	\label{table:simulations}}
	\begin{ruledtabular}
		\begin{tabular}{cccc | cccc }
			\toprule
		    q$^{-1}$			&   $t^\text{alig}_\text{in}$					&	$t^\text{alig}_\text{fin}$ &	$t_\text{merg}$			& q$^{-1}$			&   $t^\text{alig}_\text{in}$					&	$t^\text{alig}_\text{fin}$ &	$t_\text{merg}$	\\
			\midrule
			\addlinespace[.02cm]
			\hline
			\addlinespace[.1cm]
				1	& 5107	& 6911 & 9517	&	6   & 2971 & 4254 & 6000 \\
				2	& 5406	& 7078 & 9384   &   7   & 776 & 2083 & 4142\\
				3   & 3940	& 5532 & 7858   & 8 & 2652 & 3918& 5956\\
				4   & 3479	& 4975 & 7200	& 9 & 513 & 1732 & 3692\\
				5   & 4206	& 5641 & 7864   & 10  &  587 & 1771 & 3691	\\

		\end{tabular}
	\end{ruledtabular}
\end{table}

\begin{figure}
	\includegraphics[width=\columnwidth]{figuresSMR/ejq1.pdf}
	\includegraphics[width=\columnwidth]{figuresSMR/ejq10.pdf}
	\caption{\textbf{SMR vs PN binding energies:} we compare the difference $\Delta E$ in binding energy from NR for our SMR Hamiltonians versus angular momentum $l$. We compare it to similar results for PN models up to third order, in both PS and DJS gauges. The \textit{estimated} NR error is shown in grey.}
	\label{fig:compEJ}
\end{figure}

\begin{figure}
	\includegraphics[width=\columnwidth]{figuresSMR/eoq1.pdf}
	\includegraphics[width=\columnwidth]{figuresSMR/eoq10.pdf}
	\caption{\textbf{SMR vs PN binding energies:} we compare the difference $\Delta E$ in binding energy from NR for our SMR Hamiltonians versus frequency $(M\Omega)$. We compare it to similar results for PN models up to third order, in both PS and DJS gauges. The \textit{estimated} NR error is shown in grey.}
	\label{fig:compEO}
\end{figure}

Here we study the energetics of the $H_{\text{SMR}}^{\rm EOB,PS}$ and $H_{\text{SMR-3PN}}^{\rm EOB,PS}$ models and the PN EOB models in both gauges via comparisons of their binding energies against NR predictions. The EOB Hamiltonians evolved and their notation are summarized in Table~\ref{table:models}.
The (quasi) gauge-invariant relations between the dimensionless circular orbit binding energy $E \equiv (H-M)/\mu$ and angular momentum $ l \equiv \hat{p}_\phi=p_\phi/(M\mu)$ (and orbital frequency $\hat\Omega$) are used to draw comparisons against NR.
This type of comparisons is useful to understand how information of the real two-body motion is resummed into the conservative dynamics~\cite{Antonelli:2019ytb}.
In contrast to Ref.~\cite{Antonelli:2019ytb} and Sec.~\ref{sec:EOBSMR} of this paper, where the binding energy is calculated in the circular-orbit limit, the binding energies appearing in this section are obtained evolving the EOB Hamiltonians along quasi-circular orbits. This more closely matches the procedure used to extract the binding energy from NR simulations of quasi-circular inspirals, providing clearer comparisons \cite{Ossokine:2017dge}.
Finally, we calculate the dephasing $\Delta\phi_{22} \equiv  \phi_\text{NR}-\phi_\text{EOB}$ of the ($\ell,m$)=(2,2) modes of the $H_{\text{SMR}}^{\rm EOB,PS}$ and $H_{\text{SMR-3PN}}^{\rm EOB,PS}$ models against NR results. While more thorough comparisons aimed at using the models for LIGO inference studies would need a systematic calculation of the unfaithfulness (see e.g., Refs.~\cite{Pan:2011gk,Taracchini:2013rva,Bohe:2016gbl,Cotesta:2018fcv}), we find these comparisons illustrative to contextualize the $H_{\text{SMR}}^{\rm EOB,PS}$ and $H_{\text{SMR-3PN}}^{\rm EOB,PS}$ models in this paper.

We employ a set of ten non-spinning NR simulations from the Simulating eXtreme Spacetimes (SXS) collaboration \cite{Mroue:2013xna,Chu:2015kft}, with mass ratios $1/10\leq q\leq 1$. We summarize the details of these simulations in Table~\ref{table:setofNR}.
A description of how the $E(l)$ and $E(\hat\Omega)$ curves were calculated for a subset of these simulations can be found in Ref.~\cite{Ossokine:2017dge}.

We evolve EOB Hamiltonians with PN information up to third order, since 3PN is the order at which PS-gauge Hamiltonians can be uniquely derived for generic orbits (see the Appendix of Ref.~\cite{Antonelli:2019ytb} for more details). It is worthwhile to mention that the $H_{\rm 3PN}^{\rm EOB}$ Hamiltonian has better energetics and phases performances against NR than both $H_{\rm 4PN}^{\rm EOB}$ and the SEOBNR Hamiltonian used as a baseline for the current generation of EOB waveform models (defined, e.g., in the Appendix of Ref.~\cite{Steinhoff:2016rfi}), when calibration and NQC parameters are turned off. Restricting ourselves to comparisons with $H_{\rm 3PN}^{\rm EOB}$ only, we are therefore not running the risk to overestimate the performance of SMR models when comparing them to PN results.

Let us begin comparing the $E(l)$ and $E(\hat\Omega)$ curves. The difference $\Delta E \equiv |E_\text{NR}-E_\text{EOB}|$ is plotted for a variety of EOB models in Figs.~\ref{fig:compEJ} and \ref{fig:compEO}.
Considering the $E(l)$ relations first and focusing on the SMR models, it is seen that for $q=1/10$ both $H^{\rm EOB,PS}_{\text{SMR}}$ and $H^{\rm EOB,PS}_{\text{SMR-3PN}}$ perform better against NR than the 3PN model in the same gauge, e.g., $H^{\rm EOB,PS}_{{\rm 3PN}}$.
The $H^{\rm EOB,PS}_{\text{SMR-3PN}}$ model also performs better than both in the comparable-mass case.
A similar finding is obtained investigating the $E(\hat\Omega)$ curves, see Fig.~\ref{fig:compEO}.
Taken together, these results highlight the importance of SMR results to improve the modeling of both equal- and unequal-mass systems within the EOB approach.
It is also seen that, for both mass ratios considered and for both $E(l)$ and $E(\hat\Omega)$ curves, $H^{\rm EOB,PS}_{\text{SMR-3PN}}$ improves the predictions of $H^{\rm EOB,PS}_{\text{SMR}}$, suggesting that generic orbit terms are important when considering quasi-circular orbit binding energies (especially in the equal-mass-ratio case).

PN Hamiltonians in the PS gauge generically perform worse in binding energy comparisons than Hamiltonians in the DJS gauge, as found out in the adiabatic approximation already in Ref.~\cite{Antonelli:2019ytb}.
This finding suggests that, notwithstanding the already good agreement between SMR models and NR simulations for both mass ratios, a better description for the EOB dynamics than the one provided by the PS gauge could be pursued in order to maximize the performance of evolutions from both PN and SMR EOB models.

\begin{figure*}[t]
	\includegraphics[width=\columnwidth]{figuresSMR/wavephaseq1.pdf}
	\hspace{12pt}
	\includegraphics[width=\columnwidth]{figuresSMR/wavephaseq10.pdf}
	\caption{\textbf{Dephasing of EOB models:} in the top panels, the real parts $\mathfrak{R}(h_{22})$ of the ($\ell,m$)=(2,2) mode EOB waveform for the SMR, SMR-3PN models are shown and compared to the NR waveforms (in dashed-black, overlapping with the EOB waveforms up to few GW cycles to merger).
	In the lower panels, the dephasing of SMR and PN EOB models from the NR simulations is calculated.
	Also shown are the times corresponding to 8, 4 and 2 GW cycles before NR merger.  }
\label{fig:wavephase}
\end{figure*}
\begin{figure*}[t]
	\includegraphics[width=\columnwidth]{figuresSMR/dphivsqwerror_N4.pdf}
	\hspace{12pt}
	\includegraphics[width=\columnwidth]{figuresSMR/dphivsqwerror_N2.pdf}
	\caption{\textbf{Dephasing vs mass ratio:} we compare the dephasing of $H_{\text{SMR}}^{\rm EOB,PS}$, $H_{\text{SMR-PN}}^{\rm EOB,PS}$ and $H_{\text{3PN}}^{\rm EOB}$ after they have been aligned with the NR simulations from Table~\ref{table:simulations}. For each $q$, we snapshot the dephasing of the EOB models and the NR simulation at a time corresponding to 4 and 2 orbits before the merger of the binary system in the NR simulation. }
	\label{fig:dphivsq}
\end{figure*}

We complete our comparison study with the dephasing $\Delta\phi_{22}$ of the ($\ell,m$)=(2,2) modes from the EOB models and the NR simulations.
For a proper comparison, the EOB and NR waveforms must be aligned for each $q$. Here we use the alignment procedure outlined in Ref.~\cite{Pan:2011gk}, which amounts to minimizing the function:
\begin{equation}
\Xi(\Delta t,\Delta \phi)=\int_{t^\text{alig}_1}^{t^\text{alig}_2}[\phi_\text{NR}(t)-\phi_\text{EOB}(t+\Delta t)-\Delta\phi]^2 dt\\,
\end{equation}
over the time and phase shifts, $\Delta t$ and $\Delta\phi$.
The integrating interval [$t^\text{alig}_1,t^\text{alig}_2$] defines the time-domain window in which the alignment is performed: conservatively, it must be chosen in the inspiral of the NR simulation, large enough to average out the numerical noise and such as to avoid junk radiation at the beginning of the NR simulation \cite{Pan:2011gk}. From the alignment procedure described above, one can obtain the phase and amplitude time-shift to be applied to the EOB model to align it with the NR waveforms, i.e., the aligned waveforms are:
\begin{align}
&h^\text{NR}_{22}=A_\text{NR}(t)e^{i\phi_\text{NR}(t)}\,,\label{h22nr}\\
&h^\text{EOB}_{22}=A_\text{EOB}(t+\Delta t)e^{i[\phi_\text{EOB}(t+\Delta t)+\Delta\phi]}\label{h22eob}\,.
\end{align}

Our choices for the time-windows are reported in Table~\ref{table:setofNR}.
In Fig.~\ref{fig:wavephase}, we show the results of our phase comparisons for $q=1$ and $q=1/10$ up to merger. For clarity, the upper panels only include the $H_{\text{SMR}}^{\rm EOB,PS}$ and $H_{\text{SMR-3PN}}^{\rm EOB,PS}$ models and the NR simulations. They show the real parts of Eqs.~\eqref{h22nr} and ~\eqref{h22eob}, from which we infer that the SMR models do not accumulate a significant amount of dephasing. Overall, they are in very good agreement with NR for both $q=1$ and $q=1/10$.
It is important to place the above results in context. In the lower panel, the dephasing of SMR models from NR is compared to that of 3PN models\footnote{In this comparison we do not include 2PN models, which we find to have much larger dephasing than the 3PN models shown.}.
Interestingly, even in the equal-mass-ratio case $H_{\text{SMR}}^{\rm EOB,PS}$ and $H_{\text{SMR-3PN}}^{\rm EOB,PS}$ compare much better than the 3PN model in the same gauge, e.g., $H_{\rm 3PN}^{\rm EOB,PS}$. Their dephasing is comparable to  $H_{\rm 3PN}^{\rm EOB}$. In the $q=1/10$ case,
they have a smaller dephasing than any other PN model considered in this study.
In Table~\ref{table:simulations0}, we report the dephasing that the $H_{\text{SMR}}^{\rm EOB,PS}$, $H_{\text{SMR-3PN}}^{\rm EOB,PS}$, $H_{\text{3PN}}^{\rm EOB}$ and $H_{\text{3PN}}^{\rm EOB,PS}$ models accumulate up to 8 and 4 GW cycles before merger for all mass ratios (with the corresponding estimated NR error)\footnote{We have checked that shifting the time-windows by $\Delta t=\pm 100\,M$,  our $\Delta \phi$'s only change by a few hundredths of a radian.}.

Next, we want to study how the dephasing of the above models varies as a function of $q$. It would be tempting to compare the $\Delta \phi$'s reported in Table~\ref{table:simulations0} at a fixed number of cycles before merger. While this remains a valid possibility, such a comparison would neither take into account the different lengths of the NR simulations used in this set, nor the different number of GW cycles encompassed by the time-windows of Table~\ref{table:setofNR}.
To keep both parameters under control, we realign our models with alternative time-windows that are dictated by the number of GW cycles to merger $\Delta N_\text{GW}(t) \equiv N_\text{GW}(t)-N_\text{GW}^\text{merg}
$ of the NR simulations. That is, for each mass ratio we fix a different time-window [$t^\text{alig}_1,t^\text{alig}_2$], corresponding to the \textit{same} interval of cycles to merger [$\Delta N_\text{GW}(t^\text{alig}_1)$, $\Delta N_\text{GW}(t^\text{alig}_2)$].
The benefits of this choice are two-fold. To begin with, the alignment windows thus calculated depends on the position of the NR merger (peak of $h_{22}^{\rm NR}$), which is a quantifiable feature of every NR simulation.
Moreover, this choice allows us to assess trends across the mass ratios fairly, since the waveforms thus aligned are compared in the same range of GW cycles. A caveat for this alignment method is that the GW cycles of evolutions with smaller $q$ lie in a regime of stronger gravity.

We choose to align the EOB models to NR in an interval of $N_\text{GW}$ such that
$[\Delta N_\text{GW}(t^\text{alig}_1), \Delta N_\text{GW}(t^\text{alig}_2)]=[-34, -24]$, corresponding to the time-windows reported in Table~\ref{table:simulations}.
This choice stems from the length of the shortest NR simulation, e.g., $q=1/9$, which counts $N_\text{GW}^\text{merg}=37.96$ GW cycles at merger (the first ~3GW cycles of this simulation are neglected in order to avoid junk radiation). In Fig.~\ref{fig:dphivsq}, we plot the dephasing for the three models that perform best in Fig.~\ref{fig:wavephase}: that is, $H_{\rm 3PN}^{\rm EOB}$, $H_{\text{SMR}}^{\rm EOB,PS}$ and $H_{\text{SMR-3PN}}^{\rm EOB,PS}$ and study the trends across $q$. For every simulation, we calculate the dephasing 8 and 4 GW cycles before merger to show the robustness of the trends\footnote{We have also checked that the trends are unaffected by variations in the number of orbital cycles in the alignment window.}.
Noticeably, the 3PN EOB waveform in the DJS gauge starts degrading in accuracy as the mass ratio is increased, while the SMR and SMR-3PN ones improve: remarkably, for most $q$'s, the SMR-3PN model only dephases by a few hundredths of a radian
up to a 4 GW cycles before merger.
Moreover, we notice that SMR models start performing better than $H_{\rm 3PN}^{\rm EOB}$ for $q \lesssim 1/3$, hinting again to the
fact that SMR information, when reorganized in the EOB framework, could be used to model systems that are very close to the equal-mass-ratio regime \cite{LeTiec:2011bk,Barausse:2011dq}.

The picture emerging from Fig.~\ref{fig:dphivsq} is that the SMR-3PN model is the most consistent of the two models with SMR information, corroborating the findings for $q=1$ and $q=1/10$ in the binding energy comparisons. The small dephasing of the SMR-3PN model suggests that the Hamiltonian upon which it is based is a possible starting point to develop a new generation of EOB waveform models able to tackle the currently challenging intermediate-mass-ratio regime.

\section{Discussion and Conclusions}
\label{sec:conclusions}

The complete EOB Hamiltonian at linear order in SMR from Ref.~\cite{LeTiec:2011dp} suffers from a coordinate singularity at the LR radius in the deformed Schwarzschild background. Building on Refs.~\cite{Akcay:2012ea,Damour:2017zjx}, we have constructed two Hamiltonians in the post-Schwarzschild (PS) reformulation of the EOB approach \cite{Damour:2017zjx,Antonelli:2019ytb} (both with the SMR correction to the Detweiler redshift and with mixed SMR-3PN information), and checked
that they are not affected by poles  at the LR radius (and related unphysical features) by studying plunging trajectories.

We have then explored the merits of the SMR and mixed SMR-3PN Hamiltonians via comparisons of their waveforms and binding energies, and those of PN Hamiltonians in different gauges, against NR predictions. Ultimately, we find that:
\begin{enumerate}
	\item For both $q=1$ and $q=1/10$, the binding energies of SMR and SMR-3PN EOB models (see Figs.~\ref{fig:compEJ} and~\ref{fig:compEO}) generally compare better against NR than the binding energy of the PS Hamiltonian with 3PN information.
	\item The generic orbit 3PN information in the SMR-3PN EOB Hamiltonian improves the binding energy and phase comparisons of SMR EOB models.

	\item PN Hamiltonians in the EOB-PS gauge have binding energies that compare worse than those from PN Hamiltonians in the standard EOB gauge, confirming the findings of Ref.~\cite{Antonelli:2019ytb} and extending their validity to non-adiabatic evolutions.
	\item The SMR-3PN EOB model agrees remarkably well against NR simulations, see Fig.~\ref{fig:dphivsq}. The dephasing up to 4 GW cycles before merger is a few hundredths of a radian for $q \lesssim 1/3$ and a tenth of a radian for $q>1/3$. The only EOB PN model with comparable dephasing is the 3PN EOB Hamiltonian in the DJS gauge for $q \gtrsim 1/3$.
\end{enumerate}

The construction of the SMR EOB Hamiltonian in this paper depends on a number of choices. First of all, we chose to fix the coordinate freedom in the effective Hamiltonian using the PS gauge. This was chosen because of its relative simplicity, while allowing a natural path towards avoiding singularities at the light ring. However, there may exist different choices that are equally (or more) effective. Second, while the EOB Hamiltonian in principle applies to generic orbits, we fix the linear-in-$\nu$ part only by comparing to the circular-orbit binding energy. Consequently, there is considerable freedom in the ``non-circular-orbit'' part of the Hamiltonian. In practice, we fix this freedom by choosing the specific functional dependence of the effective Hamiltonian on $\hat{H}_{\text{S}}$ given by Eq.~\eqref{HSMR}. This choice is in part restricted by the requirement that the Hamiltonian be analytic, but other options are available. Third and finally, SMR data for the binding energy extends only to the light ring. The Hamiltonian in the region $u>\tfrac{1}{3}$ therefore depends only on the analytic extension of the redshift data. Given that this data is known only to finite numerical precision, there is some freedom in the choice of the exact analytical form of its fit. This choice can also affect the relative size of the different coefficient functions in Eq.~\eqref{HSMR}.

Our investigation opens up further avenues of research. To begin with, one can study whether it is possible to uniquely fix other EOB gauges that could accommodate the Detweiler redshift (without introducing a LR-divergence) and study their merits via comparisons against NR. As discussed already in Ref.~\cite{Akcay:2012ea}, to solve the LR-divergence arising in this context the non-geodesic function $\hat{Q}$ needs a term proportional to $p_{\phi}^3$, possibly resummed in another quantity (as done in the PS gauge using $\hat{H}_\text{S}$). It would be quite interesting to see whether other gauges that allow solving the LR-divergence also improve the comparisons against NR predictions. One concrete example of different resummation that was shown to improve the comparisons of the conservative sector of post-Minkowskian Hamiltonians in PS form has been given in the Appendix of Ref.~\cite{Antonelli:2019ytb}. It is worthwhile to study whether a similar choice  could work for the SMR and SMR-3PN models herein presented.
The hope is that using different resummations, and including information from the second order in the SMR, one could obtain a considerably improved EOB Hamiltonian that, after further calibration to NR, would be very useful for LIGO/Virgo analyses in the near-future.

Further research endeavours could be directed towards informing the EOB with different SMR quantities than the circular orbit Detweiler redshift. An example of a quantity that still needs to be fully exploited is the generalized redshift \cite{Barack:2011ed,Akcay:2015pza}, which includes information for arbitrarily eccentric orbits.
We envision using EOB Hamiltonians at linear and higher orders in the
mass ratio for inference studies in the future detectors' era, when
precise models will be needed to properly characterize high
signal-to-noise systems, possibly having rather small mass ratios. In
order for this program to be achieved, not only should the
conservative sector be optimized with both results at second order in
$q$ and (potentially) a better resummation, but information from other
crucial physical quantities should also be incorporated: notably missing
features in our analysis are the spin and eccentricity. Furthermore, a
more comprehensive study of the dissipative sector must be pursued. It
would be desirable, for instance, to include more self-force
information in the flux. Lastly, we would also need to build the full
inspiral, merger and ringdown waveforms, and calibrate them to NR simulations.
We leave these important investigations to future work.

\begin{acknowledgments}
It is a pleasure to thank Sergei Ossokine for providing us with the NR data used for the comparisons, and Tanja Hinderer for the original version of the EOB evolution code in Mathematica used in this paper. A.A. would like to further thank Sergei Ossokine and Roberto Cotesta for many fruitful discussions over the preparation of this paper. MvdM was supported by European Union's Horizon 2020 research and innovation programme under grant agreement No.~705229.
This work makes use of the Black-Hole Perturbation Toolkit \cite{BHPToolkit}.
\end{acknowledgments}