Hostname: page-component-78c5997874-lj6df Total loading time: 0 Render date: 2024-11-05T13:57:01.713Z Has data issue: false hasContentIssue false

Optimal coupling of jumpy Brownian motion on the circle

Published online by Cambridge University Press:  04 July 2023

Stephen B. Connor*
Affiliation:
University of York
Roberta Merli*
Affiliation:
University of York
*
*Postal address: Department of Mathematics, University of York, York, YO10 5DD, UK.
*Postal address: Department of Mathematics, University of York, York, YO10 5DD, UK.
Rights & Permissions [Opens in a new window]

Abstract

Consider a Brownian motion on the circumference of the unit circle, which jumps to the opposite point of the circumference at incident times of an independent Poisson process of rate $\lambda$. We examine the problem of coupling two copies of this ‘jumpy Brownian motion’ started from different locations, so as to optimise certain functions of the coupling time. We describe two intuitive co-adapted couplings (‘Mirror’ and ‘Synchronous’) which differ only when the two processes are directly opposite one another, and show that the question of which strategy is best depends upon the jump rate $\lambda$ in a non-trivial way. We also provide an explicit description of a (non-co-adapted) maximal coupling for any jump rate in the case that the two jumpy Brownian motions begin at antipodal points of the circle.

Type
Original Article
Copyright
© The Author(s), 2023. Published by Cambridge University Press on behalf of Applied Probability Trust

1. Introduction

Consider a continuous-time stochastic process $X=(X_t)_{t\ge 0}$ on $\left( { - \pi ,\pi } \right]$ given by

(1.1) \begin{equation} X_t= X_0 + \frac{1}{2}B_t+\pi N^\lambda_t \ (\text{mod}\ 2\pi) , \end{equation}

where $B_t$ is a standard $\mathbb{R}$ -valued Brownian motion and $N^\lambda$ is an independent Poisson process of rate $\lambda\ge0$ . (The factor of $\frac12$ is introduced purely for algebraic convenience in what follows.) X is a Lévy process which can be viewed as a Brownian motion on the circumference of the unit circle which jumps to the antipodal point of the circle at incident times of $N^\lambda$ . For this reason, we will refer to the process X as a jumpy Brownian motion of rate $\lambda$ (JBM( $\lambda$ )).

In this paper we are interested in couplings of two JBM( $\lambda$ ) processes started from different points of the circle. That is, we are interested in processes $(X,\hat X)$ on ${{ ({-}\pi,\pi]}}^2$ such that, viewed marginally, X and $\hat X$ each behave as a copy of JBM( $\lambda$ ). Given a coupling $(X,\hat X)$ , we define the coupling time by ${T} = \inf\{t\ge 0 \,:\, X_s = \hat X_s \text{ for all } s\ge t \}$ . Recall that the tail distribution of any coupling time provides an upper bound on the total variation distance between the laws of $X_t$ and $\hat X_t$ via the coupling inequality [Reference Aldous1]

(1.2) \begin{equation} {\Vert{\mathcal{L}(X_t) - \mathcal{L}(\hat X_t)}\Vert}_{\text{TV}} \le {\mathbb{P}}({T} \gt t) ,\end{equation}

where

(1.3) \begin{equation} {\Vert{\mathcal{L}(X_t) - \mathcal{L}(\hat X_t)}\Vert}_{\text{TV}} = \sup\nolimits_A \{ {\mathbb{P}}(X_t\in A) - {\mathbb{P}}(\hat{X}_t\in A)\} \end{equation}

and where the supremum is taken over all Borel subsets of $\left( { - \pi ,\pi } \right]$ . A coupling is called successful if ${\mathbb{P}}({T} \gt t)\to 0$ as $t\to \infty$ ; it is called maximal if it achieves equality in (1.2) for all times t.

It is well known that a maximal coupling of two random processes exists under extremely mild conditions: see [Reference Goldstein6,Reference Griffeath7,Reference Pitman16] for discrete-time processes, and [Reference Sverchkov and Smirnov18] for the case of càdlàg processes with Polish state-space. However, in most cases, explicit construction of a maximal coupling is extremely difficult, and it is natural for attention to focus on classes of couplings which are more readily realisable. One such class is that of co-adapted couplings.

Definition 1.1. A coupling $(X,\hat X)$ is called co-adapted if there exists a filtration $(\mathcal{F}_t)_{t\ge 0}$ such that X and $\hat X$ are both adapted to $(\mathcal{F}_t)_{t\ge 0}$ and, for any $0\le s\le t$ , $\mathcal{L}(X_t\mid\mathcal{F}_s) = \mathcal{L}(X_t\mid X_s)$ and $\mathcal{L}(\hat X_t\mid\mathcal{F}_s) = \mathcal{L}(\hat X_t\mid\hat X_s)$ .

In other words, both X and $\hat X$ are Markov with respect to the filtration $(\mathcal{F}_t)_{t\ge0}$ . Kendall [Reference Kendall11] refers to co-adapted couplings as immersed, since the condition of Definition 1.1 is equivalent to demanding that the natural filtrations of X and $\hat X$ are both immersed in a common filtration.

Even though maximal couplings certainly need not be co-adapted, there are a few processes for which maximal co-adapted couplings have been shown to exist. Probably the simplest of these is the reflection/mirror coupling of Euclidean Brownian motions, in which the path of one process, until the coupling time, is obtained by reflecting the other in the hyperplane bisecting the line joining their starting points [Reference Lindvall and Rogers13]; indeed, this is the unique maximal co-adapted coupling [Reference Hsu and Sturm8], a result which holds more generally for Brownian motion on a Riemannian manifold [Reference Kuwada12]. More recently, [Reference Banerjee and Kendall2] showed that a maximal co-adapted coupling for smooth elliptic diffusions on a complete Riemannian manifold can only possibly exist if the underlying space is a sphere, Euclidean space, or hyperbolic space.

Of particular relevance to the work of the current paper is the analysis of a symmetric random walk on the hypercube in [Reference Connor and Jacka4], subsequently generalised in [Reference Connor3]. They considered the class of co-adapted couplings for two such random walks, and showed that there exists a stochastically optimal coupling within this class. In other words, they exhibited a coupling whose coupling time $T^*$ satisfies, simultaneously for all $t\ge 0$ , ${\mathbb{P}}(T^*>t) = \min\{{\mathbb{P}}(T>t)\,:\,\text{co-adapted coupling times T}\}$ . Furthermore, they showed that this optimal co-adapted coupling does not achieve equality in (1.2), thus demonstrating that there does not exist a maximal co-adapted coupling for this random walk. A result in the same vein was proved in [Reference Kendall11], which showed that there exists a stochastically optimal co-adapted coupling for the two-dimensional process consisting of Brownian motion together with its local time at 0; numerical evidence indicates that this coupling is once again not maximal. We also highlight here [Reference Jacka, Mijatović and Širaj9], which investigated (amongst other things) whether the reflection coupling for Brownian motions optimises various functions of the coupling time of the corresponding geometric Brownian motions. It showed that the Laplace transform of the coupling time is maximised by the reflection coupling, but that whether this coupling also solves the finite time horizon problem (to minimise the coupling time’s tail distribution) depends upon an underlying drift parameter.

1.1. Motivation and main results

