# Sticky nonlinear SDEs and convergence of McKean-Vlasov equations without confinement

Arnaud Guillin, Alain Durmus, Andreas Eberle, Katharina Schuh

<figure><img src="image_1.png" /></figure>

## To cite this version:

Arnaud Guillin, Alain Durmus, Andreas Eberle, Katharina Schuh. Sticky nonlinear SDEs and convergence of McKean-Vlasov equations without confinement. 2022. hal-03536620

HAL Id: hal-03536620

https://uca.hal.science/hal-03536620

Preprint submitted on 20 Jan 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Sticky nonlinear SDEs and convergence of McKean-Vlasov equations without confinement

Alain Durmus, Andreas Eberle, Arnaud Guillin, Katharina Schuh

November 19, 2021

## Abstract

We develop a new approach to study the long time behaviour of solutions to non-linear stochastic differential equations in the sense of McKean, as well as propagation of chaos for the corresponding mean-field particle system approximations. Our approach is based on a sticky coupling between two solutions to the equation. We show that the distance process between the two copies is dominated by a solution to a one-dimensional nonlinear stochastic differential equation with a sticky boundary at zero. This new class of equations is then analyzed carefully. In particular, we show that the dominating equation has a phase transition. In the regime where the Dirac measure at zero is the only invariant probability measure, we prove exponential convergence to equilibrium both for the one-dimensional equation, and for the original nonlinear SDE. Similarly, propagation of chaos is shown by a componentwise sticky coupling and comparison with a system of one dimensional nonlinear SDEs with sticky boundaries at zero. The approach applies to equations without confinement potential and to interaction terms that are not of gradient type.

## 1 Introduction

The main objective of this paper is to study and quantify convergence to equilibrium for McKean-Vlasov type nonlinear stochastic differential equations of the form

$$d\bar{X}_t = \left[ \int_{\mathbb{R}^d} b(\bar{X}_t - x)d\bar{\mu}_t(x) \right] dt + dB_t, \quad \bar{\mu}_t = \text{Law}(\bar{X}_t), \quad (1)$$

where $(B_t)_{t \ge 0}$ is a $d$-dimensional standard Brownian motion and $b: \mathbb{R}^d \to \mathbb{R}^d$ is a Lipschitz continuous function. This nonlinear SDE is the probabilistic counterpart of the *Fokker-Planck equation*

$$\frac{\partial}{\partial t} u_t = \nabla \cdot \left[ (1/2) \nabla u_t - (b * u_t) u_t \right], \quad (2)$$
which describes the time evolution of the density $u_t$ of $\bar{\mu}_t$ with respect to the Lebesgue measure on $\mathbb{R}^d$. Moreover, we also study uniform in time propagation of chaos for the approximating mean-field interacting particle systems

$$dX_t^{i,N} = \frac{1}{N} \sum_{j=1}^{N} b(X_t^{i,N} - X_t^{j,N})dt + dB_t^i, \quad i \in \{1, \dots, N\}, \quad (3)$$

with i.i.d. initial values $X_0^{1,N}, \dots, X_0^{N,N}$, and driven by independent $d$-dimensional Brownian motions $\{(B_t^i)_{t \ge 0}\}_{i=1}^N$. Our results are based on a new probabilistic approach relying on sticky couplings and comparison with solutions to a class of nonlinear stochastic differential equations on the real interval $[0, \infty)$ with a sticky boundary at 0. The study of this type of equations carried out below might also be of independent interest.

The equations (1) and (2) have been studied in many works. Often a slightly different setup is considered, where the interaction $b$ is assumed to be of gradient type, i.e., $b = -\nabla W$ for an interaction potential function $W: \mathbb{R}^d \to \mathbb{R}$, and an additional confinement potential function $V: \mathbb{R}^d \to \mathbb{R}$ satisfying $\lim_{|x| \to \infty} V(x) = \infty$ is included in the equations. The corresponding Fokker-Planck equation

$$\frac{\partial}{\partial t} u_t = \nabla \cdot \left[ (1/2) \nabla u_t + (\nabla V + \nabla W * u_t) u_t \right], \quad (4)$$

occurs for example in the modelling of granular media, see [58, 6] and the references therein. Existence and uniqueness of solutions to (1), (2) and (4) have been studied intensively. Introductions to this topic can be found for example in [27, 44, 45, 56], while recent results have been established in [48, 30]. Under appropriate conditions, it can be shown that the solutions converge to a unique stationary distribution at some given rate, see e.g. [13, 14, 10, 22, 19, 28]. In the case without confinement considered here, convergence to equilibrium of $(\bar{\mu}_t)_{t \ge 0}$ defined by (1) can only be expected for centered solutions, or after recentering around the center of mass of $\bar{\mu}_t$. It has first been analyzed in [13, 14] by an analytic approach and under the assumption that $b = -\nabla W$ for a convex function $W$. In particular, exponential convergence to equilibrium has been established under the strong convexity assumption Hess($W$) $\ge \rho$ Id for some $\rho > 0$, and polynomial convergence in the case where $W$ is only degenerately strictly convex. Similar results and some extensions have been derived in [42, 15] using a probabilistic approach.

Our first contribution aims at complementing these results, and extending them to non-convex interaction potentials and interaction functions that are not of gradient type. More precisely, suppose that

$$b(x) = -Lx + \gamma(x), \quad x \in \mathbb{R}^d, \quad (5)$$

where $L \in (0, \infty)$ is a positive real constant, and $\gamma: \mathbb{R}^d \to \mathbb{R}^d$ is a bounded function. Then we give conditions on $\gamma$ ensuring exponential convergence of centered solutions to (1) to a
unique stationary distribution in the standard L¹ Wasserstein metric. More generally, we show in Theorem 1 that under these conditions there exist constants $M, c \in (0, \infty)$ that depend only on $L$ and $\gamma$ such that if $(\bar{\mu}_t)_{t \ge 0}$ and $(\bar{\nu}_t)_{t \ge 0}$ are the marginal distributions of two solutions of (1), then for all $t \ge 0$,

$$W_1(\bar{\mu}_t, \bar{\nu}_t) \le Me^{-ct}W_1(\bar{\mu}_0, \bar{\nu}_0).$$

Using a coupling approach, related results have been derived in the previous works [22, 19] for the case where an additional confinement term is included in the equations. However, the arguments in these works rely on treating the equation with confinement and interaction term as a perturbation of the corresponding equation without interaction term, which has good ergodic properties. In the unconfined case this approach does not work, since the equation without interaction is transient and hence does not admit an invariant probability measure. Therefore, we have to develop a new approach for analyzing the equation without confinement.

Our approach is based on sticky couplings, an idea first developed in [23] to control the total variation distance between the marginal distributions of two non degenerate diffusion processes with identical noise but different drift coefficients. Since two solutions of (1) differ only in their drifts, we can indeed couple them using a sticky coupling in the sense of [23]. It can then be shown that the coupling distance process can be controlled by the solution $(r_t)_{t \ge 0}$ of a nonlinear SDE on $[0, \infty)$ with a sticky boundary at 0 of the form

$$dr_t = [\tilde{b}(r_t) + a\mathbb{P}(r_t > 0)]dt + 2\mathbb{1}_{(0,\infty)}(r_t)dW_t, \quad (6)$$

Here $\tilde{b}$ is a real-valued function on $[0, \infty)$ satisfying $\tilde{b}(0) = 0$, $a$ is a positive constant, and $(W_t)_{t \ge 0}$ is a one-dimensional standard Brownian motion. Solutions to SDEs with diffusion coefficient $r \mapsto \mathbb{1}_{(0,\infty)}(r)$, as in (6), have a sticky boundary at 0, i.e., if the drift at 0 is strictly positive, then the set of all time points $t \in [0, \infty)$ such that $r_t = 0$ is a fractal set with strictly positive Lebesgue measure that does not contain any open interval. Sticky SDEs have attracted wide interest, starting from [25, 26] in the one-dimensional case. Multivariate extensions have been considered in [34, 60, 61] building upon results obtained in [43, 52, 53], while corresponding martingale problems have been investigated in [54]. Note that in general no strong solution for this class of SDEs exists as illustrated in [16]. We refer to [24, 3] and the references therein for recent contributions on this topic. Note, however, that in contrast to standard sticky SDEs, the equation (6) is nonlinear in the sense of McKean. We are not aware of previous studies of such nonlinear sticky equations, which seems to be a very interesting topic on its own.

Intuitively, one would expect that as time evolves, more mass gets stuck at 0, i.e., $\mathbb{P}(r_t > 0)$ decreases. As a consequence, the drift at 0 in Equation (6) decreases, which again forces even more mass to get stuck at 0. Therefore, one could hope that $\mathbb{P}(r_t = 0)$ converges to 1 as $t \to \infty$. On the other hand, if $a$ is too large then the drift at 0 might be too strong so that not all of the mass gets stuck at 0 eventually. This indicates that there might be a phase
transition for the nonlinear sticky SDE depending on the size of the constant a compared to $\tilde{b}$. In Section 3, we prove rigorously that this intuition is correct. Under appropriate conditions on $\tilde{b}$, we show at first that existence and uniqueness in law holds for solutions of (6). Then we prove that for a sufficiently small, the Dirac measure at 0 is the unique invariant probability measure, and geometric ergodicity holds. As a consequence, under corresponding assumptions, the sticky coupling approach yields exponential convergence to equilibrium for the original nonlinear SDE (1). On the other hand, we prove the existence of multiple invariant probability measures for (6) if the smallness condition on a is not satisfied. Our results for (1) can also be adapted to deal with nonlinear SDEs over the torus $T = \mathbb{R}/(2\pi\mathbb{Z})$, as considered in [18]. As an example, we discuss the application to the Kuramoto model for which a more explicit analysis is available [1, 7, 8, 12].

Finally, in addition to studying the long-time behaviour of the nonlinear SDE (1), we are also interested in establishing propagation of chaos for the mean-field particle system approximation (3). The propagation of chaos phenomenon first introduced by Kac [37] describes the convergence of the empirical measure of the mean-field particle system (3) to the solution (1). More precisely, in [56, 45] it has been shown under weak assumptions on W that for i.i.d. initial laws, the random variables $X_t^{i,N}$, $i \in \{1, \dots, N\}$, become asymptotically independent as $N \to \infty$, and the common law $\mu_t^N$ of each of these random variables converges to $\bar{\mu}_t$. However, the original results are only valid uniformly over a finite time horizon. Quantifying the convergence uniformly for all times $t \in \mathbb{R}_+$ is an important issue. The case with a confinement potential has been studied for example in [19], see also the references therein. Again, the case when there is only interaction is more difficult. Malrieu [42] seems the first to consider the case without confinement. By applying a synchronous coupling, he proved uniform in time propagation of chaos for strongly convex interaction potentials. Later on, assuming that the interaction potential is losing strict convexity only in a finite number of points (e.g., $W(x) = |x|^3$), Cattiaux, Guillin and Malrieu [15] have shown uniform in time propagation of chaos with a rate getting worse with the degeneracy in convexity. In a very recent work, Delarue and Tse [17] prove uniform in time weak propagation of chaos (i.e., observable by observable) on the torus via Lions derivative methods. Remarkably, their results are not limited to the unique invariant measure case.

Our contribution is in the same vein using probabilistic tools in place of analytic ones. We endow the space $\mathbb{R}^N$ consisting of N particle configurations $x = (x^i)_{i=1}^N$ with the semi-metric $l^1 \circ \pi$, where

$$l^1(x, y) = \frac{1}{N} \sum_{i=1}^{N} |x^i - y^i| \qquad (7)$$

is a normalized $l^1$-distance between configurations $x, y \in \mathbb{R}^N$, and

$$\pi(x, y) = \left( \left( x^i - \frac{1}{N} \sum_{j=1}^{N} x^j \right)_{i=1}^{N}, \left( y^i - \frac{1}{N} \sum_{j=1}^{N} y^j \right)_{i=1}^{N} \right), \qquad (8)$$
is a projection from $\mathbb{R}^{Nd} \times \mathbb{R}^{Nd}$ to the subspace $\mathbb{H}_N \times \mathbb{H}_N$, where

$$ \mathbb{H}_N = \{x \in \mathbb{R}^{Nd} : \sum_{i=1}^{N} x^i = 0\}. \qquad (9) $$

Let $W_{l^{1 \circ \pi}}$ denote the $L^1$ Wasserstein semimetric on probability measures on $\mathbb{R}^{Nd}$ corresponding to the cost function $l^1 \circ \pi$. Then under assumptions stated below, we prove uniform in time propagation of chaos for the mean-field particle system in the following sense: Suppose that $(X_t^{1,N}, \dots, X_t^{N,N})_{t \ge 0}$ is a solution of (3) such that $X_0^{1,N}, \dots, X_0^{N,N}$ are i.i.d. with distribution $\bar{\mu}_0$ having finite second moment. Let $\nu_t^N$ denote the joint law of the random variables $X_t^{i,N}$, $i \in \{1, \dots, N\}$, and let $\bar{\mu}_t$ denote the law of the solution of (1) with initial law $\bar{\mu}_0$. Then there exists a constant $C \in [0, \infty)$ such that for any $N \in \mathbb{N}$,

$$ \sup_{t \ge 0} W_{l^{1 \circ \pi}}(\bar{\mu}_t^{\otimes N}, \nu_t^N) \le CN^{-1/2}. \qquad (10) $$

The proof is based on a componentwise sticky coupling, and a comparison of the coupling difference process with a system of one-dimensional sticky nonlinear SDEs.

The paper is organised as follows. In Section 2, we state our main results regarding the long-time behaviour of (1). The main results on one-dimensional nonlinear SDEs with a sticky boundary at zero are stated in Section 3. Sections 4 and 5 contain the corresponding results on uniform (in time) propagation of chaos and mean-field systems of sticky SDEs. All the proofs are given in Section 6. In Appendix A, we carry the results over to nonlinear sticky SDEs over $\mathbb{T}$ and consider the application to the Kuramoto model.

**Notation** The Euclidean norm on $\mathbb{R}^d$ is denoted by $| \cdot |$. For $x \in \mathbb{R}$, we write $x_+ = \max(0, x)$. For some space $\mathcal{X}$, which here is either $\mathbb{R}^d$, $\mathbb{R}^{Nd}$ or $\mathbb{R}_+$, we denote its Borel $\sigma$-algebra by $\mathcal{B}(\mathcal{X})$. The space of all probability measures on $(\mathcal{X}, \mathcal{B}(\mathcal{X}))$ is denoted by $\mathcal{P}(\mathcal{X})$. Let $\mu, \nu \in \mathcal{P}(\mathcal{X})$. A coupling $\xi$ of $\mu$ and $\nu$ is a probability measure on $(\mathcal{X} \times \mathcal{X}, \mathcal{B}(\mathcal{X}) \otimes \mathcal{B}(\mathcal{X}))$ with marginals $\mu$ and $\nu$. $\Gamma(\mu, \nu)$ denotes the set of all couplings of $\mu$ and $\nu$. The $L^1$ Wasserstein distance with respect to a distance function $d: \mathcal{X} \times \mathcal{X} \to \mathbb{R}_+$ is defined by

$$ W_d(\mu, \nu) = \inf_{\xi \in \Gamma(\mu, \nu)} \int_{\mathcal{X} \times \mathcal{X}} d(x, y)\xi(dx dy). $$

We write $W_1$ if the underlying distance function is the Euclidean distance.

We denote by $\mathcal{C}(\mathbb{R}_+, \mathcal{X})$ the set of continuous functions from $\mathbb{R}_+$ to $\mathcal{X}$, and by $\mathcal{C}^2(\mathbb{R}_+, \mathcal{X})$ the set of twice continuously differentiable functions.

Consider a probability space $(\Omega, \mathcal{A}, P)$ and a measurable function $r: \Omega \to \mathcal{C}(\mathbb{R}_+, \mathcal{X})$. Then $P = P \circ r^{-1}$ denotes the law on $\mathcal{C}(\mathbb{R}_+, \mathcal{X})$, and $P_t = P \circ r_t^{-1}$ the marginal law on $\mathcal{X}$ at time $t$.
## 2 Long-time behaviour of McKean-Vlasov diffusions

We establish our results regarding (1) and (3) under the following assumption on b.

**B 1.** The function $b: \mathbb{R}^d \to \mathbb{R}^d$ is Lipschitz continuous and anti-symmetric, i.e., $b(z) = -b(-z)$, and there exist $L \in (0, \infty)$, a function $\gamma: \mathbb{R}^d \to \mathbb{R}^d$ and a Lipschitz continuous function $\kappa: [0, \infty) \to \mathbb{R}$ such that

$$b(z) = -Lz + \gamma(z) \quad \text{for all } z \in \mathbb{R}^d, \qquad (11)$$

and the following conditions are satisfied for all $x, \tilde{x}, y \in \mathbb{R}^d$:

$$\langle x - y, \gamma(x - \tilde{x}) - \gamma(y - \tilde{x}) \rangle \le \kappa(|x - y|)|x - y|^2, \qquad (12)$$

and

$$\limsup_{r \to \infty} (\kappa(r) - L) < 0. \qquad (13)$$

Let $\bar{b}(r) = (\kappa(r) - L)r$. If (13) holds, then there exist $R_0, R_1 \ge 0$ such that for

$$\bar{b}(r) < 0, \qquad \text{for any } r \ge R_0, \qquad (14)$$

$$\bar{b}(r)/r \le -4/[R_1(R_1 - R_0)], \qquad \text{for any } r \ge R_1. \qquad (15)$$

In addition, we assume

**B 2.**

$$\|\gamma\|_{\infty} \le (4 \int_{0}^{R_1} \exp \left( \frac{1}{2} \int_{0}^{s} \bar{b}(r) + dr \right) ds)^{-1}.$$

We consider the following condition on the initial distribution.

**B 3.** The initial distribution $\mu_0$ satisfies $\int_{\mathbb{R}^d} \|x\|^4 \mu_0(dx) < +\infty$ and $\int_{\mathbb{R}^d} x \mu_0(dx) = 0$.

Note that under conditions **B1** and **B3**, unique strong solutions $(\bar{X}_t)_{t \ge 0}$ and $(\{X_t^{i,N}\}_{i=1}^N)_{t \ge 0}$ exist for (1) and (3), see e.g. [15, Theorem 2.6]. In addition, note that since $b$ is assumed to be anti-symmetric, by an easy localisation argument, we get that $d\mathbb{E}[\bar{X}_t]/dt = \mathbb{E}[b * \mu_t(\bar{X}_t)] = 0$ and $d\mathbb{E}[N^{-1} \sum_{i=1}^N X_t^{i,N}]/dt = 0$. Thus, if $\bar{X}_0$ and $\{X_0^{i,N}\}_{i=1}^N$ have distribution $\mu_0$ and $\mu_0^{\otimes N}$, respectively, with $\mu_0$ satisfying **B 3**, then it holds $\mathbb{E}[\bar{X}_t] = 0$ and $\mathbb{E}[N^{-1} \sum_{i=1}^N X_t^{i,N}] = 0$ for all $t \ge 0$.

Suppose $f: \mathbb{R}_+ \to \mathbb{R}_+$ is an increasing, concave function vanishing at zero. Then $d(x, y) = f(|x - y|)$ defines a distance. The corresponding $L^1$ Wasserstein distance is denoted by $W_f$. Note that in the case $f(t) = t$ for any $t \ge 0$, $W_f$ is simply $W_1$.
**Theorem 1** (Contraction for nonlinear SDE). Assume **B**1. Let $\bar{\mu}_0, \bar{\nu}_0$ be probability measures on $(\mathbb{R}^d, \mathcal{B}(\mathbb{R}^d))$ satisfying **B**3. For any $t \ge 0$, let $\bar{\mu}_t$ and $\bar{\nu}_t$ denote the laws of $\bar{X}_t$ and $\bar{Y}_t$ where $(\bar{X}_s)_{s \ge 0}$ and $(\bar{Y}_s)_{s \ge 0}$ are solutions of (1) with initial distribution $\bar{\mu}_0$ and $\bar{\nu}_0$, respectively. Then, for all $t \ge 0$,

$$ W_f(\bar{\mu}_t, \bar{\nu}_t) \le e^{-\bar{c}t} W_f(\bar{\mu}_0, \bar{\nu}_0) \quad \text{and} \quad W_1(\bar{\mu}_t, \bar{\nu}_t) \le M_1 e^{-\bar{c}t} W_1(\bar{\mu}_0, \bar{\nu}_0), \quad (16) $$

where the function $f$ is defined by (37) and the constants $\bar{c}$ and $M_1$ are given by

$$ \bar{c}^{-1} = 2 \int_{0}^{R_1} \int_{0}^{s} \exp \left( \frac{1}{2} \int_{r}^{s} \bar{b}(u)_{+} \, du \right) \mathrm{d}r \, \mathrm{d}s, \quad (17) $$

$$ M_1 = 2 \exp \left( \frac{1}{2} \int_{0}^{R_0} \bar{b}(s)_{+} \mathrm{d}s \right). \quad (18) $$

*Proof.* The proof is postponed to Section 6.2.2. □

The construction and definition of the underlying distance function $f(|x-y|)$ mentioned in Theorem 1 is based on the one introduced by [20].

To prove Theorem 1 we use a coupling $(\bar{X}_t, \bar{Y}_t)_{t \ge 0}$ of two copies of solutions to the nonlinear stochastic differential equation (1) with different initial conditions. The coupling $(\bar{X}_t, \bar{Y}_t)_{t \ge 0}$ will be defined as the weak limit of a family of couplings $(\bar{X}_t^\delta, \bar{Y}_t^\delta)_{t \ge 0}$, parametrized by $\delta > 0$. Roughly, this family is mixture of synchronous and reflection couplings and can be described as follows. For $\delta > 0$, $(\bar{X}_t^\delta, \bar{Y}_t^\delta)_{t \ge 0}$ behaves like a reflection coupling if $|\bar{X}_t^\delta - \bar{Y}_t^\delta| \ge \delta$, and like a synchronous coupling if $|\bar{X}_t^\delta - \bar{Y}_t^\delta| = 0$. For $|\bar{X}_t^\delta - \bar{Y}_t^\delta| \in (0, \delta)$ we take an interpolation of synchronous and reflection coupling. We argue that the family of couplings $\{(\bar{X}_t^\delta, \bar{Y}_t^\delta)_{t \ge 0} : \delta > 0\}$ is tight and that a subsequence $\{(\bar{X}_t^{\delta_n}, \bar{Y}_t^{\delta_n})_{t \ge 0} : n \in \mathbb{N}\}$ converges to a limit $(\bar{X}_t, \bar{Y}_t)_{t \ge 0}$. This limit is a coupling which we call the sticky coupling associated to (1).

To carry out the construction rigorously, we take two Lipschitz continuous functions $\mathrm{rc}^\delta, \mathrm{sc}^\delta : \mathbb{R}_+ \to [0, 1]$ for $\delta > 0$ such that

$$ \mathrm{rc}^{\delta}(0) = 0, \quad \mathrm{rc}^{\delta}(r) = 1 \text{ for } r \ge \delta, \quad \mathrm{rc}^{\delta}(r) > 0 \text{ for } r > 0 \text{ and } \mathrm{rc}^{\delta}(r)^2 + \mathrm{sc}^{\delta}(r)^2 = 1 \text{ for } r \ge 0. \quad (19) $$

Further, we assume that there exists $\epsilon_0 > 0$ such that for any $\delta \le \epsilon_0$, $\mathrm{rc}^\delta$ satisfies

$$ \mathrm{rc}^{\delta}(r) \ge \frac{\|\gamma\|_{\mathrm{Lip}}}{2\|\gamma\|_{\infty}} r \quad \text{for any } r \in (0, \delta), \quad (20) $$

where $\|\gamma\|_{\mathrm{Lip}} < \infty$ denotes the Lipschitz norm of $\gamma$. This assumption is satisfied for example if $\mathrm{rc}^\delta(r) = \sin((\pi/2\delta)r)1_{r<\delta} + 1_{r \ge \delta}$ and $\mathrm{sc}^\delta(r) = \cos((\pi/2\delta)r)1_{r<\delta}$ with $\delta \le \epsilon_0 = 2\|\gamma\|_{\infty}/\|\gamma\|_{\mathrm{Lip}}$.

Let $(B_t^1)_{t \ge 0}$ and $(B_t^2)_{t \ge 0}$ be two $d$-dimensional Brownian motions. We define the coupling $(\bar{X}_t^\delta, \bar{Y}_t^\delta)_{t \ge 0}$ as a process in $\mathbb{R}^{2d}$ satisfying the following nonlinear stochastic differential
equation

$$ \begin{aligned} d\bar{X}_t^\delta &= b * \bar{\mu}_t^\delta(\bar{X}_t^\delta)dt + \mathrm{rc}^\delta(\bar{r}_t^\delta)dB_t^1 + \mathrm{sc}^\delta(\bar{r}_t^\delta)dB_t^2, & \bar{\mu}_t^\delta &= \mathrm{Law}(\bar{X}_t^\delta), \\ d\bar{Y}_t^\delta &= b * \bar{\nu}_t^\delta(\bar{Y}_t^\delta)dt + \mathrm{rc}^\delta(\bar{r}_t^\delta)(\mathrm{Id} - 2\bar{e}_t^\delta(\bar{e}_t^\delta)^T)dB_t^1 + \mathrm{sc}^\delta(\bar{r}_t^\delta)dB_t^2, & \bar{\nu}_t^\delta &= \mathrm{Law}(\bar{Y}_t^\delta) \end{aligned} \quad (21) $$

with initial condition $(\bar{X}_0^\delta, \bar{Y}_0^\delta) = (x_0, y_0)$. Here we set $\bar{Z}_t^\delta = \bar{X}_t^\delta - \bar{Y}_t^\delta$, $\bar{r}_t^\delta = |\bar{Z}_t^\delta|$ and $\bar{e}_t^\delta = \bar{Z}_t^\delta / \bar{r}_t^\delta$ if $\bar{r}_t^\delta \neq 0$. For $\bar{r}_t^\delta = 0$, $\bar{e}_t^\delta$ is some arbitrary unit vector, whose exact choice is irrelevant since $\mathrm{rc}^\delta(0) = 0$.

**Theorem 2.** Assume B1. Let $\bar{\mu}_0$ and $\bar{\nu}_0$ be probability measures on $(\mathbb{R}^d, \mathcal{B}(\mathbb{R}^d))$ satisfying B3. Then, $(\bar{X}_t, \bar{Y}_t)_{t \ge 0}$ is the subsequential limit as $\delta \to 0$ of $\{(\bar{X}_t^\delta, \bar{Y}_t^\delta)_{t \ge 0} : \delta > 0\}$ where $(\bar{X}_t)_{t \ge 0}$ and $(\bar{Y}_t)_{t \ge 0}$ are solutions of (1) with initial distribution $\bar{\mu}_0$ and $\bar{\nu}_0$. Further, there exists a process $(r_t)_{t \ge 0}$ defined on the same probability space as $(\bar{X}_t, \bar{Y}_t)_{t \ge 0}$ satisfying for any $t \ge 0$, $|\bar{X}_t - \bar{Y}_t| \le r_t$ almost surely and which is a weak solution of

$$ dr_t = (\bar{b}(r_t) + 2\|\gamma\|_{\infty}\mathbb{P}(r_t > 0))dt + 2\mathbb{1}_{(0,\infty)}(r_t)d\bar{W}_t, \quad (22) $$

