View source for Busemann function

{{CS1 config|mode=cs2}}
{{multiple issues|
{{technical|date=June 2020}}
{{more citations|date=June 2020}}
}}

In [[geometric topology]], '''Busemann functions''' are used to study the large-scale geometry of geodesics in [[Hadamard space]]s and in particular [[Hadamard manifold]]s ([[simply connected]] complete [[Riemannian manifold]]s of nonpositive curvature). They are named after [[Herbert Busemann]], who introduced them; he gave an extensive treatment of the topic in his 1955 book "The geometry of geodesics".

== Definition and elementary properties ==

Let <math>(X,d)</math> be a [[metric space]]. A [[Ray (geometry)|geodesic ray]] is a path <math>\gamma : [0,\infty) \to X</math> which minimizes distance everywhere along its length. i.e., for all <math>t,t' \in [0,\infty)</math>, 
<math display=block>d\big(\gamma(t), \gamma(t') \big) = \big| t - t' \big|.</math> 
Equivalently, a ray is an isometry from the "canonical ray" (the set <math>[0,\infty)</math> equipped with the Euclidean metric) into the metric space ''X''.

Given a ray ''γ'', the Busemann function <math>B_\gamma : X \to \mathbb R</math> is defined by

<math display=block>B_\gamma(x)=\lim_{t\to\infty}\big(d\big( \gamma(t), x \big)  - t \big)</math>

Thus, when ''t'' is very large, the distance <math>d\big( \gamma(t), x \big)</math> is approximately equal to <math> B_\gamma(x) + t </math>. Given a ray ''γ'', its Busemann function is always well-defined: indeed the right hand side above <math>F_t(x) \stackrel{\text{def}}{=} d\big( \gamma(t), x \big)  - t  </math>, tends pointwise to the left hand side on compacta, since <math>t-d( \gamma(t), x )=d(\gamma(t),\gamma(0)) - d(\gamma(t),x)</math> is bounded above by <math> d(\gamma(0),x)</math> and non-increasing since, if <math>s\le t</math>,

<math display=block>t-s + d(x,\gamma(s)) - d(x,\gamma(t)) \ge t -s -d (\gamma(s),\gamma(t)) = 0.</math>

It is immediate from the triangle inequality that

<math display=block>|B_\gamma(x) - B_\gamma(y)|\le d(x,y),</math>

so that <math>B_\gamma</math> is uniformly continuous. More specifically, the above estimate above shows that

*'''Busemann functions are [[Lipschitz function]]s with constant 1'''.<ref>{{harvnb|Busemann|1955|page=131}}</ref>

By [[Dini's theorem]], the functions <math>F_t(x)=d(x,\gamma(t)) -t </math> tend to <math>B_\gamma(x)</math> uniformly on compact sets as ''t'' tends to infinity.

== Example: Poincaré disk ==

Let <math>D</math> be the unit disk in the complex plane with the [[poincare disk model|Poincaré metric]]

<math display=block> ds^2 = {4 \,|dz|^2 \over (1 - |z|^2)^2}.</math>

Then, for <math>|z| < 1</math> and <math>|\zeta| = 1</math>, the Busemann function is given by<ref>{{harvnb|Bridson|Haefliger|1999|page=273}}</ref>

<math display=block> B_\zeta(z) = - \log\left({1 - |z|^2\over |z - \zeta|^2}\right),</math>

where the term in brackets on the right hand side is the [[Poisson kernel]] for the unit disk and <math>\zeta</math> corresponds to the radial geodesic <math>\gamma</math> from the origin towards <math>\zeta</math>,
<math>\gamma(t) = \zeta \tanh(t/2)</math>. The computation of <math>d(x,y)</math> can be reduced to that of <math>d(z,0) = d(|z|,0) = 2 \operatorname{artanh}(|z|) = \log \left( \tfrac{1+|z|}{1-|z|} \right) </math>, since the metric is invariant under [[Möbius transformation]]s in <math>SU(1,1)</math>; the geodesics through <math>0</math> have the form <math>\zeta g_t(0)</math> where <math>g_t</math> is the 1-parameter subgroup of <math>SU(1,1)</math>,

<math display=block>g_t=\begin{pmatrix}\cosh (t/2) & \sinh (t/2) \\ \sinh (t/2) & \cosh (t/2)\end{pmatrix}</math>

The formula above also completely determines the Busemann function by Möbius invariance.

==Busemann functions on a Hadamard space==

In a [[Hadamard space]], where any two points are joined by a unique geodesic segment, the function <math>F = F_t</math> is ''convex'', i.e. convex on geodesic segments <math>[x,y]</math>. Explicitly this means that if <math>z(s)</math> is the point which divides <math>[x,y]</math> in the ratio {{math|1=''s'' : (1 − ''s'') }}, then <math>F(z(s)) \leq s F(x) + (1 - s) F(y)</math>.  For fixed <math>a</math> the function <math>d(x,a)</math> is convex and hence so are its translates; in particular, if <math>\gamma</math> is a geodesic ray in <math>X</math>, then <math>F_t</math> is convex. Since the Busemann function <math>B_\gamma</math> is the pointwise limit of <math>F_t</math>,

*'''Busemann functions are convex on Hadamard spaces'''.<ref name="BGS">{{harvnb|Ballmann|Gromov|Schroeder|1985}}</ref>
*'''On a Hadamard space, the functions <math>F_t(y)=d(y,\gamma(t)) -t </math> converge uniformly to <math>B_\gamma</math> uniformly on any bounded subset of <math>X</math>.'''<ref>{{harvnb|Bridson|Haefliger|1999|pages=268–269}}</ref><ref>{{harvnb|Lurie|2010|page=13}}</ref>

Let {{math|1=''h''(''t'') = ''d''(''y'',γ(''t'')) − ''t'' = ''F''<sub>''t''</sub>(''y'')}}. Since <math>\gamma (t)</math> is parametrised by arclength, Alexandrov's first comparison theorem for Hadamard spaces implies that the function {{math|1=''g''(''t'') = ''d''(''y'',γ(''t''))<sup>2</sup> − ''t''<sup>2</sup>}} is convex. Hence for {{math|0< ''s'' < ''t''}}

<math display=block>g(s) \le (1-{s\over t})g(0) + {s\over t} g(t).</math>

Thus

<math display=block>2sh(s)\le (h(s) +s)^2 - s^2 =g(s)  \le (1-{s\over t})d(x,y)^2 + {s\over t} ( 2th(t) +h(t)^2 )\le d(x,y)^2 + 2s h(t) +{s\over t} d(x,y)^2,</math>

so that

<math display=block>|F_s(y) - F_t(y)|= |h(s) - h(t)| \le {1\over 2} (s^{-1} + t^{-1}) d(x,y)^2.</math>

Letting ''t'' tend to ∞, it follows that

<math display=block>|F_s(y) - B_\gamma(y)| \le {d(x,y)^2\over 2s},</math>

so convergence is uniform on bounded sets.

Note that the inequality above for <math>|F_s(y) - F_t(y)| </math> (together with its proof) also holds for geodesic segments: if {{math|&gamma;(''t'')}} is a geodesic segment starting at {{math|''x''}} and parametrised by arclength then

<math display=block>|d(y,\Gamma(s)) -s -d(y,\Gamma(t)) +t| \le (s^{-1} + t^{-1}) d(x,y)^2.</math>

Next suppose that {{math|''x'', ''y''}}  are points in a Hadamard space, and let {{math|1=&delta;(''s'')}} be the geodesic through {{math|''x''}} with {{math|1=&delta;(0) = ''y''}} and {{math|1=&delta;(''t'') = ''x''}}, where {{math|1=''t'' = ''d''(''x'',''y'')}}. This geodesic cuts the boundary of the closed ball {{math|{{overline|''B''}}(''y'',''r'')}} at the point {{math|&delta;(''r'')}}. Thus if {{math|1=''d''(''x'',''y'') > ''r''}}, there is a point {{math|''v''}} with {{math|1=''d''(''y'',''v'') = ''r''}} such that {{math|1=''d''(''x'',''v'')  = ''d''(''x'',''y'') − ''r''}}.

This condition persists for Busemann functions. The statement and proof of the property for Busemann functions relies on a fundamental theorem on closed convex subsets of a Hadamard space, which generalises [[orthogonal projection]] in a [[Hilbert space]]: if {{math|''C''}} is a closed convex set in a Hadamard space {{math|''X''}}, then every point {{math|''x''}} in {{math|''X''}} has a unique closest point {{math|1=''P''(''x'') ≡ ''P''<sub>''C''</sub>(''x'')}} in {{math|''C''}} and {{math|1=''d''(''P''(''x''),''P''(''y'')) ≤ ''d''(''x'',''y'')}}; moreover {{math|1=''a'' = ''P''(''x'')}} is uniquely determined by the property that, for {{math|''y''}} in {{math|''C''}},

<math display=block>d(x,y)^2 \ge d(x,a) ^2 + d(a,y)^2,</math>

so that the angle at {{math|''a''}} in the Euclidean [[comparison triangle]] for {{math|''a'',''x'',''y''}} is greater than or equal to {{math|''&pi;''/2}}.

*'''If {{math|h}} is a Busemann function on a Hadamard space, then, given {{math|y}} in {{math|X}} and {{math|r > 0}}, there is a unique point {{math|v}} with {{math|1=d(y,v) = r}} such that {{math|1=h(v) =  h(y) − r}}. For fixed {{math|r > 0}}, the point {{math|v}} is the closest point of {{math|y}} to the closed convex {{math|C}} set of points {{math|u}} such that {{math|h(u) &le; h(y) − r}} and therefore depends continuously on {{math|y}}.'''<ref>{{harvnb|Bridson|Haefliger|1999|pages=271–272}}</ref>

Let {{math|''v''}} be the closest point to {{math|''y''}} in ''{{math|C}}''. Then {{math|1=''h''(''v'') = ''h''(''y'') − ''r''}} and so ''{{math|h}}'' is minimised by ''{{math|v}}'' in {{math|{{overline|''B''}}(''y'',''R'')}} where {{math|1=''R'' = ''d''(''y'',''v'') and ''v''}} is the unique point where ''{{math|h}}'' is minimised. By the Lipschitz condition {{math|1=''r'' = {{!}}''h''(''y'') − ''h''(''v''){{!}} ≤ ''R''}}. To prove the assertion, it suffices to show that {{math|1=''R'' = ''r''}}, i.e. {{math|1=''d''(''y'',''v'') = ''r''}}. On the other hand, ''{{math|h}}'' is the uniform limit on any closed ball of functions {{math|''h''<sub>''n''</sub>}}. On {{math|{{overline|''B''}}(''y'',''r'')}}, these are minimised by points {{math|''v''<sub>''n''</sub>}} with {{math|1=''h''<sub>''n''</sub>(''v''<sub>''n''</sub>) = ''h''<sub>''n''</sub>(''y'') − ''r''}}. Hence the infimum of ''{{math|h}}'' on {{math|{{overline|''B''}}(''y'',''r'')}} is {{math|''h''(''y'') − ''r''}} and {{math|''h''(''v''<sub>''n''</sub>)}} tends to {{math|''h''(''y'') − ''r''}}. Thus {{math|1=''h''(''v''<sub>''n''</sub>) = ''h''(''y'') − ''r''<sub>''n''</sub>}} with {{math|''r''<sub>''n''</sub> &le; ''r''}} and {{math|''r''<sub>''n''</sub>}} tending towards ''{{math|r}}''. Let {{math|''u''<sub>''n''</sub>}} be the closest point to {{math|''y''}} with {{math|''h''(''u''<sub>''n''</sub>) &le; ''h''(''y'') − ''r''<sub>''n''</sub>}}. Let {{math|1=''R''<sub>''n''</sub> = ''d''(''y'',''u''<sub>''n''</sub>) ≤ ''r''}}. Then {{math|1=''h''(''u''<sub>''n''</sub>) =  ''h''(''y'') − ''r''<sub>''n''</sub>}}, and, by the Lipschitz condition on {{math|''h''}}, {{math|''R''<sub>''n''</sub> &ge; ''r''<sub>''n''</sub>}}. In particular {{math|''R''<sub>''n''</sub>}} tends to ''{{math|r}}''. Passing to a subsequence if necessary it can be assumed that {{math|''r''<sub>''n''</sub>}} and {{math|''R''<sub>''n''</sub>}} are both increasing (to ''{{math|r}}''). The inequality for convex optimisation implies that for {{math|''n'' > ''m''}}.

<math display=block> d(u_n,u_m)^2 \le R_n^2 - R_m^2\le 2r|R_n -R_m|,</math>

so that {{math|''u''<sub>''n''</sub>}} is a Cauchy sequence. If ''{{math|u}}'' is its limit, then {{math|1=''d''(''y'',''u'') = ''r''}} and {{math|1=''h''(''u'') = ''h''(''y'') − ''r''}}. By uniqueness it follows that {{math|1=''u'' = ''v''}} and hence {{math|1=''d''(''y'',''v'') = ''r''}}, as required.

'''Uniform limits.''' The above argument proves more generally that if {{math|''d''(''x''<sub>''n''</sub>,''x''<sub>0</sub>)}} tends to infinity and the functions {{math|1=''h''<sub>''n''</sub>(''x'') = ''d''(''x'',''x''<sub>''n''</sub>) – ''d''(''x''<sub>''n''</sub>,''x''<sub>0</sub>)}} tend uniformly on bounded sets to {{math|''h''(''x'')}}, then ''{{math|h}}'' is convex, Lipschitz with Lipschitz constant 1 and, given ''{{math|y}}'' in ''{{math|X}}'' and {{math|''r'' > 0}}, there is a unique point ''{{math|v}}'' with {{math|1=''d''(''y'',''v'') = ''r''}} such that {{math|1=''h''(''v'') = ''h''(''y'') − ''r''}}. If on the other hand the sequence {{math|(''x''<sub>''n''</sub>)}} is bounded, then the terms all lie in some closed ball and uniform convergence there implies that {{math|(''x''<sub>''n''</sub>)}} is a Cauchy sequence so converges to some {{math|''x''<sub>&infin;</sub>}} in ''{{math|X}}''. So {{math|''h''<sub>''n''</sub>}} tends uniformly to {{math|1=''h''<sub>&infin;</sub>(''x'') = ''d''(''x'',''x''<sub>&infin;</sub>) – ''d''(''x''<sub>&infin;</sub>,''x''<sub>0</sub>)}}, a function of the same form. The same argument also shows that the class of functions which satisfy the same three conditions (being convex, Lipschitz and having minima on closed balls) is closed under taking uniform limits on bounded sets.

'''Comment.''' Note that, since any closed convex subset of a Hadamard subset of a Hadamard space is also a Hadamard space, any closed ball in a Hadamard space is a Hadamard space. In particular it need not be the case that every geodesic segment is contained in a geodesic defined on the whole of '''{{math|R}}''' or even a semi-infinite interval {{math|[0,&infin;)}}. The closed unit ball of a Hilbert space gives an explicit example which is not a proper metric space.

*'''If {{math|h}} is a convex function, Lipschitz with constant 1 and {{math|h}} assumes its minimum on any closed ball centred on {{math|y}} with radius {{math|r}} at a unique point {{math|v}} on the boundary with {{math|1=h(v) = h(y) − r}}, then for each {{math|y}} in {{math|X}} there is a unique geodesic ray {{math|&delta;}} such that {{math|1=&delta;(0) = y}} and {{math|&delta;}} cuts each closed convex set {{math|h &le; h(y) – r}} with {{math|r > 0}} at {{math|&delta;(r)}}, so that {{math|1=h(&delta;(t)) = h(y) – t}}. In particular this holds for each Busemann function.'''<ref name="Bridson 1999 pages=271–272">{{harvnb|Bridson|Haefliger|1999|pages=271–272}}</ref>

The third condition implies that ''{{math|v}}'' is the closest point to ''{{math|1=y}}'' in the closed convex set {{math|1=''C''<sub>''r''</sub>}} of points ''<math>u</math>'' such that {{math|''h''(''u'') &le; ''h''(''y'') – ''r''}}. Let {{math|&delta;(''t'')}} for {{math|1=0 &le; ''t'' &le; ''r''}} be the geodesic joining ''{{math|y}}'' to ''{{math|v}}''. Then {{math|1=''k''(''t'') = ''h''(&delta;(''t'')) - ''h''(''y'')}} is a convex Lipschitz function on {{math|[0,''r'']}} with Lipschitz constant 1 satisfying {{math|''k''(''t'') ≤ – ''t''}} and {{math|1=''k''(0) = 0}} and {{math|1=''k''(''r'') = –''r''}}. So ''{{math|k}}'' vanishes everywhere, since if {{math|0 < ''s'' < ''r'', ''k''(''s'') ≤ –''s''}} and {{math|{{!}}''k''(s){{!}} ≤ ''s''}}. Hence {{math|1=''h''(&delta;(''t'')) = ''h''(''y'') – ''t''}}. By uniqueness it follows that {{math|&delta;(''t'')}} is the closest point to {{math|''y''}} in {{math|''C''<sub>''t''</sub>}} and that it is the unique point minimising ''{{math|h}}'' in {{math|1={{overline|''B''}}(''y'',''t'')}}. Uniqueness implies that these geodesics segments coincide for arbitrary ''{{math|r}}'' and therefore that {{math|&delta;}} extends to a geodesic ray with the stated property.

*'''If {{math|1=h = h<sub>γ</sub>}}, then the geodesic ray {{math|&delta;}} starting at {{math|y}} satisfies''' <math>\sup d(\gamma(t),\delta(t)) < \infty</math>. '''If {{math|&delta;<sub>1</sub>}} is another ray starting at {{math|y}} with''' <math>\sup d(\delta(t),\delta_1(t)) < \infty</math> '''then {{math|1=&delta;<sub>1</sub> = &delta;}}.'''

To prove the first assertion, it is enough to check this for ''{{math|t}}'' sufficiently large. In that case {{math|γ(''t'')}} and {{math|&delta;(''t'' − ''h''(''y''))}} are the projections of ''{{math|x}}'' and ''{{math|y}}'' onto the closed convex set {{math|1=''h'' ≤ −''t''}}. Therefore, {{math|''d''(&gamma;(''t''),&delta;(''t'' − ''h''(''y''))) &le; ''d''(''x'',''y'')}}. Hence  {{math|''d''(&gamma;(''t''),&delta;(''t'')) &le; ''d''(&gamma;(''t''),&delta;(''t'' − ''h''(''y''))) + ''d''(&delta;(''t'' − ''h''(''y'')),&delta;(''t'')) &le; ''d''(''x'',''y'') + {{!}}''h''(''y''){{!}}}}. The second assertion follows because {{math|''d''(&delta;<sub>1</sub>(''t''),&delta;(''t''))}} is convex and bounded on {{math|[0,&infin;)}}, so, if it vanishes at {{math|1=''t'' = 0}}, must vanish everywhere.

*'''Suppose that {{math|h}} is a continuous convex function and for each {{math|y}} in {{math|X}} there is a unique geodesic ray {{math|&delta;}} such that {{math|1=&delta;(0) = y}} and {{math|&delta;}} cuts each closed convex set {{math|h &le; h(y) – r with r > 0}} at {{math|&delta;(r)}}, so that {{math|1=h(&delta;(t)) = h(y) – t}}; then {{math|h}} is a Busemann function. {{math|h − h<sub>&delta;</sub>}} is a constant function.'''<ref name="Bridson 1999 pages=271–272"/>

Let {{math|''C''<sub>''r''</sub>}} be the closed convex set of points ''{{math|z}}'' with {{math|''h''(''z'') &le; −''r''}}. Since ''{{math|X}}'' is a Hadamard space for every point ''{{math|y}}'' in ''{{math|X}}'' there is a unique closest point {{math|''P''<sub>''r''</sub>(''y'')}} to {{math|''y''}} in {{math|''C''<sub>''r''</sub>}}. It depends continuously on ''{{math|y}}'' and if ''{{math|y}}'' lies outside {{math|''C''<sub>''r''</sub>}}, then {{math|''P''<sub>''r''</sub>(''y'')}} lies on the hypersurface {{math|1=''h''(''z'') = − ''r''}}—the boundary ∂{{math|''C''<sub>''r''</sub>}} of {{math|''C''<sub>''r''</sub>}}—and {{math|''P''<sub>''r''</sub>(''y'')}} satisfies the inequality of convex optimisation. Let {{math|&delta;(''s'')}} be the geodesic ray starting at ''{{math|y}}''.

Fix ''{{math|x}}'' in ''{{math|X}}''. Let {{math|&gamma;(''s'')}} be the geodesic ray starting at ''{{math|x}}''. Let {{math|1=''g''(''z'') = ''h''<sub>γ</sub>(''z'')}}, the Busemann function for {{math|&gamma;}} with base point ''{{math|x}}''. In particular {{math|1=''g''(''x'') = 0}}. It suffices to show that {{math|1=''g'' = ''h'' – ''h''(''y'')1}}. Now take ''{{math|y}}'' with {{math|1=''h''(''x'') = ''h''(''y'')}} and let {{math|&delta;(''t'')}} be the geodesic ray starting at ''{{math|y}}'' corresponding to ''{{math|h}}''. Then

<math display=block>d(x,y) \ge d(\gamma(t),\delta(t)),\,\,\, d(x,\delta(t))^2 \ge d(x,\gamma(t))^2+d(\gamma(t),\delta(t))^2,\,\,\,d(y,\gamma(t))^2 \ge d(y,\delta(t))^2+d(\gamma(t),\delta(t))^2.</math>

On the other hand, for any four points ''{{math|a}}'', ''{{math|b}}'', ''{{math|c}}'', ''{{math|d}}'' in a Hadamard space, the following quadrilateral inequality of [[Yurii Reshetnyak|Reshetnyak]] holds:

<math display=block> |d(a,c)^2 + d(b,d)^2 - d(a,d)^2 - d(b,d)^2|\le 2 d(a,b)d(c,d).</math>

Setting {{math|1=''a'' = ''x''}}, {{math|1=''b'' = ''y''}}, {{math|1=''c'' = &gamma;(''t'')}}, {{math|1=''d'' =  &delta;(''t'')}}, it follows that

<math display=block>|d(y,\gamma(t))^2 - d(x,\delta(t))^2|\le 2d(x,y)^2,</math>

so that

<math display=block>|d(y,\gamma(t)) - d(x,\gamma(t))|\le  2 {d(x,y)^2\over d(y,\gamma(t)) + d(x,\gamma(t))}\le {d(x,y)^2\over t}.</math>

Hence {{math|1=''h''<sub>&gamma;</sub>(''y'') = 0}}. Similarly {{math|1=''h''<sub>&delta;</sub>(''x'') = 0}}. Hence {{math|1=''h''<sub>&gamma;</sub>(''y'') = 0}} on the level surface of ''{{math|h}}'' containing ''{{math|x}}''. Now for {{math|''t'' ≥ 0}} and ''{{math|z}}'' in ''{{math|X}}'', let {{math|1=α<sub>''t''</sub>(''z'') = γ<sub>1</sub>(''t'')}} the geodesic ray starting at ''{{math|z}}''. Then {{math|1=α<sub>''s'' + ''t''</sub> = α<sub>''s''</sub> ∘  α<sub>''t''</sub>}} and {{math|1=''h'' ∘ α<sub>''t''</sub> = ''h''  − ''t''}}. Moreover, by boundedness, {{math|1=''d''(α<sub>''t''</sub>(''u''),α<sub>''t''</sub>(''v'')) ≤ ''d''(''u'',''v'')}}. The flow {{math|α<sub>''s''</sub>}} can be used to transport this result to all the level surfaces of ''{{math|h}}''. For general {{math|''y''<sub>1</sub>}}, if {{math|''h''(''y''<sub>1</sub>) < ''h''(''x'')}}, take {{math|''s'' > 0}} such that {{math|1=''h''(α<sub>''s''</sub>(''x'')) = ''h''(''y''<sub>1</sub>)}} and set {{math|1=''x''<sub>1</sub> = α<sub>''s''</sub>(''x'')}}. Then {{math|1=''h''<sub>γ<sub>1</sub></sub>(''y''<sub>1</sub>) = 0}}, where {{math|1=γ<sub>1</sub>(''t'') = α<sub>''t''</sub>(''x''<sub>1</sub>) = γ(''s'' + ''t'')}}. But then {{math|1=''h''<sub>γ<sub>1</sub></sub> = ''h''<sub>γ</sub> – ''s''}}, so that {{math|1=''h''<sub>γ</sub>(''y''<sub>1</sub>) = ''s''}}. Hence {{math|1=''g''(''y''<sub>1</sub>) = ''s'' = ''h''((α<sub>''s''</sub>(''x'')) – ''h''(''x'') = ''h''(''y''<sub>1</sub>) – ''h''(''x'')}}, as required. Similarly if {{math|''h''(''y''<sub>1</sub>) > ''h''(''x'')}}, take
{{math|''s'' > 0}} such that {{math|1=''h''(α<sub>''s''</sub>(''y''<sub>1</sub>)) = ''h''(''x'')}}. Let {{math|1=''y'' = α<sub>''s''</sub>(''y''<sub>1</sub>)}}. Then {{math|1=''h''<sub>γ</sub>(''y'') = 0}}, so {{math|1=''h''<sub>γ</sub>(''y''<sub>1</sub>) = –''s''}}. Hence {{math|1=''g''(''y''<sub>1</sub>) = –''s'' =  ''h''(''y''<sub>1</sub>) – ''h''(''x'')}}, as required.

Finally there are necessary and sufficient conditions for two geodesics to define the same Busemann function up to constant:

*'''On a Hadamard space, the Busemann functions of two geodesic rays <math>\gamma_1</math> and <math>\gamma_2</math> differ by a constant if and only if <math>\sup_{t\ge 0} d(\gamma_1(t),\gamma_2(t)) < \infty</math>.'''<ref>{{harvnb|Dal'bo|Peigné|Sambusetti|2012|pages=94–96}}</ref>

Suppose firstly that {{math|γ}} and {{math|δ}} are two geodesic rays with Busemann functions differing by a constant. Shifting the argument of one of the geodesics by a constant, it may be assumed that {{math|1=''B''<sub>γ</sub> = ''B''<sub>δ</sub> = ''B''}}, say. Let ''{{math|C}}'' be the closed convex set on which {{math|''B''(''x'') ≤ −''r''}}. Then {{math|1=''B''(γ(''t'')) = ''B''<sub>γ</sub>(γ(''t'')) = −''t''}} and similarly {{math|1=''B''(δ(''t'')) = − ''t''}}. Then for {{math|''s'' ≤ ''r''}}, the points {{math|γ(''s'')}} and  {{math|δ(''s'')}} have closest points {{math|γ(''r'')}} and  {{math|δ(''r'')}} in ''{{math|C}}'', so that {{math|''d''(γ(''r''), δ(''r'')) ≤ ''d''(γ(''s''), δ(''s''))}}. Hence {{math|1=sup<sub>''t'' ≥ 0</sub> ''d''(γ(''t''), δ(''t'')) < ∞}}.

Now suppose that {{math|1=sup<sub>''t'' ≥ 0</sub> ''d''(γ<sub>1</sub>(''t''), γ<sub>2</sub>(''t'')) < ∞}}. Let {{math|δ<sub>''i''</sub>(''t'')}} be the geodesic ray starting at ''{{math|y}}'' associated with {{math|''h''<sub>γ<sub>''i''</sub></sub>}}. Then {{math|1=sup<sub>''t'' ≥ 0</sub> ''d''(γ<sub>''i''</sub>(''t''), δ<sub>''i''</sub>(''t'')) < ∞}}. Hence {{math|1=sup<sub>''t'' ≥ 0</sub> ''d''(δ<sub>1</sub>(''t''), δ<sub>2</sub>(''t'')) < ∞}}. Since {{math|δ<sub>1</sub>}} and {{math|δ<sub>2</sub>}} both start at ''{{math|y}}'', it follows that {{math|1=δ<sub>1</sub>(''t'') ≡ δ<sub>2</sub>(''t'')}}. By the previous result {{math|''h''<sub>γ<sub>''i''</sub></sub>}} and {{math|''h''<sub>δ<sub>''i''</sub></sub>}} differ by a constant; so {{math|''h''<sub>γ<sub>1</sub></sub>}} and {{math|''h''<sub>γ<sub>2</sub></sub>}} differ by a constant.

To summarise, the above results give the following characterisation of Busemann functions on a Hadamard space:<ref name="Bridson 1999 pages=271–272"/>

'''THEOREM.''' '''On a Hadamard space, the following conditions on a function {{math|f}} are equivalent:'''

*'''{{math|h}} is a Busemann function.'''
*'''{{math|h}} is a convex function, Lipschitz with constant {{math|1=1}} and {{math|h}} assumes its minimum on any closed ball centred on {{math|y}} with radius {{math|r}} at a unique point {{math|v}} on the boundary with {{math|1=h(v) = h(y) − r}}.'''
*'''{{math|h}} is a continuous convex function and for each {{math|y}} in {{math|X}} there is a unique geodesic ray {{math|δ}} such that {{math|1=δ(0) = y}} and, for any {{math|r > 0}}, the ray {{math|δ}} cuts each closed convex set {{math|h ≤ h(y) – r}} at {{math|δ(r)}}.'''

==Bordification of a Hadamard space==
In the previous section it was shown that if {{math|''X''}} is a Hadamard space and {{math|''x''<sub>0</sub>}} is a fixed point in {{math|''X''}} then the union of the space of Busemann functions vanishing at {{math|''x''<sub>0</sub>}} and the space of functions {{math|1=''h''<sub>''y''</sub>(''x'') = ''d''(''x'',''y'') − ''d''(''x''<sub>0</sub>,''y'')}} is closed under taking uniform limits on bounded sets. This result can be formalised in the notion of '''bordification''' of {{math|''X''}}.<ref>{{harvnb|Bridson|Haefliger|1999|pages=260–276}}</ref> In this topology, the points {{math|''x''<sub>''n''</sub>}} tend to a geodesic ray {{math|γ}} starting at {{math|''x''<sub>0</sub>}} if and only if {{math|''d''(''x''<sub>0</sub>,''x''<sub>''n''</sub>)}} tends to {{math|∞}} and for {{math|''t''  > 0}} arbitrarily large the sequence obtained by taking the point on each segment {{math|[''x''<sub>0</sub>,''x''<sub>''n''</sub>]}} at a distance {{math|''t''}} from {{math|''x''<sub>0</sub>}} tends to {{math|γ(''t'')}}.

If {{math|''X''}} is a metric space, Gromov's bordification can be defined as follows. Fix a point {{math|''x''<sub>0</sub>}} in {{math|''X''}} and let {{math|1=''X''<sub>''N''</sub> = {{overline|''B''}}(''x''<sub>0</sub>,''N'')}}. Let {{math|1=''Y'' = ''C''(''X'')}} be the space of Lipschitz continuous functions on {{math|''X''}}, i.e. those for which {{math|{{!}} ''f''(''x'')  – ''f''(''y'') {{!}} ≤ ''A'' ''d''(''x'',''y'')}} for some constant {{math|''A'' > 0}}. The space {{math|''Y''}} can be topologised by the seminorms {{math|1={{norm| ''f'' }}<sub>''N''</sub> = sup<sub>''X''<sub>''N''</sub></sub> {{!}} ''f'' {{!}}}}, the topology of uniform convergence on bounded sets. The seminorms are finite by the Lipschitz conditions. This is the topology induced by the natural map of {{math|''C''(''X'')}} into the direct product of the Banach spaces {{math|''C''<sub>''b''</sub>(''X''<sub>''N''</sub>)}} of continuous bounded functions on {{math|''X''<sub>''N''</sub>}}. It is give by the metric {{math|1=''D''(''f'',''g'') = Σ 2<sup>−''N''</sup> {{norm| ''f'' − ''g'' }}<sub>''N''</sub>(1 +{{norm| ''f'' − ''g'' }}<sub>''N''</sub>)<sup>−1</sup>}}.

The space {{math|''X''}} is embedded into {{math|''Y''}} by sending {{math|''x''}} to the function {{math|1=''f''<sub>''x''</sub>(''y'') = ''d''(''y'',''x'') – ''d''(''x''<sub>0</sub>,''x'')}}. Let {{math|{{overline|''X''}}}} be the closure of {{math|''X''}} in {{math|''Y''}}. Then {{math|{{overline|''X''}}}} is metrisable, since {{math|''Y''}} is, and contains {{math|''X''}} as an open subset; moreover bordifications arising from different choices of basepoint are naturally homeomorphic.  Let {{math|1=''h''(''x'') = (''d''(''x'',''x''<sub>0</sub>) + 1)<sup>−1</sup>}}. Then {{math|''h''}} lies in {{math|''C''<sub>0</sub>(''X'')}}. It is non-zero on {{math|''X''}} and vanishes only at {{math|∞}}. Hence it extends to a continuous function on {{math|{{overline|''X''}}}} with zero set {{math|{{overline|''X''}} \ ''X''}}. It follows that {{math|{{overline|''X''}} \ ''X''}} is closed in {{math|{{overline|''X''}}}}, as required. To check that {{math|1={{overline|''X''}} = {{overline|''X''}}(''x''<sub>0</sub>)}} is independent of the basepoint, it suffices to show that {{math|1=''k''(''x'') = ''d''(''x'',''x''<sub>0</sub>) − ''d''(''x'',''x''<sub>1</sub>)}} extends to a continuous function on {{math|{{overline|''X''}}}}. But {{math|1=''k''(''x'') = ''f''<sub>''x''</sub>(''x''<sub>1</sub>)}}, so, for {{math|''g''}} in {{math|{{overline|''X''}}}}, {{math|1=''k''(''g'') = ''g''(''x''<sub>1</sub>)}}. Hence the correspondence between the compactifications for {{math|''x''<sub>0</sub>}} and {{math|''x''<sub>1</sub>}} is given by sending {{math|''g''}} in {{math|{{overline|''X''}}(''x''<sub>0</sub>)}} to {{math|''g'' + ''g''(''x''<sub>1</sub>)1}} in {{math|{{overline|''X''}}(''x''<sub>1</sub>)}}.

When {{math|''X''}} is a Hadamard space, Gromov's ideal boundary {{math|1=∂''X'' = {{overline|''X''}} \ ''X''}} can be realised explicitly as "asymptotic limits" of geodesic rays using Busemann functions. If {{math|''x''<sub>''n''</sub>}} is an unbounded sequence in {{math|''X''}} with {{math|1=''h''<sub>''n''</sub>(''x'') = ''d''(''x'',''x''<sub>''n''</sub>) − ''d''(''x''<sub>''n''</sub>,''x''<sub>0</sub>)}} tending to {{math|''h''}} in {{math|''Y''}}, then {{math|''h''}} vanishes at {{math|''x''<sub>0</sub>}}, is convex, Lipschitz with Lipschitz constant {{math|1}} and has minimum {{math|''h''(''y'') − ''r''}} on any closed ball {{math|{{overline|''B''}}(''y'',''r'')}}. Hence {{math|''h''}} is a Busemann function {{math|''B''<sub>γ</sub>}} corresponding to a unique geodesic ray {{math|γ}} starting at {{math|''x''<sub>0</sub>}}.

On the other hand, {{math|''h''<sub>''n''</sub>}} tends to {{math|''B''<sub>γ</sub>}} uniformly on bounded sets if and only if {{math|''d''(''x''<sub>0</sub>,''x''<sub>''n''</sub>)}} tends to {{math|∞}} and for {{math|''t''  > 0}} arbitrarily large the sequence obtained by taking the point on each segment {{math|[''x''<sub>0</sub>,''x''<sub>''n''</sub>]}} at a distance {{math|''t''}} from {{math|''x''<sub>0</sub>}} tends to {{math|γ(''t'')}}. For {{math|''d''(''x''<sub>0</sub>,''x''<sub>''n''</sub>) ≥ ''t''}}, let {{math|''x''<sub>''n''</sub>(''t'')}} be the point in {{math|[''x''<sub>0</sub>,''x''<sub>''n''</sub>]}} with {{math|1=''d''(''x''<sub>0</sub>,''x''<sub>''n''</sub>(''t'')) = ''t''}}. Suppose first that {{math|''h''<sub>''n''</sub>}} tends to {{math|''B''<sub>γ</sub>}} uniformly on {{math|{{overline|''B''}}(''x''<sub>0</sub>,''R'')}}. Then for {{math|''t'' ≤ ''R''}},
{{math|1={{!}}''h''<sub>''n''</sub>(γ(''t'')) – ''B''<sub>γ</sub>(γ(''t'')){{!}}=''d''(γ(''t''),''x''<sub>''n''</sub>) – ''d''(''x''<sub>''n''</sub>,''x''<sub>0</sub>) + ''t''}}. This is a convex function. It vanishes as {{math|1=''t'' = 0}} and hence is increasing. So it is maximised at {{math|1=''t'' = ''R''}}. So for each {{math|''t''}}, {{math|{{!}}''d''(γ(''t''),''x''<sub>''n''</sub>) – ''d''(''x''<sub>''n''</sub>,''x''<sub>0</sub>) – ''t''{{!}}}} tends towards 0. Let {{math|1=''a'' = ''X''<sub>0</sub>}}, {{math|1=''b'' = γ(''t'')}} and {{math|1=''c'' = ''x''<sub>''n''</sub>}}. Then {{math|''d''(''c'',''a'') – ''d''(''c'',''b'')}} is close to {{math|''d''(''a'',''b'')}} with {{math|''d''(''c'',''a'')}} large. Hence in the Euclidean comparison triangle {{math|''CA'' - ''CB''}} is close to {{math|''AB''}} with {{math|''CA''}} large. So the angle at {{math|''A''}} is small. So the point {{math|''D''}} on {{math|''AC''}} at the same distance as {{math|''AB''}} lies close to {{math|''B''}}. Hence, by the first comparison theorem for geodesic triangles, {{math|''d''(''x''<sub>''n''</sub>(''t''),γ(''t''))}} is small. Conversely suppose that for fixed {{math|''t''}} and {{math|''n''}} sufficiently large {{math|''d''(''x''<sub>''n''</sub>(''t''),γ(''t''))}} tends to 0. Then from the above {{math|1=''F''<sub>''s''</sub>(''y'') = ''d''(''y'',γ(''s'')) – ''s''}} satisfies

<math display=block>|F_s(y) - B_\gamma(y)| \le {d(x_0,y)^2\over 2s},</math>

so it suffices show that on any bounded set {{math|1=''h''<sub>''n''</sub>(''y'') = ''d''(''y'',''x''<sub>''n''</sub>) – ''d''(''x''<sub>0</sub>,''x''<sub>''n''</sub>)}} is uniformly close to {{math|''F''<sub>''s''</sub>(''y'')}} for {{math|''n''}} sufficiently large.<ref>{{harvnb|Ballmann|1995|pages=27–30}}</ref>

For a fixed ball {{math|{{overline|''B''}}(''x''<sub>0</sub>,''R'')}}, fix {{math|''s''}} so that {{math|''R''<sup>2</sup>/''s'' ≤ ε}}.  The claim is then an immediate consequence of the inequality for geodesic segments in a Hadamard space, since

<math display=block>|d(y,x_n) -d(y,x_0) - d(y,x_n(s)) + s |\le  {d(x_0,y)^2\over s} \le \varepsilon.</math>

Hence, if {{math|''y''}} in {{math|{{overline|''B''}}(''x''<sub>0</sub>,''R'')}} and {{math|''n''}} is sufficiently large that {{math|''d''(''x''<sub>''n''</sub>(''s''),γ(''s'')) ≤ ε}}, then

<math display=block>|h_n(y)-B_\gamma(y)|=|d(y,x_n) -d(y,x_0) - B_\gamma(y)| \le |d(y,x_n) -d(y,x_0) - d(y,x_n(s)) + s | + d(x_n(s),\gamma(s)) + |F_s(y) - B_\gamma(y)| \le 3\varepsilon.</math>

== Busemann functions on a Hadamard manifold ==

Suppose that {{math|''x'', ''y''}}  are points in a Hadamard manifold and let {{math|''γ''(''s'')}} be the geodesic through {{math|''x''}} with {{math|1=''γ''(0) = ''y''}}. This geodesic cuts the boundary of the closed ball {{math|{{overline|''B''}}(''y'',''r'')}} at the two points {{math|γ(±''r'')}}. Thus if {{math|''d''(''x'',''y'') > ''r''}}, there are points {{math|''u'', ''v''}} with {{math|1=''d''(''y'',''u'') = ''d''(''y'',''v'') = ''r''}} such that {{math|1={{!}} ''d''(''x'',''u'') − ''d''(''x'',''v'') {{!}} = 2''r''}}. By continuity this condition persists for Busemann functions:

*'''If {{math|h}} is a Busemann function on a Hadamard manifold, then, given {{math|y}} in {{math|X}} and {{math|r > 0}}, there are unique points {{math|u}}, {{math|v}} with {{math|1=d(y,u) = d(y,v) = r}} such that {{math|1=h(u) = h(y) + r}} and  {{math|1=h(v) =  h(y) − r}}. For fixed {{math|''r'' > 0}}, the points {{math|''u''}} and {{math|''v''}} depend continuously on {{math|''y''}}.'''<ref name="BGS"/>

Taking a sequence {{math|''t''<sub>''n''</sub>}} tending to {{math|∞}} and {{math|1=''h''<sub>''n''</sub> = ''F''<sub>''t''<sub>''n''</sub></sub>}}, there are points {{math|''u''<sub>''n''</sub>}} and {{math|''v''<sub>''n''</sub>}} which satisfy these conditions for {{math|''h''<sub>''n''</sub>}} for {{math|''n''}} sufficiently large. Passing to a subsequence if necessary, it can be assumed that {{math|''u''<sub>''n''</sub>}} and {{math|''v''<sub>''n''</sub>}} tend to {{math|''u''}} and {{math|''v''}}. By continuity these points satisfy the conditions for {{math|''h''}}. To prove uniqueness, note that by compactness {{math|''h''}} assumes its maximum and minimum on {{math|{{overline|''B''}}(''y'',''r'')}}. The Lipschitz condition shows that the values of {{math|''h''}} there differ by at most {{math|2''r''}}. Hence {{math|''h''}} is minimized at {{math|''v''}} and maximized at {{math|''u''}}. On the other hand, {{math|1=''d''(''u'',''v'') = 2''r''}} and for {{math|''u''}} and {{math|''v''}} the points {{math|''v''}} and {{math|''u''}} are the unique points in {{math|{{overline|''B''}}(''y'',''r'')}} maximizing this distance. The Lipschitz condition on {{math|''h''}} then immediately implies {{math|''u''}} and {{math|''v''}} must be the unique points in {{math|{{overline|''B''}}(''y'',''r'')}} maximizing and minimizing {{math|''h''}}. Now suppose that {{math|''y''<sub>''n''</sub>}} tends to {{math|''y''}}. Then the corresponding points {{math|''u''<sub>''n''</sub>}} and {{math|''v''<sub>''n''</sub>}} lie in a closed ball so admit convergent subsequences. But by uniqueness of {{math|''u''}} and {{math|''v''}} any such subsequences must tend to {{math|''u''}} and {{math|''v''}}, so that {{math|''u''<sub>''n''</sub>}} and {{math|''v''<sub>''n''</sub>}} must tend to {{math|''u''}} and {{math|''v''}}, establishing continuity.

The above result holds more generally in a Hadamard space.<ref>{{harvnb|Bridson|Haefliger|1999|pages=271–272}}</ref>

*'''If {{math|h}} is a Busemann function on a Hadamard manifold, then {{math|h}} is continuously differentiable with {{math|1={{norm| dh(y) }} = 1}} for all {{math|y}}.'''<ref name="BGS"/>

From the previous properties of {{math|''h''}}, for each {{math|''y''}} there is a unique geodesic γ(''t'') parametrised by arclength with {{math|1=γ(0) = ''y''}} such that  {{math|1=''h'' ∘ γ(''t'') = ''h''(''y'') + ''t''}}. It has the property that it cuts {{math|∂''B''(''y'',''r'')}} at {{math|1=''t'' = ±''r''}}: in the previous notation {{math|1=γ(''r'') = ''u''}} and  {{math|1=γ(–''r'') = ''v''}}. The vector field {{math|''V''<sub>''h''</sub>}} defined by the unit vector <math>\dot{\gamma}(0)</math> at {{math|''y''}} is continuous, because {{math|''u''}} is a continuous function of {{math|''y''}} and the map  sending {{math|(''x'',''v'')}}  to {{math|(''x'',exp<sub>''x''</sub> ''v'')}} is a diffeomorphism from {{math|''TX''}} onto {{math|''X'' × ''X''}} by the [[Cartan-Hadamard theorem]]. Let {{math|δ(''s'')}} be another geodesic parametrised by arclength through {{math|''y''}} with {{math|1=δ(0) = ''y''}}. Then {{math|1=''dh'' ∘ δ  (0)/ ''ds'' =}} <math>(\dot{\delta}(0),\dot{\gamma}(0))</math>. Indeed, let {{math|1=''H''(''x'') = ''h''(''x'') − ''h''(''y'')}}, so that {{math|1=''H''(''y'') = 0}}. Then

<math display=block> |H(\delta(s)) - H(x)| \le d(\delta(s),x).</math>

Applying this with {{math|1=''x'' = ''u''}} and {{math|''v''}}, it follows that for {{math|''s'' > 0}}

<math display=block>(r- d(\delta(s),u))/s \le (h(\delta(s)) - h(y))/s \le (d(\delta(s),v) -r)/s.</math>

The outer terms tend to <math>(\dot{\delta}(0),\dot{\gamma}(0))</math> as {{math|''s''}} tends to 0, so the middle term has the same limit, as claimed. A similar argument applies for {{math|''s'' < 0}}.

The assertion on the outer terms follows from the first variation formula for arclength, but can be deduced directly as follows. Let 
<math>a=\dot{\delta}(0)</math> and <math>b=\dot{\gamma}(0)</math>, both unit vectors. Then for tangent vectors {{math|''p''}} and {{math|''q''}} at {{math|''y''}} in the unit ball<ref>In [[geodesic normal coordinates]], the metric {{math|1=''g''(''x'') = I + ε {{norm| ''x'' }}}}. By geodesic convexity, a geodesic from {{math|''p''}} to {{math|''q''}} lies in the ball of radius {{math|1=''r'' = max {{norm| ''p'' }}, {{norm| ''q'' }}}}. The straight line segment gives an upper estimate for {{math|''d''(''p'',''q'')}} of the stated form. To obtain a similar lower estimate, observe that if {{math|''c''(''t'')}} is a smooth path from {{math|''p''}} to {{math|''q''}}, then {{math|''L''(''c'') ≥  (1 − ε ''r'') ⋅ ∫ {{norm| ''c'' }} ''dt'' ≥  (1 − ε ''r'') ⋅ {{norm| ''p'' − ''q'' }}}}. (Note that these inequalities can be improved using the sharper estimate {{math|1=''g''(''x'') = I + ε {{norm|''x''}}<sup>2</sup>}}).</ref>

<math display=block>d(\exp_y p,\exp_y q) = \|p-q\| + \varepsilon \max \|p\|^2, \|q\|^2</math>
 
with {{math|ε}} uniformly bounded. Let {{math|1=''s'' = ''t''<sup>3</sup>}} and {{math|1=''r'' = ''t''<sup>2</sup>}}. Then

<math display=block>(d(\delta(s),v) -r)/s= (d(\exp_y(t^3 a),\exp_y(-t^2 b)) - t^2)/t^3 = (\|t^3a +t^2b\| - t^2)/t^3 + \varepsilon |t|= (\|ta +b\| -1)/t + \varepsilon |t|.</math>

The right hand side here tends to {{math|(''a'',''b'')}} as {{math|''t''}} tends to 0 since

<math display=block> {d\over dt} \|b + t a\||_{t=0 }= {1\over 2} \,{d\over dt} \|b + t a\|^2|_{t=0} = (a,b).</math>

The same method works for the other terms.

Hence it follows that {{math|''h''}} is a {{math|C<sup>1</sup>}} function with {{math|''dh''}} dual to the vector field {{math|''V''<sub>''h''</sub>}}, so that {{math|1={{norm| ''dh''(''y'') }} = 1}}.  The vector field {{math|''V''<sub>''h''</sub>}} is thus the [[gradient vector field]] for {{math|''h''}}. The geodesics through any point are the flow lines for the flow {{math|α<sub>''t''</sub>}} for {{math|''V''<sub>''h''</sub>}}, so that {{math|α<sub>''t''</sub>}} is the [[gradient flow]] for {{math|''h''}}.

'''THEOREM.''' '''On a Hadamard manifold {{math|X}} the following conditions on a continuous function {{math|h}} are equivalent:'''<ref name="BGS"/>
# '''{{math|h}} is a Busemann function.'''
# '''{{math|h}} is a convex, Lipschitz function with constant 1, and for each {{math|y}} in {{math|X}} there are points {{math|u<sub>±</sub>}} at a distance {{math|r}} from {{math|y}} such that {{math|1=h(u<sub>±</sub>) = h(y) ± r}}.'''
# '''{{math|h}} is a convex {{math|C<sup>1</sup>}} function with {{math|1={{norm| dh(x) }} ≡ 1}}.'''

It has already been proved that (1) implies (2).

The arguments above show ''[[mutatis mutandi]]'' that (2) implies (3).

It therefore remains to show that (3) implies (1). Fix {{math|''x''}} in {{math|''X''}}. Let {{math|α<sub>''t''</sub>}} be the gradient flow for {{math|''h''}}. It follows that {{math|1=''h'' ∘ α<sub>''t''</sub> (''x'') = ''h''(''x'') + ''t''}} and that {{math|1=γ(''t'') = α<sub>''t''</sub>(''x'')}} is a geodesic through {{math|''x''}} parametrised by arclength with {{math|1=γ(0) = ''x''}}. Indeed, if {{math|''s'' < ''t''}}, then

<math display=block>|s-t|=|h(\alpha_s(x))-h(\alpha_t(x))| \le d(\alpha_s(x),\alpha_t(x))\le \int_{s}^t \|d\alpha_\tau(x)/d\tau\|  \, d\tau =\int_s^t \|dh(\alpha_\tau(x))\| \, d\tau = |s-t|,</math>

so that {{math|1=''d''(γ(''s''),γ(''t'')) = {{!}} ''s'' − ''t'' {{!}}}}. Let {{math|1=''g''(''y'') = ''h''<sub>γ</sub>(''y'')}}, the Busemann function for {{math|γ}} with base point {{math|''x''}}. In particular {{math|1=''g''(''x'') = 0}}. To prove (1), it suffices to show that {{math|1=''g'' = ''h'' – ''h''(''x'')1}}.

Let {{math|''C''(−''r'')}} be the closed convex set of points {{math|''z''}} with {{math|''h''(''z'') ≤ −''r''}}. Since {{math|''X''}} is a Hadamard space for every point {{math|''y''}} in {{math|''X''}} there is a unique closest point {{math|''P''<sub>''r''</sub>(''y'')}} to {{math|''y''}} in {{math|''C''(-''r'')}}. It depends continuously on {{math|''y''}} and if {{math|''y''}} lies outside {{math|''C''(-''r'')}}, then {{math|''P''<sub>''r''</sub>(''y'')}} lies on the hypersurface {{math|1=''h''(''z'') = − ''r''}}—the boundary {{math|∂''C''(–''r'')}} of {{math|''C''(–''r'')}}—and the geodesic from {{math|''y''}} to {{math|''P''<sub>''r''</sub>(''y'')}} is orthogonal to {{math|∂''C''(–''r'')}}. In this case the geodesic is just {{math|α<sub>''t''</sub>(''y'')}}. Indeed, the fact that {{math|α<sub>''t''</sub>}} is the gradient flow of {{math|''h''}} and the conditions {{math|{{norm| ''dh''(''y'') }} ≡ 1}} imply that the flow lines  {{math|α<sub>''t''</sub>(''y'')}} are geodesics parametrised by arclength and cut the level curves of {{math|''h''}} orthogonally. Taking {{math|''y''}} with {{math|1=''h''(''y'') = ''h''(''x'')}} and  {{math|''t'' > 0}},

<math display=block>d(x,y) \ge d(\alpha_t(x),\alpha_t(y)),\,\,\, d(x,\alpha_t(y))^2 \ge d(x,\alpha_t(x))^2+d(\alpha_t(x),\alpha_t(y))^2,\,\,\,d(y,\alpha_t(x))^2 \ge d(y,\alpha_t(y))^2+d(\alpha_t(x),\alpha_t(y))^2.</math>

On the other hand, for any four points {{math|''a''}}, {{math|''b''}}, {{math|''c''}}, {{math|''d''}} in a Hadamard space, the following quadrilateral inequality of [[Yurii Reshetnyak|Reshetnyak]] holds:

<math display=block> |d(a,c)^2 + d(b,d)^2 - d(a,d)^2 - d(b,d)^2|\le 2 d(a,b)d(c,d).</math>

Setting {{math|1=''a'' = ''x''}}, {{math|1=''b'' = ''y''}}, {{math|1=''c'' = α<sub>''t''</sub>(''x'')}}, {{math|1=''d'' =  α<sub>''t''</sub>(''y'')}}, it follows that

<math display=block>|d(y,\alpha_t(x))^2 - d(x,\alpha_t(x))^2|\le 2d(x,y)^2,</math>

so that

<math display=block>|d(y,\alpha_t(x)) - d(x,\alpha_t(x))|\le  2 {d(x,y)^2\over d(y,\alpha_t(x)) + d(x,\alpha_t(x))}\le {d(x,y)^2\over t}.</math>

Hence {{math|1=''h''<sub>γ</sub>(''y'') = 0}} on the level surface of {{math|''h''}} containing {{math|''x''}}. The flow {{math|α<sub>''s''</sub>}} can be used to transport this result to all the level surfaces of {{math|''h''}}. For general {{math|''y''<sub>1</sub>}} take {{math|''s''}} such that {{math|1=''h''(α<sub>''s''</sub>(''x'')) = ''h''(''y''<sub>1</sub>)}} and set {{math|1=''x''<sub>1</sub> = α<sub>''s''</sub>(''x'')}}. Then {{math|1=''h''<sub>γ<sub>1</sub></sub>(''y''<sub>1</sub>) = 0}}, where {{math|1=γ<sub>1</sub>(''t'') = α<sub>''t''</sub>(''x''<sub>1</sub>) = γ(''s'' + ''t'')}}. But then {{math|1=''h''<sub>γ<sub>1</sub></sub> = ''h''<sub>γ</sub> – ''s''}}, so that {{math|1=''h''<sub>γ</sub>(''y''<sub>1</sub>) = ''s''}}. Hence {{math|1=''g''(''y''<sub>1</sub>) = ''s'' = ''h''((α<sub>''s''</sub>(''x'')) – ''h''(''x'') = ''h''(''y''<sub>1</sub>) – ''h''(''x'')}}, as required.

Note that this argument could be shortened using the fact that two Busemann functions {{math|''h''<sub>γ</sub>}} and {{math|''h''<sub>δ</sub>}} differ by a constant if and only if the corresponding geodesic rays satisfy {{math|sup<sub>''t'' ≥ 0</sub> ''d''(γ(''t''),δ(''t'')) < ∞}}. Indeed, all the geodesics defined by the flow {{math|α<sub>''t''</sub>}} satisfy the latter condition, so differ by constants. Since along any of these geodesics {{math|''h''}} is linear with derivative 1, {{math|''h''}} must differ from these Busemann functions by constants.

==Compactification of a proper Hadamard space==
{{harvtxt|Eberlein|O'Neill|1973}} defined a compactification of a [[Hadamard manifold]] {{math|''X''}} which uses Busemann functions. Their construction, which can be extended more generally to proper (i.e. locally compact) [[Hadamard space]]s, gives an explicit geometric realisation of a compactification defined by Gromov—by adding an "ideal boundary"—for the more general class of [[proper metric space]]s {{math|''X''}}, those for which every closed ball is compact. Note that, since any Cauchy sequence is contained in a closed ball, any proper metric space is automatically complete.<ref>Note that a metric space {{math|''X''}}  which is complete and locally compact need not be proper, for example {{math|'''R'''}} with the metric {{math|1=''d''(''x'',''y'') = {{!}} ''x'' – ''y'' {{!}} /(1 + {{!}} ''x'' – ''y'' {{!}})}}. On the other hand, by the [[Hopf–Rinow theorem]] for metric spaces, if {{math|''X''}} is complete, locally compact and geodesic—every two points {{math|''x''}} and {{math|''y''}} are joined by a geodesic parametrised by arclength—then {{math|''X''}} is proper (see {{harvnb|Bridson|Haefliger|1999|pages=35–36}}). Indeed if not, there is a point {{math|''x''}} in {{math|''X''}} and a closed ball {{math|1=''K'' = {{overline|''B''}}(''x'',''r'')}} maximal subject to being compact; then, since by hypothesis {{math|{{overline|''B''}}(''x'',''R'')}} is non-compact for each {{math|''R'' > ''r''}}, a diagonal argument shows that there is a sequence {{math|(''x''<sub>''n''</sub>)}} with {{math|''d''(''x'',''x''<sub>''n''</sub>)}} decreasing to {{math|''r''}} but with no convergent subsequence; on the other hand taking {{math|''y''<sub>''n''</sub>}} on a geodesic joining {{math|''x''}} and {{math|''x''<sub>''n''</sub>}}, with {{math|1=''d''(''x'',''y''<sub>''n''</sub>) = ''r''}}, compactness of  {{math|''K''}} implies {{math|(''y''<sub>''n''</sub>)}}, and hence {{math|(''x''<sub>''n''</sub>)}}, has a convergent subsequence, a contradiction.</ref> The ideal boundary is a special case of the ideal boundary for a metric space. In the case of Hadamard spaces, this agrees with the space of geodesic rays emanating from any fixed point described using Busemann functions in the bordification of the space.

If {{math|''X''}} is a proper metric space, Gromov's compactification can be defined as follows. Fix a point {{math|''x''<sub>0</sub>}} in {{math|''X''}} and let {{math|1=''X''<sub>''N''</sub> = {{overline|''B''}}(''x''<sub>0</sub>,''N'')}}. Let {{math|1=''Y'' = ''C''(''X'')}} be the space of Lipschitz continuous functions on {{math|''X''}}, .e. those for which {{math|{{!}} ''f''(''x'')  – ''f''(''y'') {{!}} ≤ ''A'' ''d''(''x'',''y'')}} for some constant {{math|''A'' > 0}}. The space {{math|''Y''}} can be topologised by the seminorms {{math|1={{norm| ''f'' }}<sub>''N''</sub> = sup<sub>''X''<sub>''N''</sub></sub> {{!}} ''f'' {{!}}}}, the topology of uniform convergence on compacta. This is the topology induced by the natural map of ''C''(''X'') into the direct product of the Banach spaces {{math|''C''(''X''<sub>''N''</sub>)}}. It is give by the metric {{math|1=''D''(''f'',''g'') = Σ 2<sup>−''N''</sup> {{norm| ''f'' − ''g'' }} <sub>''N''</sub>(1 + {{norm| ''f'' − ''g'' }} <sub>''N''</sub>)<sup>−1</sup>}}.

The space {{math|''X''}} is embedded into {{math|''Y''}} by sending {{math|''x''}} to the function {{math|1=''f''<sub>''x''</sub>(''y'') = ''d''(''y'',''x'') – ''d''(''x''<sub>0</sub>,''x'')}}. Let {{math|{{overline|''X''}}}} be the closure of {{math|''X''}} in {{math|''Y''}}. Then {{math|{{overline|''X''}}}} is compact (metrisable) and contains {{math|''X''}} as an open subset; moreover compactifications arising from different choices of basepoint are naturally homeomorphic. Compactness follows from the [[Arzelà–Ascoli theorem]] since the image in {{math|''C''(''X''<sub>''N''</sub>)}} is [[equicontinuous]] and uniformly bounded in norm by {{math|''N''}}. Let {{math|''x''<sub>''n''</sub>}} be a sequence in {{math|''X'' ⊂ {{overline|''X''}}}} tending to {{math|''y''}} in {{math|{{overline|''X''}} \ ''X''}}. Then all but finitely many terms must lie outside {{math|''X''<sub>''N''</sub>}} since {{math|''X''<sub>''N''</sub>}} is compact, so that any subsequence would converge to a point in {{math|''X''<sub>''N''</sub>}}; so the sequence {{math|''x''<sub>''n''</sub>}} must be unbounded in {{math|''X''}}. Let {{math|1=''h''(''x'') = (''d''(''x'',''x''<sub>0</sub>) + 1)<sup>−1</sup>}}. Then {{math|''h''}} lies in {{math|''C''<sub>0</sub>(''X'')}}. It is non-zero on {{math|''X''}} and vanishes only at {{math|∞}}. Hence it extends to a continuous function on {{math|{{overline|''X''}}}} with zero set  {{math|{{overline|''X''}} \ ''X''}}. It follows that {{math|{{overline|''X''}} \ ''X''}} is closed in {{math|{{overline|''X''}}}}, as required. To check that the compactification  {{math|1={{overline|''X''}} = {{overline|''X''}}(''x''<sub>0</sub>)}} is independent of the basepoint, it suffices to show that {{math|1=''k''(''x'') = ''d''(''x'',''x''<sub>0</sub>) − ''d''(''x'',''x''<sub>1</sub>)}} extends to a continuous function on {{math|{{overline|''X''}}}}. But {{math|1=''k''(''x'') = ''f''<sub>''x''</sub>(''x''<sub>1</sub>)}}, so, for {{math|''g''}} in {{math|{{overline|''X''}}}}, {{math|1=''k''(''g'') = ''g''(''x''<sub>1</sub>)}}. Hence the correspondence between the compactifications for {{math|''x''<sub>0</sub>}} and {{math|''x''<sub>1</sub>}} is given by sending {{math|''g''}} in {{math|{{overline|''X''}}(''x''<sub>0</sub>)}} to {{math|''g'' + ''g''(''x''<sub>1</sub>)1}} in {{math|{{overline|''X''}}(''x''<sub>1</sub>)}}.

When {{math|''X''}} is a Hadamard manifold (or more generally a proper Hadamard space), Gromov's ideal boundary {{math|1=∂''X'' = {{overline|''X''}} \ ''X''}} can be realised explicitly as "asymptotic limits" of geodesics by using Busemann functions. Fixing a base point {{math|''x''<sub>0</sub>}}, there is a unique geodesic {{math|γ(''t'')}} parametrised by arclength such that {{math|1=γ(0) = ''x''<sub>0</sub>}} and <math>\dot{\gamma}(0)</math> is a given unit vector. If {{math|''B''<sub>γ</sub>}} is the corresponding Busemann function, then
{{math|''B''<sub>γ</sub>}} lies in {{math|∂''X''(''x''<sub>0</sub>)}} and induces a homeomorphism of the unit {{math|(''n'' − 1)}}-sphere onto {{math|∂''X''(''x''<sub>0</sub>)}}, sending <math>\dot{\gamma}(0)</math> to {{math|''B''<sub>γ</sub>}}.

== Quasigeodesics in the Poincaré disk, CAT(-1) and hyperbolic spaces==
===Morse–Mostow lemma===
In the case of spaces of negative curvature, such as the Poincaré disk, CAT(-1) and hyperbolic spaces, there is a metric structure on their Gromov boundary. This structure is preserved by the group of quasi-isometries which carry geodesics rays to quasigeodesic rays. Quasigeodesics were first studied for negatively curved surfaces—in particular the hyperbolic upper halfplane and unit disk—by [[Marston Morse|Morse]] and generalised to negatively curved [[symmetric space]]s by [[George Mostow|Mostow]], for his work on the [[Mostow rigidity theorem|rigidity of discrete groups]]. The basic result is the '''Morse–Mostow lemma''' on the stability of geodesics.<ref name="Bourdon 1995">{{harvnb|Bourdon|1995}}</ref><ref name="Buyalo 2007">{{harvnb|Buyalo|Schroeder|2007}}</ref><ref>{{harvnb|Mostow|1973}}</ref><ref>{{harvnb|Roe|2003}}</ref>

By definition a '''[[quasigeodesic]]''' Γ defined on an interval {{math|[''a'',''b'']}}  with {{math|−∞  ≤ ''a'' < ''b'' ≤ ∞}} is a map {{math|Γ(''t'')}} into a metric space, not necessarily continuous, for which there are constants {{math|λ ≥ 1}} and {{math|ε > 0}} such that for all {{math|''s''}} and {{math|''t''}}:

<math display=block>\lambda^{-1} |s-t| - \varepsilon \le d(\Gamma(s),\Gamma(t)) \le \lambda |s-t| + \varepsilon.</math>

The following result is essentially due to [[Marston Morse]] (1924).

'''Morse's lemma on stability of geodesics.''' In the hyperbolic disk there is a constant {{math|''R''}} depending on {{math|λ}} and {{math|ε}} such that any quasigeodesic segment {{math|Γ}} defined on a finite interval {{math|[''a'',''b'']}} is within a [[Hausdorff distance]] {{math|''R''}} of the geodesic segment {{math|[Γ(''a''),Γ(''b'')]}}.<ref>{{harvnb|Buyalo|Schroeder|2007|pages=1–6}}</ref><ref>{{harvnb|Bridson|Haefliger|1999|pages=399–405}}</ref>

===Classical proof for Poincaré disk===

The classical proof of Morse's lemma for the Poincaré unit disk or upper halfplane proceeds more directly by using orthogonal projection onto the geodesic segment.<ref>{{harvnb|Kapovich|2001|pages=51–52}}</ref><ref>{{harvnb|Morse|1924}}</ref><ref>{{harvnb|Ratcliffe|2006|pages=580–599}}</ref>

*It can be assumed that Γ satisfies the stronger "pseudo-geodesic" condition:<ref>{{harvnb|Kapovich|2001|page=51}}</ref>

<math display=block>\lambda^{-1}|s-t| -\varepsilon \le d(\Gamma(s),\Gamma(t)) \le \lambda |s-t|.</math>

{{math|Γ}} can be replaced by a continuous piecewise geodesic curve Δ with the same endpoints lying at a finite Hausdorff distance from {{math|Γ}} less than {{math|1=''c'' = (2λ<sup>2</sup> + 1)ε}}: break up the interval on which {{math|Γ}} is defined into equal subintervals of length {{math|2λε}} and take the geodesics between the images under {{math|Γ}} of the endpoints of the subintervals. Since {{math|Δ}} is piecewise geodesic, {{math|Δ}} is Lipschitz continuous with constant {{math|λ<sub>1</sub>}}, {{math|''d''(Δ(''s''),Δ(''t'')) ≤ λ<sub>1</sub> {{!}} ''s'' – ''t'' {{!}}}}, where  {{math|λ<sub>1</sub> ≤ λ + ε}}. The lower bound is automatic at the endpoints of intervals. By construction the other values differ from these by a uniformly bounded depending only on {{math|λ}} and {{math|ε}}; the lower bound inequality holds by increasing ε by adding on twice this uniform bound.

*If {{math|γ}} is a piecewise smooth curve segment lying outside an {{math|''s''}}-neighbourhood of a geodesic line and {{math|''P''}} is the orthogonal projection onto the geodesic line then:<ref>{{harvnb|Ratcliffe|2006|page=583}}, Lemma 4</ref>

<math display=block>\ell(P \circ \gamma) \le {\ell(\gamma)\over \cosh s}.</math>

Applying an isometry in the upper half plane, it may be assumed that the geodesic line is the positive imaginary axis in which case the orthogonal projection onto it is given by {{math|1=''P''(''z'') = ''i'' {{!}} ''z'' {{!}}}} and {{math|1={{!}} ''z'' {{!}} / Im z = cosh ''d''(''z'',''Pz'')}}. Hence the hypothesis implies {{math|{{!}} γ(''t'') {{!}} ≥ cosh(''s'') Im γ(''t'')}}, so that

<math display=block> \ell(P\circ\gamma) =  \int_a^b {|d\gamma|\over |\gamma|} \le \int_a^b {|d\gamma|\over \cosh(s)\, \Im \gamma}={\ell(\gamma)\over \cosh(s)}.</math>

*There is a constant {{math|''h'' > 0}} depending only on {{math|λ}} and {{math|ε}} such that  {{math|Γ[''a'',''b'']}} lies within an {{math|''h''}}-neighbourhood of the geodesic segment {{math|[Γ(''a''),Γ(''b'')]}}.<ref>{{harvnb|Ratcliffe|2006|pages=584–586}}, Lemmas 5–6</ref>

Let {{math|γ(''t'')}} be the geodesic line containing the geodesic segment {{math|[Γ(''a''),Γ(''b'')]}}. Then there is a constant {{math|''h'' > 0}} depending only on {{math|λ}} and {{math|ε}} such that  {{math|''h''}}-neighbourhood {{math|Γ[''a'',''b'']}} lies within an {{math|''h''}}-neighbourhood of {{math|γ('''R''')}}. Indeed for any {{math|''s'' > 0}}, the subset of {{math|[''a'',''b'']}} for which {{math|Γ(''t'')}} lies outside the closure of the {{math|''s''}}-neighbourhood of {{math|γ('''R''')}} is open, so a countable union of open intervals {{math|(''c'',''d'')}}. Then

<math display=block>\ell(\Gamma|_{[c,d]}) \le s_1 \equiv \lambda^2 (2s +\varepsilon) \left(1 -{\lambda^2\over \cosh(s)}\right),</math>

:since the left hand side is less than or equal to {{math|λ {{!}} ''c'' – ''d'' {{!}}}} and

<math display=block>{|c-d|\over \lambda} -\varepsilon \le d(\Gamma(c),\Gamma(d)) \le 2s + d(P \circ \Gamma|_{[c,d]}) \le 2s + {\lambda |c-d|\over\cosh(s)}.</math>

:Hence every point lies at a distance less than or equal to {{math|''s'' + ''s''<sub>1</sub>}} of {{math|γ('''R''')}}. To deduce the assertion, note that the subset of {{math|[''a'',''b'']}} for which {{math|Γ(''t'')}} lies outside the closure of the {{math|''s''}}-neighbourhood of {{math|[Γ(''a''),Γ(''b'')] ⊂ γ('''R''')}} is open, so a union of intervals {{math|(''c'',''d'')}} with {{math|Γ(''c'')}} and {{math|Γ(''d'')}} both at a distance {{math|''s'' + ''s''<sub>1</sub>}} from either {{math|Γ(''a'')}} or {{math|Γ(''b'')}}. Then

<math display=block>\ell(\Gamma|_{[c,d]}) \le s_2 \equiv \lambda^2(2(s + s_1) + \varepsilon),</math>

:since

<math display=block>{|c-d|\over \lambda} -\varepsilon \le d(\Gamma(c),\Gamma(d)) \le 2(s+s_1).</math>

:Hence the assertion follows taking any {{math|''h''}} greater than {{math|''s'' +''s''<sub>1</sub> + ''s''<sub>2</sub>}}.

*There is a constant {{math|''h'' > 0}} depending only on {{math|λ}} and {{math|ε}} such that the geodesic segment {{math|[Γ(''a''),Γ(''b'')]}} lies within an {{math|''h''}}-neighbourhood of {{math|Γ[''a'',''b'']}}.<ref>{{harvnb|Kapovich|2001|page=52}}</ref>

Every point of {{math|Γ[''a'',''b'']}} lies within a distance {{math|''h''}} of {{math|[Γ(''a''),Γ(''b'')]}}. Thus orthogonal projection {{math|''P''}} carries each point of {{math|Γ[''a'',''b'']}} onto a point in the closed convex set {{math|[Γ(''a''),Γ(''b'')]}} at a distance less than {{math|''h''}}. Since {{math|''P''}} is continuous and {{math|Γ[''a'',''b'']}} connected, the map {{math|''P''}} must be onto since the image contains the endpoints of {{math|[Γ(''a''),Γ(''b'')]}}. But then every point of {{math|[Γ(''a''),Γ(''b'')]}} is within a distance {{math|''h''}} of a point of {{math|Γ[''a'',''b'']}}.

===Gromov's proof for Poincaré disk ===
The generalisation of Morse's lemma to CAT(-1) spaces is often referred to as the Morse–Mostow lemma and can be proved by a straightforward generalisation of the classical proof. There is also a generalisation for the more general class of [[hyperbolic metric space]]s due to Gromov. Gromov's proof is given below for the Poincaré unit disk; the properties of hyperbolic metric spaces are developed in the course of the proof, so that it applies [[mutatis mutandi]] to CAT(-1) or hyperbolic metric spaces.<ref name="Bourdon 1995"/><ref name="Buyalo 2007"/>

Since this is a large-scale phenomenon, it is enough to check that any maps {{math|Δ}} from {{math|{{mset|0, 1, 2, ..., ''N''}}}} for any {{math|''N'' > 0}} to the disk satisfying the inequalities is within a Hausdorff distance {{math|''R''<sub>1</sub>}} of the geodesic segment {{math|[Δ(0),Δ(''N'')]}}. For then translating it may be assumed without loss of generality {{math|Γ}} is defined on {{math|[0,''r'']}} with {{math|''r'' > 1}} and then, taking {{math|1=''N'' = [''r'']}} (the integer part of {{math|''r''}}), the result can be applied to {{math|Δ}} defined by {{math|1=Δ(''i'') = Γ(''i'')}}. The Hausdorff distance between the images of {{math|Γ}} and {{math|Δ}} is evidently bounded by a constant {{math|''R''<sub>2</sub>}} depending only on {{math|λ}} and {{math|ε}}.

:Now the [[incircle]] of a geodesic triangle has diameter less than {{math|δ}} where {{math|1=δ = 2 log 3}}; indeed it is strictly maximised by that of an ideal triangle where it equals {{math|2 log 3}}. In particular, since the incircle breaks the triangle breaks the triangle into three isosceles triangles with the third side opposite the vertex of the original triangle having length less than {{math|δ}}, it follows that every side of a geodesic triangle is contained in a {{math|δ}}-neighbourhood of the other two sides. A simple induction argument shows that a geodesic polygon with {{math|2<sup>''k''</sup> + 2}} vertices for {{math|''k'' ≥ 0}} has each side within a {{math|1=(''k'' + 1)δ}} neighbourhood of the other sides (such a polygon is made by combining two geodesic polygons with {{math|1=2<sup>''k''−1</sup> + 1}} sides along a common side). Hence if {{math|1=''M'' ≤ 2<sup>''k''</sup> + 2}}, the same estimate holds for a polygon with {{math|''M''}} sides.

:For {{math|1=''y''<sub>''i''</sub> = Δ(''i'')}} let {{math|1=''f''(''x'') = min ''d''(''x'',''y''<sub>''i''</sub>)}}, the largest radius for a closed ball centred on {{math|''x''}} which contains no {{math|''y''<sub>''i''</sub>}} in its interior. This is a continuous function non-zero on {{math|[Δ(0),Δ(''N'')]}} so attains its maximum {{math|''h''}} at some point {{math|''x''}} in this segment. Then {{math|[Δ(0),Δ(''N'')]}} lies within an {{math|''h''<sub>1</sub>}}-neighbourhood of the image of {{math|Δ}} for any {{math|''h''<sub>1</sub> > ''h''}}. It therefore suffices to find an upper bound for {{math|''h''}} independent of {{math|''N''}}.

:Choose {{math|''y''}} and {{math|''z''}} in the segment {{math|[Δ(0),Δ(''N'')]}} before and after {{math|''x''}} with {{math|1=''d''(''x'',''y'') = 2''h''}} and {{math|1=''d''(''x'',''z'') = 2''h''}} (or an endpoint if it within a distance of less than {{math|2''h''}} from {{math|''x''}}). Then there are {{math|''i'', ''j''}} with {{math|''d''(''y'',Δ(''i''))}}, {{math|''d''(''z'',Δ(''j'')) ≤ ''h''}}. Hence {{math|''d''(Δ(''i''),Δ(''j'')) ≤ 6''h''}}, so that {{math|{{!}} ''i'' − ''j'' {{!}} ≤ 6λ''h'' + λε}}. By the triangle inequality all points on the segments {{math|[''y'',Δ(''i'')]}} and {{math|[''z'',Δ(''j'')]}} are at a distance {{math|≥ ''h''}} from {{math|''x''}}. Thus there is a finite sequence of points starting at {{math|''y''}} and ending at {{math|''z''}}, lying first on the segment {{math|[''y'',Δ(''i'')]}}, then proceeding through the points {{math|Δ(''i''), Δ(''i''+1), ..., Δ(''j'')}}, before taking the segment {{math|[Δ(''j''),''z'']}}. The successive points {{math|Δ(''i''), Δ(''i''+1), ..., Δ(''j'')}} are separated by a distance no greater than {{math|λ + ε}} and successive points on the geodesic segments can also be chosen to satisfy this condition. The minimum number {{math|''K''}} of points in such a sequence satisfies {{math|''K'' ≤ {{!}} ''i'' - ''j'' {{!}} + 3 + 2(λ + ε)<sup>–1</sup>''h''}}. These points form a geodesic polygon, with {{math|[''y'',''z'']}} as one of the sides. Take {{math|1=''L'' = [''h''/δ]}}, so that the {{math|(''L'' − 1)δ}}-neighbourhood of {{math|[''y'',''z'']}} does not contain all the other sides of the polygon. Hence, from the result above, it follows that {{math|''K'' > 2<sup>''L'' − 1</sup> + 2}}. Hence

<math display=block>3 + 2(\lambda+\varepsilon)^{-1}h + 6\lambda h + \varepsilon> 2^{h/\delta} + 2.</math>

:This inequality implies that {{math|''h''}} is uniformly bounded, independently of {{math|''N''}}, as claimed.

:If all points {{math|Δ(''i'')}} lie within {{math|''h''<sub>1</sub>}} of the {{math|[Δ(0),Δ(''N'')]}}, the result follows. Otherwise the points which do not fall into maximal subsets {{math|1=''S'' = {{mset|''r'', ..., ''s''}}}} with {{math|''r'' < ''s''}}. Thus points in {{math|[Δ(0),Δ(''N'')]}} have a point  {{math|Δ(''i'')}} with {{math|''i''}} in the complement of {{math|''S''}} within a distance of {{math|''h''<sub>1</sub>}}. But the complement of {{math|1=''S'' = ''S''<sub>1</sub>  ∐  ''S''<sub>2</sub>}}, a disjoint union with {{math|1=''S''<sub>1</sub> = {{mset|0, ..., ''r'' − 1}}}} and {{math|1=''S''<sub>2</sub> = {{mset|''s'' + 1, ..., ''N''}}}}. [[connected space|Connectivity]] of {{math|[Δ(0),Δ(''N'')]}} implies there is a point {{math|''x''}} in the segment which is within a distance {{math|''h''<sub>1</sub>}} of points {{math|Δ(''i'')}} and {{math|Δ(''j'')}} with {{math|''i'' < ''r''}} and {{math|''j'' > ''s''}}. But then {{math|''d''(Δ(''i''),Δ(''j'')) < 2 ''h''<sub>1</sub>}}, so {{math|{{!}} ''i'' − ''j'' {{!}} ≤ 2λ''h''<sub>1</sub> + λε}}. Hence the points {{math|Δ(''k'')}} for {{math|''k''}} in {{math|''S''}} lie within a distance from  {{math|[Δ(0),Δ(''N'')]}} of less than {{math|1=''h''<sub>1</sub> + λ {{!}} ''i'' − ''j'' {{!}} + ε ≤ ''h''<sub>1</sub> + λ (2λ''h''<sub>1</sub> + λε) + ε ≡ ''h''<sub>2</sub>}}.

===Extension to quasigeodesic rays and lines===
Recall that in a Hadamard space if {{math|[''a''<sub>1</sub>,''b''<sub>1</sub>]}} and {{math|[''a''<sub>2</sub>,''b''<sub>2</sub>]}} are two geodesic segments and the intermediate points {{math|''c''<sub>1</sub>(''t'')}} and {{math|''c''<sub>2</sub>(''t'')}} divide them in the ratio {{math|''t'':(1 – ''t'')}}, then {{math|''d''(''c''<sub>1</sub>(''t''),''c''<sub>2</sub>(''t''))}} is a convex function of {{math|''t''}}. In particular if {{math|Γ<sub>1</sub>(''t'')}} and {{math|Γ<sub>2</sub>(''t'')}} are geodesic segments of unit speed defined on {{math|[0,''R'']}} starting at the same point then

<math display=block>d(\Gamma_1(t),\Gamma_2(t)) \le {t\over R} d(\Gamma_1(R),\Gamma_2(R)).</math>

In particular this implies the following:

*'''In a CAT(–1) space {{math|''X''}}, there is a constant {{math|h > 0}} depending only on {{math|λ}} and {{math|ε}} such that any quasi-geodesic ray is within a bounded Hausdorff distance {{math|h}} of a geodesic ray. A similar result holds for quasigeodesic and geodesic lines.'''

If {{math|Γ(''t'')}} is a geodesic say with constant {{math|λ}} and {{math|ε}}, let {{math|Γ<sub>''N''</sub>(''t'')}} be the unit speed geodesic for the segment {{math|[Γ(0),Γ(''N'')]}}. The estimate above shows that for fixed {{math|''R'' > 0}} and {{math|''N''}} sufficiently large, {{math|(Γ<sub>''N''</sub>)}} is a Cauchy sequence in {{math|''C''([0,''R''],''X'')}} with the uniform metric. Thus {{math|Γ<sub>''N''</sub>}} tends to a geodesic ray {{math|γ}} uniformly on compacta the bound on the Hausdorff distances between {{math|Γ}} and the segments {{math|Γ<sub>''N''</sub>}} applies also to the limiting geodesic {{math|γ}}. The assertion for quasigeodesic lines follows by taking {{math|Γ<sub>''N''</sub>}} corresponding to the geodesic segment {{math|[Γ(–''N''),Γ(''N'')]}}.

== Efremovich–Tikhomirova theorem ==
Before discussing CAT(-1) spaces, this section will describe the '''Efremovich–Tikhomirova theorem''' for the unit disk {{math|''D''}} with the Poincaré metric. It asserts that quasi-isometries of {{math|''D''}} extend to quasi-Möbius homeomorphisms of the unit disk with the Euclidean metric. The theorem forms the prototype for the more general theory of CAT(-1) spaces. Their original theorem was proved in a slightly less general and less precise form in {{harvtxt|Efremovich|Tikhomirova|1964}} and applied to bi-Lipschitz homeomorphisms of the unit disk for the Poincaré metric;<ref>Bi-Lipschitz homeomorphisms are those for which they and their inverses are Lipschitz continuous</ref> earlier, in the posthumous paper {{harvtxt|Mori|1957}}, the Japanese mathematician Akira Mori had proved a related result within [[Teichmüller theory]] assuring that every [[quasiconformal mapping|quasiconformal homeomorphism]] of the disk is [[Hölder continuous]] and therefore extends continuously to a homeomorphism of the unit circle (it is known that this extension is quasi-Möbius).<ref>See:
*{{harvnb|Ahlfors|1966}}
*{{harvnb|Lehto|1987}}
</ref>

===Extension of quasi-isometries to boundary ===
If {{math|''X''}} is the Poincaré unit disk, or more generally a CAT(-1) space, the Morse lemma on stability of quasigeodesics implies that every quasi-isometry of {{math|''X''}} extends uniquely to the boundary. By definition two self-mappings {{math|''f'', ''g''}} of {{math|''X''}} are quasi-equivalent if {{math|sup<sub>''X''</sub> ''d''(''f''(''x''),''g''(''x'')) < ∞}}, so that corresponding points are at a uniformly bounded distance of each other. A quasi-isometry {{math|''f''<sub>1</sub>}} of {{math|''X''}} is a self-mapping of {{math|''X''}}, not necessarily continuous, which has a quasi-inverse {{math|''f''<sub>2</sub>}} such that {{math|''f''<sub>1</sub> ∘ ''f''<sub>2</sub>}} and {{math|''f''<sub>2</sub> ∘ ''f''<sub>1</sub>}} are quasi-equivalent to the appropriate identity maps and such that there are constants {{math|λ ≥ 1}} and {{math|ε > 0}} such that for all {{math|''x'', ''y''}} in {{math|''X''}} and both mappings

<math display=block> \lambda^{-1} d(x,y) - \varepsilon \le d(f_k(x),f_k(y)) \le \lambda d(x,y) +\varepsilon.</math>

Note that quasi-inverses are unique up to quasi-equivalence; that equivalent definition could be given using possibly different right and left-quasi inverses, but they would necessarily be quasi-equivalent; that quasi-isometries are closed under composition which up to quasi-equivalence depends only the quasi-equivalence classes; and that, modulo quasi-equivalence, the quasi-isometries form a group.<ref>See:
*{{harvnb|Bourdon|1995}}
*{{harvnb|Bridson|Haefliger|1999}}
*{{harvnb|Roe|2003}}
*{{harvnb|Buyalo|Schroeder|2007}}
*{{harvnb|Bourdon|2009}}
</ref>

Fixing a point {{math|''x''}} in {{math|''X''}}, given a geodesic ray {{math|γ}} starting at {{math|''x''}}, the image {{math|''f'' ∘ γ}} under a quasi-isometry {{math|''f''}} is a quasi-geodesic ray. By the Morse-Mostow lemma it is within a bounded distance of a unique geodesic ray {{math|δ}} starting at {{math|''x''}}. This defines a mapping {{math|∂''f''}} on the boundary {{math|∂''X''}} of {{math|''X''}}, independent of the quasi-equivalence class of {{math|''f''}}, such that {{math|1=∂(''f'' ∘ ''g'') = ∂''f'' ∘ ∂''g''}}. Thus there is a homomorphism of the group of quasi-isometries into the group of self-mappings of {{math|∂''X''}}.

To check that {{math|∂''f''}} is continuous, note that if {{math|γ<sub>1</sub>}} and {{math|γ<sub>2</sub>}} are geodesic rays that are uniformly close on {{math|[0,''R'']}}, within a distance {{math|η}}, then {{math|''f'' ∘ γ<sub>1</sub>}} and {{math|''f'' ∘ γ<sub>2</sub>}} lie within a distance {{math|λη + ε}}  on {{math|[0,''R'']}}, so that {{math|δ<sub>1</sub>}} and {{math|δ<sub>2</sub>}} lie within a distance {{math|λη + ε + 2''h''(λ,ε)}}; hence on a smaller interval {{math|[0,''r'']}}, {{math|δ<sub>1</sub>}} and {{math|δ<sub>2</sub>}} lie within a distance {{math|(''r''/''R'')⋅[λη + ε + 2''h''(λ,ε)]}} by convexity.<ref>{{harvnb|Bridson|Haefliger|1999|pages=430–431}}</ref>

On CAT(-1) spaces, a finer version of continuity asserts that {{math|∂''f''}} is a quasi-Möbius mapping with respect to a natural class of metric on {{math|∂''X''}}, the "visual metrics" generalising the Euclidean metric on the unit circle and its transforms under the Möbius group. These visual metrics can be defined in terms of Busemann functions.<ref>See:
*{{harvnb|Bourdon|1995}}
*{{harvnb|Buyalo|Schroeder|2007}}
*{{harvnb|Bourdon|2009}}
</ref>

In the case of the unit disk, Teichmüller theory implies that the homomorphism carries quasiconformal homeomorphisms of the disk onto the group of quasi-Möbius homeomorphisms of the circle (using for example the Ahlfors–Beurling or [[Douady–Earle extension]]): it follows that the homomorphism from the quasi-isometry group into the quasi-Möbius group is surjective.

In the other direction, it is straightforward to prove that the homomorphism is injective.<ref>{{harvnb|Roe|2003|page=113}}</ref> Suppose that {{math|''f''}} is a quasi-isometry of the unit disk such that {{math|∂''f''}} is the identity. The assumption and the Morse lemma implies that if {{math|γ('''R''')}} is a geodesic line, then {{math|''f''(γ('''R'''))}} lies in an {{math|''h''}}-neighbourhood of {{math|γ('''R''')}}. Now take a second geodesic line {{math|δ}} such that {{math|δ}} and {{math|γ}} intersect orthogonally at a given point in {{math|''a''}}. Then {{math|''f''(''a'')}} lies in the intersection of {{math|''h''}}-neighbourhoods of {{math|δ}} and {{math|γ}}. Applying a Möbius transformation, it can be assumed that {{math|''a''}} is at the origin of the unit disk and the geodesics are the real and imaginary axes. By convexity, the {{math|''h''}}-neighbourhoods of these axes intersect in a {{math|3''h''}}-neighbourhood of the origin: if {{math|''z''}} lies in both neighbourhoods, let {{math|''x''}} and {{math|''y''}} be the orthogonal projections of {{math|''z''}} onto the {{math|''x''}}- and {{math|''y''}}-axes; then {{math|1=''d''(''z'',''x'') ≤ ''h''}} so taking projections onto the {{math|''y''}}-axis, {{math|1=''d''(0,''y'') ≤ ''h''}}; hence  {{math|1=''d''(''z'',0) ≤ ''d''(''z'',''y'') + ''d''(''y'',0) ≤ 2''h''}}.  Hence {{math|1=''d''(''a'',''f''(''a'')) ≤ 2''h''}}, so that {{math|''f''}} is quasi-equivalent to the identity, as claimed.

===Cross ratio and distance between non-intersecting geodesic lines===
Given two distinct points {{math|''z'', ''w''}} on the unit circle or real axis there is a unique hyperbolic geodesic {{math|[''z'',''w'']}} joining them. It is given by the circle (or straight line) which cuts the unit circle unit circle or real axis orthogonally at those two points. Given four distinct points {{math|''a'', ''b'', ''c'', ''d''}} in the extended complex plane their [[cross ratio]] is defined by

<math display=block>(a,b;c,d) ={(a-c)(b-d)\over (a-d)(b-c)}.</math>

If {{math|''g''}} is a complex [[Möbius transformation]] then it leaves the cross ratio invariant: {{math|1=(''g''(''a''),''g''(''b'');''g''(''c''),''g''(''d'')) = (''a'',''b'':''c'',''d'')}}. Since the Möbius group acts simply transitively on triples of points, the cross ratio can alternatively be described as the complex number {{math|''z''}} in {{math|'''C'''\{{mset|0,1}}}} such that {{math|1=''g''(''a'') = 0, ''g''(''b'') = 1, ''g''(''c'') = λ, ''g''(''d'') = ∞}} for a Möbius transformation {{math|''g''}}.

Since {{math|''a''}}, {{math|''b''}}, {{math|''c''}} and {{math|''d''}} all appear in the numerator defining the cross ratio, to understand the behaviour of the cross ratio under permutations of {{math|''a''}}, {{math|''b''}}, {{math|''c''}} and {{math|''d''}}, it suffices to consider permutations that fix {{math|''d''}}, so only permute {{math|''a''}}, {{math|''b''}} and {{math|''c''}}. The cross ratio transforms according to the [[anharmonic group]] of order 6 generated by the Möbius transformations sending {{math|λ}} to {{math|1 – λ}} and {{math|λ<sup>−1</sup>}}. The other three transformations send {{math|λ}} to {{math|1 – λ<sup>−1</sup>}}, to {{math|λ(λ – 1)<sup>−1</sup>}} and to {{math|(1 – λ)<sup>−1</sup>}}.<ref>{{harvnb|Beardon|1983|pages=75–78}} Note that there is a natural homomorphism of {{math|''S''<sub>4</sub>}} onto {{math|''S''<sub>3</sub>}}, acting by conjugation on {{math|(''a'',''b'')(''c'',''d''), (''a'',''c'')(''b'',''d'')}} and {{math|(''a'',''d'')(''b'',''c'')}}. Indeed these permutations together with the identity form a normal Abelian subgroup equal to its own centraliser: the action of {{math|''S''<sub>4</sub>}} by conjugation on the non-trivial elements yields the homomorphism onto {{math|''S''<sub>3</sub>}}.</ref>

Now let {{math|''a'', ''b'', ''c'', ''d''}} be points on the unit circle or real axis in that order. Then the geodesics {{math|[''a'',''b'']}} and {{math|[''c'',''d'']}} do not intersect and the distance between these geodesics is well defined: there is a unique geodesic line cutting these two geodesics orthogonally and the distance is given by the length of the geodesic segment between them. It is evidently invariant under real Möbius transformations. To compare the cross ratio and the distance between geodesics, Möbius invariance allows the calculation to be reduced to a symmetric configuration. For {{math|0 < ''r'' < ''R''}}, take {{math|1=''a'' = –''R'', ''b'' = −''r'', ''c'' = ''r'', ''d'' = ''R''}}. Then
{{math|1=λ = (''a'',''b'';''c'',''d'') = (''R'' + ''r'')<sup>2</sup>/4''rR'' = (''t'' + 1)<sup>2</sup>/4''t''}} where {{math|1=''t'' = ''R''/''r'' > 1}}. On the other hand, the geodesics {{math|[''a'',''d'']}} and {{math|[''b'',''c'']}} are the semicircles in the upper half plane of radius {{math|''r''}} and {{math|''R''}}. The geodesic which cuts them orthogonally is the positive imaginary axis, so the distance between them is the hyperbolic distance between {{math|''ir''}} and {{math|''iR''}}, {{math|1=''d''(''ir'',''iR'') = log ''R''/''r'' = log ''t''}}. Let {{math|1=''s'' = log ''t''}}, then {{math|1=λ = cosh<sup>2</sup>(''s''/2)}}, so that there is a constant {{math|''C'' > 0}} such that, if {{math|(''a'',''b'';''c'',''d'') > 1}}, then

<math display=block>d([a,d];[b,c]) - C \le \log (a,b;c,d) \le d([a,d];[b,c]) + C,</math>

since {{math|1=log[cosh(''x'')/exp''x'')] = log (1 + exp(–2''x''))/2}} is bounded above and below in {{math|''x'' ≥ 0}}. Note that {{math|''a'', ''b'', ''c'', ''d''}} are in order around the unit circle if and only if {{math|(''a'',''b'';''c'',''d'') > 1}}.

A more general and precise geometric interpretation of the cross ratio can be given using projections of ideal points on to a geodesic line; it does not depend on the order of the points on the circle and therefore whether or not geodesic lines intersect.<ref>See: 
*{{harvnb|Paulin|1996}}
*{{harvnb|Buyalo|Schroeder|2007}}</ref>
*'''If {{math|p}} and {{math|q}} are the feet of the perpendiculars from {{math|c}} and {{math|d}} to the geodesic line {{math|ab}}, then {{math|1=d(p,q) = {{!}} log {{!}} (a,b;c,d) {{!}}{{!}}}}.'''

Since both sides are invariant under Möbius transformations, it suffices to check this in the case that {{math|1=''a'' = 0}}, {{math|1=''b'' = ∞}}, {{math|1=''c'' = ''x''}} and {{math|1=''d'' = 1}}. In this case the geodesic line is the positive imaginary axis, right hand side equals {{math|{{!}} log {{!}} ''x'' {{!}}{{!}}}}, {{math|1=''p'' = {{!}} ''x'' {{!}} ''i''}} and {{math|1=''q'' = ''i''}}. So the left hand side equals {{math|{{!}} log {{!}} ''x'' {{!}}{{!}}}}. Note that {{math|''p''}} and {{math|''q''}} are also the points where the incircles of the ideal triangles {{math|''abc''}} and {{math|''abd''}} touch {{math|''ab''}}.

===Proof of theorem===
A homeomorphism {{math|''F''}} of the circle is [[quasisymmetric map|''quasisymmetric'']] if there are constants {{math|''a'', ''b'' > 0}} such that

<math display=block>\displaystyle{{|F(z_1)-F(z_2)|\over |F(z_1)-F(z_3)|} \le a {|z_1-z_2|^b\over |z_1-z_3|^b}.}</math>

It is ''quasi-Möbius'' is there are constants {{math|''c'', ''d'' > 0}} such that

<math display=block>\displaystyle{|(F(z_1),F(z_2);F(z_3),F(z_4))| \le c |(z_1,z_2;z_3,z_4)|^d,}</math>

where

<math display=block>\displaystyle{ (z_1,z_2;z_3,z_4)={(z_1-z_3)(z_2-z_4)\over(z_2-z_3)(z_1-z_4)}}</math>

denotes the [[cross-ratio]].

It is immediate that quasisymmetric and quasi-Möbius homeomorphisms are closed under the operations of inversion and composition.

If {{math|''F''}} is quasisymmetric then it is also quasi-Möbius, with {{math|1=''c'' = ''a''<sup>2</sup>}} and {{math|1=''d'' = ''b''}}: this follows by multiplying the first inequality for {{math|(''z''<sub>1</sub>,''z''<sub>3</sub>,''z''<sub>4</sub>)}} and {{math|(''z''<sub>2</sub>,''z''<sub>4</sub>,''z''<sub>3</sub>)}}. Conversely any quasi-Möbius homeomorphism {{math|''F''}} is quasisymmetric. To see this, it can be first be checked that {{math|''F''}} (and hence {{math|''F''<sup>−1</sup>}}) is [[Hölder continuous]]. Let {{math|''S''}} be the set of cube roots of unity, so that if {{math|''a'' ≠ ''b''}} in {{math|''S''}}, then {{math|1={{!}} ''a'' − ''b'' {{!}} = 2 sin {{pi}}/3 = {{radic|3}}}}. To prove a Hölder estimate, it can be assumed that {{math|''x'' – ''y''}} is uniformly small. Then both {{math|''x''}} and {{math|''y''}} are greater than a fixed distance away from {{math|''a'', ''b''}} in {{math|''S''}} with {{math|''a'' ≠ ''b''}}, so the estimate follows by applying the quasi-Möbius inequality to {{math|''x'', ''a'', ''y'', ''b''}}. To verify that {{math|''F''}} is quasisymmetric, it suffices to find a uniform upper bound for {{math|{{!}} ''F''(''x'') − ''F''(''y'') {{!}} / {{!}} ''F''(''x'') − ''F''(''z'') {{!}}}} in the case of a triple with {{math|1={{!}} ''x'' − ''z'' {{!}} = {{!}} ''x'' − ''y'' {{!}}}}, uniformly small. In this case there is a point {{math|''w''}} at a distance greater than 1 from {{math|''x''}}, {{math|''y''}} and {{math|''z''}}. Applying the quasi-Möbius inequality to {{math|''x''}}, {{math|''w''}}, {{math|''y''}} and {{math|''z''}} yields the required upper bound. To summarise:

* '''A homeomorphism of the circle is quasi-Möbius if and only if it is quasisymmetric. In this case it and its inverse are Hölder continuous. The quasi-Möbius homeomorphisms form a group under composition.'''<ref>{{harvnb|Väisälä|1984}}</ref>

To prove the theorem it suffices to prove that if {{math|1=''F'' = ∂''f''}} then there are constants {{math|''A'', ''B'' > 0}} such that for {{math|''a'', ''b'', ''c'', ''d''}} distinct points on the unit circle<ref>{{harvnb|Bourdon|2009}}</ref>

<math display=block>\displaystyle{|(F(a),F(b);F(c),F(d))| \le A |(a,b;c,d)|^B.}</math>

It has already been checked that {{math|''F''}} (and is inverse) are continuous. Composing {{math|''f''}}, and hence {{math|''F''}}, with complex conjugation if necessary, it can further be assumed that {{math|''F''}} preserves the orientation of the circle. In this case, if {{math|''a'',''b'', ''c'',''d''}} are in order on the circle, so too are there images under {{math|''F''}}; hence both {{math|(''a'',''b'';''c'',''d'')}} and {{math|(''F''(''a''),''F''(''b'');''F''(''c''),''F''(''d''))}} are real and greater than one. In this case

<math display=block>(F(a),F(b);F(c),F(d)) \le A (a,b;c,d)^B.</math>

To prove this, it suffices to show that {{math|1=log (''F''(''a''),''F''(''b'');''F''(''c''),''F''(''d'')) ≤ ''B'' log (''a'',''b'';''c'',''d'') + ''C''}}. From the previous section it suffices show {{math|1=''d''([''F''(''a''),''F''(''b'')],[''F''(''c''),''F''(''d'')]) ≤ ''P'' ''d''([''a'',''b''],[''c'',''d'']) + ''Q''}}. This follows from the fact that the images under {{math|''f''}} of {{math|[''a'',''b'']}} and {{math|[''c'',''d'']}} lie within {{math|''h''}}-neighbourhoods of {{math|[''F''(''a''),''F''(''b'')]}} and {{math|[''F''(''c''),''F''(''d'')]}}; the minimal distance can be estimated using the quasi-isometry constants for {{math|''f''}} applied to the points on {{math|[''a'',''b'']}} and {{math|[''c'',''d'']}}
realising {{math|''d''([''a'',''b''],[''c'',''d''])}}.

Adjusting {{math|''A''}} and {{math|''B''}} if necessary, the inequality above applies also to {{math|''F''<sup>−1</sup>}}. Replacing {{math|''a''}}, {{math|''b''}}, {{math|''c''}} and {{math|''d''}} by their images under {{math|''F''}}, it follows that

<math display=block>A^{-1} |(a,b;c,d)|^{-B} \le |(F(a),F(b);F(c),F(d))| \le A |(a,b;c,d)|^B</math>

if {{math|''a''}}, {{math|''b''}}, {{math|''c''}} and {{math|''d''}} are in order on the unit circle. Hence the same inequalities are valid for the three cyclic of the quadruple {{math|''a'', ''b'', ''c'', ''d''}}. If {{math|''a''}} and {{math|''b''}} are switched then the cross ratios are sent to their inverses, so lie between 0 and 1; similarly if {{math|''c''}} and {{math|''d''}} are switched. If both pairs are switched, the cross ratio remains unaltered.  Hence the inequalities are also valid in this case. Finally if {{math|''b''}} and {{math|''c''}} are interchanged, the cross ratio changes from {{math|λ}} to {{math|1=λ<sup>–1</sup>
(λ – 1) = 1 – λ<sup>–1</sup>}}, which lies between 0 and 1. Hence again the same inequalities are valid. It is easy to check that using these transformations the inequalities are valid for all possible permutations of {{math|''a''}}, {{math|''b''}}, {{math|''c''}} and {{math|''d''}}, so that {{math|''F''}} and its inverse are quasi-Möbius homeomorphisms.

== Busemann functions and visual metrics for CAT(-1) spaces==

Busemann functions can be used to determine special visual metrics on the class of CAT(-1) spaces. These are complete geodesic metric spaces in which the distances between points on the boundary of a geodesic triangle are less than or equal to the comparison triangle in the hyperbolic upper half plane or equivalently the unit disk with the Poincaré metric. In the case of the unit disk the chordal metric can be recovered directly using Busemann functions {{math|''B''<sub>γ</sub>}} and the special theory for the disk generalises completely to any proper CAT(-1) space {{math|''X''}}. The hyperbolic upper half plane is a CAT(0) space, as lengths in a hyperbolic geodesic triangle are less than lengths in the Euclidean comparison triangle: in particular a CAT(-1) space is a CAT(0) space, so the theory of Busemann functions and the Gromov boundary applies. From the theory of the hyperbolic disk, it follows in particular that every geodesic ray in a CAT(-1) space extends to a geodesic line and given two points of the boundary there is a unique geodesic {{math|γ}} such that has these points as the limits {{math|γ(±∞)}}. The theory applies equally well to any CAT({{math|−κ}}) space with {{math|κ > 0}} since these arise by scaling the metric on a CAT(-1) space by {{math|κ<sup>−1/2</sup>}}. On the hyperbolic unit disk {{math|''D''}} quasi-isometries of {{math|''D''}} induce quasi-Möbius homeomorphisms of the boundary in a functorial way. There is a more general theory of Gromov hyperbolic spaces, a similar statement holds, but with less precise control on the homeomorphisms of the boundary.<ref name="Bourdon 1995"/><ref name="Buyalo 2007"/>

===Example: Poincaré disk===

==Applications in percolation theory==
{{see|percolation theory}}
More recently Busemann functions have been used by [[probability theory|probabilists]] to study asymptotic properties in models of [[first passage percolation|first-passage percolation]]<ref>{{harvnb|Hoffman|2005}}</ref><ref>{{harvnb|Damron|Hanson|2014}}</ref> and directed last-passage percolation.<ref>{{harvnb|Georgiou|Rassoul-Agha|Seppäläinen|2016}}</ref>

== Notes ==
{{reflist|40em}}

==References==
*{{citation|last=Ahlfors|first=Lars V.|author-link=Lars Ahlfors|title=Lectures on quasiconformal mappings|publisher=Van Nostrand|year=1966}}
*{{citation|last1=Ballmann|first1=Werner|last2= Gromov|first2=Mikhael|last3= Schroeder|first3= Viktor|title=Manifolds of nonpositive curvature|series=Progress in Mathematics|volume= 61|publisher=[[Birkhäuser]]|year= 1985|isbn= 0-8176-3181-X }}
*{{citation|last=Ballmann|first=Werner|title=Lectures on spaces of nonpositive curvature|series=DMV Seminar|volume= 25|publisher=Birkhäuser|year= 1995|isbn= 3-7643-5242-6}}
*{{citation|first=Alan F.|last= Beardon|title=The Geometry of Discrete Groups|year=1983|publisher=Springer-Verlag|isbn=0-387-90788-2}}
*{{citation|last=Bourdon|first= Marc|language=fr|title=Structure conforme au bord et flot géodésique d'un CAT(−1)-espace|journal=Enseign. Math.|volume= 41|year=1995|pages=63–102}}
*{{citation|last=Bourdon|first=Marc|chapter=Quasi-conformal geometry and Mostow rigidity|title=Géométries à courbure négative ou nulle, groupes discrets et rigidités|pages= 201–212|series=Sémin. Congr.|volume=18|publisher= Soc. Math. France|year= 2009}}
*{{citation | last1=Bridson | first1=Martin R. | last2=Haefliger | first2=André | title=Metric spaces of non-positive curvature | year=1999 | publisher=Springer}}
*{{citation|last=Busemann|first=Herbert|title=The geometry of geodesics|publisher= Academic Press|year=1955}}
*{{citation|last1=Buyalo|first1=Sergei|last2=Schroeder|first2=Viktor|title=Elements of asymptotic geometry|series=EMS Monographs in Mathematics|publisher=[[European Mathematical Society]]|year= 2007|isbn= 978-3-03719-036-4}}
*{{citation|last1=Dal'bo|first1= Françoise|last2=Peigné|first2=Marc|last3= Sambusetti|first3= Andrea|title=On the horoboundary and the geometry of rays of negatively curved manifolds|journal=Pacific J. Math.|volume= 259|year=2012|pages=55–100|url=http://msp.org/pjm/2012/259-1/pjm-v259-n1-p03-s.pdf|doi=10.2140/pjm.2012.259.55|arxiv=1010.6028|s2cid= 7309250}}, Appendix
*{{citation|last1=Damron|first1=Michael|last2=Hanson|first2= Jack|title=Busemann functions and infinite geodesics in two-dimensional first-passage percolation|journal=Comm. Math. Phys.|volume= 325|year=2014|issue=3|pages=917–963|bibcode=2014CMaPh.325..917D|doi=10.1007/s00220-013-1875-y|arxiv=1209.3036|s2cid=119589291}}
*{{citation|last1=Eberlein|first1=P.|last2=O'Neill|first2=B.|title=Visibility manifolds|journal=Pacific J. Math.|volume= 46|year=1973|pages=45–109|doi=10.2140/pjm.1973.46.45|doi-access=free}}
*{{citation|first1=V. A.|last1=Efremovich|author-link=Efremovich|first2= E. S.|last2=Tikhomirova|title=Equimorphisms of hyperbolic spaces|journal=Izv. Akad. Nauk SSSR Ser. Mat.|volume= 28|year=1964|pages=1139–1144|url=http://www.mathnet.ru/php/getFT.phtml?jrnid=im&paperid=3046&what=fullt&option_lang=eng|language=ru}}
*{{citation|last1=Georgiou|first1= Nicos|last2= Rassoul-Agha|first2= Firas|last3= Seppäläinen|first3=Timo|title=Variational formulas and cocycle solutions for directed polymer and percolation models|journal=Comm. Math. Phys. |volume=346 |year=2016|issue= 2|pages= 741–779|bibcode=2016CMaPh.346..741G|doi=10.1007/s00220-016-2613-z|arxiv=1311.3016|s2cid= 5887311}}
*{{citation|last=Hoffman|first=Christopher|title=Coexistence for Richardson type competing spatial growth models|journal=Ann. Appl. Probab.|volume= 15 |year=2005|pages= 739–747|doi=10.1214/105051604000000729|arxiv=math/0405377|s2cid=15113728}}
*{{citation|last=Kapovich|first=Michael|title=Hyperbolic manifolds and discrete groups|series=Progress in Mathematics|volume=183|publisher=[[Birkhäuser]] |year=2001|isbn= 0-8176-3904-7 }}
*{{citation|last=Lehto|first= Olli|title=Univalent functions and Teichmüller spaces|series=Graduate Texts in Mathematics|volume= 109|publisher= Springer-Verlag|year= 1987|isbn=0-387-96310-3}}
*{{citation |last=Lurie |first=J. |url=http://www.math.harvard.edu/~lurie/papers/hadamard.pdf |archive-url=https://web.archive.org/web/20060708074249/http://www.math.harvard.edu/~lurie/papers/hadamard.pdf |url-status=dead |archive-date=July 8, 2006 |title=Notes on the theory of Hadamard spaces |year=2010 |publisher=[[Harvard University]]}}
*{{citation|last=Mori|first= Akira|title=On quasi-conformality and pseudo-analyticity|journal=Trans. Amer. Math. Soc.|volume= 84 |year=1957|pages= 56–77|url=https://www.ams.org/journals/tran/1957-084-01/S0002-9947-1957-0083024-5/S0002-9947-1957-0083024-5.pdf|doi=10.1090/s0002-9947-1957-0083024-5|doi-access=free}}
*{{citation|first=H. M.|last= Morse|author-link=Marston Morse|title= A fundamental class of geodesics on any closed surface of genus greater than one|journal=Trans. Amer. Math. Soc.|volume=26|year=1924|pages= 25–60|url=https://www.ams.org/journals/tran/1924-026-01/S0002-9947-1924-1501263-9/S0002-9947-1924-1501263-9.pdf|doi=10.1090/s0002-9947-1924-1501263-9|doi-access=free}}
*{{cite book |last1=Mostow |first1=G. Daniel |title=Strong Rigidity of Locally Symmetric Spaces |series=Annals of Mathematics Studies |volume=78 |date=1973 |publisher=Princeton University Press |isbn=978-0-6910-8136-6 |jstor=j.ctt1bd6kr9}}
*{{citation|last=Papadopoulos|first= Athanase|title=Metric spaces, convexity and non-positive curvature|edition=Second|series= IRMA Lectures in Mathematics and Theoretical Physics|volume= 6|publisher=[[European Mathematical Society]]|year=2014|isbn=978-3-03719-132-3}}
*{{citation|last=Paulin|first= Frédéric|title=Un groupe hyperbolique est déterminé par son bord|language=fr|journal=J. London Math. Soc. |volume=54 |year=1996|pages= 50–74|doi=10.1112/jlms/54.1.50}}
*{{citation|last=Ratcliffe|first= John G.|title=Foundations of hyperbolic manifolds|edition=Second |series= Graduate Texts in Mathematics|volume= 149|publisher=Springer|year=2006|isbn=978-0387-33197-3}}
*{{citation|last=Roe|first=John|title=Lectures on coarse geometry|series=University Lecture Series|volume= 31|publisher=[[American Mathematical Society]]|year= 2003|isbn=0-8218-3332-4}}
*{{citation|last=Shiohama|first= Katsuhiro|chapter=Topology of complete noncompact manifolds|title=Geometry of geodesics and related topics |pages= 423–450|series=Adv. Stud. Pure Math.|volume= 3|publisher=[[North-Holland Publishing Company|North-Holland]]|year= 1984}}
*{{SpringerEOM|id=Busemann_function&oldid=13719|title=Busemann function|first=T.|last= Shioya}}
*{{citation|last=Väisälä|first= Jussi|title=Quasi-Möbius maps|journal=[[Journal d'Analyse Mathématique]] |volume=44|year=1984|pages= 218–234|doi=10.1007/bf02790198|doi-access=free|s2cid= 189767039}}

[[Category:Geometry]]