Hostname: page-component-78c5997874-94fs2 Total loading time: 0 Render date: 2024-11-18T02:19:08.126Z Has data issue: false hasContentIssue false

STABILITY ANALYSIS FOR STOCHASTIC MCKEAN–VLASOV EQUATION

Published online by Cambridge University Press:  06 November 2024

C. SHI
Affiliation:
Department of Mathematics, Nanjing University of Science and Technology, Nanjing, China; e-mail: njueduscg@163.com
W. WANG*
Affiliation:
Department of Mathematics, Nanjing University, Nanjing, China
Rights & Permissions [Opens in a new window]

Abstract

The pth ($p\geq 1$) moment exponential stability, almost surely exponential stability and stability in distribution for stochastic McKean–Vlasov equation are derived based on some distribution-dependent Lyapunov function techniques.

Type
Research Article
Copyright
© The Author(s), 2024. Published by Cambridge University Press on behalf of Australian Mathematical Publishing Association Inc.

1. Introduction

McKean–Vlasov stochastic differential equations (SDEs), originating from the seminal works [Reference McKean15, Reference Vlasov18], are also known as mean-field SDEs or distribution-dependent SDEs which are used to study the interacting particle system and mean-field games. There are numerous works on the well-posedness, ergodicity and large deviations [Reference Hong, Li and Liu10, Reference Liu and Ma13, Reference Reis, Salkeld and Tugaut17, Reference Wang19]. Moreover, there are also several works on the stability of the McKean–Vlasov SDEs. Recently, Ding and Qiao [Reference Ding and Qiao5] considered the stability for the McKean–Vlasov SDEs with non-Lipschitz coefficients

(1.1) $$ \begin{align} dX(t)&=b(X(t),\mathcal{L}(X(t)))\,dt+\sigma(X(t),\mathcal{L}(X(t)))\,dW(t),\nonumber\\ X(0)&=x_{0}, \end{align} $$

where $\mathcal {L}(X(t))$ is the distribution of $X(t)$ , $W_{\cdot }=(W^{1}_{\cdot },W^{2}_{\cdot },\dots ,W^{l}_{\cdot })$ is a $\mathcal {F}_{t}$ -adapted standard Brownian motion and the coefficients $b: \mathbb {R}^{d}\times \mathcal {M}_{\lambda ^{2}}(\mathbb {R}^{d})\to \mathbb {R}^{d}$ and $\sigma : \mathbb {R}^{d}\times \mathcal {M}_{\lambda ^{2}}(\mathbb {R}^{d})\to \mathbb {R}^{d}\times \mathbb {R}^{l}$ are Borel measurable functions. The definition of $\mathcal {M}_{\lambda ^{2}}(\mathbb {R}^{d})$ is defined in the next section. Sufficient conditions are given for the exponential stability of the second moments for their solutions in terms of a Lyapunov function [Reference Khasminskii12]. Furthermore, the almost surely (a.s.) asymptotic stability of their solutions is also discussed. Lv and Shan [Reference Lv and Shan14] considered the long time behaviour of stochastic McKean–Vlasov equations, and the exponential and logarithmic decay are discussed. Bahlali et al. [Reference Bahlali, Mezerdi and Mezerdi1] discussed the existence and uniqueness of solutions under a non-Lipschitz condition and derived various stability properties with respect to initial data, coefficients and driving processes. Wu et al. [Reference Wu, Hu, Gao and Yuan20] studied the stability of solutions of McKean–Vlasov SDEs via feedback control based on discrete-time state observation and derived the $\mathrm {H}_{\infty }$ stability, asymptotic stability and exponential stability in mean square for the controlled systems.

In this paper, we first provide a sufficient condition for the pth moment exponential stability and a.s. exponential stability for (1.1) (Theorem 3.2) by using the classical Lyapunov function method. Furthermore, asymptotic stability in distribution is derived by introducing a distribution-dependent operator, together with a similar discussion as that for SDE with Markovian switching [Reference Yuan and Mao21].

There are many recent works on the stability in distribution for stochastic differential equations with distribution-independent coefficients. Yuan et al. [Reference Yuan, Zou and Mao22] discussed the stochastic differential equations with Markovian switching and investigated the stability in distribution of the equations. Further, Du et al. [Reference Du, Dang and Dieu6] improved the result of Yuan et al. [Reference Yuan, Zou and Mao22] by giving a new sufficient condition for stability in distribution. Bao et al. [Reference Bao, Hou and Yuan2] considered a neutral stochastic differential delay equation with Markovian switching and obtained sufficient conditions for stability in distribution. Fei et al. [Reference Fei, Fei, Deng and Mao8] considered the stability in distribution for a highly nonlinear stochastic differential equation driven by G-Brownian motion [Reference Peng16].

The rest of the paper is organized as follows. In Section 2, we recall some preliminary knowledge. The pth moment exponential stability and a.s. exponential stability is presented in Section 3.1. The stability in distribution is established in Section 3.2.

2. Preliminary and main result

Let $C(\mathbb {R}^{d})$ be the collection of continuous functions on $\mathbb {R}^{d}$ . For convenience, we denote the norm of vectors and matrices by $|\cdot |$ and $\|\cdot \|$ , respectively. Furthermore, let $\langle \cdot ,\cdot \rangle $ denote the scalar product in $\mathbb {R}^{d}$ . Let $\mathcal {B}(\mathbb {R}^{d})$ be the Borel $\sigma $ -algebra on $\mathbb {R}^{d}$ and $\mathcal {P}(\mathbb {R}^{d})$ denote the space of all probability measures defined on $\mathcal {B}(\mathbb {R}^{d})$ with the topology of weak convergence.

For $\lambda (x)=1+|x|$ , $x\in \mathbb {R}^d$ , define the Banach space

$$ \begin{align*} C_{\lambda}(\mathbb{R}^{d})=\bigg\{\phi\in C(\mathbb{R}^{d})\bigg| \|\phi\|_{C_{\lambda}(\mathbb{R}^{d})}=\sup_{x\in \mathbb{R}^{d}}\frac{|\phi(x)|}{\lambda^2(x)}+\sup_{x\neq y}\frac{|\phi(x)-\phi(y)|}{|x-y|}<\infty\bigg\}. \end{align*} $$

Let $\mathcal {M}_{\lambda ^{p}}^{s}(\mathbb {R}^{d})$ be the space of signed measures m on $\mathcal {B}(\mathbb {R}^{d})$ satisfying

$$ \begin{align*} \|m\|_{\lambda^{p}}^{p}=\int_{\mathbb{R}^{d}}\lambda^p(x)|m|(dx)<\infty,\quad p\geq2, \end{align*} $$

where $|m|=m^{+}+m^{-}$ , and $m=m^{+}-m^{-}$ is the Jordan decomposition of m. Let $\mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d})=\mathcal {M}_{\lambda ^{p}}^{s}(\mathbb {R}^{d})\cap \mathcal {P}(\mathbb {R}^{d})$ be the set of probability measures on $\mathcal {B}(\mathbb {R}^{d})$ with finite pth moments equipped with the metric,