where $(\bar{W}_t)_{t \ge 0}$ is a one-dimensional Brownian motion.

*Proof.* The proof is postponed to Section 6.2.1. □

Therefore, next we study sticky nonlinear SDEs given by (6).

## 3 Nonlinear SDEs with sticky boundaries

Consider nonlinear SDEs with a sticky boundary at 0 of the form

$$ dr_t = (\tilde{b}(r_t) + P_t(g))dt + 2\mathbb{1}_{(0,\infty)}(r_t)dW_t, \quad P_t = \mathrm{Law}(r_t), \quad (23) $$

where $\tilde{b}: [0, \infty) \to \mathbb{R}$ is some continuous function and $P_t(g) = \int_{\mathbb{R}_+} g(r)P_t(dr)$ for some measurable function $g: [0, \infty) \to \mathbb{R}$.

In this section we establish existence, uniqueness in law and comparison results for solutions of (6). Consider a filtered probability space $(\Omega, \mathcal{A}, (\mathcal{F}_t), P)$ and a probability measure $\mu$ on $\mathbb{R}_+$. We call an $(\mathcal{F}_t)$ adapted process $(r_t, W_t)_{t \ge 0}$ a weak solution of (23) with initial distribution $\mu$ if the following holds: $\mu = P \circ r_0^{-1}$, the process $(W_t)_{t \ge 0}$ is a one-dimensional $(\mathcal{F}_t)$ Brownian motion w.r.t. $P$, the process $(r_t)_{t \ge 0}$ is non-negative and continuous, and satisfies almost-surely

$$ r_t - r_0 = \int_0^t (\tilde{b}(r_s) + P_s(g))ds + \int_0^t 2 \cdot \mathbb{1}_{(0,\infty)}(r_s)dW_s, \quad \text{for } t \in \mathbb{R}_+. $$

Note that the sticky nonlinear SDE given in (6) is a special case of (23) with $g(r) = a\mathbb{1}_{(0,\infty)}(r)$ since $\mathbb{P}(r_t > 0) = \int_{\mathbb{R}_+} \mathbb{1}_{(0,\infty)}(y)dP_t(y)$ with $P_t = P \circ r_t^{-1}$.
## 3.1 Existence, uniqueness in law, and a comparison result

Let $\mathbb{W} = C(\mathbb{R}_+, \mathbb{R})$ be the space of continuous functions endowed with the topology of uniform convergence on compact sets, and let $\mathcal{B}(\mathbb{W})$ be the corresponding Borel $\sigma$-algebra. Suppose $(r_t, W_t)_{t \ge 0}$ is a solution of (23) on $(\Omega, \mathcal{A}, P)$, then we denote by $\mathbb{P} = P \circ r^{-1}$ its law on $(\mathbb{W}, \mathcal{B}(\mathbb{W}))$. We say that uniqueness in law holds for (23) if for any two solutions $(r_t^1)_{t \ge 0}$ and $(r_t^2)_{t \ge 0}$ of (23) with the same initial law, the distributions of $(r_t^1)_{t \ge 0}$ and $(r_t^2)_{t \ge 0}$ on $(\mathbb{W}, \mathcal{B}(\mathbb{W}))$ are equal.

We impose the following assumptions on $\tilde{b}, g$ and the initial condition $\mu$:

**H 1.** $\tilde{b}$ is a Lipschitz continuous function with Lipschitz constant $\tilde{L}$ and $\tilde{b}(0) = 0$.

**H 2.** $g$ is a left-continuous, non-negative, non-decreasing and bounded function.

**H 3.** There exists $p > 2$ such that the $p$-th order moment of the law $\mu$ is finite.

Note that for (6), the condition **H 2** is satisfied if $a$ is a positive constant. It follows from **H 1** and **H 2** that there is a constant $C < \infty$ such that for all $r \in \mathbb{R}_+$, the following linear growth condition holds,

$$ \tilde{b}(r) + \sup_{p \in \mathcal{P}(\mathbb{R}_+)} p(g) \le C(1 + |r|). \quad (24) $$

In order to get a solution to (23) on $\mathbb{R}_+$ we extend the function $\tilde{b}$ to $\mathbb{R}$ by setting $\tilde{b}(r) = 0$ for $r < 0$. Note that any solution $(r_t)_{t \ge 0}$ with initial distribution supported on $\mathbb{R}_+$ satisfies almost surely $r_t \ge 0$ for all $t \ge 0$. This follows from Ito-Tanaka formula applied to $F(r) = 1_{(-\infty,0)}(r)r$, cf. [51, Chapter 6, Theorem 1.1]. Indeed

$$ \begin{aligned} 1_{(-\infty,0)}(r_t)r_t &= \int_0^t 1_{(-\infty,0)}(r_s)dr_s - \frac{1}{2}\ell_t^0(r) \\ &= \int_0^t 1_{(-\infty,0)}(r_s)(\tilde{b}(r_s) + P_s(g))ds + \int_0^t 1_{(-\infty,0)}21_{(0,\infty)}(r_s)dW_s - \frac{1}{2}\ell_t^0(r) \\ &= \int_0^t 1_{(-\infty,0)}P_s(g)ds > 0, \end{aligned} $$

where $\ell_t^0(r)$ is the local time at 0, which vanishes, since $d[r]_s = 1_{(0,\infty)}(r_s)ds$.

Existence and uniqueness in law of (23) is a direct consequence of a stronger result that we now introduce. To study existence and uniqueness and to compare two solutions of (23) with different drifts, we establish existence of a synchronous coupling of two copies of (23),

$$ \begin{aligned} dr_t &= (\tilde{b}(r_t) + P_t(g))dt + 21_{(0,\infty)}(r_t)dW_t, \\ ds_t &= (\hat{b}(s_t) + \hat{P}_t(h))dt + 21_{(0,\infty)}(s_t)dW_t, \quad \text{Law}(r_0, s_0) = \eta, \end{aligned} \quad (25) $$

where $P_t = P \circ r_t^{-1}$, $\hat{P}_t = P \circ s_t^{-1}$, $(W_t)_{t \ge 0}$ is a Brownian motion and where $\eta \in \Gamma(\mu, \nu)$ for $\mu, \nu \in \mathcal{P}(\mathbb{R}_+)$.
**Theorem 3.** Suppose that $(\tilde{b}, g)$ and $(\hat{b}, h)$ satisfy $\mathbf{H}1$ and $\mathbf{H}2$. Let $\eta \in \Gamma(\mu, \nu)$ where the probability measures $\mu$ and $\nu$ on $\mathbb{R}_+$ satisfy $\mathbf{H}3$. Then there exists a weak solution $(r_t, s_t)_{t \ge 0}$ of the sticky stochastic differential equation (25) with initial distribution $\eta$ defined on a probability space $(\Omega, \mathcal{A}, P)$ with values in $(\mathbb{W} \times \mathbb{W}, \mathcal{B}(\mathbb{W}) \otimes \mathcal{B}(\mathbb{W}))$. If additionally,

$$ \tilde{b}(r) \le \hat{b}(r) \quad \text{and} \quad g(r) \le h(r) \quad \text{for any } r \in \mathbb{R}_+, \text{ and} \\ P[r_0 \le s_0] = 1, $$

then $P[r_t \le s_t \text{ for all } t \ge 0] = 1$.

*Proof.* The proof is postponed to Section 6.3.1. □

**Remark 4.** We note that by the comparison result we can deduce uniqueness in law for the solution of (23).

## 3.2 Invariant measures and phase transition

Under the following conditions on the drift function $\tilde{b}$ we exhibit a phase transition phenomenon for the model (6).

**Theorem 5.** Suppose $\mathbf{H}1$ holds and $\limsup_{r \to \infty} (r^{-1}\tilde{b}(r)) < 0$. Then, the Dirac measure at 0, $\delta_0$, is an invariant probability measure for (6). If there exists $p \in (0, 1)$ solving

$$ (2/a) = (1-p)I(a,p) \quad (26) $$

with

$$ I(a,p) = \int_{0}^{\infty} \exp \left( \frac{1}{2} apx + \frac{1}{2} \int_{0}^{x} \tilde{b}(r) dr \right) dx, \quad (27) $$

then the probability measure $\pi$ on $[0, \infty)$ given by

$$ \pi(dx) \propto \frac{1}{Z} \left( \frac{2}{ap} \delta_0(dx) + \exp \left( \frac{1}{2} apx + \frac{1}{2} \int_0^x \tilde{b}(r) dr \right) \lambda_{(0,\infty)}(dx) \right) \quad (28) $$

is another invariant probability measure for (6).

*Proof.* The proof is postponed to Section 6.3.2. □

In our next result we specify a necessary and sufficient condition for the existence of a solution of (26).

**Proposition 6.** Suppose that $\tilde{b}(r)$ in (6) is of the form $\tilde{b}(r) = -\tilde{L}r$ with constant $a \tilde{L} > 0$. If $a/\sqrt{\tilde{L}} > 2/\sqrt{\pi}$, then there exists a unique $\hat{p}$ solving (27). In particular, the Dirac measure $\delta_0$ and the measure $\pi$ given in (28) with $\hat{p}$ are invariant measures for (6). On the other hand, if $a/\sqrt{\tilde{L}} \le 2/\sqrt{\pi}$, then there exists no $\hat{p}$ solving (27).

*Proof.* The proof is postponed to Section 6.3.2. □
### 3.3 Convergence for sticky nonlinear SDEs

Under $H1$ and the following additional assumption we establish geometric convergence in Wasserstein distance for the marginal law of the solution $r_t$ of (6) to the Dirac measure at 0:

**H4.** It holds $\limsup_{r \to \infty} (r^{-1}\tilde{b}(r)) < 0$ and $a \le (2 \int_0^{\tilde{R}_1} \exp(\frac{1}{2} \int_0^s \tilde{b}(u)_+ du) ds)^{-1}$ with $\tilde{R}_0, \tilde{R}_1$ defined by

$$ \tilde{R}_0 = \inf\{s \in \mathbb{R}_+ : \tilde{b}(r) \le 0 \ \forall r \ge s\} \quad \text{and} \quad (29) $$

$$ \tilde{R}_1 = \inf\{s \ge \tilde{R}_0 : -\frac{s}{r}(s - \tilde{R}_0)\tilde{b}(r) \ge 4 \ \forall r \ge s\}. \quad (30) $$

**Theorem 7.** Suppose $H1$ and $H4$ holds. Then, the Dirac measure at 0, $\delta_0$, is the unique invariant probability measure of (6). Moreover if $(r_s)_{s \ge 0}$ is a solution of (6) with $r_0$ distributed with respect to an arbitrary probability measure $\mu$ on $(\mathbb{R}_+, \mathcal{B}(\mathbb{R}_+))$, it holds for all $t \ge 0$,

$$ \mathbb{E}[f(r_t)] \le e^{-ct}\mathbb{E}[f(r_0)], \quad (31) $$

where $f$ and $c$ are given by (37) and (36) with $a$ and $\tilde{b}$ given in (6) and $\tilde{R}_0$ and $\tilde{R}_1$ given in (29) and (30).

*Proof.* The proof is postponed to Section 6.3.3. □

## 4 Uniform in time propagation of chaos

To prove uniform in time propagation of chaos, we consider the $L^1$ Wasserstein distance with respect to the cost function $\bar{f}_N \circ \pi : \mathbb{R}^{Nd} \times \mathbb{R}^{Nd} \to \mathbb{R}_+$ with $\pi$ given in (8), and $\bar{f}_N$ given by

$$ \bar{f}_N((x^{i,N})_{i=1}^N, (y^{i,N})_{i=1}^N) = \frac{1}{N} \sum_{i=1}^N f(|x^i - y^i|), \quad (32) $$

with $f: \mathbb{R}_+ \to \mathbb{R}_+$ defined in (37). This distance is denoted by $W_{f,N}$. Note that $\bar{f}_N$ is equivalent to $l^1$ defined in (7).

We note that since $\pi$ defines a projection from $\mathbb{R}^{Nd}$ to the hyperplane $H_N \subset \mathbb{R}^{Nd}$ given in (9), for $\hat{\mu}$ and $\hat{\nu}$ on $H_N$, $W_{f,N}(\hat{\mu}, \hat{\nu})$ coincides with the Wasserstein distance given by

$$ \hat{W}_{f,N}(\hat{\mu}, \hat{\nu}) = \inf_{\xi \in \Gamma(\hat{\mu}, \hat{\nu})} \int_{H_N \times H_N} \bar{f}_N(x, y) \xi(dx dy) \quad (33) $$

and $W_{l^1 \circ \pi}(\hat{\mu}, \hat{\nu}) = \hat{W}_{l^1}(\hat{\mu}, \hat{\nu})$, where $\bar{f}_N$ and $l^1$ are given in (32) and (7), respectively, and where $\hat{W}_{l^1}(\hat{\mu}, \hat{\nu})$ is defined as in (33) with respect to the distance $l^1$.
**Theorem 8** (Uniform in time propagation of chaos). Let $N \in \mathbb{N}$ and assume **B**1. Let $\bar{\mu}_0$ and $\nu_0$ be probability measures on $(\mathbb{R}^d, \mathcal{B}(\mathbb{R}^d))$ satisfying **B**3. For $t \ge 0$, denote by $\bar{\mu}_t$ and $\nu_t^N$ the law of $\bar{X}_t$ and $\{X_t^{i,N}\}_{i=1}^N$ where $(\bar{X}_s)_{s \ge 0}$ and $(\{X_s^{i,N}\}_{i=1}^N)_{s \ge 0}$ are solutions of (1) and (3), respectively, with initial distributions $\bar{\mu}_0$ and $\nu_0^{\otimes N}$. Then for all $t \ge 0$,

$$
\begin{align*}
W_{f,N}(\bar{\mu}_t^{\otimes N}, \nu_t^N) &\le e^{-\bar{c}t} W_{f,N}(\bar{\mu}_0^{\otimes N}, \nu_0^{\otimes N}) + \tilde{C}\tilde{c}^{-1}N^{-1/2}, \\
W_{l^{1\circ\pi}}(\bar{\mu}_t^{\otimes N}, \nu_t^N) &\le M_1 e^{-\bar{c}t} W_{l^{1\circ\pi}}(\bar{\mu}_0^{\otimes N}, \nu_0^{\otimes N}) + M_1 \tilde{C}\tilde{c}^{-1}N^{-1/2},
\end{align*}
$$

where $f$ is defined by (37), $M_1$ by (18), $\tilde{c}$ by (17) and $\tilde{C}$ is a finite constant depending on $\|\gamma\|_{\infty}$, $L$ and the second moment of $\bar{\mu}_0$ and given in (58).

*Proof.* The proof is postponed to Section 6.2.3. □

**Remark 9.** Denote by $\mu_t^N$ and $\nu_t^N$ the distribution of $\{X_t^{i,N}\}_{i=1}^N$ and $\{Y_t^{i,N}\}_{i=1}^N$ where the two processes $(\{X_s^{i,N}\}_{s=1}^N)_{s \ge 0}$ and $(\{Y_s^{i,N}\}_{s=1}^N)_{s \ge 0}$ are solutions of (3) with initial probability distributions $\mu_0^N$, $\nu_0^N \in \mathcal{P}(\mathbb{R}^{Nd})$, respectively, with finite forth moment. An easy inspection and adaptation of the proof of Theorem 8 show that if **B**1 holds, then

$$
W_{f,N}(\mu_t^N, \nu_t^N) \le e^{-\bar{c}t} W_{f,N}(\mu_0^{\otimes N}, \nu_0^{\otimes N}), \quad W_{l^{1\circ\pi}}(\mu_t^N, \nu_t^N) \le 2M_1 e^{-\bar{c}t} W_{l^{1\circ\pi}}(\mu_0^{\otimes N}, \nu_0^{\otimes N}),
$$

where $f$, $\tilde{c}$ and $M_1$ are defined as in Theorem 8.

**5 System of N sticky SDEs**

Consider a system of $N$ one-dimensional SDEs with sticky boundaries at 0 given by

$$
\mathrm{d}r_t^i = \left( \tilde{b}(r_t^i) + \frac{1}{N} \sum_{j=1}^N g(r_t^j) \right) \mathrm{d}t + 2 \mathbb{1}_{(0,\infty)}(r_t^i) \mathrm{d}W_t^i, \quad i = 1, \dots, N. \quad (34)
$$

The results on existence, uniqueness and the comparison theorem for solutions of sticky
nonlinear SDEs mostly carry directly over to a solution of (34) and are applied to prove
propagation of chaos in Theorem 8.

Let $\mu$ be a probability distribution on $\mathbb{R}_+$. For $N \in \mathbb{N}$, $(\{r_t^i, W_t^i\}_{i=1}^N)_{t \ge 0}$ is a weak solution on the filtered probability space $(\Omega, \mathcal{A}, (\mathcal{F}_t), P)$ of (34) with initial distribution $\mu^{\otimes N}$ if the following hold: $\mu^{\otimes N} = P \circ (\{r_0\}_{i=1}^N)^{-1}$, $(\{W_t\}_{t \ge 0})_{t \ge 0}$ is a $N$-dimensional $(\mathcal{F}_t)$ Brownian motion w.r.t. $P$, the process $(r_t^i)_{t \ge 0}$ is non-negative, continuous and satisfies almost surely for any $i \in \{1, \dots, N\}$ and $t \in \mathbb{R}_+$,

$$
r_t^i - r_0^i = \int_0^t (\tilde{b}(r_s^i) + \frac{1}{N} \sum_{j=1}^N g(r_s^j)) \mathrm{d}s + \int_0^t 2 \mathbb{1}_{(0,\infty)}(r_s^i) \mathrm{d}W_s^i.
$$

To show existence and uniqueness in law of a weak solution $(\{r_t^i, W_t^i\}_{i=1}^N)_{t \ge 0}$, we suppose $\mathbf{H}1$ and $\mathbf{H}2$ for $\tilde{b}$ and $g$.
It follows that there exists a constant $C < \infty$ such that for all $\{r^i\}_{i=1}^N \in \mathbb{R}_+^N$, it holds $\sum_{i=1}^N |\tilde{b}(r^i)| + |g(r^i)| \le C(1 + \sum_{i=1}^N |r^i|)$, and a possible solution $(\{r_t^i\}_{i=1}^N)_{t \ge 0}$ is non-explosive. If the initial distribution is supported on $\mathbb{R}_+^N$, then in the same line as for the nonlinear SDE in Section 3.1, the solution $(\{r_t^i\}_{i=1}^N)_{t \ge 0}$ satisfies $r_t^i > 0$ almost surely for any $i = 1, \dots, N$ and $t \ge 0$ by $\mathbf{H}1$ and $\mathbf{H}2$.

Existence and uniqueness in law of (34) is a direct consequence of a stronger result that we now introduce. To study existence and uniqueness and to compare two solutions of (34) with different drifts, we establish existence of a synchronous coupling of two copies of (34),

$$ 
\begin{align} 
dr_t^i &= (\tilde{b}(r_t^i) + \frac{1}{N} \sum_{j=1}^N g(r_t^j)) dt + 2 \mathbb{1}_{(0,\infty)}(r_t^i) dW_t^i, \notag \\ 
ds_t^i &= (\hat{b}(s_t^i) + \frac{1}{N} \sum_{j=1}^N h(s_t^j)) dt + 2 \mathbb{1}_{(0,\infty)}(s_t^i) dW_t^i, \tag{35} \\ 
\text{Law}(r_0^i, s_0^i) &= \eta, \notag 
\end{align} 
$$

where $(\{W_t^i\}_{i=1}^N)_{t \ge 0}$ are $N$ i.i.d.1-dimensional Brownian motions and where $\eta \in \Gamma(\mu, \nu)$ for $\mu, \nu \in \mathcal{P}(\mathbb{R}_+)$. Let $\mathbb{W}^N = \mathcal{C}(\mathbb{R}_+, \mathbb{R}^N)$ be the space of continuous functions from $\mathbb{R}_+$ to $\mathbb{R}^N$ endowed with the topology of uniform convergence on compact sets, and let $\mathcal{B}(\mathbb{W}^N)$ denote its Borel $\sigma$-Algebra.

**Theorem 10.** Assume that $(\tilde{b}, g)$ and $(\hat{b}, h)$ satisfy $\mathbf{H}1$ and $\mathbf{H}2$. Let $\eta \in \Gamma(\mu, \nu)$ where $\mu$ and $\nu$ are the probability measure on $\mathbb{R}_+$ satisfying $\mathbf{H}3$. Then there exists a weak solution $(\{r_t^i, s_t^i\}_{i=1}^N)_{t \ge 0}$ of the sticky stochastic differential equation (35) with initial distribution $\eta^{\otimes N}$ defined on a probability space $(\Omega, \mathcal{A}, \mathcal{P})$ with values in $\mathbb{W}^N \times \mathbb{W}^N$. If additionally,

$$ 
\begin{gathered} 
\tilde{b}(r) \le \hat{b}(r) \quad \text{and} \quad g(r) \le h(r), \qquad \text{for any } r \in \mathbb{R}_+, \\ 
P[r_0^i \le s_0^i \text{ for all } i = 1, \dots, N] = 1, 
\end{gathered} 
$$

then $P[r_t^i \le s_t^i \text{ for all } t \ge 0 \text{ and } i = 1, \dots, N] = 1$.

*Proof.* The proof is postponed to Section 6.4. $\square$

**Remark 11.** We note that by the comparison result we can deduce uniqueness in law for the solution of (34).

# 6 Proofs

## 6.1 Definition of the metrics

In Theorem 1, Theorem 7 and Theorem 8 we consider Wasserstein distances based on carefully designed concave function $f: \mathbb{R}_+ \to \mathbb{R}_+$ that we now define. In addition we
derive useful properties of this function that will be used in our proofs of Theorem 1, Theorem 8 and Theorem 7. Let $a \in \mathbb{R}_+$ and $\tilde{b}: \mathbb{R}_+ \to \mathbb{R}$ be such that **H**4 is satisfied with $\tilde{R}_0$ and $\tilde{R}_1$ defined in (29). We define

$$ \varphi(r) = \exp\left(-\int_{0}^{r} \{\tilde{b}(s)_{+}/2\}ds\right), \quad \Phi(r) = \int_{0}^{r} \varphi(s)ds, \quad \text{and} $$

$$ g(r) = 1 - \frac{c}{2} \int_{0}^{r \wedge \tilde{R}_{1}} \{\Phi(s)/\varphi(s)\}ds - \frac{a}{2} \int_{0}^{r \wedge \tilde{R}_{1}} \{1/\varphi(s)\}ds $$

where

$$ c = \left( 2 \int_{0}^{\tilde{R}_{1}} \{\Phi(s)/\varphi(s)\}ds \right)^{-1}, \quad (36) $$

and $\tilde{R}_1$ is given in (30). It holds $\varphi(r) = \varphi(\tilde{R}_0)$ for $r \ge \tilde{R}_0$ with $\tilde{R}_0$ given in (29), $g(r) = g(\tilde{R}_1) \in [1/2, 3/4]$ for $r \ge \tilde{R}_1$ and $g(r) \in [1/2, 1]$ for all $r \in \mathbb{R}_+$ by (36) and **H**4. We define the increasing function $f: [0, \infty) \to [0, \infty)$ by

$$ f(t) = \int_{0}^{t} \varphi(r)g(r)dr. \quad (37) $$

Note that $f$ is concave, since $\varphi$ and $g$ are decreasing. Since for all $r \in \mathbb{R}_+$

$$ \varphi(\tilde{R}_0)r/2 \le \Phi(r)/2 \le f(r) \le \Phi(r) \le r, \quad (38) $$

$(x, y) \mapsto f(|x - y|)$ defines a distance on $\mathbb{R}^d$ equivalent to the Euclidean distance on $\mathbb{R}^d$. Moreover, $f$ satisfies

$$ 2f''(0) = -\tilde{b}(0)_{+} - a = -a, \quad (39) $$

and

$$ 2f''(r) \le 2f''(0) - f'(r)\tilde{b}(r) - cf(r), \quad \text{for all } r \in \mathbb{R}_{+} \setminus \{\tilde{R}_1\}. \quad (40) $$

Indeed by construction of $f$, $f''(r) = -\tilde{b}(r)_{+}f'(r)/2 - c\Phi(r)/2 - a/2$ for $0 \le r < \tilde{R}_1$ and so (40) holds for $0 \le r < \tilde{R}_1$ by (38). To show (40) for $r > \tilde{R}_1$ note that $f''(r) = 0$ and $f'(r) \ge \varphi(\tilde{R}_0)/2$ hold for $r > \tilde{R}_1$. Hence, by the definition (30) of $\tilde{R}_1$, for $r > \tilde{R}_1$,

$$ f''(r) + f'(r)\tilde{b}(r)/2 \le \varphi(\tilde{R}_0)\tilde{b}(r)/4 \le -(\tilde{R}_1(\tilde{R}_1 - \tilde{R}_0))^{-1}\varphi(\tilde{R}_0)r. \quad (41) $$

Since $\varphi(r) = \varphi(\tilde{R}_0)$ for $r \ge \tilde{R}_0$, it holds $\Phi(r) = \Phi(\tilde{R}_0) + (r - \tilde{R}_0)\varphi(\tilde{R}_0)$ for $r \ge \tilde{R}_0$. Further, it holds $\Phi(R_0) \ge \tilde{R}_0\varphi(\tilde{R}_0)$ since $\varphi$ is decreasing for $r \le \tilde{R}_0$. Hence,

$$ \frac{r}{\tilde{R}_1} = \frac{(r - \tilde{R}_1)(\Phi(\tilde{R}_0) + (\tilde{R}_1 - \tilde{R}_0)\varphi(\tilde{R}_0))}{\tilde{R}_1\Phi(\tilde{R}_1)} + 1 \ge \frac{(r - \tilde{R}_1)\tilde{R}_1\varphi(\tilde{R}_0)}{\tilde{R}_1\Phi(\tilde{R}_1)} + 1 = \frac{\Phi(r)}{\Phi(\tilde{R}_1)}. \quad (42) $$
Furthermore, we have

$$ 
\begin{aligned}
\int_{\tilde{R}_0}^{\tilde{R}_1} \{\Phi(s)/\varphi(s)\}ds &= \int_{\tilde{R}_0}^{\tilde{R}_1} \frac{\Phi(\tilde{R}_0) + (s - \tilde{R}_0)\varphi(\tilde{R}_0)}{\varphi(\tilde{R}_0)}ds \\
&= (\tilde{R}_1 - \tilde{R}_0)\frac{\Phi(\tilde{R}_0)}{\varphi(\tilde{R}_0)} + \frac{1}{2}(\tilde{R}_1 - \tilde{R}_0)^2 \ge \frac{1}{2}(\tilde{R}_1 - \tilde{R}_0)\frac{\Phi(\tilde{R}_1)}{\varphi(\tilde{R}_0)}.
\end{aligned} 
\quad (43) $$

We insert (42) and (43) in (41) and use (36) to obtain

$$ f''(r) + f'(r)\tilde{b}(r)/2 \le -\Phi(r)\Phi(\tilde{R}_1)^{-1}(\tilde{R}_1 - \tilde{R}_0)^{-1}\varphi(\tilde{R}_0) \quad (44) $$

$$ \le -\frac{\Phi(r)}{2\int_{\tilde{R}_0}^{\tilde{R}_1} \{\Phi(s)/\varphi(s)\}ds} \le -\frac{cf(r)}{2} - \frac{c\Phi(r)}{2}. \quad (45) $$