By way of motivation for the study of JBM( $\lambda$ ), consider the following problem. Suppose that X and $\hat X$ start from opposite points on the circle (i.e. $|X_0-\hat X_0| = \pi$ , where $|\cdot|$ denotes the shortest distance between two points on the circumference of the unit circle) and that we wish to couple their evolution so as to minimise (some non-decreasing function of) the corresponding coupling time. Write $D_t = |X_t-\hat X_t| \in [0,\pi]$ . Consider the following two cases:

  1. 1. If $\lambda = 0$ then there are no jumps; in this case the fastest coupling is achieved by setting $\hat X_t = \pi - X_t\ (\text{mod}\ 2\pi)$ , i.e. having the two processes reflect about either of the two points at distance $\pi/2$ from both $X_0$ and $\hat X_0$ . The resulting process $D_t$ is a reflected Brownian motion on $[0,\pi]$ , and the corresponding coupling time is given by ${T} = \inf\{t\,:\, D_t = 0\} = \tau$ , where $\tau$ is the first hitting time of the set $\{-\pi/2,\pi/2\}$ by the Brownian motion $\tfrac12 B_t$ making up the diffusion component of X.

  2. 2. Alternatively, if $\lambda$ is large then we expect jumps to be occurring quickly; in this case it may be better to synchronise the driving Brownian motions, and to let the Poisson processes N and $\hat N$ be independent until the first time that we see an incident on either one. In other words, we maintain $D_t=\pi$ for all $t<J$ , where $J \sim \mbox{Exp}(2\lambda)$ is the first time that either X or $\hat X$ jumps, and then the coupling time is given by ${T}=J$ .

Note that the ways in which the driving Brownian motions are coupled in these two scenarios are as different as can be; the first coupling proceeds by maximising the volatility of $B_t-\hat B_t$ for all $t<{T}$ , whilst the second minimises it (setting it to zero, in fact). Taking this intuition further, it is not unreasonable to suppose that there might be a critical value of $\lambda$ with the property that the second of the two couplings sketched above is better than the other if and only if the jump rate exceeds this critical value.

Furthermore, it is clear that when $\lambda=0$ the mirror coupling described in case 1 above will be maximal, but that this will not be true of the coupling strategy outlined in case 2: if we choose to synchronise the driving Brownian motions and couple at time J, then there is a positive probability that coupling would have occurred faster if the Brownian motions had instead been reflected. Indeed, we show in Section 2 that when $|X_0-\hat X_0| = \pi$ the total variation distance between the distributions of $X_t$ and $\hat X_t$ satisfies

\[{\Vert{\mathcal{L}(X_t\mid X_0=0) - \mathcal{L}(\hat X_t\mid\hat X_0=\pi)}\Vert}_{\text{TV}} = {\mathbb{P}}(\min\{J, \tau\} \gt t),\]

where J and $\tau$ are independent. We shall use this observation to explicitly describe a maximal, but non-co-adapted, coupling of two JBM( $\lambda$ ) processes started from antipodal points of the circle.

In Section 3 we restrict our attention to the class of co-adapted couplings of two JBM( $\lambda$ ) processes with arbitrary starting points, and consider two couplings motivated by the cases outlined above. We refer to these as Mirror and Synchronous couplings, but these names relate solely to the way in which the driving Brownian motions are coupled when $D_t = \pi$ . A formal definition will be given later (Definition 3.1), but in order to state our main results we provide an informal description here.

Definition 1.2. (Informal description of Mirror and Synchronous couplings.) The Mirror and Synchronous couplings treat the jump and diffusion components of X and $\hat X$ as follows, until the coupling time T.

Jumps:

  • if $D_t \in(0,\pi/2]$ , synchronise the driving Poisson processes (i.e. make $\hat X$ jump at time t if and only if X sees a jump at time t);

  • if $D_t \in(\pi/2,\pi]$ , make the driving Poisson processes independent.

Diffusion:

  • if $D_t \in(0,\pi)$ , reflect the driving Brownian motions (i.e. set $\hat B_t = - B_t$ );

  • if $D_t=\pi$ , then:

    • for the Mirror coupling, reflect the driving Brownian motions (i.e. continue to set $\hat B_t = - B_t$ );

    • for the Synchronous coupling, synchronise the driving Brownian motions (i.e. set $\hat B_t = B_t$ ).

Trivially, the Mirror coupling is successful for all jump rates $\lambda\ge 0$ , while the Synchronous coupling is successful for any $\lambda>0$ . We note that our Synchronous coupling has a similar two-stage approach to the ‘reflection/synchronised coupling’ of Brownian motion together with local time described in [Reference Kendall11]. In both cases the Brownian motions are reflected until some stopping time $T_1$ (in our case, the first hitting time of $D_t$ on the set $\{0,\pi\}$ ); if coupling has not occurred at time $T_1$ then the Brownian motions are subsequently synchronised until the coupling time, T.

In Section 3 we calculate the Laplace transform of the associated coupling times, ${{{T_{\text{M}}}}}$ (Mirror) and ${{{T_{\text{S}}}}}$ (Synchronous), and then use ideas from stochastic control to prove that the Synchronous coupling is faster than the Mirror coupling if and only if $\lambda\ge \lambda^\star$ , where $\lambda^\star=0.083\,37\dots$ is the unique solution to the equation

(1.4) \begin{equation} \text{cosech}(\pi\sqrt{\lambda}) + 2\pi\sqrt{\lambda} = 2\coth(\pi\sqrt{\lambda}).\end{equation}

More precisely, for this range of $\lambda$ the Synchronous coupling turns out to uniquely maximise the Laplace transform of the coupling time within the class of all co-adapted couplings.

Theorem 1.1. For any $\lambda\ge \lambda^\star$ the Synchronous coupling uniquely maximises, within the class of co-adapted couplings, the Laplace transform of the corresponding coupling time. That is, for any $\lambda\ge\lambda^\star$ and any $x\in(0,\pi]$ , the coupling time ${{{T_{\text{S}}}}}$ maximises ${\mathbb{E}}\left[{\text{e}^{-\gamma T}\mid D_0=x}\right]$ simultaneously for all $\gamma\ge 0$ .

It does not appear that there is any co-adapted coupling which simultaneously maximises the Laplace transform for all values of $\gamma$ when $\lambda<\lambda^\star$ (see Remark 3.2). However, the Mirror coupling does turn out to minimise the mean coupling time for this set of jump rates.

Theorem 1.2. For any $\lambda\in[0,\lambda^\star)$ , the Mirror coupling minimises the mean coupling time within the class of co-adapted couplings. That is, for any $\lambda\in[0,\lambda^\star)$ and any $x\in(0,\pi]$ , the coupling time ${{{T_{\text{M}}}}}$ minimises ${\mathbb{E}}\left[{T\mid D_0=x}\right]$ within the class of co-adapted couplings.

In the following we use the terms ‘LT-optimal’ and ‘mean-optimal’, and usually refrain from writing ‘within the class of co-adapted couplings’. Note that when $\lambda=0$ (no jumps) the Mirror coupling reduces to the reflection coupling for Brownian motion and is therefore maximal. Similarly, when $\lambda\to \infty$ the Synchronous coupling makes jumps from $x\mapsto\pi-x$ occur immediately for any $x>\pi/2$ , and so in the limit the difference process D behaves like reflected Brownian motion on $[0,\pi/2]$ ; it is simple to see that this will once again be maximal. Table 1 summarises the non-trivial way in which the various optimality properties of our coupling strategies depend upon the jump rate.

Table 1. Summary of optimal coupling properties for different values of the jump rate $\lambda$ .

JBM( $\lambda$ ) therefore makes an interesting addition to the relatively small number of examples in the literature of processes for which precise results about optimal co-adapted couplings have been established. In particular, to the best of the authors’ knowledge, our results are the first of their kind for a Lévy process with both continuous and jump components.