$$ \begin{align*} \rho(\mu,\nu)\triangleq\sup_{\|\phi\|_{C_{\lambda}(\mathbb{R}^{d})}\leq1}\bigg|\!\int_{\mathbb{R}^{d}}\phi(x)\mu(dx)-\int_{\mathbb{R}^{d}}\phi(x)\nu(dx)\bigg|. \end{align*} $$

Then, $(\mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}),\rho )$ is a complete metric space.

Given a complete filtered probability space $(\Omega , \mathcal {F}, \{\mathcal {F}_{t}\}_{t\in [0,\infty )},\mathbb {P})$ , we recall the definition of the derivative for a function with respect to a probability measure [Reference Cardaliaguet4]. A function $f: \mathcal {M}_{\lambda ^{p}}^{s}(\mathbb {R}^{d})\to \mathbb {R}$ is differential at $\mu \in \mathcal {M}_{\lambda ^{p}}^{s}(\mathbb {R}^{d})$ if for $\tilde {f}(\xi )\triangleq f(\mathbb {P}_{\xi }), \xi \in L^{p}(\Omega; \mathbb {R}^{d})$ , there exists some $\zeta \in L^{p}(\Omega; \mathbb {R}^{d})$ with $\mathbb {P}_{\zeta }=\mu $ such that $\tilde {f}$ is the Fréchet differential at $\zeta $ , that is, there exists a linear continuous mapping $D\tilde {f}(\zeta ): L^{p}(\Omega; \mathbb {R}^{d})\to \mathbb {R}$ such that for all $\eta \in L^{p}(\Omega; \mathbb {R}^{d})$ ,

$$ \begin{align*} \tilde{f}(\zeta+\eta)-\tilde{f}(\zeta)=D\tilde{f}(\zeta)(\eta)+o(|\eta|_{L^{p}}),\quad |\eta|_{L^{p}}\to0. \end{align*} $$

Since $D\tilde {f}(\zeta )\in L(L^{p}(\Omega; \mathbb {R}^{d}),\mathbb {R})$ , by the Riesz representation theorem [Reference Brezis3], there exists a $\mathbb {P}$ -a.s. unique random variable $\mathcal {\theta }\in L^{p}(\Omega; \mathbb {R}^{d})$ such that for $\eta \in L^{p}(\Omega; \mathbb {R}^{d})$ ,

$$ \begin{align*} D\tilde{f}(\zeta)(\eta)=(\mathcal{\theta}, \eta)_{L^{p}}=\mathbb{E}[\mathcal{\theta}\cdot\eta]. \end{align*} $$

Thus, there exists a Borel measurable function $h: \mathbb {R}^{d}\to \mathbb {R}^{d}$ which depends on the distribution $\mathbb {P}_{\zeta }$ rather than $\zeta $ itself such that $\mathcal {\theta }=h(\zeta )$ , and for $\xi \in L^{2}(\Omega; \mathbb {R}^{d})$ ,

$$ \begin{align*} f(\mathbb{P}_{\xi})-f(\mathbb{P}_{\zeta})=\mathbb{E}[h(\zeta)(\xi-\zeta)]+o(|\xi-\zeta|_{L^{2}}). \end{align*} $$

We call $\partial _{\mu }f(\mathbb {P}_{\zeta })(y)\triangleq h(y), y\in \mathbb {R}^{d}$ as the derivative of $f: \mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d})\to \mathbb {R}$ at $\mathbb {P}_{\zeta }, \zeta \in L^{p}(\Omega; \mathbb {R}^{d})$ .

Definition 2.1. Function f is said to be in $C^{1}(\mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}))$ if for each $\xi \in L^{2}(\Omega; \mathbb {R}^{d})$ , there exists a $\mathbb {P}_{\xi }$ -modification of $\partial _{\mu }f(\mathbb {P}_{\zeta })(\cdot )$ which is denoted by $\partial _{\mu }f(\mathbb {P}_{\zeta })(\cdot )$ again, such that $\partial _{\mu }f: \mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d})\times \mathbb {R}^{d}\to \mathbb {R}^{d}$ is continuous and we identify the function $\partial _{\mu }f$ with the derivative of f.

Definition 2.2. A function f belongs to $C_{b}^{1,1}(\mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}))$ if $f\in C^{1}(\mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}))$ , and $\partial _{\mu }f$ is bounded and Lipschitz continuous, that is there exists a real number $C>0$ such that:

  1. (i) $\partial _{\mu }f(\mu )(x)|\leq C,\ \mu \in \mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d})$ ;

  2. (ii) $\partial _{\mu }f(\mu )(x)-\partial _{\mu }f(\nu )(y)|\leq C(\rho (\mu ,\nu )+|x-y|), \ \mu ,\nu \in \mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}) $

for $x,y\in \mathbb {R}^{d}.$

Definition 2.3. The function f is said to be in $C^{2}(\mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}))$ if for every $\mu \in \mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d})$ , $f\in C^{1}(\mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}))$ and $\partial _{\mu }f(\mathbb {P}_{\xi })(\cdot )$ is differentiable and its derivative $\partial _{y}\partial _{\mu }f: \mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d})\times \mathbb {R}^{d}\to \mathbb {R}^{d}\times \mathbb {R}^{d}$ is continuous.

Definition 2.4. The function f is said to be in $C_{b}^{2,1}(\mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}))$ if $f\in C^{2}(\mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}))\cap C_{b}^{1,1}(\mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}))$ and its derivative $\partial _{y}\partial _{\mu }f$ is bounded and continuous.

Definition 2.5. The function $\Phi \in C_{b}^{2,2,1}(\mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d})\times \mathbb {R}^{d})$ if:

  1. (i) $\Phi $ is bi-continuous with respect to $(x,\mu )$ ;

  2. (ii) for any x, $\Phi (x,\cdot )\in C_{b}^{2,1}(\mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}))$ and for any $\mu \in \mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d})$ , $\Phi (\cdot ,\mu )\in C^{2}(\mathbb {R}^{d})$ .

If $\Phi \in C_{b}^{2,2,1}(\mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}\times \mathbb {R}^{d}))$ and $\Phi \geq 0$ , then we say $\Phi \in C_{b+}^{2,2,1}(\mathbb {R}^{d}\times \mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}))$ .

Definition 2.6. The function $\Phi \in \mathcal {C}(\mathbb {R}^{d}\times \mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}))$ , if $\Phi \in C^{2,2}(\mathbb {R}^{d}\times \mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}))$ and for every compact set $K\subseteq \mathbb {R}^{d}\times \mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d})$ ,

$$ \begin{align*} \sup_{(x,\mu)\in K}\int_{\mathbb{R}^{d}}(\|\partial_{y}\partial_{\mu}\Phi(x,\mu)(y)\|^{2}+|\partial_{\mu}\Phi(x,\mu)(y)|^{2})\mu(dy)<\infty. \end{align*} $$