By H4 and (36), we get

$$ -\frac{c\Phi(r)}{2} \le -\frac{\Phi(\tilde{R}_1)}{4\int_0^{\tilde{R}_1} \{\Phi(s)/\varphi(s)\}ds} \le -\frac{1}{4\int_0^{\tilde{R}_1} \{1/\varphi(s)\}ds} \le -\frac{a}{2} = f''(0). $$

Combining this estimate with (44) gives (40) for $r > \tilde{R}_1$.

## 6.2 Proof of Section 2

### 6.2.1 Proof of Theorem 2

Note that the nonlinear SDE (21) has Lipschitz continuous coefficients. The existence and the uniqueness of the coupling $(\bar{X}_t^\delta, \bar{Y}_t^\delta)_{t \ge 0}$ follows from [45, Theorem 2.2]. By Levy's characterization, $(\bar{X}_t^\delta, \bar{Y}_t^\delta)_{t \ge 0}$ is indeed a coupling of two copies of solutions of (1). Further, we remark that $W_t^\delta = \int_0^t (\bar{e}_s^\delta)^T \cdot dB_s^1$ is a one-dimensional Brownian motion. In the next step, we analyse $|\bar{X}_t^\delta - \bar{Y}_t^\delta|$.

**Lemma 12.** Suppose that the conditions B1 and B3 are satisfied. Then, it holds for any $\epsilon < \epsilon_0$, where $\epsilon_0$ is given by (20), setting $\bar{r}_t^\delta = |\bar{X}_t^\delta - \bar{Y}_t^\delta|$

$$ d\bar{r}_t^\delta = (-L\bar{r}_t^\delta + \langle \bar{e}_t^\delta, \mathbb{E}_{x\sim\bar{\mu}_t^\delta, y\sim\bar{\nu}_t^\delta}[\gamma(\bar{X}_t^\delta - x) - \gamma(\bar{Y}_t^\delta - y)]\rangle)dt + 2rc^\delta(\bar{r}_t^\delta)dW_t^\delta \quad (46) $$

$$ \le (\bar{b}(\bar{r}_t^\delta) + 2\|\gamma\|_\infty \int_{\mathbb{R}^d} \int_{\mathbb{R}^d} \mathrm{rc}^\epsilon(|x-y|)\bar{\mu}_t^\delta(dx)\bar{\nu}_t^\delta(dy))dt + 2\mathrm{rc}^\delta(\bar{r}_t^\delta)dW_t^\delta, \quad (47) $$

almost surely for all $t \ge 0$, where $\bar{\mu}_t^\delta$ and $\bar{\nu}_t^\delta$ are the laws of $\bar{X}_t^\delta$ and $\bar{Y}_t^\delta$, respectively.

*Proof.* Using (21), B1 and B3, the stochastic differential equation of the process $((\bar{r}_t^\delta)^2)_{t \ge 0}$ is given by

$$ d((\bar{r}_t^\delta)^2) = 2\langle Z_t^\delta, -LZ_t^\delta + \int_{\mathbb{R}^d} \int_{\mathbb{R}^d} \gamma(\bar{X}_t^\delta - x) - \gamma(\bar{Y}_t^\delta - y)\bar{\mu}_t^\delta(dx)\bar{\nu}_t^\delta(dy) \rangle dt \\ + 4\mathrm{rc}^\delta(\bar{r}_t^\delta)^2 dt + 4\mathrm{rc}^\delta(\bar{r}_t^\delta)\langle Z_t^\delta, e_t^\delta \rangle dW_t^\delta. $$
For $\epsilon > 0$ we define as in [23, Lemma 8] a $C^2$ approximation of the square root by

$$S_{\epsilon}(r) = \begin{cases} (-1/8)\epsilon^{-3/2}r^2 + (3/4)\epsilon^{-1/2}r + (3/8)\epsilon^{1/2} & \text{for } r < \epsilon \\ \sqrt{r} & \text{otherwise.} \end{cases}$$

Then, by Ito's formula,

$$\begin{aligned} dS_{\epsilon}((\bar{r}_t^{\delta})^2) &= S'_{\epsilon}((\bar{r}_t^{\delta})^2)d(\bar{r}_t^{\delta})^2 + \frac{1}{2}S''_{\epsilon}((\bar{r}_t^{\delta})^2)d[(\bar{r}^{\delta})^2]_t \\ &= 2S'_{\epsilon}((\bar{r}_t^{\delta})^2)\langle Z_t^{\delta}, -LZ_t^{\delta} + \int_{\mathbb{R}^d} \int_{\mathbb{R}^d} \gamma(\bar{X}_t^{\delta} - x) - \gamma(\bar{Y}_t^{\delta} - y)\bar{\mu}_t^{\delta}(dx)\bar{\nu}_t^{\delta}(dy) \rangle dt \\ &\quad + S'_{\epsilon}((\bar{r}_t^{\delta})^2)4\mathrm{rc}^{\delta}(\bar{r}_t^{\delta})^2dt + S'_{\epsilon}((\bar{r}_t^{\delta})^2)4\mathrm{rc}^{\delta}(\bar{r}_t^{\delta})\langle Z_t^{\delta}, e_t^{\delta} \rangle dW_t^{\delta} + 8S''_{\epsilon}((\bar{r}_t^{\delta})^2)(\mathrm{rc}^{\delta}(\bar{r}_t^{\delta}))^2(\bar{r}_t^{\delta})^2dt. \end{aligned}$$

We take the limit $\epsilon \to 0$. Then for $r > 0$, $\lim_{\epsilon \to 0} S'_{\epsilon}(r) = (1/2)r^{-1/2}$ and $\lim_{\epsilon \to 0} S''_{\epsilon}(r) = -(1/4)r^{-3/2}$. Since $\sup_{0 \le r \le \epsilon} |S'_{\epsilon}(r)| \lesssim \epsilon^{-1/2}$, $\sup_{0 \le r \le \epsilon} |S''_{\epsilon}(r)| \lesssim \epsilon^{-3/2}$ and $\mathrm{rc}^{\delta}$ is Lipschitz continuous with $\mathrm{rc}^{\delta}(0) = 0$, we apply Lebesgue's dominated convergence theorem to show convergence for the integrals with respect to time $t$. More precisely, we note that the integrand $(4S'_{\epsilon}((\bar{r}_t^{\delta})^2) + 8S''_{\epsilon}((\bar{r}_t^{\delta})^2))\mathrm{rc}^{\delta}(\bar{r}_t^{\delta}))^2(\bar{r}_t^{\delta})^2$ is dominated by $3\epsilon^{1/2}\|\mathrm{rc}^{\delta}\|_{\mathrm{Lip}}$. For any $\epsilon < \epsilon_0$ for fixed $\epsilon_0 > 0$, the integrand $2S'_{\epsilon}((\bar{r}_t^{\delta})^2)\langle Z_t^{\delta}, -LZ_t^{\delta} + \int \int \gamma(\bar{X}_t^{\delta} - x) - \gamma(\bar{Y}_t^{\delta} - y)d\bar{\mu}_t^{\delta}(x)d\bar{\nu}_t^{\delta}(y) \rangle$ is dominated by $(3/2)(L \max(\epsilon_0^{1/2}, \bar{r}_t^{\delta}) + 2\|\gamma\|_{\infty})$.

For the stochastic integral it holds $|S'_{\epsilon}((\bar{r}_t^{\delta})^2)4\mathrm{rc}^{\delta}(\bar{r}_t^{\delta})\bar{r}_t^{\delta}| \le 3$. Hence, the stochastic integral converges along a subsequence almost surely, to $\int_0^t 2\mathrm{rc}^{\delta}(\bar{r}_s^{\delta})dW_s^{\delta}$, see [51, Chapter 4, Theorem 2.12]. Hence, we obtain (46). By B1 and (20), we obtain for $\epsilon < \epsilon_0$

$$\begin{aligned} \langle \bar{e}_t^{\delta}, \mathbb{E}[\gamma(\bar{X}_t^{\delta} - \hat{X}_t^{\delta}) - \gamma(\bar{Y}_t^{\delta} - \hat{Y}_t^{\delta})] \rangle \\ \le \langle \bar{e}_t^{\delta}, \mathbb{E}_{x \sim \bar{\mu}_t^{\delta}, y \sim \bar{\nu}_t^{\delta}}[\gamma(\bar{X}_t^{\delta} - x) - \gamma(\bar{Y}_t^{\delta} - x) + \gamma(\bar{Y}_t^{\delta} - x) - \gamma(\bar{Y}_t^{\delta} - y)] \rangle \\ \le \kappa(\bar{r}_t^{\delta})\bar{r}_t^{\delta} + \mathbb{E}_{x \sim \bar{\mu}_t^{\delta}, y \sim \bar{\nu}_t^{\delta}}[2\|\gamma\|_{\infty}\mathrm{rc}^{\epsilon}(|x - y|)], \end{aligned}$$

and hence (47) holds. □

We define a one-dimensional process $(r_t^{\delta, \epsilon})_{t \ge 0}$ by

$$dr_t^{\delta, \epsilon} = (\bar{b}(r_t^{\delta, \epsilon}) + 2\|\gamma\|_{\infty} \int \mathrm{rc}^{\epsilon}(u)dP_t^{\delta, \epsilon}(u))dt + 2 \cdot \mathrm{rc}^{\delta}(r_t^{\delta, \epsilon})dW_t^{\delta} \quad (48)$$

with initial condition $r_0^{\delta, \epsilon} = \bar{r}_0^{\delta}$, $P_t^{\delta, \epsilon} = \mathrm{Law}(r_t^{\delta, \epsilon})$ and $W_t^{\delta} = \int_0^t (\bar{e}_s^{\delta})^T dB_s^1$. This process will allow us to control the distance of $\bar{X}_t^{\delta}$ and $\bar{Y}_t^{\delta}$.

By [45, Theorem 2.2], under B1 and B3, $(U_t^{\delta, \epsilon})_{t \ge 0} = (\bar{X}_t^{\delta}, \bar{Y}_t^{\delta}, r_t^{\delta, \epsilon})_{t \ge 0}$ exists and is unique, where $(\bar{X}_t^{\delta}, \bar{Y}_t^{\delta})_{t \ge 0}$ solves uniquely (21), $(\bar{r}_t^{\delta})_{t \ge 0}$ and $(r_t^{\delta, \epsilon})_{t \ge 0}$ solve uniquely (46) and (48), respectively, with $W_t^{\delta} = \int_0^t (\bar{e}_s^{\delta})^T dB_s^1$.

**Lemma 13.** Assume B1 and B3. Then, $|\bar{X}_t^{\delta} - \bar{Y}_t^{\delta}| = \bar{r}_t^{\delta} \le r_t^{\delta, \epsilon}$, almost surely for all $t$ and $\epsilon < \epsilon_0$.
Proof. Note that $(\bar{r}_t^\delta)_{t \ge 0}$ and $(r_t^{\delta, \epsilon})_{t \ge 0}$ have the same initial distribution and are driven by the same noise. Since the drift of $(\bar{r}_t^\delta)_{t \ge 0}$ is smaller than the drift of $(r_t^{\delta, \epsilon})_{t \ge 0}$ for $\epsilon < \epsilon_0$, the result follows by Lemma 18.

<figure><img src="image_3.png" /></figure>

Proof of Theorem 2. We consider the nonlinear process $(U_t^{\delta, \epsilon})_{t \ge 0} = (\bar{X}_t^\delta, \bar{Y}_t^\delta, r_t^{\delta, \epsilon})_{t \ge 0}$ on $\mathbb{R}^{2d+1}$ for each $\epsilon, \delta > 0$. We denote by $\mathbb{P}^{\delta, \epsilon}$ the law of $U^{\delta, \epsilon}$ on the space $C(\mathbb{R}_+, \mathbb{R}^{2d+1})$. We define by $\mathbf{X}, \mathbf{Y}: C(\mathbb{R}_+, \mathbb{R}^{2d+1}) \to C(\mathbb{R}_+, \mathbb{R}^d)$ and $\mathbf{r}: C(\mathbb{R}_+, \mathbb{R}^{2d+1}) \to C(\mathbb{R}_+, \mathbb{R})$ the canonical projections onto the first $d$ components, onto the second $d$ components and onto the last component, respectively. By **B1** and **B3** following the same line as the proof of Lemma 19, see (65), it holds for each $T > 0$

$$E[|U_{t_2}^{\delta, \epsilon} - U_{t_1}^{\delta, \epsilon}|^4] \le C|t_2 - t_1|^2 \quad \text{for } t_1, t_2 \in [0, T], \qquad (49)$$

for some constant $C$ depending on $T$, $L$, $\|\gamma\|_{\text{Lip}}$, $\|\gamma\|_{\infty}$ and on the fourth moment of $\mu_0$ and $\nu_0$. As in Lemma 19 the law $\mathbb{P}_T^{\delta, \epsilon}$ of $(U_t^{\delta, \epsilon})_{0 \le t \le T}$ on $C([0, T], \mathbb{R}^{2d+1})$ is tight for each $T > 0$ by [38, Corollary 14.9] and for each $\epsilon > 0$ there exists a subsequence $\delta_n \to 0$ such that $(\mathbb{P}_T^{\delta, \epsilon})_{n \in \mathbb{N}}$ on $C([0, T], \mathbb{R}^{2d+1})$ converge to a measure $\mathbb{P}_T^\epsilon$ on $C([0, T], \mathbb{R}^{2d+1})$. By a diagonalization argument and since $\{\mathbb{P}_T^\epsilon : T \ge 0\}$ is a consistent family, cf. [38, Theorem 5.16], there exists a probability measure $\mathbb{P}^\epsilon$ on $C(\mathbb{R}_+, \mathbb{R}^{2d+1})$ such that for all $\epsilon$ there exists a subsequence $\delta_n$ such that $(\mathbb{P}^{\delta, \epsilon})_{n \in \mathbb{N}}$ converges along this subsequence to $\mathbb{P}^\epsilon$, i.e., the family $(\mathbb{P}^{\delta, \epsilon})_{\delta > 0}$ is tight. As in the proof of Lemma 20 we repeat this argument for the family of measures $(\mathbb{P}^\epsilon)_{\epsilon > 0}$. Hence, there exists a subsequence $\epsilon_n \to 0$ such that $(\mathbb{P}^{\epsilon_n})_{n \in \mathbb{N}}$ converges to a measure $\mathbb{P}$. Let $(\bar{X}_t, \bar{Y}_t, r_t)_{t \ge 0}$ be some process on $\mathbb{R}^{2d+1}$ with distribution $\mathbb{P}$ on $(\bar{\Omega}, \bar{\mathcal{F}}, \bar{P})$.

Since $(\bar{X}_t^\delta)_{t \ge 0}$ and $(\bar{Y}_t^\delta)_{t \ge 0}$ are solutions of (1) which are unique in law, we have that for any $\epsilon, \delta > 0$, $\mathbb{P}^{\delta, \epsilon} \circ (\mathbf{X})^{-1} = \mathbb{P} \circ (\mathbf{X})^{-1}$ and $\mathbb{P}^{\delta, \epsilon} \circ (\mathbf{Y})^{-1} = \mathbb{P} \circ (\mathbf{Y})^{-1}$. And therefore $(\bar{X}_t)_{t \ge 0}$ and $(\bar{Y}_t)_{t \ge 0}$ are solutions of (1) as well with the same initial condition. Hence $\mathbb{P} \circ (\mathbf{X}, \mathbf{Y})^{-1}$ is a coupling of two copies of (1).

In addition, the statement of Lemma 13 carries over to the limiting process $(r_t)_{t \ge 0}$, i.e., $|\bar{X}_t - \bar{Y}_t| \le r_t$ for all $t \ge 0$, since by the weak convergence and the Portmanteau theorem, $\bar{P}(|\bar{X}_t - \bar{Y}_t| \le r_t) \ge \limsup_{\epsilon \to 0} \limsup_{\delta \to 0} P(|\bar{X}_t^\delta - \bar{Y}_t^\delta| \le r_t^{\delta, \epsilon}) = 1$.

Finally, similarly to the proof of Lemma 19 and Lemma 20 there exist an extended probability space and a one-dimensional Brownian motion $(W_t)_{t \ge 0}$ such that $(r_t, W_t)_{t \ge 0}$ is a solution to

$$dr_t = (\bar{b}(r_t) + 2\|\gamma\|_{\infty}\mathbb{P}(r_t > 0))dt + 2\mathbb{1}_{(0,\infty)}(r_t)dW_t.$$

<figure><img src="image_4.png" /></figure>
### 6.2.2 Proof of Theorem 1

*Proof of Theorem 1.* We consider the process $(\bar{X}_t, \bar{Y}_t, r_t)_{t \ge 0}$ defined in Theorem 2 and satisfying $|\bar{X}_t - \bar{Y}_t| \le r_t$ for any $t \ge 0$, and $(r_t)_{t \ge 0}$ is a weak solution of (22). Set $a = 2\|\gamma\|_{\infty}$ and $\tilde{b}(r) = \bar{b}(r)$. With this notation, **B1** and **B2** imply **H4** and $\tilde{R}_0 = R_0$ and $\tilde{R}_1 = R_1$ by (14), (15), (29) and (30). By Ito-Tanaka formula, cf. [51, Chapter 6, Theorem 1.1], using that $f'$ is absolutely continuous, we have,

$$d f(r_t) \le f'(r_t)(\bar{b}(r_t) + 2\|\gamma\|_{\infty}\mathbb{P}(r_t > 0))dt + 2f''(r_t)1_{(0,\infty)}(r_t)dt + f'(r_t)21_{(0,\infty)}(r_t)dW_t.$$

Taking expectation we obtain by (39) and (40)