Remark 1.1. Our original proofs of Theorems 1.1 and 1.2 were based on excursion theory of Brownian motion, using calculations similar to those in [Reference Rogers and Williams17, Section VI.56]. We are grateful to an anonymous referee who suggested the more direct method used in Section 3. However, we still find the excursion approach very appealing, and the corresponding calculations can be found in the original arXiv version of this paper [Reference Connor and Merli5]. Further details can also be found in the PhD thesis of the second author [Reference Merli14], alongside supporting evidence for the correctness of our results obtained via direct simulation of the two couplings.

2. Maximal coupling

In this section we briefly describe the construction of a maximal, but non-co-adapted, coupling in the case that $D_0 = \pi$ (i.e. the two coupled JBM( $\lambda$ ) processes begin at opposite points of the circle). Let $(X_0,\hat X_0) = (0,\pi)$ . For any set $A\subset {{{({-}\pi,\pi]}}} $ ,

\[ {\mathbb{P}}(X_t \in A) = \sum_{k=-\infty}^\infty {\mathbb{P}}(\tfrac12 B_t \in A + 2k\pi,\,N^\lambda(t) \text{ even}) + {\mathbb{P}}(\tfrac12 B_t \in A + (2k+1)\pi, \, N^\lambda(t) \text{ odd}).\]

Since $\mathbb{P}(N^\lambda(t) \text{ even}) = \frac12({1+\text{e}^{-2\lambda t}}) \ge \mathbb{P}(N^\lambda(t) \text{ odd})$ , it is clear that the set on which $X_t$ has greater density than $\hat X_t$ is, for all $t\ge 0$ , the half-circle $A^\star=({-}\pi/2,\pi/2)$ centred at $X_0$ . As above, let $\tau$ denote the first hitting time of the set $\{-\pi/2,\pi/2\}$ by the Brownian motion $\tfrac12 B_t$ . By the definition of total variation distance (1.3), we see that

(2.1) \begin{align} {\Vert{\mathcal{L}(X_t\mid X_0=0) - \mathcal{L}(\hat X_t\mid \hat X_0 = \pi)}\Vert}_{\text{TV}} & = {\mathbb{P}}(X_t\in A^\star) - \mathbb P({\hat X_t\in A^\star}) \nonumber \\ & = 2{\mathbb{P}}(X_t\in A^\star) - 1 \nonumber \\ & = \text{e}^{-2\lambda t} \big(2{\mathbb{P}}(\tfrac12 B_t \in A^\star + 2\pi\mathbb Z) - 1 \big) \nonumber \\ & = \text{e}^{-2\lambda t} {\mathbb{P}}(\tau \gt t), \end{align}

where the final equality follows by the reflection principle. Thus we see that when X and $\hat X$ begin from antipodal points, the total variation distance between their laws is given by the tail distribution of the random variable $\min\{J, \tau\}$ , where $J\sim\mbox{Exp}(2\lambda)$ .

This observation leads to an explicit construction of a maximal coupling, as follows. Let $N^{2\lambda}$ be a marked Poisson process of rate $2\lambda$ , whose incident times are denoted by $J_k$ , and with marks $Y_k{{{\stackrel{i.i.d.}{\sim}}}}\mbox{Bernoulli}(1/2)$ . Let B be a standard Brownian motion started at 0, with hitting time $\tau$ defined as above. We first of all use these to define the JBM( $\lambda$ ) process X: the diffusion component of X is given by $\tfrac12 B$ , and X makes a jump of size $\pi$ at incident times $J_k$ if and only if the corresponding mark $Y_k$ equals 1.

We now define $\hat X$ , started at $\pi$ , in such a way that $\hat X$ and X couple almost surely at time $T^\star=\min\{J_1,\tau\}$ :

  • On the event $\{\tau< J_1\}$ , we let the diffusion component of $\hat X$ equal $-\tfrac12 B$ , so that $\hat X$ is the reflection of X about the points $\pm\pi/2$ . In this case $\hat X$ and X will meet at time $\tau$ .

  • On the event $\{J_1<\tau\}$ , we let the diffusion component of $\hat X$ equal $\tfrac12 B$ (i.e. $\hat X$ and X move synchronously, remaining at distance $\pi$ from each other) until time $J_1$ . We then make $\hat X$ jump by $\pi$ at time $J_1$ if and only if the mark $Y_1$ equals 0. That is, $\hat X$ jumps at time $J_1$ if and only if X does not jump at that time, and since their diffusion components were synchronised, $\hat X$ and X will meet exactly when one of them jumps at time $J_1$ .

  • From time $\min\{J_1,\tau\}$ onwards, the diffusion and jump components of $\hat X$ equal those of X.

This is a valid coupling (the process $\hat X$ , viewed marginally, really is a JBM( $\lambda$ ) process started from $\pi$ ). Furthermore, it is clearly maximal but not co-adapted: the evolution of X is adapted to the natural filtration generated by $\{N^{2\lambda},B\}$ , whereas the evolution of $\hat X$ until the coupling time depends on knowledge of which of the times $J_1$ and $\tau$ occurs first.

Remark 2.1. For general starting states ( $D_0=x<\pi$ ) such an explicit description of a maximal coupling is significantly more challenging. The main complication is that it is no longer true that the set on which $X_t$ has greater density than $\hat X_t$ is independent of t. Rather, with $(X_0,\hat X_0) = ({-}x/2,x/2)$ , the set $A^\star_t$ on which $X_t$ has greater density takes the form $A^\star_t = ({-}r_t,0)\cup(r_t,\pi)$ , where $r_t=r_t(x)$ is the unique point in $(0,\pi)$ at which the densities of $X_t$ and $\hat X_t$ agree. Depending upon the value of $\lambda$ , $r_t$ tends to either $\pi/2$ or $\pi$ as $t\to\infty$ , but not necessarily monotonically. The only exception to this rule is in the limit as $\lambda\to \infty$ , when $r_t \to \pi/2$ for all values of t (and for all $x<\pi$ ); in this case $A^\star_t$ does not depend upon t, and an explicit maximal coupling is once again straightforward to describe. (See the comment before Table 1.)

3. Co-adapted couplings

We begin the search for optimal co-adapted couplings by adopting the perspective of stochastic control. We wish to work with a pair of JBM( $\lambda$ ) processes which are adapted to a common filtration of $\sigma$ -algebras $(\mathcal{F}_t)$ ; it is more convenient to work initially with jumpy Brownian motions $(X,\hat X)$ on the real line (i.e. processes satisfying (1.1) without the (mod $2\pi$ )) and then project to the circle. Let $(\mathcal{F}_t)_{t\ge 0}$ be any filtration to which the following independent random processes are all adapted: two standard Brownian motions, B and $\tilde B$ ; two marked Poisson processes $N^\lambda$ and $\tilde N^\lambda$ of rate $\lambda\ge 0$ , whose marks ( $U_t$ and $\tilde U_t$ , respectively) are independent and identically distributed as Uniform $[0,1]$ .

First, note that any co-adapted coupling of two Brownian motions B and $\hat B$ can be represented by the stochastic differential equation ${{{\text{d}}}} \hat B_t = \theta_t \,{{{\text{d}}}} B_t + \sqrt{1-\theta_t^2}\, {{{\text{d}}}} \tilde B_t$ , where $\theta_t$ is a predictable random process taking values in $[{-}1,1]$ [Reference Kendall10, Lemma 6]. Similarly, as explained in [Reference Connor and Jacka4], any co-adapted coupling of two Poisson processes $(N^\lambda,\hat N^\lambda)$ can be written as