If $\Phi \in \mathcal {C}(\mathbb {R}^{d}\times \mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}))$ and $\Phi \geq 0$ , then we say that $\Phi \in \mathcal {C}_{+}(\mathbb {R}^{d}\times \mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d}))$ .

For (1.1), we make the following assumptions.

Assumption 2.7. Functions $b,\sigma $ are continuous with respect to $(x,\mu )\kern1.4pt{\in}\kern1.4pt \mathbb {R}^{d}\kern1.4pt{\times}\kern1.4pt \mathcal {M}_{\lambda ^{2}}(\mathbb {R}^{d})$ , and there is a constant $L_1>0$ such that

$$ \begin{align*} |b(x,\mu)|^{2}+\|\sigma(x,\mu)\|^{2}\leq L_{1}(1+|x|^{2}+\|\mu\|_{\lambda^{2}}^{2}). \end{align*} $$

Assumption 2.8. There is a constant $L_2>0$ such that

$$ \begin{align*} 2\langle x-y, b(x,\mu)-b(y,\mu)\rangle+\|\sigma(x,\mu)-\sigma(y,\nu)\|^{2}\leq L_{2}(|x-y|^{2}+\rho(\mu,\nu)^{2}). \end{align*} $$

Assumption 2.9. There exists a function $v(\cdot ,\cdot ): \mathbb {R}^{d}\times \mathcal {M}_{\lambda ^{p}}(\mathbb {R}^{d})\to \mathbb {R}$ such that:

  1. (i) $v\in \mathcal {C}_{+}(\mathbb {R}^{d}\times \mathcal {M}_{\lambda ^{2}}(\mathbb {R}^{d}))$ ;

  2. (ii) $\int _{\mathbb {R}^{d}}(L^{\mu }v(x,\mu )+\gamma v(x,\mu ))\mu (dx)\leq 0$ ;

  3. (iii) for some $p\geq 1$ , $a_{1}\int _{\mathbb {R}^{d}}|x|^{p}\mu (dx)\leq \int _{\mathbb {R}^{d}}v(x,\mu )\mu (dx)\leq a_{2}\int _{\mathbb {R}^{d}}|x|^{p}\mu (dx)$ .

By the classical result of Wang [Reference Wang19], under Assumptions 2.72.9, there exists a unique strong solution $X_t^{x_0}$ , with initial value $x_0$ , to (1.1), and for $p\geq 2$ , $\mathbb {E}\sup _{0\leq t\leq T}|X^{x_0}_{t}|^{p}<\infty $ . Let $C^{2}(\mathbb {R}^{d}\times \mathcal {M}_{\lambda ^{2}}(\mathbb {R}^{d});\mathbb {R}^{+})$ denote the space of nonnegative functions which are continuous and twice differentiable. For $V\in C^{2}(\mathbb {R}^{d}\times \mathcal {M}_{\lambda ^{2}}(\mathbb {R}^{d});\mathbb {R}^{+})$ , we have the following generator of (1.1):

$$ \begin{align*} L^{\mu}V(x,\mu)&=(b^{i}\partial_{x_{i}})(x,\mu)+\frac{1}{2}((\sigma\sigma^{*})^{ij}\partial^{2}_{x_{i}x_{j}})(x,\mu)+\int_{\mathbb{R}^{d}}b^{i}(y,\mu)(\partial_{\mu}V)_{i}(x,\mu)(y)\mu(dy)\nonumber\\&\quad+\frac{1}{2}\int_{\mathbb{R}^{d}}(\sigma\sigma^{*})^{ij}(y,\mu)\partial_{y_{i}}(\partial_{\mu}V)_{j}(y)\mu(dy)\,. \end{align*} $$

Remark 2.10. In fact, the strong solution to (1.1) defines a Markov process [Reference Wang19]. Let $p(t,x_{0},dz)$ be the transition probability distribution to process $X_{t}^{x_{0}}$ and $p(t,x_{0},\Gamma )$ be the probability for the event $\{X^{x_0}_{t}\in \Gamma \}$ with the initial value $x_{0}$ , that is,

$$ \begin{align*} p(t,x_{0},\Gamma)=\int_{\Gamma}p(t,x_{0},dz),\quad \Gamma\in\mathcal{B}(\mathbb{R}^{d}). \end{align*} $$

Definition 2.11. The process $X_{t}^{x_{0}}$ with initial value $x_{0}$ is called stable in distribution if there exists a probability measure $\Pi (\cdot )$ such that for any initial value $x_{0}$ , its transition probability $p(t,x_{0},dz)$ weakly converges to $\Pi (\cdot )$ , as $t\to \infty $ . Equation (1.1) is said to be stable in distribution if $X_{t}^{x_{0}}$ is stable in distribution.

To study the stability in distribution, we need the following assumption. First, for a given function $U\in C^{2}(\mathbb {R}^{d}; \mathbb {R})$ , we define the operator

$$ \begin{align*} L(x,y,\mu,\nu)U(x-y)&=[b(x,\mu)-b(y,\nu)]U_{x}(x-y)\nonumber\\ &\quad +\tfrac{1}{2}Tr[(\sigma(x,\mu)-\sigma(y,\nu))^{*}U_{xx}(x-y)(\sigma(x,\mu)-\sigma(y,\nu))]. \end{align*} $$

Assumption 2.12. There exist a function $U\in C^{2}(\mathbb {R}^{d};\mathbb {R}_{+})$ and a constant $K>0$ , such that for any two solutions $(X_{t}^{x_{0}})_{t\geq 0}$ and $(X_{t}^{y_{0}})_{t\geq 0}$ with its distributions $\mathcal {L}(X_{t}^{x_{0}})=\mu _{t}$ and $\mathcal {L}(X_{t}^{y_{0}})=\nu _{t}$ , and for all couplings $\pi \in \Pi (\mu _{t},\nu _{t})$ ,

$$ \begin{align*} \int_{\mathbb{R}^{d}\times\mathbb{R}^{d}}L(x,y,\mu_{t},\nu_{t})U(x-y)\pi(dx,dy)\leq -K\int_{\mathbb{R}^{d}\times\mathbb{R}^{d}}U(x-y)\pi(dx,dy). \end{align*} $$

3. Stability analysis

3.1. Exponential stability

This section gives the exponential stability in the pth moment and in a.s. sense.

Definition 3.1. Let $p\geq 1$ . The solution $X_{t}^{x_{0}}$ of (1.1) is said to be pth moment exponentially stable if there is a pair of constants $\gamma>0$ and $C>0$ such that

$$ \begin{align*} \mathbb{E}|X_{t}^{x_{0}}|^{p}\leq C|x_{0}|^{p}e^{-\gamma t},\quad t\geq0. \end{align*} $$

Further, it is said to be a.s. exponentially stable if

$$ \begin{align*} \limsup_{t\to\infty}\frac{\mathrm{\log}|X_{t}^{x_{0}}|}{t}\leq -\gamma,\quad \text{a.s}. \end{align*} $$

For this, we further assume that for some constant $N>0$ and $p\geq 1$ ,