$$\frac{d}{dt} \mathbb{E}[f(r_t)] \le \mathbb{E}[f'(r_t)\tilde{b}(r_t)_+ + 2(f''(r_t) - f''(0))] + \mathbb{E}[(a + 2f''(0))1_{r_t>0}] \le -\tilde{c}\mathbb{E}[f(r_t)],$$

where $\tilde{c}$ is given by (17). Therefore by Gronwall's lemma,

$$\mathbb{E}[f(|\bar{X}_t - \bar{Y}_t|)] \le \mathbb{E}[f(r_t)] \le e^{-\tilde{c}t} \mathbb{E}[f(r_0)] = e^{-\tilde{c}t} \mathbb{E}[f(|\bar{X}_0 - \bar{Y}_0|)].$$

Hence, it holds

$$W_f(\bar{\mu}_t, \bar{\nu}_t) \le \mathbb{E}[f(|\bar{X}_t - \bar{Y}_t|)] \le e^{-\tilde{c}t} \int f(|x - y|)\xi(dx dy)$$

for an arbitrary coupling $\xi \in \Gamma(\mu_0, \nu_0)$. Taking the infimum over all couplings $\xi \in \Gamma(\mu_0, \nu_0)$, we obtain the first inequality of (16). By (38), we get the second inequality of (16). $\square$

### 6.2.3 Proof of Theorem 8

To show Theorem 8, we first give a uniform in time bound for the second moment of the process $(\bar{X}_t)_{t \ge 0}$ solving (1).

**Lemma 14.** Let $(\bar{X}_t)_{t \ge 0}$ be a solution of (1) with $\mathbb{E}[|\bar{X}_0|^2] < \infty$. Assume **B1**. Then there exists $C \in (0, \infty)$ depending on $d$, $W$ and the second moment of $\bar{X}_0$ such that

$$C = \sup_{t \ge 0} \mathbb{E}[|\bar{X}_t|^2] < \infty. \qquad (50)$$

*Proof of Lemma 14.* By Ito's formula, it holds

$$\frac{1}{2} d|\bar{X}_t|^2 = \langle \bar{X}_t, b * \bar{\mu}_t(\bar{X}_t) \rangle dt + \bar{X}_t^T dB_t + \frac{1}{2} d dt.$$
Taking expectation and using symmetry, we get

$$
\begin{align*}
\frac{d}{dt} \mathbb{E}[|\bar{X}_t|^2] &= \mathbb{E}[\langle \bar{X}_t - \tilde{X}_t, b(\bar{X}_t - \tilde{X}_t) \rangle + d \\
&= -\mathbb{E}[\langle \bar{X}_t - \tilde{X}_t, L(\bar{X}_t - \tilde{X}_t) - \gamma(\bar{X}_t - \tilde{X}_t) \rangle 1_{|\bar{X}_t - \tilde{X}_t| > R_0}] \\
&\quad - \mathbb{E}[\langle \bar{X}_t - \tilde{X}_t, L(\bar{X}_t - \tilde{X}_t) - \gamma(\bar{X}_t - \tilde{X}_t) \rangle 1_{|\bar{X}_t - \tilde{X}_t| \le R_0}] + d \\
&\le \mathbb{E}[|\bar{X}_t|^2 (-2L + \kappa(|\bar{X}_t - \tilde{X}_t|) 1_{|\bar{X}_t - \tilde{X}_t| > R_0})] + \|\gamma\|_{\infty} R_0 + d.
\end{align*}
$$

Hence by definition (14) of $R_0$ and by Gronwall's lemma we obtain the result (50). $\square$

Let $N \in \mathbb{N}$. We construct a sticky coupling of $N$ i.i.d. realizations of solutions
($\{\bar{X}_t^i\}_{i=1}^N)_{t \ge 0}$ to (1) and of the solution ($\{Y_t^i\}_{i=1}^N)_{t \ge 0}$ to the mean field particle system (3).
Then, we consider a weak limit for $\delta \to 0$ of Markovian couplings which are constructed
similar as in Section 2. Let $\text{rc}^\delta$, $\text{sc}^\delta$ satisfy (19) and (20). The coupling ($\{\bar{X}_t^{i,\delta}, Y_t^{i,\delta}\}_{i=1}^N)_{t \ge 0}$
is defined as process in $\mathbb{R}^{2Nd}$ satisfying a system of SDEs given by

$$
\begin{align}
d\bar{X}_t^{i,\delta} &= b * \bar{\mu}_t^\delta(\bar{X}_t^{i,\delta})dt + \text{rc}^\delta(\tilde{r}_t^{i,\delta})dB_t^{i,1} + \text{sc}^\delta(\tilde{r}_t^{i,\delta})dB_t^{i,2} \nonumber \\
dY_t^{i,\delta} &= \frac{1}{N} \sum_{j=1}^N (b(Y_t^{i,\delta} - Y_t^{j,\delta})dt + \text{rc}^\delta(\tilde{r}_t^{i,\delta})(\text{Id} - 2\tilde{e}_t^{i,\delta}(\tilde{e}_t^{i,\delta})^T)dB_t^{i,1} + \text{sc}^\delta(\tilde{r}_t^{i,\delta})dB_t^{i,2}, \tag{51}
\end{align}
$$

where Law($\{\bar{X}_0^{i,\delta}, Y_0^{i,0}\}_{i=1}^N$) = $\bar{\mu}_0^{\otimes N} \otimes \nu_0^{\otimes N}$, and where $(\{B_t^{i,1}\}_{i=1}^N)_{t \ge 0}$, $(\{B_t^{i,2}\}_{i=1}^N)_{t \ge 0}$ are i.i.d. $d$-dimensional standard Brownian motions. We set $\tilde{X}_t^{i,\delta} = \bar{X}_t^{i,\delta} - \frac{1}{N} \sum_{j=1}^N \bar{X}_t^{j,\delta}$, $\tilde{Y}_t^{i,\delta} = Y_t^{i,\delta} - \frac{1}{N} \sum_{j=1}^N Y_t^{j,\delta}$, $\tilde{Z}_t^{i,\delta} = \tilde{X}_t^{i,\delta} - \tilde{Y}_t^{i,\delta}$, $\tilde{r}_t^{i,\delta} = |\tilde{Z}_t^{i,\delta}|$ and $\tilde{e}_t^{i,\delta} = \tilde{Z}_t^{i,\delta}/\tilde{r}_t^{i,\delta}$ for $\tilde{r}_t^{i,\delta} \neq 0$. The value $\tilde{e}_t^{i,\delta}$ for $\tilde{r}_t^{i,\delta} = 0$ is irrelevant as $\text{rc}^{i,\delta}(0) = 0$. By Levy's characterization $(\{\bar{X}_t^{i,\delta}, Y_t^{i,\delta}\}_{i=1}^N)_{t \ge 0}$ is indeed a coupling of (1) and (3). Existence and uniqueness of the coupling given in (51) hold by [45, Theorem 2.2]. In the next step we analyse $\tilde{r}_t^{i,\delta}$.

**Lemma 15.** Assume **B** 1 holds. Then, for $\epsilon < \epsilon_0$, where $\epsilon_0$ is given in (20), and for any $i \in \{1, \dots, N\}$, it holds almost surely,

$$
\begin{align}
d\tilde{r}_t^{i,\delta} &= -L\tilde{r}_t^{i,\delta}dt + \langle \tilde{e}_t^{i,\delta}, \frac{1}{N}\sum_{j=1}^N \gamma(\tilde{X}_t^{i,\delta} - \tilde{X}_t^{j,\delta}) - \gamma(\tilde{Y}_t^{i,\delta} - \tilde{Y}_t^{j,\delta}) \rangle dt \nonumber \\
&\quad + 2\sqrt{1+\frac{1}{N}}\text{rc}^\delta(\tilde{r}_t^{i,\delta})dW_t^{i,\delta} + \langle \tilde{e}_t^{i,\delta}, \Theta_t^{i,\delta} + \frac{1}{N}\sum_{k=1}^N \Theta_t^{k,\delta} \rangle dt \tag{52} \\
&\le \left(\bar{b}(\tilde{r}_t^{i,\delta}) + 2\|\gamma\|_\infty \frac{1}{N}\sum_{j=1}^N \text{rc}^\epsilon(\tilde{r}_t^{j,\delta})\right)dt + 2\sqrt{1+\frac{1}{N}}\text{rc}^\delta(\tilde{r}_t^{i,\delta})dW_t^{i,\delta} \nonumber \\
&\quad + \left(A_t^{i,\delta} + \frac{1}{N}\sum_{k=1}^N A_t^{k,\delta}\right)dt. \nonumber
\end{align}
$$
with $\Theta_t^{i,\delta} = b * \bar{\mu}_t^{\delta}(\bar{X}_t^{i,\delta}) - \frac{1}{N} \sum_{j=1}^{N} b(\bar{X}_t^{i,\delta} - \bar{X}_t^{j,\delta})$ and

$$A_t^{i,\delta} = |\Theta_t^{i,\delta}| = |b * \bar{\mu}_t^{\delta}(\bar{X}_t^{i,\delta}) - \frac{1}{N} \sum_{j=1}^{N} b(\bar{X}_t^{i,\delta} - \bar{X}_t^{j,\delta})| \quad (53)$$

and where $(\{W_t^{i,\delta}\}_{i=1}^N)_{t \ge 0}$ are $N$ one-dimensional Brownian motions given by

$$W_t^{i,\delta} = \sqrt{\frac{N}{N+1}} \left( \int_0^t (\bar{e}_s^{i,\delta})^T dB_s^{i,1} + \frac{1}{N} \sum_{j=1}^N \int_0^t (\bar{e}_s^{j,\delta})^T dB_s^{j,1} \right), \quad i = 1, \dots, N. \quad (54)$$

*Proof.* By (51) and since $\gamma$ is anti-symmetric, it holds by Ito's formula for any $i \in \{1, \dots, N\}$,

$$\begin{aligned} d(\tilde{r}_t^{i,\delta})^2 = & -2L(\tilde{r}_t^{i,\delta})^2 dt + 2\langle \tilde{Z}_t^{i,\delta}, \frac{1}{N} \sum_{j=1}^N \gamma(\tilde{X}_t^{i,\delta} - \tilde{X}_t^{j,\delta}) - \gamma(\tilde{Y}_t^{i,\delta} - \tilde{Y}_t^{j,\delta}) \rangle dt \\ & + 4\left(1 + \frac{1}{N}\right) \mathrm{rc}^{\delta}(\tilde{r}_t^{i,\delta})^2 dt + 4\sqrt{1 + \frac{1}{N}} \mathrm{rc}^{\delta}(\tilde{r}_t^{i,\delta}) \langle \tilde{Z}_t^{i,\delta}, \tilde{e}_t^{i,\delta} \rangle dW_t^{i,\delta} \\ & + 2\langle \tilde{Z}_t^{i,\delta}, b * \bar{\mu}_t^{\delta}(\bar{X}_t^{i,\delta}) - \frac{1}{N} \sum_{j=1}^N b(\bar{X}_t^{i,\delta} - \bar{X}_t^{j,\delta}) \rangle dt \\ & + 2\langle \tilde{Z}_t^{i,\delta}, -\frac{1}{N} \sum_{k=1}^N \left( b * \bar{\mu}_t^{\delta}(\bar{X}_t^{k,\delta}) - \frac{1}{N} \sum_{j=1}^N b(\bar{X}_t^{k,\delta} - \bar{X}_t^{j,\delta}) \right) \rangle dt. \end{aligned}$$

where $(\{W_t^i\}_{i=1}^N)_{t \ge 0}$ are $N$ i.i.d.one-dimensional Brownian motions given by (54). Note that the prefactor $(N/(N+1))^{1/2}$ ensures that the quadratic variation satisfies $[W^i]_t = t$ for $t \ge 0$, and hence $(\{W_t^i\}_{i=1}^N)_{t \ge 0}$ are Brownian motions. This definition of $(\{W_t^i\}_{i=1}^N)_{t \ge 0}$ leads to $(1+1/N)^{1/2}$ in the diffusion term of the SDE. Applying the $C^2$ approximation of the square root used in the proof of Lemma 12 and taking $\epsilon \to 0$ in the approximation yields the stochastic differential equations of $(\{\tilde{r}_t^{i,\delta}\}_{i=1}^N)_{t \ge 0}$. We obtain its upper bound for $\epsilon < \epsilon_0$ by B1 and (20) similarly to the proof of Lemma 12. $\square$

Next, we state a bound for (53).

**Lemma 16.** Under the same assumption as in Lemma 17, it holds for any $i = 1, \dots, N$

$$E[|A_t^{i,\delta}|^2] \le C_1 N^{-1} \quad \text{and} \quad E[A_t^{i,\delta}] \le C_2 N^{-1/2},$$

where $A_t^{i,\delta}$ is given in (53) and $C_1$ and $C_2$ are constants depending on $\|\gamma\|_{\infty}$, $L$ and $C$ given in Lemma 14.
Proof. By **B**3, it holds $\mathbb{E}[|\bar{X}_0^{i,\delta}|^2] < \infty$ for $i = 1, \dots, N$. Note that given $\bar{X}_t^{i,\delta}$, $\bar{X}_t^{j,\delta}$ are i.i.d. with law $\bar{\mu}_t^\delta$ for all $j \neq i$. Hence,

$$ \mathbb{E}[b(\bar{X}_t^{i,\delta} - \bar{X}_t^{j,\delta}) | \bar{X}_t^{i,\delta}] = b * \bar{\mu}_t^\delta(\bar{X}_t^{i,\delta}). $$

Since $\gamma$ is anti-symmetric, $b(0) = 0$, and we have

$$ \begin{aligned} & \mathbb{E}[|b * \bar{\mu}_t^\delta(\bar{X}_t^{i,\delta}) - \frac{1}{N-1} \sum_{j=1}^N b(\bar{X}_t^{i,\delta} - \bar{X}_t^{j,\delta})|^2 | \bar{X}_t^{i,\delta}] \\ &= \mathbb{E}[\left| \frac{1}{N-1} \sum_{j=1}^N \mathbb{E}[b(\bar{X}_t^{i,\delta} - \bar{X}_t^{j,\delta}) | \bar{X}_t^{i,\delta}] - \frac{1}{N-1} \sum_{j=1}^N b(\bar{X}_t^{i,\delta} - \bar{X}_t^{j,\delta})|^2 \right| \bar{X}_t^{i,\delta}] \\ &= \frac{1}{N-1} \mathrm{Var}_{\bar{\mu}_t^\delta}(b(\bar{X}_t^{i,\delta} - \cdot)). \end{aligned} $$

By (11), **B**1, **B**3 and Lemma 14, we obtain

$$ \begin{aligned} \mathrm{Var}_{\bar{\mu}_t^\delta}(b(\bar{X}_t^{i,\delta} - \cdot)) &= \int_{\mathbb{R}^d} \left| \left( -L(\bar{X}_t^{i,\delta} - x) + \int_{\mathbb{R}^d} L(\bar{X}_t^{i,\delta} - \tilde{x}) \bar{\mu}_t^\delta(d\tilde{x}) \right) \right. \\ & \qquad \left. + \left( \gamma(\bar{X}_t^{i,\delta} - x) - \int_{\mathbb{R}^d} \gamma(\tilde{X}_t^{i,\delta} - \tilde{x}) \bar{\mu}_t^\delta(d\tilde{x}) \right) \right|^2 \bar{\mu}_t^\delta(dx) \\ &= \int_{\mathbb{R}^d} |Lx + \left( \gamma(\bar{X}_t^{i,\delta} - x) - \int_{\mathbb{R}^d} \gamma(\tilde{X}_t^{i,\delta} - \tilde{x}) \bar{\mu}_t^\delta(d\tilde{x}) \right)|^2 \bar{\mu}_t^\delta(dx) \\ &\le 2L^2 \int_{\mathbb{R}^d} |x|^2 \bar{\mu}_t^\delta(dx) + 8 \|\gamma\|_{\infty}^2 \le 2L^2 C^2 + 8 \|\gamma\|_{\infty}^2. \end{aligned} $$

By the Cauchy-Schwarz inequality, we have

$$ \begin{aligned} \mathbb{E}[(A_t^{i,\delta})^2] &\le 2\mathbb{E}[|b * \bar{\mu}_t(\bar{X}_t^{i,\delta}) - \frac{1}{N-1} \sum_{j=1}^N b(\bar{X}_t^{i,\delta} - \bar{X}_t^{j,\delta})|^2] \\ &\quad + 2\left(\frac{1}{N-1} - \frac{1}{N}\right)^2 \mathbb{E}\left[\left|\sum_{j=1}^N b(\bar{X}_t^{i,\delta} - \bar{X}_t^{j,\delta})\right|^2\right] \\ &\le 2\frac{1}{N-1}\mathbb{E}[\mathrm{Var}_{\bar{\mu}_t^\delta}(b(\bar{X}_t^{i,\delta} - \cdot))] + \frac{1}{N^2(N-1)}\mathbb{E}\left[\sum_{j=1}^N |b(X_t^{i,\delta} - X_t^{j,\delta})|^2\right] \\ &\le \frac{4L^2}{N-1}C + \frac{16\|\gamma\|_{\infty}^2}{N-1} + \frac{1}{N^2}(8CL^2 + 4\|\gamma\|_{\infty}^2) \\ &\le N^{-1}C_1 < \infty, \end{aligned} $$
where $C_1$ depends on $\|\gamma\|_{\infty}$, $L$ and the second moment bound $C$. Similarly, it holds

$$
\begin{aligned}
\mathbb{E}[A_t^{i, \delta}] &\le \mathbb{E} \left[ |b * \bar{\mu}_t^{\delta}(\bar{X}_t^{i, \delta}) - \frac{1}{N-1} \sum_{j=1}^{N} b(\bar{X}_t^{i, \delta} - \bar{X}_t^{j, \delta})| \right] \\
&\quad + \left( \frac{1}{N-1} - \frac{1}{N} \right) \sum_{j=1}^{N} \mathbb{E} \left[ |b(\bar{X}_t^{i, \delta} - \bar{X}_t^{j, \delta})| \right] \\
&\le \frac{\sqrt{2L}}{\sqrt{N-1}} C^{1/2} + \frac{\sqrt{8} \|\gamma\|_{\infty}}{\sqrt{N-1}} + \frac{1}{N} \left( \sqrt{2C^{1/2}} L + \|\gamma\|_{\infty} \right) \\
&\le N^{-1/2} C_2 < \infty,
\end{aligned}
$$

$$ \text{where } C_2 = 2LC^{1/2} + 4\|\gamma\|_{\infty} + (\sqrt{2}C^{1/2} + \|\gamma\|_{\infty}). \quad \square $$

To control $(\{\tilde{r}_t^{i, \delta}\}_{i=1}^N)_{t \ge 0}$, we consider $(\{r_t^{i, \delta, \epsilon}\}_{i=1}^N)_{t \ge 0}$ given as solution of

$$
\begin{aligned}
dr_t^{i, \delta, \epsilon} = \bar{b}(r_t^{i, \delta, \epsilon})dt + \frac{1}{N} \sum_{j=1}^{N} 2\|\gamma\|_{\infty} \mathrm{rc}^{\epsilon}(r_t^{j, \delta, \epsilon})dt + (A_t^{i, \delta} + \frac{1}{N} \sum_{k=1}^{N} A_t^{k, \delta})dt \\
+ 2\sqrt{1 + \frac{1}{N}} \mathrm{rc}^{\delta}(r_t^{i, \delta, \epsilon})dW_t^{i, \delta}
\end{aligned}
\qquad (55)
$$

with initial condition $r_0^{i, \delta, \epsilon} = \tilde{r}_0^{i, \delta}$ for all $i = 1, \dots, N$, $A_t^{i, \delta}$ given in (53) and $W_t^{i, \delta}$ given in (54).

By [45, Theorem 2.2], under **B1** and **B3**, $(\{U_t^{i, \delta, \epsilon}\}_{i=1}^N)_{t \ge 0} = (\{\bar{X}_t^{i, \delta}, Y_t^{i, \delta}, r_t^{i, \delta, \epsilon}\}_{i=1}^N)_{t \ge 0}$ exists and is unique, where $(\{\bar{X}_t^{i, \delta}, \bar{Y}_t^{i, \delta}\}_{i=1}^N)_{t \ge 0}$ solves uniquely (51), $(\{\tilde{r}_t^{i, \delta}\}_{i=1}^N)_{t \ge 0}$ and $(\{r_t^{i, \delta, \epsilon}\}_{i=1}^N)_{t \ge 0}$ solve uniquely (52) and (55), respectively, with $(\{W_t^{i, \delta}\}_{i=1}^N)_{t \ge 0}$ given by (54).

**Lemma 17.** Assume **B1** and **B3**. Then for any $i = 1, \dots, N$, $|\bar{X}_t^{i, \delta} - Y_t^{i, \delta} - \frac{1}{N} \sum_j (\bar{X}_t^{j, \delta} - Y_t^{j, \delta})| = \tilde{r}_t^{i, \delta} \le r_t^{i, \delta, \epsilon}$, almost surely for all $t \ge 0$ and $\epsilon < \epsilon_0$.

*Proof.* Note, that both processes $(\{\tilde{r}_t^{i, \delta}\}_{i=1}^N)_{t \ge 0}$ and $(\{r_t^{i, \delta, \epsilon}\}_{i=1}^N)_{t \ge 0}$ have the same initial condition and are driven by the same noise. Since the drift for $(\{r_t^{i, \delta, \epsilon}\}_{i=1}^N)_{t \ge 0}$ is larger than the drift for $(\{\tilde{r}_t^{i, \delta}\}_{i=1}^N)_{t \ge 0}$ for $\epsilon < \epsilon_0$ by (20), we can conclude $\tilde{r}_t^{i, \delta} \le r_t^{i, \delta, \epsilon}$ almost surely for all $t \ge 0$, $\epsilon < \epsilon_0$ and $i = 1, \dots, N$ by Lemma 21. $\square$

*Proof of Theorem 8.* Consider the process $(\{U_t^{i, \delta, \epsilon}\}_{i=1}^N)_{t \ge 0} = (\{\bar{X}_t^{i, \delta}, Y_t^{i, \delta}, r_t^{i, \delta, \epsilon}\}_{i=1}^N)_{t \ge 0}$ on $\mathbb{R}^{N(2d+1)}$ for each $\epsilon, \delta > 0$. We denote by $\mathbb{P}^{\delta, \epsilon}$ the law of $\{U^{\delta, \epsilon}\}_{i=1}^N$ on $\mathcal{C}(\mathbb{R}_+, \mathbb{R}^{N(2d+1)})$. We define the canonical projections $\mathbf{X}, \mathbf{Y}, \mathbf{r}$ onto the first $Nd$, second $Nd$ and last $N$ components.

By **B1** and **B3** it holds in the same line as in the proof of Lemma 22 for each $T > 0$

$$ E[|\{U_{t_2}^{i, \delta, \epsilon} - U_{t_1}^{i, \delta, \epsilon}\}_{i=1}^N|^4] \le C|t_2 - t_1|^2 \quad \text{for } t_1, t_2 \in [0, T], \qquad (56) $$
for some constant C depending on T, L, ||γ||Lip, ||γ||∞, N and on the fourth moment of μ₀ and ν₀. Note that we used here that the additional drift terms (A_t^i,δ)_{t≥0} occurring in the SDE of ({r_t^i,δ,ε}^N_{i=1})_{t≥0} are Lipschitz continuous in ({\bar{X}_t^i,δ}^N_{i=1})_{t≥0}. Then as in the proofs of Lemma 22 and Lemma 23, P^δ,ε is tight and converges weakly along a subsequence to a measure P by Kolmogorov's continuity criterion, cf. [38, Corollary 14.9].

As in Lemma 22 the law P_T^δ,ε of ({U_t^i,δ,ε}^N_{i=1})_{0≤t≤T} on C([0, T], ℝ^{N(2d+1)}) is tight for each T > 0 by [38, Corollary 14.9] and for each ε > 0 there exists a subsequence δ_n → 0 such that (P_T^δ,ε)_{n∈N} on C([0, T], ℝ^{N(2d+1)}) converge to a measure P_T^ε on C([0, T], ℝ^{N(2d+1)}). By a diagonalization argument and since {P_T^ε : T ≥ 0} is a consistent family, cf. [38, Theorem 5.16], there exists a probability measure P^ε on C(ℝ+, ℝ^{N(2d+1)}) such that for all ε there exists a subsequence δ_n such that (P^δ,ε)_{n∈N} converges along this subsequence to P^ε, i.e., the family (P^δ,ε)_{δ>0} is tight. As in the proof of Lemma 23 we repeat this argument for the family of measures (P^ε)_{ε>0}. Hence, there exists a subsequence ε_n → 0 such that (P^ε_n)_{n∈N} converges to a measure P. Let ({\bar{X}_t^i, Y_t^i, r_t^i}^N_{i=1})_{t≥0} be some process on ℝ^{N(2d+1)} with distribution P on (Ω, F, P).

Since ({\bar{X}_t^i,δ}^N_{i=1})_{t≥0} and ({Y_t^i,δ}^N_{i=1})_{t≥0} are solutions that are unique in law, we have that for any δ, ε > 0, P^δ,ε ◦ (X)^{-1} = P ◦ (X)^{-1} and P^δ,ε ◦ (Y)^{-1} = P ◦ (Y)^{-1}. Hence, P ◦ (X, Y)^{-1} is a coupling of (1) and (3). In addition, the statement of Lemma 17 carries over to the limiting process ({r_t^i}^N_{i=1})_{t≥0}, since by the weak convergence and the Portmanteau, it is preserved under weak convergence by Portmanteau theorem, P(|\tilde{X}_t^i - \tilde{Y}_t^i| ≤ r_t^i for i = 1, ..., N) ≥ limsup_{ε→0} limsup_{δ→0} P(|\tilde{X}_t^i,δ - \tilde{Y}_t^i,δ| ≤ r_t^i,δ,ε for i = 1, ..., N) = 1, where \tilde{X}_t^i = \bar{X}_t^i - (1/N) \sum_{j=1}^N \bar{X}_t^j and \tilde{Y}_t^i = \bar{X}_t^i - (1/N) \sum_{j=1}^N \bar{Y}_t^j for all t ≥ 0 and i = 1, ..., N.

Similarly to the proof of Lemma 22 and Lemma 23 there exist an extended underlying probability space and N i.i.d.one-dimensional Brownian motion ({W_t^i}^N_{i=1})_{t≥0} such that ({r_t^i, W_t^i}^N_{i=1})_{t≥0} is a solution of

$$
\begin{align*}
dr_t^i = \bar{b}(r_t^i)dt + \frac{1}{N} \sum_{j=1}^{N} 2 \|\gamma\|_{\infty} \mathbb{1}_{(0,\infty)}(r_t^j)dt + \left(A_t^i + \frac{1}{N} \sum_{k=1}^{N} A_t^k\right)dt \\
+ 2\sqrt{1 + \frac{1}{N}} \mathbb{1}_{(0,\infty)}(r_t^i) dW_t^i,
\end{align*}
$$

$$
\text{where } A_t^i = |b * \bar{\mu}_t(\bar{X}_t^i) - \frac{1}{N} \sum_{j=1}^N b(\bar{X}_t^i - \bar{X}_t^j)|.
$$

Using Ito-Tanaka formula, c.f. [51, Chapter 6, Theorem 1.1], and f' is absolutely
continuous, we obtain for $f$ defined in (37) with $\tilde{b}(r) = (\kappa(r) - L)r$ and $a = 2\|\gamma\|_{\infty}$,

$$d\left(\frac{1}{N}\sum_{i=1}^{N} f(r_t^i)\right) = \frac{1}{N}\sum_{i=1}^{N} \left(\bar{b}(r_t^i)f'(r_t^i) + f''(r_t^i)2\frac{N+1}{N}1_{(0,\infty)}(r_t^i)\right)dt \\ + \frac{1}{N^2}\sum_{i=1}^{N}\sum_{j=1}^{N} 2f'(r_t^i)\|\gamma\|_{\infty}1_{(0,\infty)}(r_t^j)dt \\ + \frac{1}{N}\sum_{i=1}^{N} f'(r_t^i)2\sqrt{1+\frac{1}{N}}1_{(0,\infty)}(r_t^i)dW_t^i + \frac{1}{N}\sum_{i=1}^{N} f'(r_t^i)\left(A_t^i + \frac{1}{N}\sum_{k=1}^{N} A_t^k\right)dt.$$

Taking expectation, we get using $f'(r) \le 1$ for all $r \ge 0$,

$$\frac{d}{dt}\mathbb{E}\left[\frac{1}{N}\sum_{i=1}^{N} f(r_t^i)\right] \le \frac{1}{N}\sum_{i=1}^{N} \left\{ \mathbb{E}\left[\bar{b}(r_t^i)f'(r_t^i) + 2\frac{N+1}{N}(f''(r_t^i) - f''(0))\right] \right. \\ \left. + \mathbb{E}\left[2\left(\|\gamma\|_{\infty} + \frac{N+1}{N}f''(0)\right)1_{(0,\infty)}(r_t^i)\right] + \mathbb{E}\left[2A_t^i\right] \right\}. \quad (57)$$

By (39) and (40), the first two terms are bounded by $-\tilde{c}\frac{1}{N}\sum_i f(r_t^i)$ with $\tilde{c}$ given in (17). By Lemma 16 the last term in (57) is bounded by

$$2E[A_t^i] \le \tilde{C}N^{-1/2},$$

where

$$\tilde{C} = 2C_2 = 4LC^{1/2} + 8\|\gamma\|_{\infty} + 2(\sqrt{2}C^{1/2}L + \|\gamma\|_{\infty}). \quad (58)$$

Hence, we obtain

$$\frac{d}{dt}\mathbb{E}\left[\frac{1}{N}\sum_{i} f(r_t^i)\right] \le -\tilde{c}\frac{1}{N}\sum_{i} \mathbb{E}[f(r_t^i)] + \tilde{C}N^{-1/2}$$

for $t \ge 0$ which leads by Gronwall's lemma to

$$\mathbb{E}\left[\frac{1}{N}\sum_{i} f(r_t^i)\right] \le e^{-\tilde{c}t}\mathbb{E}\left[\frac{1}{N}\sum_{i} f(r_0^i)\right] + \frac{1}{\tilde{c}}\tilde{C}N^{-1/2}.$$

For an arbitrary coupling $\xi \in \Gamma(\bar{\mu}_0^{\otimes N}, \nu_0^{\otimes N})$, we have

$$W_{f,N}((\bar{\mu}_t)^{\otimes N}, \nu_t^N) \le e^{-\tilde{c}t} \int_{\mathbb{R}^{2Nd}} \frac{1}{N} \sum_{i=1}^{N} f\left(\left|x^i - y^i - \frac{1}{N}\sum_{j=1}^{N} (x^j - y^j)\right|\right) \xi(dx dy) + \frac{\tilde{C}}{\tilde{c}N^{1/2}},$$

as $\mathbb{E}[f(r_0^i)] \le \int_{\mathbb{R}^{2Nd}} \frac{1}{N} \sum_{i=1}^{N} f(|x^i - y^i - \frac{1}{N}\sum_{j=1}^{N} (x^j - y^j)|)\xi(dx dy)$. Taking the infimum over all couplings $\xi \in \Gamma(\bar{\mu}_0^{\otimes N}, \nu_0^{\otimes N})$ gives the first bound. By (38), the second bound follows. $\square$
## 6.3 Proof of Section 3

### 6.3.1 Proof of Theorem 3

We show Theorem 3 via a family of stochastic differential equations, indexed by $n, m \in \mathbb{N}$, with Lipschitz continuous coefficients,

$$ 
\begin{aligned}
dr_t^{n,m} &= (\tilde{b}(r_t^{n,m}) + P_t^{n,m}(g^m))dt + 2\theta^n(r_t^{n,m})dW_t \\
ds_t^{n,m} &= (\hat{b}(s_t^{n,m}) + \hat{P}_t^{n,m}(h^m))dt + 2\theta^n(s_t^{n,m})dW_t, \quad \text{Law}(r_0^{n,m}, s_0^{n,m}) = \eta_{n,m},
\end{aligned} 
\quad (59) $$

where $P_t^{n,m} = \text{Law}(r_t^{n,m})$, $\hat{P}_t^{n,m} = \text{Law}(s_t^{n,m})$ and $P_t^{n,m}(g^m) = \int_{\mathbb{R}_+} g^m(r)P_t^{n,m}(dx)$ and $\hat{P}_t^{n,m}(h^m) = \int_{\mathbb{R}_+} h^m(r)\hat{P}_t^{n,m}(dx)$ for some measurable functions $(g^m)_{m \in \mathbb{N}}$ and $(h^m)_{m \in \mathbb{N}}$, and where $\eta_{n,m} \in \Gamma(\mu_{n,m}, \nu_{n,m})$ for $\mu_{n,m}, \nu_{n,m} \in \mathcal{P}(\mathbb{R}_+)$. We identify the weak limit for $n \to \infty$ as solution of a family of stochastic differential equations, indexed by $m \in \mathbb{N}$, given by

$$ 
\begin{aligned}
dr_t^m &= (\tilde{b}(r_t^m) + P_t^m(g^m))dt + 2\mathbb{1}_{(0,\infty)}(r_t^m)dW_t \\
ds_t^m &= (\hat{b}(s_t^m) + \hat{P}_t^m(h^m))dt + 2\mathbb{1}_{(0,\infty)}(s_t^m)dW_t, \quad \text{Law}(r_0^m, s_0^m) = \eta_m.
\end{aligned} 
\quad (60) $$

with $P_t^m = \text{Law}(r_t^m)$ and $\hat{P}_t^m = \text{Law}(s_t^m)$, and where $\eta_m \in \Gamma(\mu_m, \nu_m)$ for $\mu_m, \nu_m \in \mathcal{P}(\mathbb{R}_+)$. Taking the limit $m \to \infty$, we show in the next step that the solution of (60) converges to a solution of (25).

We assume for $(g^m)_{m \in \mathbb{N}}$, $(h^m)_{m \in \mathbb{N}}$, $(\theta^n)_{n \in \mathbb{N}}$ and the initial distributions:

**H5.** $(g^m)_{m \in \mathbb{N}}$ and $(h^m)_{m \in \mathbb{N}}$ are sequences of non-decreasing non-negative uniformly bounded Lipschitz continuous functions such that for all $r \ge 0$, $g^m(r) \le g^{m+1}(r)$ and $h^m(r) \le h^{m+1}(r)$ and $\lim_{m \to +\infty} g^m(r) = g(r)$ and $\lim_{m \to +\infty} h^m(r) = h(r)$ where $g, h$ are left-continuous non-negative non-decreasing bounded functions. In addition, there exists $K_m < \infty$ for any $m$ such that for all $r, s \in \mathbb{R}$

$$ |g^m(r) - g^m(s)| \le K_m|r - s| \quad \text{and} \quad |h^m(r) - h^m(s)| \le K_m|r - s|. $$

**H6.** $(\theta^n)_{n \in \mathbb{N}}$ is a sequence of Lipschitz continuous functions from $\mathbb{R}_+$ to $[0, 1]$ with $\theta^n(0) = 0$, $\theta^n(r) = 1$ for all $r \ge 1/n$ and $\theta^n(r) > 0$ for all $r > 0$.

**H7.** $(\mu_{n,m})_{m,n \in \mathbb{N}}$, $(\nu_{n,m})_{m,n \in \mathbb{N}}$, $(\mu_m)_{m \in \mathbb{N}}$, $(\nu_m)_{m \in \mathbb{N}}$ are families of probability distributions on $\mathbb{R}_+$ and $(\eta_{n,m})_{n,m \in \mathbb{N}}$, $(\eta_m)_{m \in \mathbb{N}}$ families of probability distributions on $\mathbb{R}_+^2$ such that for any $n, m \in \mathbb{N}$, $\eta_{n,m} \in \Gamma(\mu_{n,m}, \nu_{n,m})$ and $\eta_m \in \Gamma(\mu_m, \nu_m)$ and for any $m \in \mathbb{N}$, $(\eta_{n,m})_{n \in \mathbb{N}}$ converges weakly to $\eta_m$ and $(\eta_m)_{m \in \mathbb{N}}$ converges weakly to $\eta$. Further, the p-th order moments of $(\mu_{n,m})_{n,m \in \mathbb{N}}$, $(\nu_{n,m})_{n,m \in \mathbb{N}}$, $(\mu_m)_{m \in \mathbb{N}}$ and $(\nu_m)_{m \in \mathbb{N}}$ are uniformly bounded for $p > 2$ given in **H3**.

Note that by **H5** for any non-decreasing sequence $(u_m)_{m \in \mathbb{N}}$, which converges to $u \in \mathbb{R}_+$, $g^m(u_m)$ and $h^m(u_m)$ converge to $g(u)$ and $h(u)$, respectively. More precisely, it holds for
for all $m \in \mathbb{N}$, $g^m(u_m) - g(u) \le 0$ and for $m \ge n$, $g^m(u_m) \ge g^m(u_n)$ and therefore, $\lim_{m \to \infty} g^m(u_m) - g(u) \ge \lim_{n \to \infty} \lim_{m \to \infty} = \lim_{n \to \infty} g(u_n) - g(u) = 0$ by left-continuity of $g$. Hence, $\lim_{m \to \infty} g^m(u_m) - g(u) = 0$ and analogously $\lim_{m \to \infty} h^m(u_m) - h(u) = 0$. By $\mathbf{H}5$, $\Gamma = \max(\|h\|_{\infty}, \|g\|_{\infty})$ is a uniform upper bound of $(g^m)_{m \in \mathbb{N}}$ and $(h^m)_{m \in \mathbb{N}}$.

Consider a probability space $(\Omega_0, \mathcal{A}_0, \mathcal{Q})$ and a one-dimensional Brownian motion $(W_t)_{t \ge 0}$. Under $\mathbf{H}5$, $\mathbf{H}6$ and $\mathbf{H}7$, for all $m, n \in \mathbb{N}$, there exists random variables $r^{n,m}, s^{n,m} : \Omega_0 \to \mathbb{W}$ for each $n, m$ such that $(r_t^{n,m}, s_t^{n,m})_{t \ge 0}$ is a unique strong solution to (59) associated to $(W_t)_{t \ge 0}$ by [45, Theorem 2.2]. We denote by $\mathbb{P}^{n,m} = \mathcal{Q} \circ (r^{n,m}, s^{n,m})^{-1}$ the corresponding distribution on $\mathbb{W} \times \mathbb{W}$.

Before studying the two limits $n, m \to \infty$ and proving Theorem 3, we state a modification of the comparison theorem by Ikeda and Watanabe to compare two solutions of (59), cf. [35, Section VI, Theorem 1.1].

**Lemma 18.** Let $(r_t^{n,m}, s_t^{n,m})_{t \ge 0}$ be a solution of (59) for fixed $n, m \in \mathbb{N}$. Assume $\mathbf{H}1$, $\mathbf{H}5$ and $\mathbf{H}6$. If $Q[r_0^{n,m} \le s_0^{n,m}] = 1$, $\tilde{b}(r) \le \hat{b}(r)$ and $g^m(r) \le h^m(r)$ for any $r \in \mathbb{R}_+$, then

$$Q[r_t^{n,m} \le s_t^{n,m} \text{ for all } t \ge 0] = 1. \qquad (61)$$

*Proof.* For simplicity, we drop the dependence on $n, m$ in $(r_t^{n,m})$ and $(s_t^{n,m})$. Denote by $\rho$ the Lipschitz constant of $\theta^n$. Let $(a_k)_{k \in \mathbb{N}}$ be a decreasing sequence, $1 > a_1 > a_2 > \dots > a_k > \dots > 0$, such that $\int_{a_1}^1 \rho^{-2}x^{-1}dx = 1$, $\int_{a_2}^{a_1} \rho^{-2}x^{-1}dx = 2, \dots$, $\int_{a_k}^{a_{k-1}} \rho^{-2}x^{-1}dx = k$. We choose a sequence $\Psi_k(u), k = 1, 2, \dots$, of continuous functions such that its support is contained in $(a_k, a_{k-1})$, $\int_{a_k}^{a_{k-1}} \Psi_k(u)du = 1$ and $0 \le \Psi_k(u) \le 2/k \cdot \rho^{-2}u^{-2}$. Such a function exists. We set

$$\varphi_k(x) = \begin{cases} \int_0^x dy \int_0^y \Psi_k(u) du & \text{if } x \ge 0, \\ 0 & \text{if } x < 0. \end{cases}$$

Note that for any $k \in \mathbb{N}$, $\varphi_k \in C^2(\mathbb{R}_+)$, $|\varphi_k'(x)| \le 1$, $\varphi_k(x) \to x_+$ as $k \uparrow \infty$ and $\varphi_k'(x) \uparrow 1_{(0,\infty)}(x)$. Applying Ito's formula to $\varphi_k(r_t - s_t)$, we obtain

$$\varphi_k(r_t - s_t) = \varphi_k(r_0 - s_0) + I_1(k) + I_2(k) + I_3(k),$$

where

$$I_1(k) = \int_0^t \varphi_k'(r_u - s_u)[\theta^n(r_u) - \theta^n(s_u)]dB_u,$$

$$I_2(k) = \int_0^t \varphi_k'(r_u - s_u)[b(r_u) - \hat{b}(s_u) + P_u(g^m) - \hat{P}_u(h^m)]du,$$

$$I_3(k) = \frac{1}{2} \int_0^t \varphi_k''(r_u - s_u)[\theta^n(r_u) - \theta^n(s_u)]^2 du,$$
with $P_u = Q \circ r_u^{-1}$ and $\hat{P}_u = Q \circ s_u^{-1}$. It holds by boundedness and Lipschitz continuity of $\theta^n$

$$ \mathbb{E}[I_1(k)] = 0, \quad \text{and} \quad \mathbb{E}[I_3(k)] \le \frac{1}{2} \mathbb{E} \left[ \int_0^t \varphi_k''(r_u - s_u) \rho^2 |r_u - s_u|^2 du \right] \le \frac{t}{k}. $$

We note that by $\mathbf{H}5$ $\mathbb{E}[(g^m(r_u) - h^m(s_u))1_{r_u-s_u<0}] \le 0$ and

$$ \begin{aligned} \mathbb{E}[(g^m(r_u) - h^m(s_u))1_{r_u-s_u \ge 0}] &\le \mathbb{E}[(g^m(r_u) - g^m(s_u) + g^m(s_u) - h^m(s_u))1_{r_u-s_u \ge 0}] \\ &\le \mathbb{E}[(g^m(r_u) - g^m(s_u))1_{r_u-s_u \ge 0}] \\ &\le K_m \mathbb{E}[|r_u - s_u|1_{r_u-s_u \ge 0}] \end{aligned} \quad (62) $$

by Lipschitz continuity of $g^m$, by $g^m(r) \le h^m(r)$ and since $g^m$ and $h^m$ are non-decreasing. Hence for $I_2$, we obtain

$$ \begin{aligned} I_2(k) &= \int_0^t \varphi'_k(r_u - s_u)[\tilde{b}(r_u) - \hat{b}(r_u) + \hat{b}(r_u) - \hat{b}(s_u)]du \\ &\quad + \int_0^t \varphi'_k(r_u - s_u)(\mathbb{E}[(g^m(r_u) - h^m(s_u))1_{r_u-s_u \ge 0}] + \mathbb{E}[(g^m(r_u) - h^m(s_u))1_{r_u-s_u < 0}])du \\ &\le \int_0^t \varphi'_k(r_u - s_u)\tilde{L}|r_u - s_u|du + \int_0^t \varphi'_k(r_u - s_u)K_m\mathbb{E}[|r_u - s_u|1_{r_u-s_u \ge 0}]du. \end{aligned} $$

Taking the limit $k \to \infty$ and using that $\mathbb{E}[r_0 - s_0] = 0$, we obtain

$$ \mathbb{E}[(r_t - s_t)_+] \le \tilde{L} \mathbb{E} \left[ \int_0^t (r_u - s_u)_+ du \right] + K_m \mathbb{E} \left[ \int_0^t 1_{(0,\infty)}(r_u - s_u) \mathbb{E}[(r_u - s_u)_+] du \right], \quad (63) $$

by the monotone convergence theorem and since $(\varphi'_k)_{k \in \mathbb{N}}$ is a monotone increasing sequence which converges pointwise to $1_{(0,\infty)}(x)$. Assume there exists $t^* = \inf\{t \ge 0 : \mathbb{E}[(r_t - s_t)_+] > 0\} < \infty$. Then, $\int_0^{t^*} \mathbb{E}[(r_u - s_u)_+]du > 0$ or $\int_0^{t^*} \mathbb{E}[1_{(0,\infty)}(r_u - s_u)]\mathbb{E}[(r_u - s_u)_+]du > 0$. By definition of $t^*$, $\mathbb{E}[(r_u - s_u)_+] = 0$ for all $u < t^*$ and hence both terms are zero. This contradicts the definition of $t^*$. Hence, (61) holds. $\square$

Next, we show that the distribution of the solution of (59) converges as $n \to \infty$.

**Lemma 19.** Assume that $\tilde{b}$, $\hat{b}$, $g$ and $h$ satisfy $\mathbf{H}1$ and $\mathbf{H}2$. Let $\eta \in \Gamma(\mu, \nu)$ where the probability measures $\mu$ and $\nu$ on $\mathbb{R}_+$ satisfy $\mathbf{H}3$. Assume that $(g^m)_{m \in \mathbb{N}}$, $(h^m)_{m \in \mathbb{N}}$, $(\theta^n)_{n \in \mathbb{N}}$, $(\mu_{n,m})_{m,n \in \mathbb{N}}$, $(\nu_{n,m})_{m,n \in \mathbb{N}}$ and $(\eta_{n,m})_{m,n \in \mathbb{N}}$ satisfy condition $\mathbf{H}5$, $\mathbf{H}7$ and $\mathbf{H}6$. Then for any $m \in \mathbb{N}$, there exists a random variable $(r^m, s^m)$ defined on some probability space $(\Omega^m, \mathcal{A}^m, P^m)$ with values in $\mathbb{W} \times \mathbb{W}$, such that $(r_t^m, s_t^m)_{t \ge 0}$ is a weak solution of the stochastic differential equation (60). More precisely, for all $m \in \mathbb{N}$ the sequence of laws $Q \circ (r^{n,m}, s^{n,m})^{-1}$ converges weakly to the distribution $P^m \circ (r^m, s^m)^{-1}$. If additionally,

$$ \begin{aligned} \tilde{b}(r) \le \hat{b}(r) \quad \text{and} \quad g^m(r) \le h^m(r), & \quad \text{for any } r \in \mathbb{R}_+ \text{ and} \\ Q[r_0^{n,m} \le s_0^{n,m}] = 1 & \quad \text{for any } n, m \in \mathbb{N}, \end{aligned} $$

then $P^m[r_t^m \le s_t^m \text{ for all } t \ge 0] = 1$.
**Proof.** Fix $m \in \mathbb{N}$. The proof is divided in three parts. First we show tightness of the sequences of probability measures. Then we identify the limit of the sequence of stochastic processes. Finally, we compare the two limiting processes.

**Tightness:** We show that the sequence of probability measures $(\mathbb{P}^{n,m})_{n \in \mathbb{N}}$ on $(\mathbb{W} \times \mathbb{W}, \mathcal{B}(\mathbb{W}) \otimes \mathcal{B}(\mathbb{W}))$ is tight by applying Kolmogorov's continuity theorem. Consider $p > 2$ such that the $p$-th moment in $\mathbf{H}3$ and $\mathbf{H}7$ are uniformly bounded. Fix $T > 0$. Then the $p$-th moment of $r_t^{n,m}$ for $t < T$ can be bounded using Ito's formula,

$$ 
\begin{aligned}
d|r_t^{n,m}|^p &\le p|r_t^{n,m}|^{p-2}\langle r_t^{n,m}, (\tilde{b}(r_t^{n,m}) + P_t^{n,m}(g^m))\rangle dt + 2\theta^n(r_t^{n,m})p|r_t^{n,m}|^{p-2}r_t^{n,m}dW_t \\
&\quad + p(p-1)|r_t^{n,m}|^{p-2}2\theta^n(r_t^{n,m})^2dt \\
&\le p(|r_t^{n,m}|^p \tilde{L} + \Gamma|r_t^{n,m}|^{p-1} + 2(p-1)|r_t^{n,m}|^{p-2})dt + 2\theta^n(r_t^{n,m})p(r_t^{n,m})^{p-1}dW_t \\
&\le p(\tilde{L} + \Gamma + 2(p-1))|r_t^{n,m}|^p dt + p(\Gamma + 2(p-1))dt + 2\theta^n(r_t^{n,m})p(r_t^{n,m})^{p-1}dW_t,
\end{aligned}
$$

where $\Gamma = \max(\|g\|_{\infty}, \|h\|_{\infty})$. Taking expectation yields

$$ \frac{d}{dt}\mathbb{E}[|r_t^{n,m}|^p] \le p(\tilde{L} + \Gamma + 2(p-1))\mathbb{E}|r_t^{n,m}|^p + p(\Gamma + 2(p-1)). $$

Then by Gronwall's lemma

$$ \sup_{t \in [0,T]} \mathbb{E}[|r_t^{n,m}|^p] \le e^{p(\tilde{L}+\Gamma+2(p-1))^T} (\mathbb{E}[|r_0^{n,m}|^p] + Tp(\Gamma+2(p-1))) < C_p < \infty, \quad (64) $$

where $C_p$ depends on $T$ and the $p$-th moment of the initial distribution, which is finite by **H**6. Similarly, it holds $\sup_{t \in [0,T]} \mathbb{E}[|s_t^{n,m}|^p] < C_p$ for $t \le T$. Using this moment bound, it holds for all $t_1, t_2 \in [0, T]$ by **H**1, **H**5 and **H**6,

$$ 
\begin{aligned}
\mathbb{E}[|r_{t_2}^{n,m} - r_{t_1}^{n,m}|^p] &\le C_1(p)\left(\mathbb{E}\left[\int_{t_1}^{t_2} \tilde{b}(r_u^{n,m}) + P_u^{n,m}(g^m)du\right]^p + \mathbb{E}\left[\int_{t_1}^{t_2} 2\theta^n(r_u^{n,m})dW_u\right]^p\right) \\
&\le C_2(p)\left(\left(\mathbb{E}\left[\frac{\tilde{L}^p}{|t_2-t_1|}\int_{t_1}^{t_2} |r_u^{n,m}|^p du\right] + \Gamma^p\right)|t_2-t_1|^p + \mathbb{E}\left[\int_{t_1}^{t_2} 2\theta^n(r_u^{n,m})du\right]^{p/2}\right) \\
&\le C_2(p)\left(\left(\frac{\tilde{L}^p}{|t_2-t_1|}\int_{t_1}^{t_2} \mathbb{E}[|r_u^{n,m}|^p]du + \Gamma^p\right)|t_2-t_1|^p + 2^{p/2}|t_2-t_1|^{p/2}\right) \\
&\le C_3(p, T, \tilde{L}, \Gamma, C_p)|t_2-t_1|^{p/2},
\end{aligned}
$$

where $C_i(\cdot)$ are constants depending on the stated argument and which are independent of $n, m$. Note that in the second step, we used Burkholder-Davis-Gundy inequality, see [50, Chapter IV, Theorem 48]. It holds similarly, $\mathbb{E}[|s_{t_2}^{n,m} - s_{t_1}^{n,m}|^p] \le C_3(p, T, \tilde{L}, \Gamma, C_p)|t_2-t_1|^{p/2}$. Hence,

$$ \mathbb{E}[|(r_{t_2}^{n,m}, s_{t_2}^{n,m}) - (r_{t_1}^{n,m}, s_{t_1}^{n,m})|^p] \le C_4(p, T, \tilde{L}, \Gamma, C_p)|t_2 - t_1|^{p/2} \quad (65) $$
for all $t_1, t_2 \in [0, T]$. Hence, by Kolmogorov's continuity criterion, cf. [38, Corollary 14.9], there exists a constant $\tilde{C}$ depending on $p$ and $\gamma$ such that

$$ \mathbb{E}[(r^{n,m}, s^{n,m})]_{\gamma}^{p} \leq \tilde{C} \cdot C_4(p, T, \tilde{L}, \Gamma, C_p), \quad (66) $$

where $[\cdot]_{\gamma}^{p}$ is defined by $[x]_{\gamma} = \sup_{t_1, t_2 \in [0, T]} \frac{|x(t_1) - x(t_2)|}{|t_1 - t_2|^{\gamma}}$ and $(r_t^{n,m}, s_t^{n,m})_{n \in \mathbb{N}, t \ge 0}$ is tight in $C([0, T], \mathbb{R}^2)$. Hence, for each $T > 0$ there exists a subsequence $n_k \to \infty$ and a probability measure $\mathbb{P}_T^m$ on $C([0, T], \mathbb{R}^2)$. Since $\{\mathbb{P}_T^m\}_T$ is a consistent family, there exists by [38, Theorem 5.16] a probability measure $\mathbb{P}^m$ on $(\mathbb{W} \times \mathbb{W}, \mathcal{B}(\mathbb{W}) \otimes \mathcal{B}(\mathbb{W}))$ such that there is a subsequence $(n_k)_{k \in \mathbb{N}}$ such that $\mathbb{P}^{n_k,m}$ converges along this subsequence to $\mathbb{P}^m$. Note that here we can take by a diagonalization argument the same subsequence $(n_k)_{k \in \mathbb{N}}$ for all $m$.

**Characterization of the limit measure:** In the following we drop for simplicity the index $k$ in the subsequence. Denote by $(\mathbf{r}_t, \mathbf{s}_t)(\omega) = \omega(t)$ the canonical process on $\mathbb{W} \times \mathbb{W}$. Since $\mathbb{P}^{n,m} \circ (\mathbf{r}_0, \mathbf{s}_0)^{-1} = \eta_{n,m}$ converges weakly to $\eta_m$ by **H7**, it holds $\mathbb{P}^m \circ (\mathbf{r}_0, \mathbf{s}_0)^{-1} = \eta_m$. We define the maps $M^{n,m}, N^{n,m}: \mathbb{W} \times \mathbb{W} \to \mathbb{W}$ by

$$ M_t^{n,m} = \mathbf{r}_t - \mathbf{r}_0 - \int_0^t (\tilde{b}(\mathbf{r}_u) + P_u^n(g^m))du \quad \text{and} \quad N_t^{n,m} = \mathbf{s}_t - \mathbf{s}_0 - \int_0^t (\hat{b}(\mathbf{s}_u) + \hat{P}_u^n(h^m))du, $$

where $P_u^n = \mathbb{P}^{n,m} \circ (\mathbf{r}_u)^{-1}$ and $\hat{P}_u^n = \mathbb{P}^{n,m} \circ (\mathbf{s}_u)^{-1}$. For each $m, n \in \mathbb{N}$, $(M_t^{n,m}, \mathcal{F}_t, \mathbb{P}^{n,m})$ and $(N_t^{n,m}, \mathcal{F}_t, \mathbb{P}^{n,m})$ are martingales with respect to the canonical filtration $\mathcal{F}_t = \sigma((\mathbf{r}_u, \mathbf{s}_u)_{0 \le u \le t})$ by Ito's formula and the moment estimate (64). Further the family $(M_t^{n,m}, \mathbb{P}^{n,m})_{n \in \mathbb{N}, t \ge 0}$ and $(N_t^{n,m}, \mathbb{P}^{n,m})_{n \in \mathbb{N}, t \ge 0}$ are uniformly integrable by Lipschitz continuity of $\tilde{b}$ and $\hat{b}$ and by boundedness of $g^m$ and $h^m$. Further, the mappings $M^{n,m}$ and $N^{n,m}$ are continuous in $\mathbb{W}$. We show that $\mathbb{P}^{n,m} \circ (\mathbf{r}, \mathbf{s}, M^{n,m}, N^{n,m})^{-1}$ converges weakly to $\mathbb{P}^m \circ (\mathbf{r}, \mathbf{s}, M^m, N^m)^{-1}$ as $n \to \infty$, where

$$ M_t^m = \mathbf{r}_t - \mathbf{r}_0 - \int_0^t (\tilde{b}(\mathbf{r}_u) + P_u(g^m))du \quad \text{and} \quad N_t^m = \mathbf{s}_t - \mathbf{s}_0 - \int_0^t (\hat{b}(\mathbf{s}_u) + \hat{P}_u(h^m))du, \quad (67) $$

with $P_u = \mathbb{P}^m \circ \mathbf{r}_u^{-1}$ and $\hat{P}_u = \mathbb{P}^m \circ \mathbf{s}_u^{-1}$. To show weak convergence to $\mathbb{P}^m \circ (\mathbf{r}, \mathbf{s}, M^m, N^m)^{-1}$, we note that $(M^m, N^m)$ is continuous in $\mathbb{W}$ and we consider for a Lipschitz continuous and bounded function $G: \mathbb{W} \to \mathbb{R}$,

$$ \left| \int_{\mathbb{W}} G(\omega) d\mathbb{P}^{n,m} \circ (M^{n,m})^{-1}(\omega) - \int_{\mathbb{W}} G(\omega) d\mathbb{P}^m \circ (M^m)^{-1}(\omega) \right| \\ \le \left| \int_{\mathbb{W}} G(\omega) d\mathbb{P}^{n,m} \circ (M^{n,m})^{-1}(\omega) - \int_{\mathbb{W}} G(\omega) d\mathbb{P}^{n,m} \circ (M^m)^{-1}(\omega) \right| \\ + \left| \int_{\mathbb{W}} G(\omega) d\mathbb{P}^{n,m} \circ (M^m)^{-1}(\omega) - \int_{\mathbb{W}} G(\omega) d\mathbb{P}^m \circ (M^m)^{-1}(\omega) \right|. $$
The second term converges to 0 as $n \to \infty$, since $(M^m)$ is continuous. For the first term it
holds

$$
\begin{align*}
\left| \int_{\mathbb{W}} G(\omega) d\mathbb{P}^{n,m} \circ (M^{n,m})^{-1}(\omega) - \int_{\mathbb{W}} G(\omega) d\mathbb{P}^{n,m} \circ (M^m)^{-1}(\omega) \right|
\end{align*}
\begin{align*}
&= \left| \int_{\mathbb{W}} (G \circ M^{n,m})(\omega) d\mathbb{P}^{n,m}(\omega) - \int_{\mathbb{W}} (G \circ M^m)(\omega) d\mathbb{P}^{n,m}(\omega) \right|
\end{align*}
\begin{align*}
&\le \|G\|_{\text{Lip}} \sup_{\omega \in \mathbb{W}} d_{\mathbb{W}}(M^{n,m}(\omega), M^m(\omega)),
\end{align*}
$$

where $d_W(f, g) = \sum_{k=1}^{\infty} \sup_{t \in [0,k]} 2^{-k} |f(t) - g(t)|$. This term converges to 0 for $n \to \infty$, since for all $T > 0$ and $\omega \in W$, for $n \to \infty$

$$
\sup_{t \in [0, T]} |M_t^{n,m}(\omega) - M_t^m(\omega)| \le \int_0^T \left| (\mathbb{P}^{n,m} \circ \mathbf{r}_s^{-1})(g^m) - (\mathbb{P}^m \circ \mathbf{r}_s^{-1})(g^m) \right| ds \to 0,
$$

by Lebesgue dominated convergence theorem, since *g* is bounded. Hence,

$$
\left| \int_{\mathbb{W}} G(\omega) d\mathbb{P}^{n,m} \circ (M^{n,m})^{-1}(\omega) - \int_{\mathbb{W}} G(\omega) d\mathbb{P}^{m} \circ (M^{m})^{-1}(\omega) \right| \to 0 \quad \text{for } n \to \infty,
$$

and similarly for $(N_t^{n,m})$, and therefore by the Portmanteau theorem [39, Theorem 13.16],
weak convergence of $\mathbb{P}^{n,m} \circ (\mathbf{r}, \mathbf{s}, M^{n,m}, N^{n,m})^{-1}$ to $\mathbb{P}^m \circ (\mathbf{r}, \mathbf{s}, M^m, N^m)^{-1}$ holds.

Let $G : \mathbb{W} \to \mathbb{R}_+$ be a $\mathcal{F}_s$-measurable, bounded, non-negative function. By uniformly integrability of $(M_t^{n,m}, \mathbb{P}^{n,m})_{n \in \mathbb{N}, t \ge 0}$, for any $s \le t$,

$$
\begin{align}
\mathbb{E}^m[G(M_t^m - M_s^m)] &= \mathbb{E}^m[G(\int_s^t (\tilde{b}(\mathbf{r}_u) + P_u(g^m))du)] \nonumber \\
&= \lim_{n \to \infty} \mathbb{E}^{n,m}[G(\int_s^t (\tilde{b}(\mathbf{r}_u) + P_u^n(g^m))du)] \tag{68} \\
&= \lim_{n \to \infty} \mathbb{E}^{n,m}[G(M_t^{n,m} - M_s^{n,m})] = 0, \nonumber
\end{align}
$$

and analogously for $(N_t^{n,m})_{t \ge 0}$ and hence, $(M_t^m, \mathcal{F}_t, \mathbb{P}^m)$ and $(N_t^m, \mathcal{F}_t, \mathbb{P}^m)$ are continuous martingales. The quadratic variation $([M^m, N^m]_t)$ exists $\mathbb{P}^m$-almost surely. To complete the identification of the limit, it suffices to note that the quadratic variation is given by

$$
\begin{align}
[M^m] &= 4 \int_0^t 1_{(0,\infty)}(\mathbf{r}_u) du \qquad \mathbb{P}^m\text{-almost surely,} \notag \\
[N^m] &= 4 \int_0^t 1_{(0,\infty)}(\mathbf{s}_u) du \qquad \mathbb{P}^m\text{-almost surely, and} \tag{69} \\
[M^m, N^m] &= 4 \int_0^t 1_{(0,\infty)}(\mathbf{r}_u) 1_{(0,\infty)}(\mathbf{s}_u) du \qquad \mathbb{P}^m\text{-almost surely,} \notag
\end{align}
$$

which holds following the computations in the proof of [23, Theorem 22]. We show that
$((M_t^m)^2 - 4 \int_0^t 1_{(0,\infty)} \mathbf{r}_u du)$ is a sub- and a supermartingale and hence a martingale using a
monotone class argument by noting first that for any bounded continuous and non-negative function $G: W \to \mathbb{R}_+$,

$$ \mathbb{E}^m[G(M_t^m)^2] = \lim_{n \to \infty} \mathbb{E}^{n,m}[G(M_t^{n,m})^2] \quad (70) $$

holds using uniform integrability of $((M_t^{n,m})^2, \mathbb{P}^{n,m})_{n \in \mathbb{N}, t \ge 0}$ which holds similarly as above. Note that

$$ \mathbb{E}^m \left[ G \int_s^t 1_{(0,\infty)}(\mathbf{r}_u) du \right] \le \lim_{\epsilon \downarrow 0} \liminf_{n \to \infty} \mathbb{E}^{n,m} \left[ G \int_s^t 1_{(\epsilon,\infty)}(\mathbf{r}_u) du \right] \quad (71) $$

holds by lower semicontinuity of $\omega \to \int_0^t 1_{(\epsilon,\infty)}(\omega_s)ds$ for each $\epsilon > 0$, Fatou's lemma and the Portmanteau theorem. For any fixed $\epsilon > 0$,

$$ \liminf_{n \to \infty} \mathbb{E}^{n,m} \left[ G \left( \int_s^t \theta^n(\mathbf{r}_u)^2 du - \int_s^t 1_{(\epsilon,\infty)}(\mathbf{r}_u) du \right) \right]. \quad (72) $$

Then by (70), (71) and (72)

$$ \begin{aligned} & \mathbb{E}^m \left[ G \left( (M_t^m)^2 - (M_s^m)^2 - 4 \int_s^t 1_{(0,\infty)}(\mathbf{r}_u) du \right) \right] \\ & \ge \lim_{\epsilon \downarrow 0} \liminf_{n \to \infty} \mathbb{E}^{n,m} \left[ G \left( (M_t^{n,m})^2 - (M_s^{n,m})^2 - 4 \int_s^t \theta^n(\mathbf{r}_u)^2 du \right) \right] = 0 \end{aligned} $$

and hence by a monotone class argument, cf. [50, Chapter 1, Theorem 8], $((M_t^m)^2 - 4 \int_0^t 1_{(0,\infty)}(\mathbf{r}_u) du, \mathbb{P}^m)$ is a submartingale. To show that it is also a supermartingale we note that $((M_t^m)^2 - 4t, \mathbb{P}^m)$ is a supermartingale by (70). By the uniqueness of the Doob-Meyer decomposition, cf. [50, Chapter 3, Theorem 8], $t \to [M^m]_t - 4t$ is $\mathbb{P}^m$-almost surely decreasing. Note further, that $(\mathbf{r}_t, \mathcal{F}_t, \mathbb{P}^m)$ is a continuous semimartingale with $[\mathbf{r}] = [M^m]$. Then by Ito-Tanaka formula, cf. [51, Chapter 6, Theorem 1.1],

$$ \int_0^t 1_{\{0\}}(\mathbf{r}_u)d[M^m]_u = \int_0^t 1_{\{0\}}(\mathbf{r}_u)d[\mathbf{r}]_u = \int_0^t 1_{\{0\}}(y)\ell_t^y(\mathbf{r})dy = 0, $$

where $\ell_t^y(\mathbf{r})$ is the local time of $\mathbf{r}$ in $y$. Therefore, for any $0 \le s < t$,

$$ [M^m]_t - [M^m]_s = \int_0^t 1_{(0,\infty)}(\mathbf{r}_u)d[M^m]_u \le 4 \int_0^t 1_{(0,\infty)}(\mathbf{r}_u)du $$

and hence, for any $\mathcal{F}_s$-measurable, bounded, non-negative function $G: W \to \mathbb{R}_+$,

$$ \mathbb{E}^m \left[ G \left( (M_t^m)^2 - (M_s^m)^2 - 4 \int_s^t 1_{(0,\infty)}(\mathbf{r}_u) du \right) \right] \le 0. $$

As before, by a monotone class argument, $((M_t^m)^2 - 4 \int_0^t 1_{(0,\infty)}(\mathbf{r}_u) du, \mathbb{P}^m)$ is a supermartingale, and hence a martingale.
Hence, we obtain the quadratic variation $[M^m]_t$ given in (69). The other characterizations in (69) follow by analogous arguments. Then by a martingale representation theorem, see [35, Chapter II, Theorem 7.1], we conclude, that there are a probability space $(\Omega^m, \mathcal{A}^m, P^m)$ and a Brownian motion motion $W$ and random variables $(r^m, s^m)$ on this space such that $P^m \circ (r^m, s^m)^{-1} = \mathbb{P}^m \circ (\mathbf{r}^m, \mathbf{s}^m)^{-1}$ and such that $(r^m, s^m, W)$ is a weak solution of (60). Finally, note that we have weak convergence of $Q \circ (r^{n,m}, s^{n,m})^{-1}$ to $P^m \circ (r^m, s^m)^{-1}$ not only along a subsequence since the characterization of the limit holds for any subsequence $(n_k)_{k \in \mathbb{N}}$.

**Comparison of two solutions:** To show $P^m[r_t^m \le s_t^m$ for all $t \ge 0] = 1$ we note that by Lemma 18, $Q[r_t^n \le s_t^n$ for all $t \ge 0] = 1$. The monotonicity carries over to the limit by the Portmanteau theorem for closed sets, since we have weak convergence of $\mathbb{P}^{n,m} \circ (\mathbf{r}, \mathbf{s})^{-1}$ to $\mathbb{P}^m \circ (\mathbf{r}, \mathbf{s})^{-1}$. $\square$

We show in the next step that the distribution of the solution of (60) converges as $m \to \infty$. For each $m \in \mathbb{N}$ let $(\Omega^m, \mathcal{A}^m, P^m)$ be a probability space and random variables $r^m, s^m: \Omega^m \to \mathbb{W}$ such that $(r_t^m, s_t^m)_{t \ge 0}$ is a solution of (60). Let $\mathbb{P}^m = P^m \circ (r^m, s^m)^{-1}$ denote the law on $\mathbb{W} \times \mathbb{W}$.

**Lemma 20.** Assume that $(\tilde{b}, g)$ and $(\hat{b}, h)$ satisfy **H** 1 and **H** 2. Let $\eta \in \Gamma(\mu, \nu)$ where the probability measures $\mu$ and $\nu$ on $\mathbb{R}_+$ satisfy **H** 3. Assume that $(g^m)_{m \in \mathbb{N}}$, $(h^m)_{m \in \mathbb{N}}$, $(\mu_m)_{m \in \mathbb{N}}$, $(\nu_m)_{m \in \mathbb{N}}$ and $(\eta_m)_{m \in \mathbb{N}}$ satisfy conditions **H** 5 and **H** 7. Then there exists a random variable $(r, s)$ defined on some probability space $(\Omega, \mathcal{A}, P)$ with values in $\mathbb{W} \times \mathbb{W}$, such that $(r_t, s_t)_{t \ge 0}$ is a weak solution of the sticky stochastic differential equation (25). Furthermore, the sequence of laws $P^m \circ (r^m, s^m)^{-1}$ converges weakly to the law $P \circ (r, s)^{-1}$. If additionally,

$$ \tilde{b}(r) \le \hat{b}(r), \quad g(r) \le h(r) \quad \text{and} \quad g^m(r) \le h^m(r) \quad \text{for any } r \in \mathbb{R}_+, \quad \text{and} \\ P^m[r_0^m \le s_0^m] = 1 \quad \text{for any } m \in \mathbb{N} $$

then $P[r_t \le s_t$ for all $t \ge 0] = 1$.

*Proof.* The proof is structured as the proof of Lemma 19. First analogously to the proof of (64) we show under **H** 1, **H** 5 and **H** 7,

$$ \sup_{t \in [0,T]} \mathbb{E}[|r_t^m|^p] < \infty. \qquad (73) $$

Tightness of the sequence of probability measures $(\mathbb{P}^m)_{m \in \mathbb{N}}$ on $(\mathbb{W} \times \mathbb{W}, \mathcal{B}(\mathbb{W}) \otimes \mathcal{B}(\mathbb{W}))$ holds adapting the steps of the proof of Lemma 19 to (60). Note that (64) and (65) hold analogously for $(r_t^m, s_t^m)_{m \in \mathbb{N}}$ by **H** 1, **H** 5 and **H** 7. Hence by Kolmogorov's continuity criterion, cf. [38, Corollary 14.9], we can deduce that there exists a probability measure $\mathbb{P}$ on $(\mathbb{W} \times \mathbb{W}, \mathcal{B}(\mathbb{W}) \otimes \mathcal{B}(\mathbb{W}))$ such that there is a subsequence $(m_k)_{k \in \mathbb{N}}$ along which $\mathbb{P}^{m_k}$ converge towards $\mathbb{P}$. To characterize the limit, we first note that by Skorokhod representation
theorem, cf. [9, Chapter 1, Theorem 6.7], without loss of generality we can assume that $(r^m, s^m)$ are defined on a common probability space $(\Omega, \mathcal{A}, P)$ with expectation $E$ and converge almost surely to $(r, s)$ with distribution $P$. By $\mathbf{H}5$, $P_t^m(g^m) = E[g^m(r_t^m)]$ and the monotone convergence theorem, $P_t^m(g^m)$ converges to $P_t(g)$ for any $t \ge 0$. Then, by Lebesgue convergence theorem it holds almost surely for all $t \ge 0$

$$ \lim_{m \to \infty} \int_0^t (\tilde{b}(r_t^m) + P_u^m(g^m)) du = \int_0^t (\tilde{b}(r_t) + P_u(g)) du, \quad (74) $$

where $P_u^m = P \circ (r_u^m)^{-1}$ and $P_u = P \circ (r_u)^{-1}$. A similar statement holds for $(s_t)_{t \ge 0}$.

Consider the mappings $M^m, N^m : \mathbb{W} \times \mathbb{W} \to \mathbb{W}$ defined by (67). Then for all $m \in \mathbb{N}$, $(M_t^m, \mathcal{F}_t, \mathbb{P}^m)$ and $(N_t^m, \mathcal{F}_t, \mathbb{P}^m)$ are martingales with respect to the canonical filtration $\mathcal{F}_t = \sigma((\mathbf{r}_u, \mathbf{s}_u)_{0 \le u \le t})$. Further the family $(M_t^m, \mathbb{P}^m)_{m \in \mathbb{N}, t \ge 0}$ and $(N_t^m, \mathbb{P}^m)_{m \in \mathbb{N}, t \ge 0}$ are uniformly integrable by (73). In the same line as in the proof of Lemma 19 and by (74), $\mathbb{P}^m \circ (\mathbf{r}, \mathbf{s}, M^m, N^m)$ converges weakly to $\mathbb{P} \circ (\mathbf{r}, \mathbf{s}, M, N)$ where

$$ M_t = \mathbf{r}_t - \mathbf{r}_0 - \int_0^t (\tilde{b}(\mathbf{r}_u) + P_u(g)) du \quad \text{and} \quad N_t = \mathbf{s}_t - \mathbf{s}_0 - \int_0^t (\hat{b}(\mathbf{s}_u) + \hat{P}_u(h)) du. $$

Let $G : \mathbb{W} \to \mathbb{R}_+$ be a $\mathcal{F}_s$-measurable bounded, non-negative function. By uniform integrability, for any $s \le t$,

$$ \begin{aligned} \mathbb{E}[G(M_t - M_s)] &= \mathbb{E}[G(\int_s^t (\tilde{b}(\mathbf{r}_u) + P_u(g)) du)] = \lim_{m \to \infty} \mathbb{E}^m[G(\int_s^t (\tilde{b}(\mathbf{r}_u) + P_u(g^m)) du)] \\ &= \lim_{m \to \infty} \mathbb{E}^m[G(M_t^m - M_s^m)] = 0, \end{aligned} $$

and analogously for $(N_t)_{t \ge 0}$. Hence, $(M_t, \mathcal{F}_t, \mathbb{P})$ and $(N_t, \mathcal{F}_t, \mathbb{P})$ are martingales. Further, the quadratic variation $([M, N])_t$ exists $\mathbb{P}$-almost surely and is given by (69) $\mathbb{P}$-almost surely, which holds following the computations in the proof of Lemma 19. As in Lemma 19, we conclude by a martingale representation theorem that there are a probability space $(\Omega, \mathcal{A}, P)$ and a Brownian motion $W$ and random variables $(r, s)$ on this space such that $P \circ (r, s)^{-1} = \mathbb{P} \circ (\mathbf{r}, \mathbf{s})^{-1}$ and such that $(r, s, W)$ is a weak solution of (25). Note that the limit identification holds for all subsequences $(m_k)_{k \in \mathbb{N}}$ and hence $P^m \circ (r^m, s^m)^{-1}$ converges weakly to $P \circ (r, s)^{-1}$ for $m \to \infty$. The monotonicity $P^m[r_t^m \le s_t^m$ for all $t \ge 0] = 1$ carries over to the limit by Portmanteau theorem, since $\mathbb{P}^m \circ (\mathbf{r}, \mathbf{s})^{-1}$ converges weakly to $\mathbb{P} \circ (\mathbf{r}, \mathbf{s})^{-1}$. $\square$

*Proof of Theorem 3.* The proof is a direct consequence of Lemma 19 and Lemma 20. $\square$

### 6.3.2 Proof of Theorem 5

*Proof of Theorem 5.* Note that the Dirac at 0, $\delta_0$, is by definition an invariant measure of $(r_t)_{t \ge 0}$ solving (6). Assume that the process starts from an invariant probability measure
$\pi$, hence $\mathbb{P}(r_t > 0) = p = \pi((0, \infty))$ for any $t \ge 0$. Note that for $p = 0$ the drift vanishes. If the initial measure is the Dirac measure in 0, $\delta_0$, then the diffusion coefficient disappears. Hence, $\text{Law}(r_t) = \delta_0$ for any $t \ge 0$. It remains to investigate the case $p \ne 0$. Here, we are in the regime of [23, Lemma 24] where an invariant measure is of the form (28). Since $p = \mathbb{P}(r_t > 0)$, the invariant measure $\pi$ satisfies additionally the necessary condition

$$p = \pi((0, \infty)) = \frac{I(a, p)}{2/(ap) + I(a, p)} \qquad (75)$$

with $I(a, p)$ given in (27). For $p \ne 0$, this expression is equivalent to (26). $\square$

**Proof of Proposition 6.** By Theorem 5, it suffices to study the solutions of (26). By (27) and since $\tilde{b}(r) = -\tilde{L}r$, it holds for $\hat{I}(a, p) = (1-p)I(a, p)$,

$$\hat{I}(a, p) = \left(\sqrt{\frac{\pi}{2}} + \int_{0}^{\frac{ap}{\sqrt{2L}}} \exp(-x^2/2)dx\right) \sqrt{\frac{2}{L}} \exp\left(\frac{a^2p^2}{4L}\right)(1-p). \qquad (76)$$

In the case $a/\sqrt{L} \le 2/\sqrt{\pi}$, $\hat{I}(a, 0) = \sqrt{\pi/\tilde{L}}$ by (76). Further, by $1+x \le e^x$ and $a/\sqrt{L} \le 2/\sqrt{\pi}$,

$$\left(\sqrt{\frac{\pi}{2}} + \int_{0}^{\frac{ap}{\sqrt{2L}}} e^{-\frac{x^2}{2}} dx\right) (1-p) e^{\frac{a^2 p^2}{4L}} \le \sqrt{\frac{\pi}{2}} \left(1 + \sqrt{\frac{2}{\pi}} \int_{0}^{\frac{ap}{\sqrt{2L}}} e^{-\frac{x^2}{2}} dx\right) e^{-p} e^{\frac{p^2}{\pi}} \\ \le \sqrt{\frac{\pi}{2}} \left(1 + \frac{2p}{\pi}\right) e^{-p} e^{\frac{p^2}{\pi}} \le \sqrt{\frac{\pi}{2}} e^{p\left(\frac{3}{\pi}-1\right)} < \sqrt{\frac{\pi}{2}}$$

for $p \in (0, 1]$. Hence, $\hat{I}(a, p) < \hat{I}(a, 0)$ by (76). Therefore, $\hat{I}(a, p) < \hat{I}(a, 0) \le \frac{2}{a}$ for all $p \in (0, 1]$ and so $\delta_0$ is the unique invariant probability measure for $a/\sqrt{L} \le 2/\sqrt{\pi}$.

To show that for $a/\sqrt{L} > 2/\sqrt{\pi}$, there exists a unique $p$ solving (26), we note that $\hat{I}(a, p)$ is continuous with $\hat{I}(a, 0) > 2/a$ and $\hat{I}(a, 1) = 0$. By the mean value theorem, there exists at least one $p \in (0, 1)$ satisfying (26). In the following we drop the dependence on $a$ in $I(a, p)$ and $\hat{I}(a, p)$. We show uniqueness of the solution $p$ by contradiction. Assume that $p_1 < p_2$ are the two smallest solutions of (26). Hence, it holds either $\hat{I}'(p_1) < 0$ or $\hat{I}'(p) = 0$ for $p_1$. Note that the derivative is given by

$$\begin{aligned} \hat{I}'(p_i) &= -I(p_i) + (1-p_i)I'(p_i) = -I(p_i) + (1-p_i)\left(p_i \frac{a^2}{2\tilde{L}} I(p_i) + \frac{a}{\tilde{L}}\right) \\ &= -\frac{2}{a(1-p_i)} + (1-p_i)\frac{a}{\tilde{L}}\left(\frac{p_i}{1-p_i} + 1\right) = -\frac{2}{a(1-p_i)} + \frac{a}{\tilde{L}}. \end{aligned} \qquad (77)$$

Then, for $p_2 > p_1$, it holds

$$\hat{I}'(p_2) = -\frac{2}{a(1-p_2)} + \frac{a}{\tilde{L}} < -\frac{2}{a(1-p_1)} + \frac{a}{\tilde{L}} = \hat{I}'(p_1) \le 0.$$
If $\hat{I}'(p_1) < 0$, it holds $\hat{I}'(p_2) < 0$ which contradicts that $p_1$ and $p_2$ are the two smallest solutions. In the second case, when $\hat{I}'(p_1) = 0$, we note that the second derivative of $\hat{I}(p)$ at $p_1$ is given by

$$ 
\begin{aligned}
\hat{I}''(p_1) &= -2I'(p_1) + (1-p_1)I''(p_1) \\
&= \left(-2 + (1-p_1)\frac{a^2p_1}{2\tilde{L}}\right)\left(I(p_1)\frac{a^2p_1}{2\tilde{L}} + \frac{a}{\tilde{L}}\right) + (1-p_1)I(p_1)\frac{a^2}{2\tilde{L}} \\
&= \left(-2 + (1-p_1)\frac{a^2p_1}{2\tilde{L}}\right)\frac{a}{\tilde{L}(1-p_1)} + \frac{a}{\tilde{L}} = -\frac{a}{\tilde{L}(1-p_1)} < 0.
\end{aligned}
$$

Hence, in this case there is a maximum at $p_1$, which contradicts that $p_1$ is the smallest solution. Thus, there exists a unique solution $p_1$ of (26) for $a/\sqrt{\tilde{L}} > 2/\sqrt{\pi}$. $\square$

### 6.3.3 Proof of Theorem 7

*Proof of Theorem 7.* To show (31) we extend the function $f$ to a concave function on $\mathbb{R}$ by setting $f(x) = x$ for $x < 0$. Note that $f$ is continuously differentiable and $f'$ is absolutely continuous and bounded. Using Ito-Tanaka formula, c.f. [51, Chapter 6, Theorem 1.1] we obtain

$$ df(r_t) = f'(r_t)(\tilde{b}(r_t) + a\mathbb{P}(r_t > 0))dt + 2f''(r_t)1_{(0,\infty)}(r_t)dt + dM_t, $$

where $M_t = 2 \int_0^t f'(r_s)1_{(0,\infty)}(r_s)dB_s$ is a martingale. Taking expectation, we get

$$ 
\begin{aligned}
\frac{d}{dt}\mathbb{E}[f(r_t)] &= \mathbb{E}[f'(r_t)(\tilde{b}(r_t) + a\mathbb{P}(r_t > 0))] + 2\mathbb{E}[f''(r_t)1_{(0,\infty)}(r_t)] \\
&= \mathbb{E}[f'(r_t)\tilde{b}(r_t) + 2(f''(r_t) - f''(0))] + \mathbb{E}[af'(r_t) + 2f''(0)]\mathbb{P}(r_t > 0) \\
&\le -c\mathbb{E}[f(r_t)],
\end{aligned}
$$

where the last step holds by (39) and (40). By applying Gronwall's lemma, we obtain (31). $\square$

## 6.4 Proof of Section 5

As for the nonlinear case we show Theorem 10 via a family of stochastic differential equations, with Lipschitz continuous coefficients,

$$ 
\begin{aligned}
dr_t^{i,n,m} &= (\tilde{b}(r_t^{i,n,m}) + \frac{1}{N}\sum_{j=1}^N g^m(r_t^{j,n,m}))dt + 2\theta^n(r_t^{i,n,m})dW_t^i \\
ds_t^{i,n,m} &= (\hat{b}(s_t^{i,n,m}) + \frac{1}{N}\sum_{j=1}^N h^m(s_t^{j,n,m}))dt + 2\theta^n(s_t^{i,n,m})dW_t^i
\end{aligned}
\qquad (78) $$

$$ \text{Law}(r_0^{i,n,m}, s_0^{i,n,m}) = \eta_{n,m}, \quad i \in \{1, \dots, N\}, $$
where $\eta_{n,m} \in \Gamma(\mu_{n,m}, \nu_{n,m})$. Under **H1**, **H2**, **H5**, **H6** and **H7** we identify the weak limit of $(\{r_t^{i,n,m}, s_t^{i,n,m}\}_{i=1,n,m \in \mathbb{N}}^N)_{t \ge 0}$ solving (78) for $n \to \infty$ by $(\{r_t^{i,m}, s_t^{i,m}\}_{i=1,m \in \mathbb{N}}^N)_{t \ge 0}$ solving the family of SDEs given by

$$
\begin{align}
\mathrm{d}r_t^{i,m} &= (\tilde{b}(r_t^{i,m}) + \frac{1}{N} \sum_{j=1}^N g^m(r_t^{j,m})) \mathrm{d}t + 2\mathbb{1}_{(0,\infty)}(r_t^{i,m}) \mathrm{d}W_t^i, \nonumber \\
\mathrm{d}s_t^{i,m} &= (\hat{b}(s_t^{i,m}) + \frac{1}{N} \sum_{j=1}^N h^m(s_t^{j,m})) \mathrm{d}t + 2\mathbb{1}_{(0,\infty)}(s_t^{i,m}) \mathrm{d}W_t^i, \tag{79} \\
\text{Law}(r_0^{i,m}, s_0^{i,m}) &= \eta_m, \quad i \in \{1, \dots, N\}, \nonumber
\end{align}
$$

where $\eta_m \in \Gamma(\mu_m, \nu_m)$.

Taking the limit $m \to \infty$, we obtain (35) as the weak limit of (79). In the case $g(r) = \mathbb{1}_{(0,\infty)}(r)$, we can choose $g^m = \theta^m$.

Consider a probability space $(\Omega_0, \mathcal{A}_0, \mathcal{Q})$ and $N$ i.i.d.1-dimensional Brownian motions $(\{W_t^i\}_{i=1}^N)_{t \ge 0}$. Note that under **H1–H7**, there are random variables $\{r_t^{i,n,m}\}_{i=1}^N$, $\{s_t^{i,n,m}\}_{i=1}^N$: $\Omega_0 \to \mathbb{W}^N$ for each $n, m$ such that $(\{r_t^{i,n,m}, s_t^{i,n,m}\}_{i=1}^N)$ is a unique solution to (78) by [45, Theorem 2.2]. We denote by $\mathbb{P}^{n,m} = \mathcal{Q} \circ (\{r_t^{i,n,m}, s_t^{i,n,m}\}_{i=1}^N)^{-1}$ the law on $\mathbb{W}^N \times \mathbb{W}^N$.

Before taking the two limits and proving Theorem 10, we introduce a modification of Ikeda and Watanabe's comparison theorem, to compare two solutions of (78), cf. [35, Section VI, Theorem 1.1].

**Lemma 21.** Suppose a solution $(\{r_t^{i,n,m}, s_t^{i,n,m}\}_{i=1}^N)_{t \ge 0}$ of (78) is given for fixed $n, m \in \mathbb{N}$. Assume **H5** for $g^m$ and $h^m$, **H1** for $\tilde{b}$ and $\hat{b}$, **H6** for $\theta^n$. If $Q[r_0^{i,n,m} \le s_0^{i,n,m}$ for all $i = 1, \dots, N] = 1$, $\tilde{b}(r) \le \hat{b}(r)$ and $g^m(r) \le h^m(r)$ for any $r \in \mathbb{R}_+$, then

$$Q[r_t^{i,n,m} \le s_t^{i,n,m} \text{ for all } t \ge 0 \text{ and } i = 1, \dots, N] = 1$$

*Proof.* The proof is similar for each component $i = 1, \dots, N$ to the proof of Lemma 18. It holds for the interaction part similarly to (62) using the properties of $g^m$ and $h^m$,

$$\frac{1}{N} \sum_{j=1}^{N} (g^m(r_t^{j,n,m}) - h^m(s_t^{j,n,m})) \le K_m \frac{1}{N} \sum_{j=1}^{N} |r_t^{j,n,m} - s_t^{j,n,m}| \mathbb{1}_{(0,\infty)}(r_t^{j,n,m} - s_t^{j,n,m}).$$

Hence, we obtain analogously to (63),

$$\mathbb{E}[(r_t^{i,n,m} - s_t^{i,n,m})_+] \le \tilde{L} \mathbb{E} \left[ \int_0^t (r_u^{i,n,m} - s_u^{i,n,m})_+ du \right] + K_m \mathbb{E} \left[ \int_0^t \frac{1}{N} \sum_{j=1}^N (r_u^{j,n,m} - s_u^{j,n,m})_+ du \right]$$

for all $i = 1, \dots, N$. Assume analogously to the nonlinear case $t^* = \inf\{t \ge 0: \mathbb{E}[(r_t^{i,n,m} - s_t^{i,n,m})_+] > 0$ for some $i\} < \infty$. Then, there exists $i \in \{1, \dots, N\}$ such that $\int_0^{t^*} \mathbb{E}[(r_u^{i,n,m} - s_u^{i,n,m})_+] du > 0$. But, by definition of $t^*$, for all $i, u < t^*$, $\mathbb{E}[(r_u^{i,n,m} - s_u^{i,n,m})_+] = 0$. This contradicts the definition of $t^*$. Hence, $Q[r_t^{i,n,m} \le s_t^{i,n,m}$ for all $i, t \ge 0] = 1$. $\square$
In the next step, we prove that the distribution of the solution of (78) converges as $n \to \infty$.

**Lemma 22.** Assume that $H1$ and $H2$ is satisfied for $(\tilde{b}, g)$ and $(\hat{b}, h)$. Further, let $(\theta^n)_{n \in \mathbb{N}}$, $(g^m)_{m \in \mathbb{N}}$, $(h^m)_{m \in \mathbb{N}}$, $(\mu_{n,m})_{n,m \in \mathbb{N}}$, $(\nu_{n,m})_{n,m \in \mathbb{N}}$ and $(\eta_{n,m})_{n,m \in \mathbb{N}}$ be such that $H5$, $H6$ and $H7$ hold. Let $m \in \mathbb{N}$. Then there exists a random variable $(\{r^{i,m}, s^{i,m}\}_{i=1}^N)$ defined on some probability space $(\Omega^m, \mathcal{A}^m, P^m)$ with values in $\mathbb{W}^N \times \mathbb{W}^N$ such that $(\{r_t^{i,m}, s_t^{i,m}\}_{i=1}^N)_{t \ge 0}$ is a weak solution of (79). Moreover, the laws $Q \circ (\{r^{i,n,m}, s^{i,n,m}\}_{i=1}^N)^{-1}$ converge weakly to $P^m \circ (\{r^{i,m}, s^{i,m}\}_{i=1}^N)^{-1}$. If in addition,

$$ \begin{array}{ll} \tilde{b}(r) \le \hat{b}(r) \quad \text{and} \quad g^m(r) \le h^m(r) & \text{for any } r \in \mathbb{R}_+, \\ Q[r_0^{i,n,m} \le s_0^{i,n,m}] = 1 & \text{for any } n \in \mathbb{N}, i = 1, \dots, N, \end{array} $$

then $P^m[r_t^{i,m} \le s_t^{i,m} \text{ for all } t \ge 0 \text{ and } i \in \{1, \dots, N\}] = 1$.

*Proof.* Fix $m \in \mathbb{N}$. The proof is divided in three parts and is similar to the proof of Lemma 19. First we show tightness of the sequences of probability measures. Then we identify the limit of the sequence of stochastic processes. Finally, we compare the two limiting processes.

**Tightness:** We show analogously as in the proof of Lemma 19 that the sequence of probability measures $(\mathbb{P}^{n,m})_{n \in \mathbb{N}}$ on $(\mathbb{W}^N \times \mathbb{W}^N, \mathcal{B}(\mathbb{W}^N) \otimes \mathcal{B}(\mathbb{W}^N))$ is tight by applying Kolmogorov's continuity theorem. We consider $p > 2$ such that the $p$-th moment in $H7$ are uniformly bounded. Fix $T > 0$. Then the $p$-th moment of $r_t^{i,n,m}$ and $s_t^{i,n,m}$ for $t < T$ is bounded using Ito's formula,

$$ \begin{aligned} \mathrm{d}|r_t^{i,n,m}|^p & \le p|r_t^{i,n,m}|^{p-2}\langle r_t^{i,n,m}, (\tilde{b}(r_t^{i,n,m}) + \frac{1}{N}\sum_{j=1}^N g^m(r_t^{j,n,m}))\rangle \mathrm{d}t \\ & \quad + 2\theta^n(r_t^{i,n,m})p|r_t^{n,m}|^{p-2}r_t^{i,n,m}\mathrm{d}W_t^i + p(p-1)|r_t^{i,n,m}|^{p-2}2\theta^n(r_t^{i,n,m})^2\mathrm{d}t \\ & \le p(|r_t^{i,n,m}|^p \tilde{L} + \Gamma|r_t^{i,n,m}|^{p-1} + 2(p-1)|r_t^{i,n,m}|^{p-2})\mathrm{d}t + 2\theta^n(r_t^{i,n,m})p(r_t^{i,n,m})^{p-1}\mathrm{d}W_t^i \\ & \le p(\tilde{L} + \Gamma + 2(p-1))|r_t^{i,n,m}|^p \mathrm{d}t + p(\Gamma + 2(p-1))\mathrm{d}t + 2\theta^n(r_t^{i,n,m})p(r_t^{n,m})^{p-1}\mathrm{d}W_t^i. \end{aligned} $$

Taking expectation yields

$$ \frac{\mathrm{d}}{\mathrm{d}t} \mathbb{E}[|r_t^{i,n,m}|^p] \le p(L + \Gamma + 2(p-1)) \mathbb{E}|r_t^{i,n,m}|^p + p(\Gamma + 2(p-1)). $$

Then by Gronwall's lemma

$$ \sup_{t \in [0,T]} \mathbb{E}[|r_t^{i,n,m}|^p] \le e^{p(L+\Gamma+2(p-1))T} (\mathbb{E}[|r_0^{i,n,m}|^p] + Tp(\Gamma+2(p-1))) < C_p < \infty, \quad (80) $$
where $C_p$ depends on $T$ and the $p$-th moment of the initial distribution, which is by assumption finite. Similarly, it holds $\sup_{t \in [0,T]} \mathbb{E}[|s_t^{i,n,m}|^p] < C_p$ for $t \le T$. Using these moment bounds, it holds for all $t_1, t_2 \in [0,T]$ by $\mathbf{H}1, \mathbf{H}5$ and $\mathbf{H}6$,

$$
\begin{align*}
& \mathbb{E}[|r_{t_2}^{i,n,m} - r_{t_1}^{i,n,m}|^p] \\
& \le C_1(p) \left( \mathbb{E}[\left| \int_{t_1}^{t_2} \tilde{b}(r_u^{i,n,m}) + \frac{1}{N} \sum_{j=1}^N g^m(r_t^{j,n,m}) du \right|^p] + \mathbb{E}[\left| \int_{t_1}^{t_2} 2\theta^n(r_u^{i,n,m}) dW_u^i \right|^p] \right) \\
& \le C_2(p) \left( \left( \mathbb{E}\left[\frac{\tilde{L}^p}{|t_2-t_1|} \int_{t_1}^{t_2} |r_u^{i,n,m}|^p du \right] + \Gamma^p \right) |t_2-t_1|^p + \mathbb{E}[\left| \int_{t_1}^{t_2} 2\theta^n(r_u^{i,n,m}) du \right|^{p/2}] \right) \\
& \le C_2(p) \left( \left( \frac{\tilde{L}^p}{|t_2-t_1|} \int_{t_1}^{t_2} \mathbb{E}[|r_u^{i,n,m}|^p] du + \Gamma^p \right) |t_2-t_1|^p + 2^{p/2} |t_2-t_1|^{p/2} \right) \\
& \le C_3(p, T, \tilde{L}, \Gamma, C_p) |t_2-t_1|^{p/2},
\end{align*}
$$

where $C_k(\cdot)$ are constants depending on the stated arguments, but independent of $n, m$. Note that in the second step, we use Burkholder-Davis-Gundy inequality, see [50, Chapter IV, Theorem 48]. It holds similarly, $\mathbb{E}[|s_{t_2}^{i,n,m} - s_{t_1}^{i,n,m}|^p] \le C_3(p, T, \tilde{L}, \Gamma, C_p)|t_2 - t_1|^{p/2}$. Hence,

$$
\begin{align*}
& \mathbb{E}[|\{(r_{t_2}^{i,n,m}, s_{t_2}^{i,n,m})\}_{i=1}^N - \{(r_{t_1}^{i,n,m}, s_{t_1}^{i,n,m})\}_{i=1}^N|^p] \\
& \le C_4(p, N) \left( \sum_{i=1}^N \left( \mathbb{E}[|r_{t_2}^{i,n,m} - r_{t_1}^{i,n,m}|^p] + \mathbb{E}[|s_{t_2}^{i,n,m} - s_{t_1}^{i,n,m}|^p] \right) \right) \\
& \le C_5(p, N, T, \tilde{L}, \Gamma, C_p) |t_2 - t_1|^{p/2}
\end{align*}
$$

for all $t_1, t_2 \in [0, T]$. Hence, by Kolmogorov's continuity criterion, cf. [38, Corollary 14.9], there exists a constant $\tilde{C}$ depending on $p$ and $\gamma$ such that

$$
\mathbb{E}\left[\left[\left(\{r^{i,n,m}, s^{i,n,m}\}_{i=1}^N\right)\right]^p\right] \le \tilde{C} \cdot C_5(p, N, T, \tilde{L}, \Gamma, C_p). \quad (81)
$$

where $[\cdot]^p_\gamma$ is defined by $[x]_\gamma = \sup_{t_1,t_2 \in [0,T]} \frac{|x(t_1)-x(t_2)|}{|t_1-t_2|^\gamma}$ and $(\{r_t^{i,n,m}, s_t^{i,n,m}\}_{i=1}^N)_{n \in \mathbb{N}, t \ge 0}$ is tight in $C([0,T], \mathbb{R}^{2N})$. Hence, for each $T > 0$ there exists a subsequence $n_k \to \infty$ and a probability measure $P_T$ on $C([0,T], \mathbb{R}^{2N})$. Since $\{\mathbb{P}_T^m\}_T$ is a consistent family, there exists by [38, Theorem 5.16] a probability measure $\mathbb{P}^m$ on $(\mathbb{W}^N \times \mathbb{W}^N, \mathcal{B}(\mathbb{W}^N) \otimes \mathcal{B}(\mathbb{W}^N))$ such that $\mathbb{P}^{n_k,m}$ converges weakly to $\mathbb{P}^m$. Note that we can take here the same subsequence $(n_k)$ for all $m$ using a diagonalization argument.

**Characterization of the limit measure:** Denote by $(\{\mathbf{r}_t^i, \mathbf{s}_t^i\}_{i=1}^N) = \omega(t)$ the canonical process on $\mathbb{W}^N \times \mathbb{W}^N$. To characterize the measure $\mathbb{P}^m$ we first note that $\mathbb{P}^m \circ (\mathbf{r}_0^i, \mathbf{s}_0^i)^{-1} = \eta_m$ for all $i \in \{1, \dots, N\}$, since $\mathbb{P}^{n,m}(\mathbf{r}_0^i, \mathbf{s}_0^i)^{-1} = \eta_{n,m}$ converges weakly to $\eta_m$ by assump-
tion. We define maps $M^{i,m}, N^{i,m} : \mathbb{W}^N \times \mathbb{W}^N \to \mathbb{W}$ by

$$ 
\begin{align} 
M_t^{i,m} &= \mathbf{r}_t^i - \mathbf{r}_0^i - \int_0^t \left( \tilde{b}(\mathbf{r}_u^i) + \frac{1}{N} \sum_{j=1}^N g^m(\mathbf{r}_u^j) \right) du, \quad \text{and} \tag{82} \\ 
N_t^{i,m} &= \mathbf{s}_t^i - \mathbf{s}_0^i - \int_0^t \left( \hat{b}(\mathbf{s}_u^i) + \frac{1}{N} \sum_{j=1}^N h^m(\mathbf{s}_u^j) \right) du. \nonumber 
\end{align} 
$$

For each $n, m \in \mathbb{N}$ and $i = 1, \dots, N$, $(M_t^{i,m}, \mathcal{F}_t, \mathbb{P}^{n,m})$ is a martingale with respect to the filtration $\mathcal{F}_t = \sigma((\mathbf{r}_u^j, \mathbf{s}_u^j) : j = 1, \dots, N, 0 \le u \le t)$. The family $(\{M_t^{i,m}\}_{i=1}^N, \mathbb{P}^{n,m})_{n \in \mathbb{N}, t \ge 0}$ and $(\{N_t^{i,m}\}_{i=1}^N, \mathbb{P}^{n,m})_{n \in \mathbb{N}, t \ge 0}$ are uniformly integrable. Since the mappings $M^{i,m}$ and $N^{i,m}$ are continuous in $\mathbb{W}, \mathbb{P}^{n,m} \circ (\{\mathbf{r}^i, \mathbf{s}^i, M^{i,m}, N^{i,m}\}_{i=1}^N)^{-1}$ converges weakly to $\mathbb{P}^m \circ (\{\mathbf{r}^i, \mathbf{s}^i, M^{i,m}, N^{i,m}\}_{i=1}^N)^{-1}$ by the continuous mapping theorem. Then applying the same argument as in (68), $(M_t^{m,i}, \mathcal{F}_t, \mathbb{P}^m)$ and $(N_t^{m,i}, \mathcal{F}_t, \mathbb{P}^m)$ are continuous martingales for all $i = 1, \dots, N$ and the quadratic variation $([\{M^{i,m}, N^{i,m}\}_{i=1}^N]_t)_{t \ge 0}$ exists $\mathbb{P}^m$-almost surely. To complete the identification of the limit, it suffices to identify the quadratic variation. Similar to the computations in the proof of Lemma 19, it holds

$$ 
\begin{align} 
[M^{i,m}] &= 4 \int_0^1 \mathbb{1}_{(0,\infty)}(\mathbf{r}_u^i) du \qquad \mathbb{P}^m\text{-almost surely,} \notag \\ 
[N^{i,m}] &= 4 \int_0^1 \mathbb{1}_{(0,\infty)}(\mathbf{s}_u^i) du \qquad \mathbb{P}^m\text{-almost surely, and} \tag{83} \\ 
[M^{i,m}, N^{i,m}] &= 4 \int_0^1 \mathbb{1}_{(0,\infty)}(\mathbf{r}_u^i) \mathbb{1}_{(0,\infty)}(\mathbf{s}_u^i) du \qquad \mathbb{P}^m\text{-almost surely,} \notag 
\end{align} 
$$

Further, $[M^{i,m}, M^{j,m}]_t = [N^{i,m}, N^{j,m}]_t = [M^{i,m}, N^{j,m}]_t = 0$ $\mathbb{P}^{n,m}$-almost surely for $i \neq j$ and $(M_t^{i,m}, M_t^{j,m}, \mathbb{P}^{n,m}), (N_t^{i,m}, N_t^{j,m}, \mathbb{P}^{n,m})$ and $(M_t^{i,m}, N_t^{j,m}, \mathbb{P}^{n,m})$ are martingales. For any bounded, continuous non-negative function $G : \mathbb{W} \to \mathbb{R}$, it holds

$$ E^m[G(M_t^{i,m} M_t^{j,m} - M_s^{i,m} M_s^{j,m})] = \lim_{n \to \infty} E^{n,m}[G(M_t^{i,m} M_t^{j,m} - M_s^{i,m} M_s^{j,m})] = 0, $$

respectively, $E^m[G(N_t^{i,m} N_t^{j,m} - N_s^{i,m} N_s^{j,m})] = 0$ and $E^m[G(M_t^{i,m} N_t^{j,m} - M_s^{i,m} N_s^{j,m})] = 0$. Then

$$ [M^{i,m}, M^{j,m}] = [N^{i,m}, N^{j,m}] = [M^{i,m}, N^{j,m}] = 0 \quad \mathbb{P}^m\text{-almost surely, for all } i \neq j. \quad (84) $$

Then by a martingale representation theorem, cf. [35, Chapter II, Theorem 7.1], there is a probability space $(\Omega^m, \mathcal{A}^m, \mathbb{P}^m)$ and a Brownian motion $\{W^i\}_{i=1}^N$ and random variables $(\{r^{i,m}, s^{i,m}\}_{i=1}^N)$ on this space, such that $\mathbb{P}^m \circ (\{r^{i,m}, s^{i,m}\}_{i=1}^N)^{-1} = \mathbb{P}^m \circ (\{\mathbf{r}^i, \mathbf{s}^i\}_{i=1}^N)^{-1}$ and such that $(\{r^{i,m}, s^{i,m}, W^i\}_{i=1}^N)$ is a weak solution of (79).

**Comparison of two solutions:** To show $\mathbb{P}^m[r_t^{i,m} \le s_t^{i,m}$ for all $t \ge 0$ and $i = 1, \dots, N] = 1$ it suffices to note that $\mathbb{P}^{n,m}[r_t^{i,n,m} \le s_t^{i,n,m}$ for all $t \ge 0$ and $i = 1, \dots, N] = 1$, which holds by Lemma 21, carries over to the limit by the Portmanteau theorem, since we have weak convergence of $\mathbb{P}^{n,m} \circ (\{\mathbf{r}^i, \mathbf{s}^i\}_{i=1}^N)^{-1}$ to $\mathbb{P}^m \circ (\{\mathbf{r}^i, \mathbf{s}^i\}_{i=1}^N)^{-1}$. $\square$
We show in the next step that the distribution of the solution of (79) converges as $m \to \infty$. Consider a probability space $(\Omega^m, \mathcal{A}^m, P^m)$ for each $m \in \mathbb{N}$ and random variables $\{r^{i,m}\}_{i=1}^N, \{s^{i,m}\}_{i=1}^N : \Omega^m \to \mathbb{W}^N$ such that $(\{r_t^{i,m}, s_t^{i,m}\}_{i=1}^N)_{t \ge 0}$ is a solution to (79). Denote by $P^m = P^m \circ (\{r^{i,m}, s^{i,m}\}_{i=1}^N)^{-1}$ the law on $\mathbb{W}^N \times \mathbb{W}^N$.

**Lemma 23.** Assume that $\mathbf{H}1$ and $\mathbf{H}2$ is satisfied for $(\tilde{b}, g)$ and $(\hat{b}, h)$. Let $\eta \in \Gamma(\mu, \nu)$ where the probability measures $\mu$ and $\nu$ on $\mathbb{R}_+$ satisfy $\mathbf{H}3$. Further, let $(g^m)_{m \in \mathbb{N}}$, $(h^m)_{m \in \mathbb{N}}$, $(\mu_m)_{m \in \mathbb{N}}$, $(\nu_m)_{m \in \mathbb{N}}$ and $(\eta_m)_{m \in \mathbb{N}}$ be such that $\mathbf{H}5$ and $\mathbf{H}7$ hold. Then there exists a random variable $(\{r^i, s^i\}_{i=1}^N)$ defined on some probability space $(\Omega, \mathcal{A}, P)$ with values in $\mathbb{W}^N \times \mathbb{W}^N$ such that $(\{r_t^i, s_t^i\}_{i=1}^N)$ is a weak solution of (35). Moreover, the laws $P^m \circ (\{r^{i,m}, s^{i,m}\}_{i=1}^N)^{-1}$ converge weakly to $P \circ (\{r^i, s^i\}_{i=1}^N)^{-1}$. If in addition,

$$ \tilde{b}(r) \le \hat{b}(r), \quad g(r) \le h(r), \quad \text{and} \quad g^m(r) \le h^m(r) \quad \text{for any } r \in \mathbb{R}_+, \text{ and} $$

$$ P^m[r_0^{i,m} \le s_0^{i,m} \text{ for all } t \ge 0 \text{ and } i \in \{1, \dots, N\}] = 1 \quad \text{for any } m \in \mathbb{N}, $$

$$ \text{then } P[r_t^i \le s_t^i \text{ for all } t \ge 0 \text{ and } i \in \{1, \dots, N\}] = 1. $$

*Proof.* The proof is structured as the proof of Lemma 22. Tightness of the sequence of probability measures $(P^m)_{m \in \mathbb{N}}$ on $(\mathbb{W}^N \times \mathbb{W}^N, \mathcal{B}(\mathbb{W}^N) \otimes \mathcal{B}(\mathbb{W}^N))$ holds adapting the steps of the proof of Lemma 22 to (79). Note that (80) and (81) hold analogously for $(\{r_t^{i,m}, s_t^{i,m}\}_{i=1}^N)$ by $\mathbf{H}1$, $\mathbf{H}5$ and $\mathbf{H}7$. Hence by Kolmogorov's continuity criterion, cf. [38, Corollary 14.9], we can deduce that there exists a probability measure $P$ on $(\mathbb{W}^N \times \mathbb{W}^N, \mathcal{B}(\mathbb{W}^N) \otimes \mathcal{B}(\mathbb{W}^N))$ such that there is a subsequence $(m_k)_{k \in \mathbb{N}}$ along which $P^{m_k}$ converge towards $P$.

To characterize the limit, we first note that by Skorokhod representation theorem, cf. [9, Chapter 1, Theorem 6.7], without loss of generality we can assume that $(\{r^{i,m}, s^{i,m}\}_{i=1}^N)$ are defined on a common probability space $(\Omega, \mathcal{A}, P)$ with expectation $E$ and converge almost surely to $(\{r^i, s^i\}_{i=1}^N)$ with distribution $P$. Then, by $\mathbf{H}5$ and Lebesgue convergence theorem it holds almost surely for all $t \ge 0$,

$$ \lim_{m \to \infty} \int_0^t \tilde{b}(r_t^{i,m}) + \frac{1}{N} \sum_{j=1}^N g^m(r_j^{i,m}) du = \int_0^t \tilde{b}(r_t^i) + \frac{1}{N} \sum_{j=1}^N g^m(r_j^i) du. \quad (85) $$

Consider the mappings $M^{i,m}, N^{i,m} : \mathbb{W}^N \times \mathbb{W}^N \times \mathcal{P}(\mathbb{W}^N \times \mathbb{W}^N) \to \mathbb{W}$ defined by (82). Then for all $m \in \mathbb{N}$ and $i = 1, \dots, N$, $(M_t^{i,m}, \mathcal{F}_t, P^m)$ and $(N_t^{i,m}, \mathcal{F}_t, P^m)$ are martingales with respect to the canonical filtration $\mathcal{F}_t = \sigma((\{\mathbf{r}_u^i, \mathbf{s}_u^i\}_{i=1}^N)_{0 \le u \le t})$. Further the family $(\{M_t^{i,m}\}_{i=1}^N, P^m)_{m \in \mathbb{N}, t \ge 0}$ and $(\{N_t^{i,m}\}_{i=1}^N, P^m)_{m \in \mathbb{N}, t \ge 0}$ are uniformly integrable. In the same line as weak convergence is shown in the proof of Lemma 19 and by (85),
$\mathbb{P}^m \circ (\{\mathbf{r}^i, \mathbf{s}^i, M^{i,m}, N^{i,m}\}_{i=1}^N)^{-1}$ converges weakly to $\mathbb{P} \circ (\{\mathbf{r}^i, \mathbf{s}^i, M^i, N^i\}_{i=1}^N)^{-1}$ where

$$M_t^i = \mathbf{r}_t^i - \mathbf{r}_0^i - \int_0^t \left( \tilde{b}(\mathbf{r}_u^i) + \frac{1}{N} \sum_{j=1}^N g(\mathbf{r}_u^j) \right) du, \quad \text{and}$$

$$N_t^i = \mathbf{s}_t^i - \mathbf{s}_0^i - \int_0^t \left( \hat{b}(\mathbf{s}_u^i) + \frac{1}{N} \sum_{j=1}^N h(\mathbf{s}_u^j) \right) du.$$

Then $(\{M_t^i\}_{i=1}^N, \mathcal{F}_t, \mathbb{P})$ and $(\{N_t^i\}_{i=1}^N, \mathcal{F}_t, \mathbb{P})$ are continuous martingales using the same argument as in (68). Further, the quadratic variation $([\{M_t^i, N_t^i\}_{i=1}^N]_t)_{t \ge 0}$ exists $\mathbb{P}$-almost surely and is given by (83) and (84) $\mathbb{P}$-almost surely, which holds following the computations in the proof of Lemma 19 and Lemma 22. As in Lemma 22, we conclude by a martingale representation theorem that there are a probability space $(\Omega, \mathcal{A}, \mathcal{P})$ and a Brownian motion $\{W^i\}_{i=1}^N$ and random variables $(\{r^i\}_{i=1}^N, \{s^i\}_{i=1}^N)$ on this space such that $\mathbb{P} \circ (\{r^i, s^i\}_{i=1}^N)^{-1} = \mathbb{P} \circ (\{\mathbf{r}^i, \mathbf{s}^i\}_{i=1}^N)^{-1}$ and such that $(\{r^i, s^i, W^i\}_{i=1}^N)$ is a weak solution of (25).

By the Portmanteau theorem the monotonicity carries over to the limit, since $\mathbb{P}^m \circ (\{\mathbf{r}^i, \mathbf{s}^i\}_{i=1}^N)^{-1}$ converges weakly to $\mathbb{P} \circ (\{\mathbf{r}^i, \mathbf{s}^i\}_{i=1}^N)^{-1}$. $\square$

*Proof of Theorem 10.* The proof is a direct consequence of Lemma 22 and Lemma 23. $\square$

# A Appendix

## A.1 Kuramoto model

Lower bounds on the contraction rate can also be shown for nonlinear SDEs on the one-dimensional torus using the same approach. Here, we consider the Kuramoto model given by

$$dX_t = -k \left[ \int_{\mathbb{T}} \sin(X_t - x) d\mu_t(x) \right] dt + dB_t \quad (86)$$

on the torus $\mathbb{T} = \mathbb{R}/(2\pi\mathbb{Z})$.

**Theorem 24.** Let $\mu_t$ and $\nu_t$ be laws of $X_t$ and $Y_t$ where $(X_s)_{s \ge 0}$ and $(Y_s)_{s \ge 0}$ are two solutions of (86) with initial distributions $\mu_0$ and $\nu_0$ on $(\mathbb{T}, \mathcal{B}(\mathbb{T}))$, respectively. If

$$4k \int_{0}^{\pi} \exp(2k - 2k \cos(r/2)) dr \le 1 \quad (87)$$

holds, then for all $t \ge 0$,

$$W_{\tilde{f}}(\mu_t, \nu_t) \le e^{-c_{\mathbb{T}}t} W_{\tilde{f}}(\mu_0, \nu_0) \quad \text{and} \quad W_1(\mu_t, \nu_t) \le 2 \exp(2k) e^{-c_{\mathbb{T}}t} W_1(\mu_0, \nu_0),$$
where

$$c_{\mathbb{T}} = 1/(2 \int_{0}^{\pi} \int_{0}^{r} \exp[2k(\cos(r/2) - \cos(s/2))]dsdr) \quad (88)$$

and $\tilde{f}$ is a concave, increasing function given in (92).

In [18, Appendix A], a contraction result is stated for a general drift using a similar approach.

We prove Theorem 24 via a sticky coupling approach. In the same line as in Section 2 the coupling $(X_t, Y_t)_{t \ge 0}$ is defined as the weak limit of Markovian couplings $\{(X_t^\delta, Y_t^\delta)_{t \ge 0} : \delta > 0\}$ on $\mathbb{T} \times \mathbb{T} = \mathbb{R}/(2\pi\mathbb{Z}) \times \mathbb{R}/(2\pi\mathbb{Z})$ given by

$$\begin{aligned} dX_t^\delta &= -k \left[ \int_{\mathbb{T}} \sin(X_t^\delta - x) d\mu_t^\delta(x) \right] dt + \mathrm{rc}^\delta(\bar{r}_t^\delta) dB_t^1 + \mathrm{sc}^\delta(\bar{r}_t^\delta) dB_t^2 \\ dY_t^\delta &= -k \left[ \int_{\mathbb{T}} \sin(Y_t^\delta - x) d\nu_t^\delta(x) \right] dt - \mathrm{rc}^\delta(\bar{r}_t^\delta) dB_t^1 + \mathrm{sc}^\delta(\bar{r}_t^\delta) dB_t^2, \end{aligned} \quad (89)$$

where $\bar{r}_t^\delta = d_{\mathbb{T}}(X_t^\delta, Y_t^\delta)$ with $d_{\mathbb{T}}(\cdot, \cdot)$ defined by

$$d_{\mathbb{T}}(x, y) = \begin{cases} (|x - y| \bmod 2\pi) & \text{if } (|x - y| \bmod 2\pi) \le \pi, \\ (2\pi - |x - y| \bmod 2\pi) & \text{otherwise.} \end{cases} \quad (90)$$

The functions $\mathrm{rc}^\delta, \mathrm{sc}^\delta$ are given by (19) and satisfy that there exists $\epsilon_0 > 0$ such that $\mathrm{rc}^\delta(r) \ge r/2$ for any $0 \le r \le \delta \le \epsilon_0$.

**Theorem 25.** Assume (87). Let $\mu_0$ and $\nu_0$ be probability measures on $(\mathbb{T}, \mathcal{B}(\mathbb{T}))$ having finite fourth moment. Then, $(X_t, Y_t)_{t \ge 0}$ is the subsequential limit as $\delta \to 0$ of $\{(X_t^\delta, Y_t^\delta)_{t \ge 0} : \delta > 0\}$, where $(X_t)_{t \ge 0}$ and $(Y_t)_{t \ge 0}$ are solutions of (86) with initial distributions $\mu_0$ and $\nu_0$, respectively. Further, there exists a process $(r_t)_{t \ge 0}$ satisfying for any $t \ge 0$, $d_{\mathbb{T}}(X_t, Y_t) \le r_t$ almost surely, and which is a weak solution of

$$d r_t = (2k \sin(r_t/2) + 2k \mathbb{P}(r_t)) dt + 2 \mathbb{1}_{(0,\pi]}(r_t) dW_t - 2 d \ell_t^\pi, \quad (91)$$

where $(W_t)_{t \ge 0}$ is a one-dimensional Brownian motion on $\mathbb{T}$ and $\ell^\pi$ is the local time at $\pi$.

*Proof.* The proof works analogously to the proof of Theorem 2 stated in Section 6.2.2. It holds similarly to Lemma 12 by Meyer-Tanaka's formula, cf. [51, Chapter 6, Theorem 1.1], and using (90),

$$\begin{aligned} \bar{r}_t^\delta - \bar{r}_0^\delta = & \int_0^t \operatorname{sgn}(X_t^\delta - Y_t^\delta)(-k)e_t \left[ \int_{\mathbb{T}} \sin(X_t^\delta - x) d\mu_t(x) - \int_{\mathbb{T}} \sin(Y_t^\delta - x) d\nu_t(x) \right] dt \\ & + \int_0^t \operatorname{sgn}(X_t^\delta - Y_t^\delta) 2\mathrm{rc}^\delta(\bar{r}_t^\delta) e_t dB_t^1 + \int_{\mathbb{R}} 2\mathrm{rc}^\delta(\bar{r}_t^\delta)^2 \ell_t^a (\delta_0 - \delta_\pi)(da), \end{aligned}$$
where sgn($x$) = 1$_{(0,π]}$(x) - 1$_{(π,2π]}$(x), ($ℓ_t^a$)$_{t≥0}$ is the local time at a associated with ($X_t^δ$ - $Y_t^δ$)$_{t≥0}$ and $e_t = (X_t^δ - Y_t^δ)/d_Π(X_t^δ, Y_t^δ)$ for $\bar{r}_t^δ \neq 0$. For $\bar{r}_t^δ = 0$, $e_t$ is some arbitrary unit vector. For any a the support of $ℓ_t^a$ as a function of t is a subset of the set of t such that $r_t = a$ [38, Theorem 19.1], hence 1$_{(0,π]}$(r_t)ℓ_t^0 = 0 almost surely. so the term involving the local time reduces to -2ℓ_t^π. Further, we note that $W_t = \int_0^t \text{sgn}(X_t^δ - Y_t^δ)e_t dB_t^1$ is a Brownian motion. As in Lemma 12, it holds for the process $(\bar{r}_t^δ)_{t≥0}$ for $\epsilon < \epsilon_0$ with $\epsilon_0$ given by (20),

$$d\bar{r}_t^\delta \leq (2k \sin(\bar{r}_t^\delta/2) + 2k \mathbb{E}_{x \sim \mu_t^\delta, y \sim \nu_t^\delta} (\text{rc}^\epsilon(d_\Pi(x, y)))) dt + 2\text{rc}^\delta(\bar{r}_t^\delta) dW_t - 2d\ell_t^\pi,$$

where we used the properties of rc$^δ$ and

$$(x - y) \cdot (\sin(x - \tilde{x}) - \sin(y - \tilde{x})) \leq 2 \sin(|x - y|/2) |x - y|$$

for any $x, y, \tilde{x} \in \mathbb{T}$. Consider $(r_t^{\delta, \epsilon})_{t \ge 0}$ given by

$$d r_t^{\delta, \epsilon} = (2k \sin(r_t^{\delta, \epsilon}/2) + 2k \int_0^{\pi} \text{rc}^{\epsilon}(u) d P_t^{\delta, \epsilon}(u)) dt + 2\text{rc}^{\delta}(r_t^{\delta, \epsilon}) d W_t - 2d \ell_t^{\pi},$$

where $P_t^{\delta, \epsilon}$ is the law of $r_t^{\delta, \epsilon}$. Then as in Lemma 13, for the processes $(\bar{r}_t^\delta)_{t \ge 0}$ and $(r_t^{\delta, \epsilon})_{t \ge 0}$ with the same initial condition and driven by the same noise it holds $\bar{r}_t^\delta \le r_t^{\delta, \epsilon}$ almost surely for every $t$ and $\epsilon < \epsilon_0$.

Consider the process $(U_t^{\delta, \epsilon})_{t \ge 0} = (X_t^\delta, Y_t^\delta, r_t^{\delta, \epsilon})_{t \ge 0}$ on $\mathbb{T}^2 \times [0, \pi]$ for each $\epsilon, \delta > 0$. We define by $\mathbf{X}, \mathbf{Y} : C(\mathbb{R}_+, \mathbb{T}^2 \times [0, \pi]) \to C(\mathbb{R}_+, \mathbb{T})$ and $\mathbf{r} : C(\mathbb{R}_+, \mathbb{T}^2 \times [0, \pi]) \to C(\mathbb{R}_+, [0, \pi])$ the canonical projections onto the first component, onto the second component and onto the last component, respectively. Analogously to the proof of Theorem 2, the law $P^{\delta, \epsilon}$ of the process $(U_t^{\delta, \epsilon})_{t \ge 0}$ converges along a subsequence $(\delta_k, \epsilon_k)_{k \in \mathbb{N}}$ to a probability measure $P$. Let $(X_t, Y_t, r_t)_{t \ge 0}$ be some process on $\mathbb{T}^2 \times [0, \pi]$ with distribution $P$ on $(\bar{\Omega}, \bar{\mathcal{F}}, \bar{P})$. Since $(X_t^\delta)_{t \ge 0}$ and $(Y_t^\delta)_{t \ge 0}$ are solutions of (86) which are unique in law, we have that for any $\epsilon, \delta > 0$, $P^{\delta, \epsilon} \circ (\mathbf{X})^{-1} = P \circ (\mathbf{X})^{-1}$ and $P^{\delta, \epsilon} \circ (\mathbf{Y})^{-1} = P \circ (\mathbf{Y})^{-1}$. And therefore $(X_t)_{t \ge 0}$ and $(Y_t)_{t \ge 0}$ are solutions of (86) as well with the same initial condition. Hence $P \circ (\mathbf{X}, \mathbf{Y})^{-1}$ is a coupling of two copies of (86).

Further, the monotonicity $\bar{r}_t^\delta \le r_t^{\delta, \epsilon}$ carries over to the limit by the Portmanteau theorem. Finally, similarly to the proof of Lemma 19 and Lemma 20 there exist an extended probability space and a one-dimensional Brownian motion $(W_t)_{t \ge 0}$ such that $(r_t, W_t)_{t \ge 0}$ is a solution to (97). □

*Proof of Theorem 24.* Similarly to (37) we consider a function $\tilde{f}$ on $[0, \pi]$ defined by

$$\tilde{f}(t) = \int_{0}^{t} \tilde{\varphi}(r) \tilde{g}(r) dr, \qquad (92)$$
where

$$
\begin{align*}
\tilde{\varphi}(r) &= \exp\{2k(\cos(r/2) - 1)\}, \quad \tilde{\Phi}(r) = \int_0^r \tilde{\varphi}(s)ds, \\
\tilde{g}(r) &= 1 - \frac{c_T}{2} \int_0^r \{\tilde{\Phi}(s)/\tilde{\varphi}(s)\}ds - k \int_0^r \{1/\tilde{\varphi}(s)\}ds, \\
c_T &= \left(2 \int_0^\pi \{\tilde{\Phi}(s)/\tilde{\varphi}(s)\}ds\right)^{-1}.
\end{align*}
$$

Then for $k$ satisfying (87), $\tilde{g}(r) \in [1/2, 1]$ and $\tilde{f}$ is a concave function satisfying similarly to (38)

$$
\exp(-2k)/2r \le \tilde{f} \le \tilde{\Phi}(r) \le r \qquad (93)
$$

and

$$
\begin{align}
\tilde{f}''(0) &= -k \nonumber \\
2(\tilde{f}''(r) - \tilde{f}''(0)) &\le -2k \sin(r/2) \tilde{f}'(r) - c_T \tilde{f}(r) \quad \text{for all } r \in [0, \pi]. \tag{94}
\end{align}
$$

By Ito's formula it holds

$$
\mathrm{d}\tilde{f}(r_t) = \tilde{f}'(r_t)(2k \sin(r/2) + 2k\mathbb{P}(r_t > 0))\mathrm{d}t + 2\tilde{f}'(r_t)\mathbb{1}_{(0,\pi]}(r_t)\mathrm{d}W_t - 2\tilde{f}'(r_t)\mathrm{d}\ell_t^\pi \\
+ 2\tilde{f}''(r_t)\mathbb{1}_{(0,\pi]}(r_t)\mathrm{d}t.
$$

Taking expectation and using that the term involving the local time is negative, we obtain

$$
\begin{align*}
\frac{\mathrm{d}}{\mathrm{d}t} \mathbb{E}[\tilde{f}(r_t)] &\le \mathbb{E}[2(\tilde{f}''(r_t) - f''(0)) + \tilde{f}'(r_t)2k \sin(r_t/2)] + (2\tilde{f}''(0) + 2k)\mathbb{P}(r_t > 0) \\
&\le -c_{\mathbb{T}}\mathbb{E}[\tilde{f}(r_t)],
\end{align*}
$$

where the last step holds by (94). Then

$$
\mathbb{E}[\tilde{f}(d_{\mathbb{T}}(\bar{X}_t, \bar{Y}_t))] \le \mathbb{E}[\tilde{f}(r_t)] \le e^{-c_{\mathbb{T}}t} \mathbb{E}[\tilde{f}(r_0)] = e^{-c_{\mathbb{T}}t} \mathbb{E}[\tilde{f}(d_{\mathbb{T}}(\bar{X}_0, \bar{Y}_0))], \quad (95)
$$

provided (87) holds. Thus

$$
\mathcal{W}_{\tilde{f}}(\mu_t, \nu_t) \le e^{-c_T t} \mathcal{W}_{\tilde{f}}(\mu_0, \nu_0),
$$

and by (93)

$$
W_1(\mu_t, \nu_t) \le 2 \exp(2k) e^{-c_T t} W_1(\mu_0, \nu_0).
$$

<figure><img src="image_5.png" /></figure>
**Remark 26.** Let us finally remark that we can relax the condition (87) and we can obtain contraction with a modified contraction rate $c_T$ for all $k < k_0$, where $k_0$ is given by

$$k_0 \int_0^\pi \exp(2k_0 - 2k_0 \cos(r/2)) dr = 1. \qquad (96)$$

More precisely, set $\zeta = 1 - k \int_0^\pi \exp(2k - 2k \cos(r/2)) dr$ and $c_T = \zeta (\int_0^\pi \{\tilde{\Phi}(s)/\tilde{\varphi}(s)\} ds)^{-1}$. Then, $\tilde{g}(r) \in [\zeta/2, 1]$ and $\zeta \exp(-2k)/2r \le \tilde{f}(r) \le r$. Following the previous computations, we obtain

$$W_1(\mu_t, \nu_t) \le 2 \exp(2k)/\zeta e^{-c_T t} W_1(\mu_0, \nu_0),$$

where for $k$ close to $k_0$, the contraction rate becomes small and the prefactor $2 \exp(2k)/\zeta$ explodes.

## A.2 Sticky nonlinear SDEs on bounded state space

In the same line as in Theorem 3, existence, uniqueness in law and comparison results hold for solutions to the sticky SDE on $[0, \pi]$ given by

$$dr_t = (\tilde{b}(r_t) + 2kP(r_t > 0))dt + 21_{(0,\pi)}(r_t)dW_t - 2d\ell_t^\pi, \qquad (97)$$

where $k \in \mathbb{R}_+$ and $\ell^\pi$ is the local time at $\pi$.

The analysis of invariant measures and phase transitions can be easily adapted to the case of the sticky SDE on $[0, \pi]$ given by (97).

**Theorem 27.** Let $(r_t)_{t \ge 0}$ be a solution of (97) with drift $\tilde{b}$ satisfying **H1**. Then, the Dirac measure at zero, $\delta_0$, is an invariant probability measure on $[0, \pi]$ for (97). If there exists $p \in (0, 1)$ solving $(1/k) = (1-p)I(k, p)$ where

$$I(k, p) = \int_0^\pi \exp \left( kpx + \frac{1}{2} \int_0^x \tilde{b}(r) dr \right) dx,$$

then the probability measure $\pi$ on $[0, \pi]$ given by

$$\pi(dx) \propto \frac{1}{kp} \delta_0(dx) + \exp \left( kpx + \frac{1}{2} \int_0^x \tilde{b}(r) dr \right) \lambda_{(0,\pi)}(dx) \qquad (98)$$

is another invariant probability measure for (97).

*Proof of Theorem 27.* The proof works analogously to the proof of Theorem 5 for sticky SDEs on $\mathbb{R}_+$. Note that here the condition (75) transforms for $p \in (0, 1]$ to

$$p = \pi((0, \pi)) = \frac{I(k, p)}{1/(kp) + I(k, p)} \Leftrightarrow (1-p)I(k, p) = 1/k.$$

□
**Example 28.** Consider a solution $(r_t)_{t \ge 0}$ of (97) with drift $\tilde{b}(r) = 2k \sin(r/2)$. Consider a solution $p \in (0, 1]$ solving $1/k = (1-p)I(k, p)$ with

$$I(k, p) = \int_{0}^{\pi} \exp \left( kpx + \int_{0}^{x} k \sin(r/2) dr \right) dx = \int_{0}^{\pi} \exp \left( kpx + 2k - 2k \cos(x/2) \right) dx.$$

Then by Theorem 27, the Dirac measure at zero, $\delta_0$ and the probability measure

$$\pi(dx) \propto \frac{1}{kp} \delta_0(dx) + \exp(kpx + 2k - 2k \cos(x/2)) \lambda_{(0,\pi)}(dx) \quad (99)$$

are invariant probability measures for (97). We specify a necessary and sufficient condition for the existence of a solution $p$ satisfying $1/k = (1-p)I(k, p)$. We define $\hat{I}(k, p) = (1-p)I(k, p)$. We first consider the case $1/k < \hat{I}(k, 0) = \int_{0}^{\pi} \exp(2k - 2k \cos(x/2))dx$. Then since $1/k > \hat{I}(k, 1) = 0$ and by the mean value theorem there exists a $p$ solving $1/k = \hat{I}(k, p)$ and therefore there exist multiple invariant distributions for (99). On the other hand, if $1/k > \hat{I}(k, 0) = \int_{0}^{\pi} \exp(2k - 2k \cos(x/2))dx$, since $\pi \le \int_{0}^{\pi} \exp(2k - 2k \cos(x/2))dx$ and for $k < 1/\pi$, it holds

$$\begin{aligned} \frac{d}{dp} \hat{I}(k, p) &= -I(k, p) + (1-p) \int_{0}^{\pi} kx \exp(kpx + 2k - 2k \cos(x/2))dx \\ &= \int_{0}^{\pi} ((1-p)kx - 1) \exp(kpx + 2k - 2k \cos(x/2))dx \le 0, \end{aligned}$$

there is no $p$ satisfying (99).

**Remark 29.** The contraction result given in Theorem 7 carries over to the sticky diffusion $(r_t)$ given by (97) on $[0, \pi]$ with $\tilde{b}(r) = 2k \sin(r/2)$. If (87) holds, then for $t \ge 0$, (31) holds with $\tilde{f}$ defined in (92) and $c_T$ defined in (88) using (95). Moreover by Remark 26, we can deduce that if (96) holds, the Dirac measure at zero, $\delta_0$, is the unique invariant measure and contraction towards $\delta_0$ holds.

## Acknowledgments

### Support:...

A. E. and K. S. have been supported by the Hausdorff Center for Mathematics. Gefördert durch die Deutsche Forschungsgemeinschaft (DFG) im Rahmen der Exzellenzstrategie des Bundes und der Länder - GZ 2047/1, Projekt-ID 390685813. The work of A.G. has been (partially) supported by the Project EFI ANR-17-CE40-0030 of the French National Research Agency
References

[1] J. Acebron, L. Bonilla, C. Pérez-Vicente, F. Farran, and R. Spigler. The kuramoto model: A simple paradigm for synchronization phenomena. *Reviews of Modern Physics*, 77, 04 2005.

[2] Julien Barré, Dan Crisan, and Thierry Goudon. Two-dimensional pseudo-gravity model: particles motion in a non-potential singular force field. *Trans. Amer. Math. Soc.*, 371(4):2923–2962, 2019.

[3] R. F. Bass. A stochastic differential equation with a sticky point. *Electron. J. Probab.*, 19:no. 32, 22, 2014.

[4] S. Benachour, B. Roynette, D. Talay, and P. Vallois. Nonlinear self-stabilizing processes. I. Existence, invariant probability, propagation of chaos. *Stochastic Process. Appl.*, 75(2):173–201, 1998.

[5] S. Benachour, B. Roynette, and P. Vallois. Nonlinear self-stabilizing processes. II. Convergence to invariant probability. *Stochastic Process. Appl.*, 75(2):203–224, 1998.

[6] D. Benedetto, E. Caglioti, J. A. Carrillo, and M. Pulvirenti. A non-Maxwellian steady distribution for one-dimensional granular media. *J. Statist. Phys.*, 91(5-6):979–990, 1998.

[7] L. Bertini, G. Giacomin, and K. Pakdaman. Dynamical aspects of mean field plane rotators and the Kuramoto model. *J. Stat. Phys.*, 138(1-3):270–290, 2010.

[8] L. Bertini, G. Giacomin, and C. Poquet. Synchronization and random long time dynamics for mean-field plane rotators. *Probab. Theory Related Fields*, 160(3-4):593–653, 2014.

[9] Patrick Billingsley. *Convergence of probability measures*. Wiley Series in Probability and Statistics: Probability and Statistics. John Wiley & Sons, Inc., New York, second edition, 1999. A Wiley-Interscience Publication.

[10] F. Bolley, I. Gentil, and A. Guillin. Uniform convergence to equilibrium for granular media. *Arch. Ration. Mech. Anal.*, 208(2):429–445, 2013.

[11] N. Bou-Rabee, A. Eberle, and R. Zimmer. Coupling and convergence for Hamiltonian Monte Carlo. *Ann. Appl. Probab.*, 30(3):1209–1250, 2020.

[12] J. A. Carrillo, R. S. Gvalani, G. A. Pavliotis, and A. Schlichting. Long-time behaviour and phase transitions for the Mckean-Vlasov equation on the torus. *Arch. Ration. Mech. Anal.*, 235(1):635–690, 2020.
[13] J. A. Carrillo, R. J. McCann, and C. Villani. Kinetic equilibration rates for granular media and related equations: entropy dissipation and mass transportation estimates. *Rev. Mat. Iberoamericana*, 19(3):971–1018, 2003.

[14] J. A. Carrillo, R. J. McCann, and C. Villani. Contractions in the 2-Wasserstein length space and thermalization of granular media. *Arch. Ration. Mech. Anal.*, 179(2):217–263, 2006.

[15] P. Cattiaux, A. Guillin, and F. Malrieu. Probabilistic approach for granular media equations in the non-uniformly convex case. *Probab. Theory Related Fields*, 140(1-2):19–40, 2008.

[16] R. Chitashvili. On the nonexistence of a strong solution in the boundary problem for a sticky Brownian motion. *Proc. A. Razmadze Math. Inst.*, 115:17–31, 1997.

[17] F. Delarue and A. Tse. Uniform in time weak propagation of chaos on the torus. *arXiv preprint arXiv:2104.14973*, 2021.

[18] M. Delgadino, R. Gvalani, and G. Pavliotis. On the diffusive-mean field limit for weakly interacting diffusions exhibiting phase transitions. *arXiv preprint arXiv:2001.03920v1*, 2020.

[19] A. Durmus, A. Eberle, A. Guillin, and R. Zimmer. An elementary approach to uniform in time propagation of chaos. *Proc. Amer. Math. Soc.*, 148(12):5387–5398, 2020.

[20] A. Eberle. Reflection couplings and contraction rates for diffusions. *Probab. Theory Related Fields*, 166(3-4):851–886, 2016.

[21] A. Eberle, A. Guillin, and R. Zimmer. Couplings and quantitative contraction rates for Langevin dynamics. *Ann. Probab.*, 47(4):1982–2010, 2019.

[22] A. Eberle, A. Guillin, and R. Zimmer. Quantitative Harris-type theorems for diffusions and McKean-Vlasov processes. *Trans. Amer. Math. Soc.*, 371(10):7135–7173, 2019.

[23] A. Eberle and R. Zimmer. Sticky couplings of multidimensional diffusions with different drifts. *Ann. Inst. Henri Poincaré Probab. Stat.*, 55(4):2370–2394, 2019.

[24] H.-J. Engelbert and G. Peskir. Stochastic differential equations for sticky Brownian motion. *Stochastics*, 86(6):993–1021, 2014.

[25] W. Feller. Diffusion processes in one dimension. *Trans. Amer. Math. Soc.*, 77:1–31, 1954.

[26] W. Feller. The general diffusion operator and positivity preserving semi-groups in one dimension. *Ann. of Math. (2)*, 60:417–436, 1954.
[27] T. Funaki. A certain class of diffusion processes associated with nonlinear parabolic equations. *Z. Wahrsch. Verw. Gebiete*, 67(3):331–348, 1984.

[28] A. Guillin, W. Liu, L. Wu, and C. Zhang. Uniform poincaré and logarithmic sobolev inequalities for mean field particles systems. *arXiv preprint arXiv:1909.07051v1*, to appear in *Annals of Applied Probability*, 2019.

[29] A. Guillin and P. Monmarché. Uniform long-time and propagation of chaos estimates for mean field kinetic particles in non-convex landscapes. *arXiv preprint arXiv:2003.00735v2*, 2020.

[30] W. Hammersley, D. Siska, and L. Szpruch. McKean-Vlasov SDEs under Measure Dependent Lyapunov Conditions. *arXiv preprint arXiv:1802.03974*, 2018.

[31] M. Hauray and S. Mischler. On Kac's chaos and related problems. *J. Funct. Anal.*, 266(10):6055–6157, 2014.

[32] M. Hofmanová and J. Seidler. On weak solutions of stochastic differential equations. *Stoch. Anal. Appl.*, 30(1):100–121, 2012.

[33] M. Hofmanová and J. Seidler. On weak solutions of stochastic differential equations II. *Stoch. Anal. Appl.*, 31(4):663–670, 2013.

[34] N. Ikeda. On the construction of two-dimensional diffusion processes satisfying Wentzell's boundary conditions and its application to boundary value problems. *Mem. Coll. Sci. Univ. Kyoto Ser. A. Math.*, 33:367–427, 1960/61.

[35] N. Ikeda and S. Watanabe. *Stochastic differential equations and diffusion processes*, volume 24 of *North-Holland Mathematical Library*. North-Holland Publishing Co., Amsterdam; Kodansha, Ltd., Tokyo, second edition, 1989.

[36] J. Jacod and A. N. Shiryaev. *Limit theorems for stochastic processes*, volume 288 of *Grundlehren der Mathematischen Wissenschaften* [Fundamental Principles of Mathematical Sciences]. Springer-Verlag, Berlin, second edition, 2003.

[37] M. Kac. Foundations of kinetic theory. In *Proceedings of the Third Berkeley Symposium on Mathematical Statistics and Probability, 1954–1955, vol. III*, pages 171–197. University of California Press, Berkeley and Los Angeles, 1956.

[38] Olav Kallenberg. *Foundations of modern probability*. Probability and its Applications (New York). Springer-Verlag, New York, second edition, 2002.

[39] A. Klenke. *Probability Theory: A Comprehensive Course*. Universitext. Springer London, 2013.
[40] D. Lacker. Hierarchies, entropy and quantitative propagation of chaos for mean field diffusionss. *arXiv preprint arXiv:2105.02983*, 2021.

[41] W. Liu, L. Wu, and C. Zhang. Long-time behaviors of mean-field interacting particle systems related to mckean-vlasov equations. *arXiv preprint arXiv:2007.09462v1*, 2020.

[42] F. Malrieu. Convergence to equilibrium for granular media equations and their Euler schemes. *Ann. Appl. Probab.*, 13(2):540–560, 2003.

[43] H. P. McKean, Jr. A. Skorohod's stochastic integral equation for a reflecting barrier diffusion. *J. Math. Kyoto Univ.*, 3:85–88, 1963.

[44] H. P. McKean, Jr. A class of Markov processes associated with nonlinear parabolic equations. *Proc. Nat. Acad. Sci. U.S.A.*, 56:1907–1911, 1966.

[45] S. Méléard. Asymptotic behaviour of some interacting particle systems; McKean-Vlasov and Boltzmann models. In *Probabilistic models for nonlinear partial differential equations (Montecatini Terme, 1995)*, volume 1627 of *Lecture Notes in Math.*, pages 42–95. Springer, Berlin, 1996.

[46] S. Mischler and C. Mouhot. Kac's program in kinetic theory. *Invent. Math.*, 193(1):1–147, 2013.

[47] S. Mischler, C. Mouhot, and B. Wennberg. A new approach to quantitative propagation of chaos for drift, diffusion and jump processes. *Probab. Theory Related Fields*, 161(1-2):1–59, 2015.

[48] S. Mishura and A.Y. Veretennikov. Existence and uniqueness theorems for solutions of mckean-vlasov stochastic equations. *arXiv preprint arXiv:1603.02212*, 2016.

[49] G. Peskir. *On boundary behaviour of one-dimensional diffusions: From Brown to Feller and beyond*. In William Feller – Selected Papers II. Springer, Berlin, 2015.

[50] Philip E. Protter. *Stochastic integration and differential equations*, volume 21 of *Applications of Mathematics (New York)*. Springer-Verlag, Berlin, second edition, 2004. Stochastic Modelling and Applied Probability.

[51] Daniel Revuz and Marc Yor. *Continuous martingales and Brownian motion*, volume 293 of *Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]*. Springer-Verlag, Berlin, third edition, 1999.

[52] A.V. Skorokhod. Stochastic equations for diffusion processes in a bounded region. *Theory Probab. Appl.*, 6(3):264–274, 1961.

[53] A.V. Skorokhod. Stochastic equations for diffusion processes in a bounded region. ii. *Theory Probab. Appl.*, 7(1):3–23, 1962.
[54] D. W. Stroock and S. R. S. Varadhan. Diffusion processes with boundary conditions. *Comm. Pure Appl. Math.*, 24:147–225, 1971.

[55] D. W. Stroock and S. R. S. Varadhan. *Multidimensional diffusion processes*, volume 233 of *Grundlehren der Mathematischen Wissenschaften* [Fundamental Principles of Mathematical Sciences]. Springer-Verlag, Berlin-New York, 1979.

[56] A.-S. Sznitman. Topics in propagation of chaos. In *École d'Été de Probabilités de Saint-Flour XIX—1989*, volume 1464 of *Lecture Notes in Math.*, pages 165–251. Springer, Berlin, 1991.

[57] J. Tugaut. Self-stabilizing processes in multi-wells landscape in $\mathbb{R}^d$-invariant probabilities. *J. Theoret. Probab.*, 27(1):57–79, 2014.

[58] C. Villani. *Topics in optimal transportation*, volume 58 of *Graduate Studies in Mathematics*. American Mathematical Society, Providence, RI, 2003.

[59] J. Warren. Branching processes, the Ray-Knight theorem, and sticky Brownian motion. In *Séminaire de Probabilités, XXXI*, volume 1655 of *Lecture Notes in Math.*, pages 1–15. Springer, Berlin, 1997.

[60] S. Watanabe. On stochastic differential equations for multi-dimensional diffusion processes with boundary conditions. *J. Math. Kyoto Univ.*, 11:169–180, 1971.

[61] S. Watanabe. On stochastic differential equations for multi-dimensional diffusion processes with boundary conditions. II. *J. Math. Kyoto Univ.*, 11:545–551, 1971.

UNIVERSITÉ PARIS-SACLAY, ENS PARIS-SACLAY, CNRS,
CENTRE BORELLI, F-91190 GIF-SUR-YVETTE, FRANCE.

*E-mail address: alain.durmus@cmla.ens-cachan.fr*

UNIVERSITÄT BONN, INSTITUT FÜR ANGEWANDTE MATHEMATIK, ENDENICHER ALLEE 60, 53115 BONN, GERMANY

*E-mail address: eberle@uni-bonn.de, katharina.schuh@uni-bonn.de*

LABORATOIRE DE MATHÉMATIQUES BLAISE PASCAL, CNRS-UMR 6620, UNIVERSITÉ CLERMONT-AUVERGNE, AVENUE DE LANDAIS, 63177 AUBIERE CEDEX, FRANCE

*E-mail address: guillin@math.univ-bpclermont.fr*