(3.1) \begin{equation} \hat N^\lambda ({{{\text{d}}}} t) = {\mathbf{1}}_{[U_t\le p_{t-}]} N^\lambda({{{\text{d}}}} t) + {\mathbf{1}}_{[\tilde U_t> p_{t-}]} \tilde N^\lambda({{{\text{d}}}} t) ,\end{equation}

where $p_t$ is a càdlàg control process adapted to $(\mathcal{F}_t)$ , taking values in [0,1]. (More accurately, [Reference Connor and Jacka4] describes any co-adapted coupling of two unit-rate random walks on the hypercube, $\mathbb Z_2^n$ , in terms of $(n+1)^2$ independent marked Poisson processes which are controlled by a doubly stochastic matrix-valued process; (3.1) is just a simplified parametrisation of their result when $n=1$ .)

Combining these two results, it is clear that the joint process $(X_t,\hat X_t)$ satisfies

(3.2) \begin{equation} \begin{pmatrix} {{{\text{d}}}} X_t\\ {{{\text{d}}}} \hat X_t \end{pmatrix} = \frac12 \begin{pmatrix} 1 & 0 \\ \theta_t &\sqrt{1-\theta_t^2} \end{pmatrix} \begin{pmatrix} {{{\text{d}}}} B_t\\ {{{\text{d}}}} \tilde B_t \end{pmatrix} + \pi \begin{pmatrix} 1 & 0 \\ {\mathbf{1}}_{[U_t\le p_{t-}]} & {\mathbf{1}}_{[\tilde U_t> p_{t-}]} \end{pmatrix} \begin{pmatrix} N^\lambda ({{{\text{d}}}} t) \\ \tilde N^\lambda ({{{\text{d}}}} t) \end{pmatrix} ,\end{equation}

where the control process $c_t = (\theta_t,p_t)$ takes values in $[{-}1,1]\times [0,1]$ and is adapted to the filtration $(\mathcal{F}_t)$ . The control process explicitly determines the dependence between X and $\hat X$ . In particular, if $\theta_t=0$ then their continuous components are independent; if $\theta_t = 1$ then they are synchronised; and if $\theta_t = -1$ then they are mirror or reflection coupled. Similarly, if $p_t=1$ then the jump components of X and $\hat X$ are synchronised, whereas if $p_t=0$ then they are independent.

In what follows, we will be interested in the difference process $Z_t = X_t-\hat X_t$ , and in particular the time it takes for this to hit a multiple of $2\pi$ (at which time the projections of the two processes onto the unit circle will meet). Using (3.2), we see that

\[ {{{\text{d}}}} Z_t = \frac12\big({(1-\theta_t)\,{{{\text{d}}}} B_t - \sqrt{1-\theta_t^2}\,{{{\text{d}}}} \tilde B_t }\big) + \pi \big({(1-{\mathbf{1}}_{[U_t\le p_{t-}]}) N^\lambda ({{{\text{d}}}} t) - {\mathbf{1}}_{[\tilde U_t \gt p_{t-}]} \tilde N^\lambda ({{{\text{d}}}} t)}\big).\]

Thus, Z has the same dynamics as a Lévy process on $\mathbb{R}$ whose continuous component has volatility $(1-\theta)/2$ and which makes jumps of size $+\pi$ and $-\pi$ each at rate $\lambda(1-p)$ . For any $z\in\mathbb{R}$ the distances of $z+\pi$ and $z-\pi$ from the set $2\pi\mathbb Z$ are equal, and so the distribution of the time taken for Z to hit the set $2\pi\mathbb Z$ is unchanged if we alter the dynamics to make all jumps of size $+\pi$ . Furthermore, the independence of N and $\tilde N$ means that

\[ (1-{\mathbf{1}}_{[U_t\le p_{t-}]}) N^\lambda ({{{\text{d}}}} t) + {\mathbf{1}}_{[\tilde U_t \gt p_{t-}]} \tilde N^\lambda ({{{\text{d}}}} t) \ {{{\stackrel{{\text{D}}}{=}}}}\ {\mathbf{1}}_{[U^{\prime}_t \gt p_{t-}]} N^{2\lambda}({{{\text{d}}}} t) ,\]

where ${{{\stackrel{{\text{D}}}{=}}}}$ denotes equality in distribution, and where $N^{2\lambda}$ is a marked Poisson process of rate $2\lambda$ with marks $U^{\prime}\sim\mbox{Uniform}[0,1]$ . Thus, it suffices, for any given adapted control process $c_t = (\theta_t,p_t)$ , to consider the hitting time on the set $2\pi\mathbb Z$ of the process given by

\[ {{{\text{d}}}} Z_t = \sqrt{\frac{1-\theta_t}{2}}\,{{{\text{d}}}} B_t+\pi {\mathbf{1}}_{[U^{\prime}_t> p_{t-}]} N^{2\lambda}({{{\text{d}}}} t); \qquad Z_0 = X_0-\hat X_0 .\]

As in Section 1, we shall work in the following with the process $D_t=|Z_t-2\pi\mathbb Z|\in[0,\pi]$ , which measures the distance between the projections of X and $\hat X$ on the circumference of the unit circle. Until the coupling time $T=\inf\{t\,:\, D_t=0\}$ , D behaves like a reflected Brownian motion (with volatility controlled by $\theta$ ), and with a jump (controlled by p) at time t having the effect of making $D_t = \pi -D_{t-}$ . (That is, a jump of size $\pi$ in Z results in D being reflected about the point $\pi/2$ .)

This allows us to view the search for an optimal co-adapted coupling as a stochastic control problem. For a given value function, we seek a control process $c_t=(\theta_t,p_t) \in [{-}1,1]\times[0,1]$ such that the time taken for the corresponding difference process D to hit zero minimises/maximises the value function, as appropriate. As noted in [Reference Øksendal and Sulem15], it suffices to restrict attention to Markov controls of the form $c_t = c(D_{t-})$ .

Using this setup we can now give a precise definition of the Mirror and Synchronous couplings for two JBM( $\lambda$ ) processes (cf. Definition 1.2).

Definition 3.1. Suppose $D_{t-}=x\in[0,\pi]$ . The Mirror coupling is the co-adapted coupling with control $c_{\text{M}} = (\theta_{\text{M}},p_{\text{M}}) $ at time t given by

\[ \text{Diffusion:} \quad\theta_{\text{M}}(x) = 2{\mathbf{1}}_{[x=0]} -1; \qquad \text{Jumps:} \quad p_{\text{M}}(x) = {\mathbf{1}}_{[x\le\pi/2]} . \]

Similarly, the Synchronous coupling uses control $c_{\text{S}} = (\theta_{\text{S}},p_{\text{S}} )$ , where:

\[ \text{Diffusion:} \quad\theta_{\text{S}}(x) = 2{\mathbf{1}}_{[x=0\text{ or }x=\pi]} -1; \qquad \text{Jumps:} \quad p_{\text{S}}(x) = {\mathbf{1}}_{[x\le\pi/2]} \,. \]

As explained above, D behaves like a Brownian motion on $(0,\pi)$ , for which we can control the speed using $\theta$ , and if D sees a jump then this simply reflects it around the point $\pi/2$ . Since we wish to minimise the time taken for D to hit 0, it seems intuitively sensible to maximise the speed of the diffusion when $D_t\in(0,\pi)$ , and to maximise/minimise the jump rate according to whether or not a jump would reduce the value of $D_t$ ; this is exactly what both Mirror and Synchronous couplings achieve. The only difference between the two couplings is in the choice of $\theta(\pi)$ : $\theta_{\text{M}}(\pi) = -1$ forces the Brownian component to reflect downwards from the barrier at $\pi$ , whereas $\theta_{\text{S}}(\pi) = 1$ means that D waits an $\mbox{Exp}(2\lambda)$ amount of time at level $\pi$ before jumping directly to 0.