(3.1) $$ \begin{align} |b(x,\mu)|^{p}+\|\sigma(x,\mu)\|^{p}\leq N\bigg(|x|^{p}+\int_{\mathbb{R}^{d}}|x|^{p}\mu(dx)\bigg). \end{align} $$

Theorem 3.2. Assume that Assumptions 2.72.9 hold. For every $x_{0}\in \mathbb {R}^{d}$ , $X_{t}^{x_{0}}$ is pth moment exponentially stable and a.s. exponentially stable. Furthermore, for every ${\epsilon>0}$ , there exists an $R>0$ such that for all $t\geq 0$ , $\mathbb {P}\{|X_{t}^{x_{0}}|\geq R\}<\epsilon $ .

Proof. For $x_{0}\in \mathbb {R}^{d}$ , let $X_{t}^{x_{0}}$ be the solution of (1.1), for positive integer k, define the stopping times

$$ \begin{align*} \rho_{k}=\inf\{t>0\mid |X_{t}^{x_{0}}|\geq k\}, \end{align*} $$

then obviously, $\rho _{k}\to \infty $ , a.s. as $k\to \infty $ . Then for the stopped processes $(X_{t\wedge \rho _{k}}^{x_{0}})_{t\geq 0}$ , function v and distribution processes $\mathcal {L}(X_{t}^{x_{0}})_{t\geq 0}$ , by It $\mathrm {\hat {o}}$ ’s formula [Reference Hammersley, Siska and Szpruch9],