3.1. LT-optimality

In this section we prove the LT-optimality of the Synchronous coupling when $\lambda\ge\lambda^\star$ , as claimed in Theorem 1.1. We begin by proving explicit formulas for the Laplace transforms of our two co-adapted coupling strategies.

Given the jump rate $\lambda$ and a constant $\gamma\ge 0$ , we write $\alpha =\sqrt{2(2\lambda+\gamma)}$ and $\beta=\sqrt{2\gamma}$ , and use these to define the following four non-negative constants (depending implicitly upon $\lambda$ and $\gamma$ ):

\begin{align*} \kappa_1 & = \frac{2 \cosh\big(\alpha\tfrac{\pi}{2}\big)\sinh\big(\beta\tfrac{\pi}{2}\big)} {2\cosh\big(\alpha\tfrac{\pi}{2}\big)\cosh\big(\beta\tfrac{\pi}{2}\big)-1} ; & \kappa_2 & = \frac{\beta}{\alpha} \kappa_1 \text{sech}\bigg(\alpha\frac{\pi}{2}\bigg); \nonumber \\ \kappa_3 & = \text{sech}\bigg(\beta\frac{\pi}{2}\bigg)\bigg(\frac{\beta}{2\alpha}\text{cosech}\bigg(\alpha\frac{\pi}{2}\bigg) + \sinh\bigg(\beta\frac{\pi}{2}\bigg)\bigg); & \kappa_4 & = \bigg({\frac{\beta}{\alpha}}\bigg)^2\text{cosech}\bigg(\alpha\frac{\pi}{2}\bigg). \end{align*}

Lemma 3.1. The Laplace transforms for the coupling times ${{{T_{\text{M}}}}}$ and ${{{T_{\text{S}}}}}$ under the Mirror and Synchronous couplings started from distance $x\in[0,\pi]$ are given by the following formulas.

Mirror coupling:

\[ \mathbb{E}_{x}\big[{\text{e}^{-\gamma {T_{\text{M}}}}}\big] = w_{\text{M}}(x) \,:\!=\, \begin{cases} \cosh(\beta x) - \kappa_1 \sinh(\beta x), & \quad 0\le x\le\tfrac{\pi}{2}; \\ w_{\text{M}}(\pi-x)-\kappa_2\sinh\big(\alpha\big(x-\tfrac{\pi}{2}\big)\big), & \quad \tfrac{\pi}{2}<x\le\pi. \end{cases} \]

Synchronous coupling:

\[ \mathbb{E}_{x}\big[{\text{e}^{-\gamma {T_{\text{S}}}}}\big] = w_{\text{S}}(x) \,:\!=\, \begin{cases} \cosh(\beta x) - \kappa_3 \sinh(\beta x), & \quad 0\le x\le\tfrac{\pi}{2}; \\ w_{\text{S}}(\pi-x) - \kappa_4\sinh\big(\alpha\big(x-\tfrac{\pi}{2}\big)\big) & \quad \tfrac{\pi}{2}<x\le\pi. \end{cases} \]

Remark 3.1. Note that when $x=\pi$ , the final formula in Lemma 3.1 simplifies to

\[ \mathbb{E}_{\pi}\big[{\textrm{e}^{-\gamma {T_{\textrm{S}}}}}\big] = 1 - \bigg({\frac{\beta}{\alpha}}\bigg)^2 = \frac{2\lambda}{2\lambda+\gamma}; \]

this is just the Laplace transform of an $\mbox{Exp}(2\lambda)$ distribution, which is consistent with the observation that under the Synchronous coupling the two processes will meet precisely when one of them jumps to the other side of the circle. Similarly, the formula for the Mirror coupling when $\lambda=0$ and $x=\pi$ reduces to $\mathbb{E}_{\pi}[{\textrm{e}^{-\gamma T_{\textrm{M}}}}] = \textrm{sech}(\pi\sqrt{2\gamma})$ , which is the Laplace transform of the time taken for a standard Brownian motion to hit $\pm\pi$ (and hence of the time taken for a Brownian motion started from, and reflected at, $\pi$ to hit 0).

Before proceeding any further, it will be helpful to quickly establish some basic properties of the function $w_{\text{S}}$ defined in Lemma 3.1.

Lemma 3.2. For any fixed $\gamma\ge 0$ , $w_{\text{S}}(x)$ is a convex function of $x\in[0,\pi]$ , and $w^{\prime\prime}_{\text{S}}(\pi)=0$ .

Proof. It is clear from the definition of $w_{\text{S}}$ that this function is continuous on $[0,\pi]$ and $C^2$ on $(0,\pi)$ , with second derivative given by

\[ w^{\prime\prime}_{\text{S}}(x) = \begin{cases} \beta^2w_{\text{S}}(x), & 0\le x\le \pi/2; \\ \beta^2w_{\text{S}}(\pi-x) - \alpha^2\kappa_4\sinh\big(\alpha\big(x-\tfrac{\pi}{2}\big)\big), & \pi/2<x\le\pi. \end{cases} \]

This is clearly non-negative for $x\le \pi/2$ , and zero for $x=\pi$ . For $x>\pi/2$ , inserting the formula for $w_{\text{S}}(\pi-x)$ and expanding shows that

(3.3) \begin{equation} w^{\prime\prime}_{\text{S}}(x) \propto s_\beta c_\beta\cosh(\beta(\pi-x)) - \bigg(\frac{\beta s_\beta}{2\alpha s_\alpha} + s_\beta^2 \bigg)\sinh(\beta(\pi-x)) - \frac{s_\beta c_\beta}{s_\alpha} \sinh\bigg(\alpha\bigg(x-\frac{\pi}{2}\bigg)\bigg), \end{equation}

where we have written $c_u = \cosh\big(u\tfrac{\pi}{2}\big)$ and $s_u = \sinh\big(u\tfrac{\pi}{2}\big)$ . Now note that $u\mapsto us_u$ is an increasing function, and that $u\mapsto \sinh\big(u\big(x-\tfrac{\pi}{2}\big)\big)/s_u$ is decreasing for any fixed $x\in[\pi/2,\pi]$ . Thus (recalling that $\beta\le \alpha$ ) the right-hand side of (3.3) can be bounded below by

\begin{align*}s_\beta c_\beta\cosh&\ (\beta(\pi-x)) - \bigg(\frac12 + s_\beta^2\bigg)\sinh(\beta(\pi-x)) - c_\beta \sinh\bigg(\beta\bigg(x-\frac{\pi}{2}\bigg)\bigg) \\& = \frac12 \sinh(\beta (\pi-x)) \ge 0 . \\[-30pt] \end{align*}

Proof of Lemma 3.1. We deal first with the Synchronous coupling. It is clear that $w_{\text{S}}(0) = 1$ , and we have already observed that $w_{\text{S}}(\pi) = 2\lambda/(2\lambda+\gamma)$ , so $w_{\text{S}}$ certainly takes the correct values at the two boundary points. Given $D_0 = x\in[0,\pi]$ , consider the (bounded) process $W_t(x) = \text{e}^{-\gamma t} w_{\text{S}}(D_t)$ . We shall show that W is a martingale until time $T_{\text{S}}$ , and then the optional stopping theorem will quickly yield that $\mathbb{E}_{x}\big[{\text{e}^{-\gamma {T_{\text{S}}}}}\big] = \mathbb{E}\big[{W_{T_{\textrm{S}}}(x)}\big] = W_0(x) = w_{\text{S}}(x)$ , as required.

The martingale property is trivial when $x=\pi$ . For $x\in(0,\pi)$ , we apply Itô’s formula to $W_t(x)$ to obtain

(3.4) \begin{align} \text{e}^{\gamma t}\,{{{\text{d}}}} W_t(x) & = \bigg({-}\gamma w_{\text{S}}(D_t) + \frac{1}{2}w^{\prime\prime}_{\text{S}}(D_t)\bigg)\,{{{\text{d}}}} t + w^{\prime}_{\text{S}}(D_t)\,{{{\text{d}}}} D_t \nonumber \\ & = \bigg({-}\gamma w_{\text{S}}(D_t) + \frac12 w^{\prime\prime}_{\text{S}}(D_t)\bigg)\,{{{\text{d}}}} t + {{{\text{d}}}} Q_t \nonumber \\ & \quad + {\mathbf{1}}_{[D_t> \pi/2]}\{w_{\text{S}}(\pi-D_t) - w_{\text{S}}(D_t)\} N^{2\lambda}({{{\text{d}}}} t), \end{align}

where $Q_t$ is a martingale. Using the formula for $w^{\prime\prime}_{\text{S}}$ established in the proof of Lemma 3.2, and recalling that $\alpha^2 = 2(2\lambda+\gamma)$ and $\beta^2 = 2\gamma$ , (3.4) simplifies to become

(3.5) \begin{equation} \text{e}^{\gamma t}\,{{{\text{d}}}} W_t(x) = {{{\text{d}}}} Q_t + {\mathbf{1}}_{[D_t> \pi/2]}\kappa_4\sinh\bigg(\alpha\bigg(D_t-\frac{\pi}{2}\bigg)\bigg) (N^{2\lambda}({{{\text{d}}}} t)-2\lambda\, {{{\text{d}}}} t). \end{equation}

The final bracketed term on the right-hand side is of course a compensated Poisson process, and we conclude that W is indeed a martingale.

The Mirror coupling case is almost identical: repeating the argument above using $w_{\text{M}}$ in place of $w_{\text{S}}$ shows that, when $x\in(0,\pi)$ , (3.5) holds with $\kappa_4$ replaced by $\kappa_2$ . It only remains to check that the appropriate boundary conditions are satisfied for $w_{\text{M}}$ : here we require $w_{\text{M}}(0)=1$ and $w^{\prime}_{\text{M}}(\pi) = 0$ , and these follow trivially from the definition of $w_{\text{M}}$ .

We are now in a position to complete the proof of Theorem 1.1.

Proof of Theorem 1.1. Let $\mathcal{C}$ denote the set of all successful co-adapted couplings. For any coupling $c\in\mathcal{C}$ , and any fixed $\gamma\ge 0$ , write $w_c(x) = \mathbb{E}_x \big[{\text{e}^{-\gamma T_c}}\big]$ for the Laplace transform of the associated coupling time when starting from $D_0=x\in[0,\pi]$ . The value function $\hat w(x) = \sup_{c\in \mathcal{C}} w_c(x)$ solves the Hamilton–Jacobi–Bellman (HJB) equation [Reference Øksendal and Sulem15, Chapter 3], which may be derived here using Itô’s formula in a similar manner to (3.4). We first of all deal with the case when $x\in[0,\pi)$ , for which the HJB equation is as follows:

(3.6) \begin{align} \sup_{\theta,\, p}\bigg({-}\gamma\hat w(x) + \frac{(1-\theta(x))}{4}\hat w^{\prime\prime}(x) + 2\lambda(1-p(x))\{\hat w(\pi-x)-\hat w(x)\}\bigg) = 0 . \end{align}

Consider the Laplace transform for the Synchronous coupling, $w_{\text{S}}(x)$ . We saw in Lemma 3.2 that this is convex on $[0,\pi]$ ; moreover, using the formula in Lemma 3.1 it may quickly be checked that $w_{\text{S}}(\pi-x)-w_{\text{S}}(x)>0$ if and only if $x\in(\pi/2,\pi]$ . Thus, if we replace $\hat w(x)$ with $w_{\text{S}}(x)$ on the left-hand side of (3.6), the corresponding supremum for $x\in[0,\pi)$ is obtained by taking $\theta(x) = -1$ and $p(x) = {\mathbf{1}}_{[x\le\pi/2]}$ ; since these values agree with the control $(\theta_{\text{S}},p_{\text{S}})$ in Definition 3.1, it also follows that this supremum is indeed zero. Therefore, the function $w_{\text{S}}$ satisfies (3.6) for $x\in[0,\pi)$ ; if we can also show that $w_c(\pi)\le w_{\text{S}}(\pi)$ for all couplings $c\in\mathcal{C}$ , then it will follow that $\hat w = w_{\text{S}}$ , i.e. that the Synchronous coupling is LT-optimal (see [Reference Øksendal and Sulem15, Theorem 3.1]).

Starting from $D_0=\pi$ , if $\theta_c(\pi)\neq 1$ then the diffusion component of D reflects downwards off the barrier at $\pi$ . This, along with the fact that $w^{\prime\prime}_{\text{S}}(\pi) =0$ (Lemma 3.2), implies that the relevant quantity to be maximised over admissible controls $(\theta_c,p_c)$ becomes

(3.7) \begin{align} -\gamma w_{\text{S}}(\pi) + 2\lambda (1-p_c(\pi)) \left\{1- w_{\text{S}}(\pi)\right\} - \sqrt{\frac{1-\theta_c(\pi)}{2}}w^{\prime}_{\text{S}}(\pi). \end{align}

We therefore need to show that this expression is maximised by taking $p_c(\pi) = p_{\text{S}}(\pi) = 0$ and $\theta_c(\pi) = \theta_{\text{S}}(\pi) = 1$ . The first of these is trivial, since $w_{\text{S}}(\pi)\in(0,1]$ for all values of $\lambda$ and $\gamma$ ; the second will follow if and only if $w^{\prime}_{\text{S}}(\pi)\ge 0$ for all values of $\gamma\ge 0$ . To make explicit the dependence of $w^{\prime}_{\text{S}}(\pi)$ on the underlying parameters $\lambda$ and $\gamma$ , let us define $h(\lambda,\gamma) \,:\!=\, w^{\prime}_{\text{S}}(\pi) = \beta\kappa_3 - \alpha\kappa_4\cosh\big(\alpha\tfrac{\pi}{2}\big)$ . The function h is continuous in both arguments, with $h(\lambda,0)=0$ and

\[ \frac{{{{\text{d}}}} h(\lambda,\gamma)}{{{{\text{d}}}}\gamma}\bigg\vert_{\gamma=0} = \pi + \frac{\text{cosech}(\pi\sqrt{\lambda})}{2\sqrt{\lambda}} - \frac{\coth(\pi\sqrt{\lambda})}{\sqrt{\lambda}}. \]

This derivative at $\gamma=0$ is a strictly increasing function of $\lambda$ , equalling zero if and only if $\lambda=\lambda^\star$ (recall the defining equation for $\lambda^\star$ , (1.4)). Furthermore, $h(\lambda,\gamma)$ is non-negative for all $\gamma$ if and only if $\lambda\ge\lambda^\star$ ; if $\lambda<\lambda^\star$ then there exist values of $\gamma>0$ for which h is negative, and others for which it is positive: see Figure 1. It follows that $w^{\prime}_{\text{S}}(\pi)$ is non-negative for all $\gamma\ge 0$ exactly when $\lambda\ge \lambda^\star$ , and we conclude that for this range of jump rates the Synchronous coupling is LT-optimal, as required.