$$ \begin{align*} e^{\gamma(t\wedge\rho_{k})}v(X_{t\wedge\rho_{k}}^{x_{0}},\ &\mathcal{L}(X_{t}^{x_{0}}))-v(x_{0},\delta_{x_{0}}) \\ &=\int_{0}^{t}\gamma e^{\gamma(s\wedge\rho_{k})}v(X_{s\wedge\rho_{k}}^{x_{0}},\mathcal{L}(X_{s}^{x_{0}}))\,ds \\ &\quad +\int_{0}^{t}e^{\gamma(s\wedge\rho_{k})}(b^{i}\partial_{x_{i}}v)(X_{s\wedge\rho_{k}}^{x_{0}},\mathcal{L}(X_{s}^{x_{0}}))\,ds\nonumber\\ &\quad+\frac{1}{2}\int_{0}^{t}e^{\gamma(s\wedge\rho_{k})}((\sigma\sigma^{*})^{ij}\partial_{x_{i}x_{j}}^{2}v)(X_{s\wedge\rho_{k}}^{x_{0}},\mathcal{L}(X_{s}^{x_{0}}))\,ds\nonumber\\ &\quad+\int_{0}^{t}e^{\gamma(s\wedge\rho_{k})}(\partial_{x_{i}}v\sigma^{ij})(X_{s\wedge\rho_{k}}^{x_{0}},\mathcal{L}(X_{s}^{x_{0}}))\,dW_{s}^{j}\nonumber\\ &\quad+\int_{0}^{t}\int_{\mathbb{R}^{d}}e^{\gamma(s\wedge\rho_{k})}((b^{i}\partial_{x_{i}}v)(X_{s\wedge\rho_{k}}^{x_{0}},\mathcal{L}(X_{s}^{x_{0}}))(y)\mathcal{L}(X_{s}^{x_{0}})(dy)\,ds\nonumber\\ &\quad+\frac{1}{2}\int_{0}^{t}\int_{\mathbb{R}^{d}}e^{\gamma(s\wedge\rho_{k})}(\sigma\sigma^{*})^{ij}(X_{s\wedge\rho_{k}}^{x_{0}},\mathcal{L}(X_{s}^{x_{0}}))(y)\mathcal{L}(X_{s}^{x_{0}})(dy)\,ds. \end{align*} $$

Then taking expectation, by Assumption 2.9,

$$ \begin{align*} \mathbb{E}[e^{\gamma(t\wedge\rho_{k})}v(X_{t\wedge\rho_{k}}^{x_{0}},\mathcal{L}(X_{t}^{x_{0}}))]-v(x_{0},\delta_{x_{0}})\leq0. \end{align*} $$

Let $k\to \infty $ and together with the Fatou lemma [Reference Ethier and Kurtz7],

$$ \begin{align*} \mathbb{E}[e^{\gamma t}v(X_{t},\mathcal{L}(X_{t}^{x_{0}}))]\leq v(x_{0},\delta_{x_{0}})\leq0. \end{align*} $$

Furthermore, by Assumption 2.7,

$$ \begin{align*} a_{1}\mathbb{E}|X_{t}^{x_{0}}|^{p}\leq \mathbb{E}v(X_{t},\mathcal{L}(X_{t}^{x_{0}}))\leq e^{-\gamma t}v(x_{0},\delta_{x_{0}})\leq a_{2}e^{-\gamma t}|x_{0}|^{p}. \end{align*} $$

Thus,

$$ \begin{align*} \mathbb{E}|X_{t}^{x_{0}}|^{p}\leq\frac{a_{2}}{a_{1}}e^{-\gamma t}|x_{0}|^{p}. \end{align*} $$

Then by Chebyshev’s inequality [Reference Ethier and Kurtz7], for $R>0$ ,

$$ \begin{align*} \mathbb{P}\{|X_{t}^{x_{0}}|\geq R\}\leq \frac{\mathbb{E}|X_{t}^{x_{0}}|^{p}}{R^{p}}. \end{align*} $$

Noticing that from (1.1),

$$ \begin{align*} X_{t+s}^{x_{0}}=X_{t}^{x_{0}}+\int_{t}^{t+s}b(X_{u}^{x_{0}},\mathcal{L}(X_{u}^{x_{0}}))\,du+\int_{t}^{t+s}\sigma(X_{u}^{x_{0}},\mathcal{L}(X_{u}^{x_{0}}))\,dW_{u}, \end{align*} $$

and for $p\geq 1$ , $\tau>0$ , by (3.1),

$$ \begin{align*} \mathbb{E}\sup_{0\leq s\leq \tau}|X_{t+s}^{x_{0}}|^{p}&\leq C_{p}\mathbb{E}|X_{t}^{x_{0}}|^{p}+C_{p}\mathbb{E}\sup_{0\leq s\leq \tau}\Big|\int_{t}^{t+s}b(X_{u}^{x_{0}},\mathcal{L}(X_{u}^{x_{0}}))\,du\Big|^{p}\nonumber\\&\quad +C_{p}\mathbb{E}\sup_{0\leq s\leq \tau}\Big|\int_{t}^{t+s}\sigma(X_{u}^{x_{0}},\mathcal{L}(X_{u}^{x_{0}}))\,dW_{u}\Big|^{p}\nonumber\\&\leq C_{p}\mathbb{E}|X_{t}^{x_{0}}|^{p}+C_{p,\tau}\mathbb{E}\int_{t}^{t+\tau}|b(X_{u}^{x_{0}},\mathcal{L}(X_{u}^{x_{0}}))|^{p}\,du \\&\quad +C_{p}\mathbb{E}\int_{t}^{t+\tau}|\sigma(X_{u}^{x_{0}},\mathcal{L}(X_{u}^{x_{0}}))|^{p}\,du\nonumber\\&\leq C_{p}\mathbb{E}|X_{t}^{x_{0}}|^{p} +C_{p,\tau,N}\mathbb{E}\int_{t}^{t+\tau}\bigg(|X_{u}^{x_{0}}|^{p}+\int_{\mathbb{R}^{d}}|x|^{p}\mathcal{L}(X_{u}^{x_{0}})(dx)\bigg)\,du\nonumber\\&\leq C_{p}\mathbb{E}|X_{t}^{x_{0}}|^{p}+C_{p,\tau,N}\int_{t}^{t+\tau}\mathbb{E}|X_{u}^{x_{0}}|^{p}\,du\nonumber\\&\leq C_{p,x_{0}}e^{-\gamma t}+C_{p,\tau,x_{0},N}\int_{t}^{t+\tau}e^{-\gamma u}\,du\leq C_{p,\tau,x_{0},\gamma,N}e^{-\gamma t}. \end{align*} $$

Then for $n=1,2, \dots $ ,

$$ \begin{align*} \mathbb{E}\sup_{n\tau\leq t\leq (n+1)\tau}|X_{t}^{x_{0}}|^{p}\leq C_{p,\tau,x_{0},\gamma,N}e^{-\gamma n\tau}, \end{align*} $$

thus, for $\epsilon \in (0,\gamma )$ and $n\in \mathbb {N}$ , by Chebychev’s inequality,

$$ \begin{align*} \mathbb{P}\bigg(\omega: \sup_{n\tau\leq t\leq (n+1)\tau}|X_{t}^{x_{0}}|^{p}>e^{-(\gamma-\epsilon)n\tau}\bigg)\leq C_{p,\tau,x_{0},\gamma,N}e^{-\epsilon n\tau}. \end{align*} $$

By the Borel–Cantelli lemma [Reference Kallenberg11], there exists a random constant $n_{0}(\omega )$ such that for almost all $\omega \in \Omega $ , for $n>n_{0}(\omega )$ ,

$$ \begin{align*} \sup_{n\tau\leq t\leq (n+1)\tau}|X_{t}^{x_{0}}|^{p}\leq C_{p,\tau,x_{0},\gamma,N}e^{-(\gamma-\epsilon) n\tau},\quad \text{a.s.} \end{align*} $$

Thus, for any $n\tau \leq t\leq (n+1)\tau $ ,

$$ \begin{align*} \frac{\mathrm{\log}|X_{t}^{x_{0}}|}{t}=\frac{\mathrm{\log}|X_{t}^{x_{0}}|^{p}}{pt}\leq \frac{\mathrm{\log}\sup_{n\tau\leq t\leq (n+1)\tau}|X_{t}^{x_{0}}|^{p}}{pn\tau}\leq -\frac{\gamma-\epsilon}{p},\quad\text{a.s.}, \end{align*} $$

and

$$ \begin{align*} \limsup_{t\to\infty}\frac{\mathrm{\log}|X_{t}^{x_{0}}|}{t}\leq -\frac{\gamma-\epsilon}{p},\quad \text{a.s.} \end{align*} $$

Now letting $\epsilon \to 0$ , the proof is complete.

Remark 3.3. From Theorem 3.2, the transition probability family $\{p(t,x_{0},dz)\mid t\geq 0\}$ is tight, that is, for $\epsilon>0$ , there exists a compact set $\mathcal {K}=\mathcal {K}(x_{0},\epsilon )$ such that

$$ \begin{align*} \mathbb{P}(t,x_{0},\mathcal{K})\geq 1-\epsilon. \end{align*} $$

3.2. Stability in distribution

Next, we consider the stability in distribution. For this, we need to consider the difference between two solutions with different initial values, that is,

$$ \begin{align*} X_{t}^{x_{0}}-X_{t}^{y_{0}}&=x_{0}-y_{0}+\int_{0}^{t}[b(X_{s}^{x_{0}},\mathcal{L}(X_{s}^{x_{0}}))-b(X_{s}^{y_{0}},\mathcal{L}(X_{s}^{y_{0}}))]\,ds\nonumber\\ &\quad +\int_{0}^{t}[\sigma(X_{s}^{x_{0}},\mathcal{L}(X_{s}^{x_{0}}))-\sigma(X_{s}^{y_{0}},\mathcal{L}(X_{s}^{y_{0}}))]\,dW_{s}. \end{align*} $$

We need two more notation. Let $\mathcal {H}$ be the set consisting of nondecreasing functions $K: \mathbb {R}_{+}\to \mathbb {R}_{+}$ such that $K(0)=0$ , and $\mathcal {H}_{\infty }$ be the set of functions $K\in \mathcal {H}$ such that $K(x)\to \infty $ as $x\to \infty $ .

Lemma 3.4. If there exists a function $U\in C^{2}(\mathbb {R}^{d};\mathbb {R}_{+})$ satisfying $U(0)=0$ and a function $\alpha _{1}\in \mathcal {H}_{\infty }$ such that

$$ \begin{align*} \alpha_{1}(|x|)\leq U(x) \quad \text{for} \;\; x\in \mathbb{R}^{d}, \end{align*} $$

then for every $\epsilon>0$ and compact set $\mathcal {K}$ on $\mathbb {R}^{d}$ , there exists $T=T(\mathcal {K},\epsilon )>0$ such that

$$ \begin{align*} \mathbb{P}\{\|X_{t}^{x_{0}}-X_{t}^{y_{0}}\|>\epsilon\}<1-\epsilon,\quad t\geq T, \; x_{0}, \;y_{0}\in\mathcal{K}. \end{align*} $$

For the convenience of presentation in the following, we rewrite (1.1) as

(3.2) $$ \begin{align} X_{t}^{x_{0}}=x_{0}+\int_{0}^{t}\tilde{b}(u,X_{u}^{x_{0}})\, du+\int_{0}^{t}\tilde{\sigma}(u,X_{u}^{x_{0}})\,dW_{u}, \quad t\geq0, \end{align} $$

where $\tilde {b}(u,X_{u}^{x_{0}})=b(X_{u}^{x_{0}},\mathcal {L}(X_{u}^{x_{0}}))$ and $\tilde {\sigma }(u,X_{u}^{x_{0}})=\sigma (X_{u}^{x_{0}},\mathcal {L}(X_{u}^{x_{0}}))$ . Since (1.1) has a unique strong solution $(X_{t}^{x_{0}})_{t\geq 0}$ with the initial distribution $\delta _{x_{0}}$ , the distribution of process $(X_{t}^{x_{0}})_{t\geq 0}$ is known, and (3.2) is a classic SDE. Then, the solution of (3.2) is a strong Markov process [Reference Ding and Qiao5, Lemma 5.3].

Proof of Lemma 3.4.

For $\epsilon>0$ , by the continuity of function U with $U(0)=0$ , we choose $\alpha \in (0,\epsilon )$ small enough such that

$$ \begin{align*} \frac{\sup_{|x|\leq \alpha}U(x)}{\mu_{1}(\epsilon)}<\frac{\epsilon}{2}. \end{align*} $$

Let $\mathcal {K}$ be a compact set on $\mathbb {R}^{d}$ and for fixed $x_{0}$ , $y_{0}\in \mathcal {K}$ and $\beta>\alpha $ , we define two stopping times

$$ \begin{align*} \tau_{\alpha}&=\inf\{t\geq0\mid |X_{t}^{x_{0}}-X_{t}^{y_{0}}|\leq \alpha\}, \\[.1cm] \tau_{\beta}&=\inf\{t\geq0\mid |X_{t}^{x_{0}}-X_{t}^{y_{0}}|\geq \beta\}. \end{align*} $$

By It $\mathrm {\hat {o}}$ ’s formula for the stopped process $U(X_{\tau _{\beta }\wedge t}^{x_{0}}-X_{\tau _{\beta }\wedge t}^{y_{0}})$ and Assumption 2.12,

$$ \begin{align*} \mathbb{E}U(X_{\tau_{\beta}\wedge t}^{x_{0}}-X_{\tau_{\beta}\wedge t}^{y_{0}})&\leq U(x_{0}-y_{0})-K\int_{0}^{\tau_{\beta}\wedge t}U(X_{s}^{x_{0}}-X_{s}^{y_{0}})\,ds\nonumber\\ &\quad +\mathbb{E}\int_{0}^{\tau_{\beta}\wedge t}U_{x}(X_{s}^{x_{0}}-X_{s}^{y_{0}})(\tilde{\sigma}(s,X_{s}^{x_{0}})-\tilde{\sigma}(u,X_{s}^{y_{0}}))\,dW_{s}\nonumber\\ &=U(x_{0}-y_{0})-K\int_{0}^{\tau_{\beta}\wedge t}U(X_{s}^{x_{0}}-X_{s}^{y_{0}})\,ds. \end{align*} $$

Then,

$$ \begin{align*} \alpha_{1}(\beta)\mathbb{P}\{\tau_{\beta}\leq t\}\leq U(x_{0}-y_{0}), \end{align*} $$

that is,

$$ \begin{align*} \mathbb{P}\{\tau_{\beta}\leq t\}\leq \frac{U(x_{0}-y_{0})}{\alpha_{1}(\beta)}. \end{align*} $$

Notice that for all $x_{0}, y_{0}\in \mathcal {K}$ and $U(x_{0}-y_{0})$ bounded, there exists $\beta =\beta (\mathcal {K},\epsilon )>0$ such that

$$ \begin{align*} \mathbb{P}\{\tau_{\beta}<\infty\}\leq \frac{\epsilon}{4}. \end{align*} $$

Fix the $\beta $ and let $t_{\alpha }=\tau _{\alpha }\wedge \tau _{\beta }\wedge t$ , then similar discussion yields

$$ \begin{align*} \mathbb{E}U(X_{t_{\alpha}}^{x_{0}}-X_{t_{\alpha}}^{y_{0}})&\leq U(x_{0}-y_{0})-K\mathbb{E}\int_{0}^{t_{\alpha}}U(X_{s}^{x_{0}}-X_{s}^{y_{0}})\,ds\nonumber\\ &\leq U(x_{0}-y_{0})-K\mathbb{E}\int_{0}^{t_{\alpha}}\alpha_{1}(|X_{s}^{x_{0}}-X_{s}^{y_{0}}|)\,ds\nonumber\\ &\leq U(x_{0}-y_{0})-K\alpha_{1}(\alpha)\mathbb{E}(\tau_{\alpha}\wedge\tau_{\beta}\wedge t), \end{align*} $$

which implies that

$$ \begin{align*} \mathbb{P}\{\tau_{\alpha}\wedge\tau_{\beta}\geq t\}\leq \frac{U(x_{0}-y_{0})}{K\alpha_{1}(\alpha)t}. \end{align*} $$

Moreover, this implies that for a given $\epsilon \in (0,1)$ , there exists $T=T(\mathcal {K},\epsilon )>0$ such that

$$ \begin{align*} \mathbb{P}\{\tau_{\alpha}\wedge\tau_{\beta}\leq T\}>1-\frac{\epsilon}{4}. \end{align*} $$

Thus,

$$ \begin{align*} 1-\frac{\epsilon}{4}<\mathbb{P}\{\tau_{\alpha}\wedge\tau_{\beta}\leq T\}&\leq \mathbb{P}\{\tau_{\alpha}\leq T\}+\mathbb{P}(\tau_{\beta}\leq T)\nonumber\\ &\leq \mathbb{P}\{\tau_{\alpha}\leq T\}+\mathbb{P}(\tau_{\beta}\leq\infty)\nonumber\\ &\leq \mathbb{P}\{\tau_{\alpha}\leq T\}+\frac{\epsilon}{4} \end{align*} $$

and

$$ \begin{align*} \mathbb{P}\{\tau_{\alpha}\leq T\}\geq 1-\frac{\epsilon}{2}. \end{align*} $$

Now we define the stopping time

$$ \begin{align*} \sigma=\inf\{t\geq\tau_{\alpha}\wedge T\mid |X_{t}^{x_{0}}-X_{t}^{y_{0}}|\geq\epsilon\}. \end{align*} $$

Let $t>T$ , then

$$ \begin{align*} \mathbb{P}({\tau_{\alpha}\leq T}\cap{\sigma\leq t})\mu_{1}(\epsilon)&\leq\mathbb{E}(I_{\tau_{\alpha}\leq T,\sigma\leq t}U(X_{t\wedge\sigma}^{x_{0}}-X_{t\wedge\sigma}^{y_{0}}))\nonumber\\ &\leq\mathbb{E}(I_{\tau_{\alpha}\leq T}U(X_{t\wedge\tau_{\alpha}}^{x_{0}}-X_{t\wedge\tau_{\alpha}}^{y_{0}}))\nonumber\\ &\leq\mathbb{E}(I_{\tau_{\alpha}\leq T}U(X_{\tau_{\alpha}}^{x_{0}}-X_{\tau_{\alpha}}^{y_{0}}))\leq\mathbb{P}({\tau_{\alpha}\leq T})\sup_{|x|\leq\alpha}U(x). \end{align*} $$

Thus,

$$ \begin{align*} \mathbb{P}(\{\tau_{\alpha}\leq T\}\cap\{\sigma\leq t\})\leq\frac{\epsilon}{2} \end{align*} $$

and

$$ \begin{align*} \mathbb{P}\{\sigma\leq t\}\leq \mathbb{P}(\{\tau_{\alpha}\leq T\}\cap\{\sigma\leq t\})+\mathbb{P}\{\tau_{\alpha>T}\}<\epsilon. \end{align*} $$

Let $t\to \infty $ , then

$$ \begin{align*} \mathbb{P}\{\sigma<\infty\}\leq \epsilon. \end{align*} $$

This indicates that for all $x_{0}$ , $y_{0}\in \mathcal {K}$ , $t\geq T$ ,

$$ \begin{align*} \mathbb{P}\{|X_{t}^{x_{0}}-X_{t}^{y_{0}}|<\epsilon\}\geq1-\epsilon. \end{align*} $$

This completes the proof.

Lemma 3.5. For every compact set $\mathcal {K}$ ,

$$ \begin{align*} \lim_{t\to\infty}\rho(p(t,x_{0},\cdot),p(t,y_{0},\cdot))=0,\;\; x_{0},\; y_{0}\in\mathcal{K}. \end{align*} $$

Proof. We only need to show that there exists $T>0$ , such that for all $\epsilon>0$ and $x_{0}$ , $y_{0}\in \mathcal {K}$ ,

$$ \begin{align*} \rho(p(t,x_{0},\cdot),p(t,y_{0},\cdot))\leq \epsilon,\quad t\geq T. \end{align*} $$

It is equivalent to show that

$$ \begin{align*} \sup_{\phi\in C_{\lambda}(\mathbb{R}^{d})}|\mathbb{E}\phi(X_{t}^{x_{0}})-\mathbb{E}\phi(X_{t}^{y_{0}})|\leq \epsilon,\quad t\geq T. \end{align*} $$

However, notice that for every $\phi \in C_{\lambda }(\mathbb {R}^{d})$ ,

$$ \begin{align*} |\mathbb{E}\phi(X_{t}^{x_{0}})-\mathbb{E}\phi(X_{t}^{y_{0}})|\leq\mathbb{E}[2\wedge|X_{t}^{x_{0}}-X_{t}^{y_{0}}|]. \end{align*} $$

By Lemma 3.4, there exists a $T_{1}>0$ such that

$$ \begin{align*} \mathbb{E}[2\wedge|X_{t}^{x_{0}}-X_{t}^{y_{0}}|]<\epsilon,\quad t\geq T_{1}. \end{align*} $$

Due to the arbitrariness of $\phi \in C_{\lambda }(\mathbb {R}^{d})$ ,

$$ \begin{align*} \sup_{\phi\in C_{\lambda}(\mathbb{R}^{d})}|\mathbb{E}\phi(X_{t}^{x_{0}})-\mathbb{E}\phi(X_{t}^{y_{0}})|\leq \epsilon,\quad t\geq T_{1}. \end{align*} $$

The proof is now complete.

Lemma 3.6. Under the assumptions of Theorem 3.2 and Lemma 3.4 , for $x_0\in \mathbb {R}^d$ , $\{p(t,x_{0},\cdot )\,, t\geq 0\}$ is a Cauchy sequence.

Proof. We need to show that for every $x_{0}\in \mathbb {R}^{d}$ and $\epsilon>0$ , there exists $T>0$ such that for $t\geq T$ and $s>0$ ,

$$ \begin{align*} \rho(p(t+s,x_{0},\cdot),p(t,x_{0},\cdot))\leq \epsilon\,, \end{align*} $$

which is equivalent that for every $\phi \in C_{\lambda }(\mathbb {R}^{d})$ ,

$$ \begin{align*} \sup_{\phi\in C_{\lambda}(\mathbb{R}^{d})}|\mathbb{E}\phi(X_{t+s}^{x_{0}})-\mathbb{E}\phi(X_{t}^{x_{0}})|\leq \epsilon,\quad t\geq T, \;s>0. \end{align*} $$

By Lemma 3.2, there exists a compact set $\mathcal {K}$ on $\mathbb {R}^{d}$ such that for $\epsilon>0$ ,

$$ \begin{align*} p(t,x_{0},\mathcal{K})>1-\frac{\epsilon}{8}. \end{align*} $$

Furthermore, by the strong Markov property of $X_{t}^{x_{0}}$ , for $\phi \in C_{\lambda }(\mathbb {R}^{d})$ and t, $s>0$ ,

$$ \begin{align*} |\mathbb{E}\phi(X_{t+s}^{x_{0}})-\mathbb{E}\phi(X_{t}^{x_{0}})|&=|\mathbb{E}[\mathbb{E}(\phi(X_{t+s}^{x_{0}})|\mathcal{F}_{s})]-\mathbb{E}\phi(X_{t}^{x_{0}})|\nonumber\\ &=|\mathbb{E}[\mathbb{E}\phi(X_{t}^{s,X_{s}^{x_{0}}})]-\mathbb{E}\phi(X_{t}^{x_{0}})|\nonumber\\ &=\Big|\!\int_{\mathbb{R}^{d}}\mathbb{E}\phi(X_{t}^{z})p(s,x_{0},dz)-\mathbb{E}\phi(X_{t}^{x_{0}})\Big|\nonumber\\ &\leq\int_{\mathbb{R}^{d}}\mathbb{E}|\phi(X_{t}^{z})-\phi(X_{t}^{x_{0}})|p(s,x_{0},dz)\nonumber\\ &=\int_{\mathcal{K}}\mathbb{E}|\phi(X_{t}^{z})-\phi(X_{t}^{x_{0}})|p(s,x_{0},dz)\nonumber\\ &\quad+\int_{\mathbb{R}^{d}-\mathcal{K}}\mathbb{E}|\phi(X_{t}^{z})-\phi(X_{t}^{x_{0}})|p(s,x_{0},dz)\nonumber\\ &\leq\int_{\mathcal{K}}\mathbb{E}|\phi(X_{t}^{z})-\phi(X_{t}^{x_{0}})|p(s,x_{0},dz)+\frac{\epsilon}{4}. \end{align*} $$

By Lemma 3.4, there exists $T>0$ such that for every $\epsilon>0$ ,

$$ \begin{align*} \mathbb{E}|\phi(X_{t}^{z})-\phi(X_{t}^{x_{0}})|\leq \mathbb{E}[2\wedge|X_{t}^{z}-X_{t}^{x_{0}}|]\leq \frac{3\epsilon}{4}, \;\;t\geq T. \end{align*} $$

Thus,

$$ \begin{align*} |\mathbb{E}\phi(X_{t+s}^{x_{0}})-\mathbb{E}\phi(X_{t}^{x_{0}})|\leq \epsilon,\quad t\geq T,\; s>0, \end{align*} $$

which completes the proof.

Theorem 3.7. Under assumptions of Theorem 3.2 and Lemma 3.4 , (1.1) is stable in distribution.

Proof. By Definition 2.11, we need to show that there exists a probability measure $\pi (\cdot )$ such that for every $x_{0}\in \mathbb {R}^{d}$ , the transition probability family $\{p(t,x_{0},\cdot ): t\geq 0\}$ weakly converges to $\pi (\cdot )$ . In fact, we show that for every $x\in \mathbb {R}^{d}$ ,

$$ \begin{align*} \lim_{t\to\infty}\rho(p(t,x_{0},\cdot),\pi(\cdot))=0. \end{align*} $$

By Lemma 3.6, $\{p(t,0,\cdot ): t\geq 0\}$ is a Cauchy sequence in $\mathcal {P}(\mathbb {R}^{d})$ with the metric $\rho $ . Since $\mathcal {P}(\mathbb {R}^{d})$ is a complete metric space, there exists a probability measure $\pi (\cdot )\in \mathcal {P}(\mathbb {R}^{d})$ such that

$$ \begin{align*} \lim_{t\to\infty}\rho(p(t,x,\cdot),\pi(\cdot))=0. \end{align*} $$

By a triangle inequality,

$$ \begin{align*} \lim_{t\to\infty}\rho(p(t,x,\cdot),\pi(\cdot))&\leq \lim_{t\to\infty}\rho(p(t,x,\cdot),p(t,0,\cdot))+\lim_{t\to\infty}\rho(p(t,0,\cdot),\pi(\cdot)) \\ &=0. \end{align*} $$

The proof is complete.

4. Conclusions

In the study, we first prove the pth moment exponential stability and a.s. exponential stability of the solution of (1.1) by using the distribution-dependent It $\mathrm {\hat {o}}$ 's formula, and then obtain the tightness of the transition probability family corresponding to the solution of (1.1). Based on this, we introduce a distribution-dependent operator, that is, Assumption 2.12, and combined with the method of Yuan and Mao [Reference Yuan and Mao21], we get that when the time is long enough, the transition probability family tends to a unique probability measure, that is, the solution of (1.1) is asymptotically stable in distribution. It would be valuable to use a similar method to analyse the long time behaviour of (1.1) with jump noise.

Acknowledgements

This research is supported by the Natural Science Foundation of Jiangsu Province (BK20230899) and the National Natural Science Foundation of China (11771207).

References

Bahlali, K., Mezerdi, M. A. and Mezerdi, B., “Stability of McKean–Vlasov stochastic differential equations and applications”, Stoch. Dyn. 20 (2020) Article ID: 2050007; doi:10.1142/S0219493720500070.CrossRefGoogle Scholar
Bao, J. H., Hou, Z. T. and Yuan, C. G., “Stability in distribution of neutral stochastic differential delay equations with Markovian switching”, Stat. Probab. Lett. 79 (2009) 16631673; doi:10.1016/j.spl.2009.04.006.CrossRefGoogle Scholar
Brezis, H., Functional analysis, Sobolev spaces and partial differential equations (Springer, New York, 2010).Google Scholar
Cardaliaguet, P., Notes on mean field games (from P. L. Lion’s lectures at College de France). https://www.ceremade.dauphine.fr/cardalia/MFG100629.pdf.Google Scholar
Ding, X. J. and Qiao, H. J., “Stability for stochastic McKean–Vlasov equations with non-Lipschitz coefficients”, SIAM J. Control Optim. 59 (2021) Article ID: 19M1289418; doi:10.1137/19M1289418.CrossRefGoogle Scholar
Du, N. H., Dang, N. H. and Dieu, N. T., “On stability in distribution of stochastic differential delay equations with Markovian switching”, Syst. Control Lett. 65 (2014) 4349; doi:10.1016/j.sysconle.2013.12.006.CrossRefGoogle Scholar
Ethier, S. N. and Kurtz, T. G., Markov processes: characterization and convergence (John Wiley & Sons, NJ, 1986).CrossRefGoogle Scholar
Fei, C., Fei, W. Y., Deng, S. N. and Mao, X. R., “Asymptotic stability in distribution of highly nonlinear stochastic differential equations with G-Brownian motion”, Qual. Theory Dyn. Syst. 22 (2023) Article ID: 57; doi:10.1007/s12346-023-00760-9.CrossRefGoogle Scholar
Hammersley, W. R. P., Siska, D. and Szpruch, L., “McKean–Vlasov SDEs under measure dependent Lyapunov conditions”, Ann. Inst. Henri Poincaré Probab. Stat. 57 (2021) 10321057; doi:10.1214/20-AIHP1106.CrossRefGoogle Scholar
Hong, W., Li, S. H. and Liu, W., “Large deviation principle for McKean–Vlasov quasilinear stochastic evolution equations”, Appl. Math. Optim. 84 (2021) 11191147; doi:10.1007/s00245-021-09796-2.CrossRefGoogle Scholar
Kallenberg, O., Foundations of modern probability (Springer, New York, 2000).Google Scholar
Khasminskii, R., Stochastic stability of differential equations (Springer, Berlin–Heidelberg, 2011).Google Scholar
Liu, Z. X. and Ma, J., “Existence, uniqueness and ergodicity for McKean–Vlasov SDEs under distribution-dependent Lyapunov conditions”, Preprint, 2023, arXiv:2309.05411.CrossRefGoogle Scholar
Lv, G. Y. and Shan, Y. Q., “Long time behavior of stochastic McKean–Vlasov equations”, Appl. Math. Lett. 128 (2022) Article ID: 107879; doi:10.1016/j.aml.2021.107879.CrossRefGoogle Scholar
McKean, H. P., “A class of Markov processes associated with nonlinear parabolic equations”, Proc. Natl. Acad. Sci. USA 56 (1966) 19071911; https://www.jstor.org/stable/57643.CrossRefGoogle ScholarPubMed
Peng, S. G., “G–Brownian motion and dynamic risk measure under volatility uncertainty”, Preprint, 2007, arXiv:0711.2834.Google Scholar
Reis, D., Salkeld, G. and Tugaut, W., “Freidlin–Wentzell LDPs in path space for McKean–Vlasov equations and the functional iterated logarithm law”, Ann. Appl. Probab. 29 (2019) 14871540; https://www.jstor.org/stable/26729307.Google Scholar
Vlasov, A. A., “The vibrational properties of an electron gas”, Sov. Phys. 10 (1968) 721733; doi:10.1070/PU1968v010n06ABEH003709.CrossRefGoogle Scholar
Wang, F. Y., “Distribution dependent SDEs for Landau type equations”, Stoch. Process. Appl. 128 (2018) 596621; doi:10.1016/j.spa.2017.05.006.CrossRefGoogle Scholar
Wu, H., Hu, J. H., Gao, S. B. and Yuan, C. G., “Stabilization of stochastic McKean–Vlasov equations with feedback control based on discrete-time state observation”, SIAM J. Control Optim. 60 (2022) 28842901; doi:10.1137/21M1454997.CrossRefGoogle Scholar
Yuan, C. G. and Mao, X. R., “Asymptotic stability in distribution of stochastic differential equations with Markovian switching”, Stoch. Process. Appl. 103 (2003) 277291; doi:10.1016/S0304-4149(02)00230-2.CrossRefGoogle Scholar
Yuan, C. G., Zou, J. Z. and Mao, X. R., “Stability in distribution of stochastic differential delay equations with Markovian switching”, Syst. Control Lett. 50 (2003) 195207; doi:10.1016/S0167-6911(03)00154-3.CrossRefGoogle Scholar