Figure 1. The function $h(\lambda,\gamma)$ for four representative values of $\lambda$ . For $\lambda<\lambda^\star$ the function takes first negative then positive values as $\gamma$ increases; for $\lambda\ge \lambda^\star$ the function is non-negative for all $\gamma$ .

Remark 3.2. For $\lambda<\lambda^\star$ , the proof given above fails only when $x=\pi$ , when it is no longer true that $h(\lambda,\gamma)$ is non-negative for all $\gamma$ . For such a $\lambda$ , if $\gamma$ is such that $h(\lambda,\gamma)<0$ , then the expression in (3.7) is clearly maximised by setting $\theta_c(\pi) =\theta_{\text{M}}(\pi) = -1$ , i.e. by using the Mirror coupling. However, any attempt to prove LT-optimality of the Mirror coupling when $\lambda<\lambda^\star$ is thwarted by the function $w^{\prime\prime}_{\text{M}}(\pi)$ taking both positive and negative values as $\gamma$ varies: as above, we arrive at the situation where for some values of $\gamma$ the optimal control is to set $\theta(\pi) =1$ , and for others it is to set $\theta(\pi)=-1$ . Figure 2 shows that when $\lambda<\lambda^\star$ , $w_{\text{M}}(\pi) \ge w_{\text{S}}(\pi)$ when $\gamma$ is small, but that this relationship is reversed for larger $\gamma$ , and so neither the Mirror nor Synchronous strategy is LT-optimal for these jump rates. This implies that when $\lambda<\lambda^\star$ there is no co-adapted coupling which maximises the Laplace transform ${\mathbb{E}}\left[{\text{e}^{-\gamma T}}\right]$ simultaneously for all values of $\gamma$ .

Figure 2. Difference between Laplace transforms corresponding to Mirror and Synchronous strategies: both graphs show $\mathbb{E}_{\pi}[\text{e}^{-\gamma{{{T_{\text{M}}}}}} - \text{e}^{-\gamma{{{T_{\text{S}}}}}}]$ as a function of $\gamma$ , for different jump rates $\lambda\le\lambda^\star$ . Neither coupling is LT-optimal when $\lambda<\lambda^\star$ .

3.2. Mean-optimality

In this final section we compare the mean coupling times under the Mirror and Synchronous couplings.

Lemma 3.3. The expectations of the coupling times under the Mirror and Synchronous couplings started from distance $x\in[0,\pi]$ are given by

(3.8) \begin{equation} {\mathbb{E}}_x [{T}]= \begin{cases} x(\pi-x)+C(\lambda)x & \hbox{if } 0\le x\le\frac{\pi}{2}, \\ x(\pi-x)+C(\lambda)(\pi-x)+C(\lambda)\frac{\sinh(\sqrt{\lambda}(2x-\pi))}{\sqrt{\lambda}} & \hbox{if } \frac{\pi}{2}<x\le \pi , \end{cases} \end{equation}

where

\begin{equation*} C(\lambda)= \begin{cases} C_{\text{M}}(\lambda) = \frac{\pi}{2\cosh(\pi\sqrt{\lambda})-1} & \quad \text{if $T={{{T_{\text{M}}}}}$ (Mirror coupling)}, \\ C_{\text{S}}(\lambda) = \frac{\text{cosech}(\pi\sqrt{\lambda})}{2\sqrt{\lambda}} & \quad \text{if $T={{{T_{\text{S}}}}}$ (Synchronous coupling).} \end{cases} \end{equation*}

Proof. This can be proved in a couple of ways. Probably the simplest, but most tedious, is to start from Lemma 3.1 and calculate ${\mathbb{E}}_x[T] = \lim_{\gamma\to0}\big({-}({{{{\text{d}}}}}/{{{{\text{d}}}}\gamma}){\mathbb{E}}_x\big[{\text{e}^{-\gamma T}}\big]\big)$ . Alternatively, the line of argument used in Lemma 3.1 can be followed: letting v(x) be given by the function in (3.8), it is straightforward to check that this has appropriate boundary conditions and that the process $V_t(x) = v(D_t) + t$ is a martingale until the time that D first hits 0.

Note that $C(\lambda) = \min\{C_{\text{M}}(\lambda),C_{\text{S}}(\lambda)\}$ , and that $C_{\text{M}}(\lambda)$ and $C_{\text{S}}(\lambda)$ are both positive, decreasing functions of $\lambda$ , which agree precisely when $\lambda = \lambda^\star$ (recall the definition of $\lambda^\star$ in (1.4)). Thus, from any starting state $x\in(0,\pi]$ , ${\mathbb{E}}_x[T_{\text{M}}]<{\mathbb{E}}_{x}[T_{\text{S}}]$ if and only if $\lambda<\lambda^\star$ .

The proof of Theorem 1.2 (mean-optimality of the Mirror coupling for $\lambda<\lambda^\star$ ) is very similar to that given above for the LT-optimality of the Synchronous coupling. In this case it suffices to show that the function $v_{\text{M}}(x) = {\mathbb{E}}_x[T_{\text{M}}]$ satisfies the HJB equation

\[ \inf_{\theta,\, p} \big(1+\tfrac14{(1-\theta(x))}v^{\prime\prime}_{\text{M}}(x) + 2\lambda(1-p(x))\left\{v_{\text{M}}(\pi-x)-v_{\text{M}}(x)\right\}\big) = 0 ;\]

this is straightforward, and the details are omitted for the sake of brevity.

Figure 3 gives an impression of how the mean coupling time formulas from Lemma 3.3 behave as functions of the jump rate and the starting distance x. The left-hand plot shows that ${\mathbb{E}}_x[T_{\text{M}}] < {\mathbb{E}}_x[T_{\text{S}}]$ for all $x\in(0,\pi]$ when $\lambda=0.05<\lambda^\star$ , with this inequality being reversed for $\lambda=0.1>\lambda^\star$ . The right-hand plot graphs ${\mathbb{E}}_x[T_{\text{M}}]$ and ${\mathbb{E}}_x[T_{\text{S}}]$ as functions of $\lambda$ when $x\in\{\pi/4,\pi/2\}$ : for both starting distances the intersection of the two mean coupling times at $\lambda=\lambda^\star$ is evident.

Figure 3. Expected coupling time under the Synchronous and Mirror couplings. The left-hand plot shows the mean coupling time as a function of $D_0=x$ for three representative values of $\lambda$ . The right-hand plot shows the mean coupling time as a function of $\lambda$ for two particular values of $D_0$ .

Now let ${{{T_{\text{min}}}}}$ denote the coupling time of the mean-optimal co-adapted coupling (i.e. Mirror for $\lambda<\lambda^\star$ , Synchronous otherwise). As expected, when $\lambda=0$ we have ${\mathbb{E}}_{x}[T_{\text{min}}] = x(2\pi-x)$ , and when $\lambda\to\infty$ we obtain ${\mathbb{E}}_{x}[T_{\text{min}}] = x(\pi-x)$ (i.e. the mean time for a Brownian motion started at x to hit $\{0,2\pi\}$ or $\{0,\pi\}$ , respectively; recall the comment immediately following Theorem 1.2). Figure 4 plots ${\mathbb{E}}_{x}[T_{\text{min}}]={\mathbb{E}}_x[T_{\text{M}}]\wedge{\mathbb{E}}_x[T_{\text{S}}]$ as a function of the initial distance $x\in[0,\pi]$ , nicely showing the monotonic dependence upon the jump rate $\lambda$ .

Figure 4. Expected coupling time from $D_0=x$ when using the mean-optimal strategy.

Finally, let us compare ${{{T_{\text{min}}}}}$ to a maximal coupling. We know that the optimal co-adapted coupling is maximal when $\lambda=0$ ; given the construction of a maximal coupling in Section 2 it should come as no surprise that this is not the case for any $\lambda>0$ . From (2.1) we can easily calculate the Laplace transform and subsequently the mean of the maximal coupling time $T^\star$ when starting from opposite sides of the circle ( $x=\pi$ ):

(3.9) \begin{equation} \mathbb{E}_{\pi}[\text{e}^{-\gamma T^\star}] = \frac{1}{2\lambda+\gamma}\bigg(2\lambda + \frac{\gamma}{\cosh(\pi\sqrt{2(2\lambda+\gamma)})}\bigg) ; \quad {\mathbb{E}}_{\pi}[{T^\star}] = \frac{1-\text{sech}(2\pi\sqrt{\lambda})}{2\lambda} .\end{equation}

On the other hand, Lemma 3.3 tells us that ${{{T_{\text{min}}}}}$ satisfies

(3.10) \begin{equation} {\mathbb{E}}_{\pi}[T_{\text{min}}] = \min\left\{C_{\text{M}}(\lambda),C_{\text{S}}(\lambda)\right\}\frac{\sinh(\pi\sqrt{\lambda})}{\sqrt{\lambda}}.\end{equation}

Thus,

\[ \frac{{\mathbb{E}}_{\pi}[{T^\star}]}{{\mathbb{E}}_{\pi}[T_{\text{min}}]} = \begin{cases} (\pi\sqrt\lambda)^{-1}(2\cosh(\pi\sqrt{\lambda})-1)\sinh(\pi\sqrt{\lambda})\text{sech}(2\pi\sqrt{\lambda}), & \quad 0\le\lambda< \lambda^\star, \\ 1-\text{sech}(2\pi\sqrt{\lambda}), & \quad\lambda^\star\le \lambda , \end{cases}\]

and we note that for all values of $\lambda>0$ this ratio is strictly less than one. (For the case $\lambda<\lambda^\star$ , this follows from the observation that

\[ \frac{(2\cosh(z)-1)\sinh(z)\ \text{sech}(2z)}{z} = \frac{\int_{z}^{2z}\cosh(y)\,{{{\text{d}}}} y}{z\cosh(2z)} ,\]

along with the convexity of $\cosh$ .) Graphs of the expressions in (3.9) and (3.10) are shown in Figure 5: we see that the difference between the mean coupling times is largest when $\lambda = \lambda^\star$ .

Figure 5. Expectation of the mean-optimal co-adapted coupling time ( ${{{T_{\text{min}}}}}$ ) and maximal coupling time ( $T^\star$ ) when $D_0 = \pi$ .

Funding information

There are no funding bodies to thank relating to the creation of this article.

Competing interests

There were no competing interests to declare which arose during the preparation or publication process of this article.

References

Aldous, D. (1983). Random walks on finite groups and rapidly mixing Markov chains. In Séminaire de Probabilités XVII 1981/82, eds. J. Azéma and M. Yor. Springer, Berlin, pp. 243297.CrossRefGoogle Scholar
Banerjee, S. and Kendall, W. S. (2017). Rigidity for Markovian maximal couplings of elliptic diffusions. Prob. Theory Relat. Fields 168, 55112.CrossRefGoogle Scholar
Connor, S. B. (2013). Optimal coadapted coupling for a random walk on the hyper-complete graph. J. Appl. Prob. 50, 11171130.CrossRefGoogle Scholar
Connor, S. B. and Jacka, S. D. (2008). Optimal co-adapted coupling for the symmetric random walk on the hypercube. J. Appl. Prob. 45, 703713.CrossRefGoogle Scholar
Connor, S. B. and Merli, R. (2022). Optimal coupling of jumpy Brownian motion on the circle. Preprint, arXiv:2203.14791.Google Scholar
Goldstein, S. (1979). Maximal coupling. Z. Wahrscheinlichkeitsth. 46, 193204.CrossRefGoogle Scholar
Griffeath, D. (1975). A maximal coupling for Markov chains. Z. Wahrscheinlichkeitsth. 31, 95106.CrossRefGoogle Scholar
Hsu, E. P. and Sturm, K. T. (2013). Maximal coupling of Euclidean Brownian motions. Commun. Math. Statist. 1, 93104.CrossRefGoogle Scholar
Jacka, S. D., Mijatović, A. and Širaj, D. (2014). Mirror and synchronous couplings of geometric Brownian motions. Stoch. Process. Appl. 124, 10551069.CrossRefGoogle Scholar
Kendall, W. S. (2009). Brownian couplings, convexity, and shy-ness. Electron. Commun. Prob. 14, 6680.CrossRefGoogle Scholar
Kendall, W. S. (2015). Coupling, local times, immersions. Bernoulli 21, 10141046.CrossRefGoogle Scholar
Kuwada, K. (2007). On uniqueness of maximal coupling for diffusion processes with a reflection. J. Theor. Prob. 20, 935957.CrossRefGoogle Scholar
Lindvall, T. and Rogers, L. C. G. (1986). Coupling of multidimensional diffusions by reflection. Ann. Prob. 14, 860872.CrossRefGoogle Scholar
Merli, R. (2021). Probabilistic coupling: Mixing times and optimality. PhD thesis, University of York.Google Scholar
Øksendal, B. K. and Sulem, A. (2007). Applied Stochastic Control of Jump Diffusions. Springer, Berlin.CrossRefGoogle Scholar
Pitman, J. W. (1976). On coupling of Markov chains. Z. Wahrscheinlichkeitsth. 35, 315322.CrossRefGoogle Scholar
Rogers, L. C. G. and Williams, D. (2000). Diffusions, Markov Processes and Martingales, Vol. 2. Cambridge University Press.Google Scholar
Sverchkov, M. Y. and Smirnov, S. N. (1990). Maximal coupling of D-valued processes. Dokl. Math. 41, 352354.Google Scholar
Figure 0

Table 1. Summary of optimal coupling properties for different values of the jump rate $\lambda$.

Figure 1

Figure 1. The function $h(\lambda,\gamma)$ for four representative values of $\lambda$. For $\lambda<\lambda^\star$ the function takes first negative then positive values as $\gamma$ increases; for $\lambda\ge \lambda^\star$ the function is non-negative for all $\gamma$.

Figure 2

Figure 2. Difference between Laplace transforms corresponding to Mirror and Synchronous strategies: both graphs show $\mathbb{E}_{\pi}[\text{e}^{-\gamma{{{T_{\text{M}}}}}} - \text{e}^{-\gamma{{{T_{\text{S}}}}}}]$ as a function of $\gamma$, for different jump rates $\lambda\le\lambda^\star$. Neither coupling is LT-optimal when $\lambda<\lambda^\star$.

Figure 3

Figure 3. Expected coupling time under the Synchronous and Mirror couplings. The left-hand plot shows the mean coupling time as a function of $D_0=x$ for three representative values of $\lambda$. The right-hand plot shows the mean coupling time as a function of $\lambda$ for two particular values of $D_0$.

Figure 4

Figure 4. Expected coupling time from $D_0=x$ when using the mean-optimal strategy.

Figure 5

Figure 5. Expectation of the mean-optimal co-adapted coupling time (${{{T_{\text{min}}}}}$) and maximal coupling time ($T^\star$) when $D_0 = \pi$.