\documentclass[10pt,a4paper]{article}
\usepackage[utf8]{inputenc}
\usepackage{amsmath}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{dsfont}
\usepackage{fancyhdr}
\usepackage{indentfirst}
\usepackage{graphicx}
\usepackage{newlfont}
\usepackage{amssymb}
\usepackage{amsmath}
\usepackage{latexsym}
\usepackage{amsthm}
\usepackage{mathtools}
\usepackage{nicefrac}
\usepackage{epstopdf}
\usepackage{caption}
\usepackage{subcaption}
\usepackage{color}
\usepackage{dsfont}
\usepackage{comment}
\usepackage[left=2cm,right=2cm,top=2cm,bottom=2cm]{geometry}
\usepackage{hyperref}
%\hypersetup{
%    colorlinks=true,
%    linkcolor=blue,
%    filecolor=magenta,      
%    urlcolor=cyan,
%    pdftitle={Overleaf Example},
%    pdfpagemode=FullScreen,
%    }


\newtheorem{theorem}{Theorem}
\newtheorem{proposition}[theorem]{Proposition}
\newtheorem{corollary}[theorem]{Corollary}
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{definition}[theorem]{Definition}
\newtheorem{assumption}[theorem]{Assumption}


\newtheorem{remark}[theorem]{Remark}
\newtheorem{example}[theorem]{Example}
\newtheorem{notation}{Notation}

\def \R {\mathbb{R}}
\def \N {\mathbb{N}}
\def \d {\mathrm{d}}
\newcommand{\red}[1]{\textcolor{red}{#1}}

\begin{document}
\section{Introduction}
TBD
\section{Preliminaries}
\begin{definition}[$bC^\alpha$ space]
$bC^\alpha$ is the space of bounded $\alpha$-Holder functions normed by
$$
||f||_{bC^\alpha}=\sup_x|f(x)|+\sup_{x\neq y}\frac{|f(x)-f(y)|}{|x-y|^\alpha}.
$$
\end{definition}
\begin{definition}[$d_{bC^\alpha}$ on $\mathcal{P}(\R^N)$]
We define on $\mathcal{P}(\R^N)$ the distance
$$
d_{bC^\alpha}(\mu,\nu)=\sup_{||f||_{bC^\alpha}\leq 1}\int f(x)\left(\mu(dx)-\nu(dx)\right).
$$
\end{definition}
\begin{definition}[MKV SDE]
A MKV SDE is an SDE of the type
$$
dX_t = B_t(X_t,[X_t])dt + \Sigma_t(X_t,[X_t])dW_t,
$$
where $W_t$ is a $N$-dimensional Brownian Motion, the coefficients are functions $B:\Omega\times[0,T]\times\R^N\times\mathcal{P}(\R^N)\rightarrow\R^N$ and $\Sigma:\Omega\times[0,T]\times\R^N\times\mathcal{P}(\R^{N})\rightarrow\R^{N\times N}$ and lastly $[X_t]$ is the distribution of $X_t$.
\end{definition}
\begin{assumption}[Structural Assumption]\label{a1}
Consider a MKV SDE where the coefficients are defined as
$$
B(t,x,\mu)=\int_{\R^N}b(t,x,y)\mu(dy),\qquad \Sigma(t,x,\mu)=\int_{\R^N}\sigma(t,x,y)\mu(dy),
$$
where $b$ and $\sigma$ are deterministic functions of $b:[0,T]\times\R^N\times\R^N\rightarrow\R^N$ and $\sigma:[0,T]\times\R^N\times\R^N\rightarrow\R^{N\times N}$.
In particular the coefficients are deterministic and the dependence on the distribution is fixed and integral-like.
\end{assumption}
\begin{assumption}[Non degeneracy]\label{a2}
The diffusion coefficient $\Sigma$ defines a uniformly parabolic operator, precisely
$$
\exists \lambda>0\ s.t.\ \frac{1}{\lambda}|v|^2\leq \left\langle \Sigma^*(t,x,\mu)\Sigma(t,x,\mu) v, v \right\rangle\leq \lambda |v|^2,\ \forall v\in \R^N, \forall (t,x,\mu)\in [0,T]\times\R^N\times\mathcal{P}(\R^N),
$$
a sufficient condition for this to occur under assupmtion \ref{a1} is if $\sigma$ is a uniformly positive definite matrix.
\end{assumption}
\begin{assumption}[Regularity]\label{a3}
Assume assumption \ref{a1}. The coefficients $b$ and $\sigma$ are functions in $L^\infty([0,T],bC^\alpha)$. Explicitly
\begin{gather*}
\exists C>0,\ s.t.\ \forall (x_1,x_2,y_1,y_2)\in \R^{4N},\ t-\mathrm{a.s.},\qquad |b(t,x,y)|+|\sigma(t,x,y)|\leq C,\\
\mathrm{and}\ s.t.\ |b(t,x_1,y_1)-b(t,x_2,y_2)| + |\sigma(t,x_1,y_1)-\sigma(t,x_2,y_2)|\leq C\left(|x_1-x_2|^\alpha+|y_1-y_2|^\alpha\right).
\end{gather*}
\end{assumption}
\section{Results}
\begin{theorem}[Non-degenerate case]\label{t1}
Consider a MKV SDE under assumptions \ref{a1}, \ref{a2} and \ref{a3} and initial distribution $\mu_0\in\mathcal{P}(\R^N)$; then we have weak existance and uniqueness of the solution of the MKV SDE.
\end{theorem}
\section{Proofs}
\begin{proof}of Theorem \ref{t1}.

The core of the proof is a contraction argument on the space of the flows of marginals, indeed if it is possible to show that the flow of marginals is unique it can be fixed in the MKV SDE to reduce it to a classical SDE, at this point weak well posedness results for classical SDEs will conclude the proof.

To begin if we consider a flow of marginals $(\mu_t)_{t\in[0,T]}\in C([0,T],\mathcal{P}(\R^N))$ and fix it inside the coefficients we are able to construct "linearized" coefficients and a classical SDE
\begin{gather*}
B^\mu(t,x)=B(t,x,\mu_t),\qquad \Sigma^\mu(t,x)=\Sigma(t,x,\mu_t),\\
dX^\mu_t = B^\mu(t,X^\mu_t)dt + \Sigma^\mu(t,X^\mu_t)dW_t\qquad X^\mu_0\sim\mu_0.
\end{gather*}
We first notice that the function $\mu\rightarrow B(t,x,\mu)\in C(\mathcal{P}(\R^N),\R^N)$ since
$$
|B(t,x,\mu)-B(t,x,\nu)|=\left| \int b(t,x,y)\left(\mu(dy)-\nu(dy)\right) \right|\leq ||b||_{bC^\alpha}d_{bC^\alpha}(\mu,\nu).
$$
Indeed the bounded $\alpha$-Holder distance metrizes the weak convergence, the proof will be postponed to the Appendix as theorem \ref{ap1}. Thus since as a function of $t$ $B(t,x,\mu)\in L^\infty$ we have that $B^\mu\in L^\infty([0,T],bC^\alpha)$. The same is true for $\Sigma^\mu$. For this reason due to Theorem 18.2.3 and 18.2.6 of \cite{Pascu_PB2} we have weak well posedness of the "linearized" SDE (with solution $X^\mu$), existance of a fundamental solution and uniform gaussian estimates. In particular $p^\mu(s,x;t,y)$ will be the fundamental solution of
\begin{align*}
(\partial_s+\mathcal{A}^\mu_s)p^\mu(s,x;t,y)&=0,\\
(\partial_t-(\mathcal{A}^\mu_t)^*)p^\mu(s,x;t,y)&=0,
\end{align*}
where $\mathcal{A}^\mu_s$ is the infinitesimal generator of the "linearized" SDE:
\begin{gather*}
\mathcal{A}^\mu_s =\frac{1}{2}\sum_{i,j=1}^N c^\mu_{i,j}(t,x)\partial_{x_ix_j} + \sum_{i=1}^NB^\mu_i(t,x)\partial_{x_i},\\
\mathrm{with}\qquad (c^\mu_{i,j})_{i,j=1,\cdots N}=\Sigma^*\Sigma.
\end{gather*}
At this point we may define respectively the forward and the backward translation operators via $(\mu_t)_{t\in[0,T]}$:
\begin{align*}
U^{t,s}_\mu u(y) = \int p^\mu(s,x;t,y)u(dx),\qquad u\in\mathcal{P}(\R^N),\\
V^{s,t}_\mu u(x) = \int p^\mu(s,x;t,y)u(dy),\qquad u\in\mathcal{P}(\R^N),
\end{align*}
and observe that the distribution of $X^\mu_t$ is absolutely continuous with density
$$
u^\mu_t(y)=U^{t,0}_\mu\mu_0.
$$
Thus we may define the application $\mathcal{L}:C([0,T], \mathcal{P}(\R^N))\rightarrow C([0,T], \mathcal{P}(\R^N))$ such that $\mathcal{L}((\mu_t)_{t\in[0,T]})=(\mathcal{L}^\mu_t)_{t\in[0,T]}=([X^\mu_t])_{t\in[0,T]}$. Less abstractly this is the application that given a flow of marginals returns the flow of marginals of the solution of the "linearized" SDE with the first flow of marginals and initial datum $\mu_0$.
\begin{remark}
Since $\mathcal{L}^\mu_0=[X^\mu_0]=\mu_0$ for any flow of marginals the image through $\mathcal{L}$ has initial law equal to $\mu_0$. 
%Thus we image of $\mathcal{L}$ is actually contained inside $$ \left\lbrace (\nu_t)_{t\in[0,T]}\ |\ \nu_0=\mu_0 \right\rbrace.$$
\end{remark}
Now we need to observe carefully the definition of $d_{bC^\alpha}([X^\mu_t],[X^\nu_t])$: if we remove the $\sup$ we get
$$
I(f)=\int f(x)\left( [X^\mu_t](dx)-[X^\nu_t](dx) \right)=\int f(x)\left( u^\mu_t(x)-u^\nu_t(x) \right)dx,
$$
which may be seen as an indicator of closeness between the two distributions. This is useful because by unraveling the definition of the densities we are able to exchange the order of integration to use the backward translation operators instead:
\begin{align*}
I(f) &= \int\int f(x)(p^\mu-p^\nu)(0,y;t,x)\mu_0(dy)dx = \int\int f(x)(p^\mu-p^\nu)(0,y;t,x)dx\ \mu_0(dy)\\
&=\int\int\int_0^t\frac{d}{ds}\left(\int p^\mu(0,y;s,z)p^\nu(s,z;t,x)f(x)dz\right)ds\ dx\ \mu_0(dy),
\intertext{Since the fundamental solutions are $C^{1,2}$ in the time and space variables locally around $(s,z)$ we may exchange the integral and the derivative}
&=\int\int\int_0^t\int\partial_{t_2}p^\mu(0,y;s,z)p^\nu(s,z;t,x)f(x)dz +\int p^\mu(0,y;s,z)\partial_{t_1}p^\nu(s,z;t,x)f(x)dz\ ds\ dx\ \mu_0(dy),
\intertext{by the fact that $p$ is the fundamental solution for both the forward and backward PDEs}
&=\int\int\int_0^t\int(\mathcal{A}_s^\mu)^* p^\mu(0,y;s,z)p^\nu(s,z;t,x)f(x)dz - \int p^\mu(0,y;s,z)\mathcal{A}_s^\nu p^\nu(s,z;t,x)f(x)dz\  ds\  dx\  \mu_0(dy)\\
&=\int \int_0^t \int(\mathcal{A}_s^\mu)^* p^\mu(0,y;s,z)V_\nu^{s,t}f(z)dz - \int (\mathcal{A}_s^\nu)^*p^\mu(0,y;s,z) V_\nu^{s,t}f(z)dz\  ds\  \mu_0(dy)\\
&=\int \int_0^t \int p^\mu(0,y;s,z)\left(\mathcal{A}_s^\mu - \mathcal{A}_s^\nu\right)V_\nu^{s,t}f(z)dz\  ds\ \mu_0(dy)\\
&=\int\int_0^t V_\mu^{0,s}\left(\mathcal{A}_s^\mu - \mathcal{A}_s^\nu\right)V_\nu^{s,t}f(y) ds\ \mu_0(dy)=\int_0^t \int V_\mu^{0,s}\left(\mathcal{A}_s^\mu - \mathcal{A}_s^\nu\right)V_\nu^{s,t}f(y) \mu_0(dy) ds.
\end{align*}
\red{Make rigorous}: Which in some sense means that instead of transporting forward the initial law via different flows and testing against $f$ we are transporting backward the test function and then testing the difference against the initial law.

At this point it's easy to see that since $\mu_0$ is a probability distribution
\begin{align*}
I(f) &= \int_0^t \int V_\mu^{0,s}\left(\mathcal{A}_s^\mu - \mathcal{A}_s^\nu\right)V_\nu^{s,t}f(y) \mu_0(dy) ds \leq \int_0^t ||V_\mu^{0,s}\left(\mathcal{A}_s^\mu - \mathcal{A}_s^\nu\right)V_\nu^{s,t}f||_{L^\infty} ds,
\intertext{due to the definition of $V^{0,s}_\mu$ and the fact that $p^\mu(0,x;s,y)dy$ is a probability measure we have}
&\leq \int_0^t ||\left(\mathcal{A}_s^\mu - \mathcal{A}_s^\nu\right)V_\nu^{s,t}f||_{L^\infty} ds\\
&\leq \sup_{s,x}\left(|B^\mu(s,x)-B^\nu(s,x)|+|c^\mu(s,x)-c^\nu(s,x)|\right)\int_0^t ||\bigtriangledown V^{s,t}_\nu f||_{L^\infty} + ||Hess\ V^{s,t}_\nu f||_{L^\infty} ds.
\end{align*}
Now by Proposition \ref{ap2} if $f\in bC^\alpha$ with $||f||_{bC^\alpha}\leq 1$ we have 
$$
||\bigtriangledown V^{s,t}_\nu f||_{L^\infty}\leq \frac{C}{\sqrt{t-s}},\qquad ||Hess\ V^{s,t}_\nu f||_{L^\infty}\leq \frac{C}{(t-s)^{1-\frac{\alpha}{2}}},
$$
which yields for $I(f)$
\begin{align*}
I(f)&\leq C\sup_{s,x} \left(|B^\mu(s,x)-B^\nu(s,x)|+|c^\mu(s,x)-c^\nu(s,x)|\right)\int_0^t \left(\frac{1}{|t-s|^{1-\frac{\alpha}{2}}} + \frac{1}{\sqrt{t-s}}\right)ds\\
&\leq C_T|t|^\frac{\alpha}{2}\sup_{s,x} \left(|B^\mu(s,x)-B^\nu(s,x)|+|c^\mu(s,x)-c^\nu(s,x)|\right).
\end{align*}
Now we observe that since the coefficients are uniformly $bC^\alpha$ we have
$$
|B^\mu(s,x)-B^\nu(s,x)|=\left| \int b(s,x,y)\mu_s(dy) - \int b(s,x,y)\nu_s(dy) \right|\leq C d_{bC^{\alpha}}(\mu_s,\nu_s),
$$
where $C=||b||_{C^{\alpha}_B}$, a priori it depends on $(s,x)$ but since $b$ uniformly $bC^\alpha$ it can be taken uniformly in $(s,x)$. It is also possible to prove that
$$
|c^\mu(s,x)-c^\nu(s,x)|\leq C ||\sigma||_{\infty}d_{bC^{\alpha}}(\mu_s,\nu_s).
$$
with these we can conclude that
$$
I(f)\leq C|t|^\frac{\alpha}{2}\sup_{s\in[0,t]} d_{C^{\alpha}_B}(\mu_s,\nu_s),
$$
for any $f\in bC^\alpha$ with $||f||_{bC^\alpha}\leq 1$ and thus
$$
\sup_{t\in[0,T]}d_{C^{\alpha}_B}(\mathcal{L}^\mu_t,\mathcal{L}^\nu_t)= \sup_{t\in[0,T]}\sup_{||f||_{bC^\alpha}\leq 1}I(f)\leq C|T|^\frac{\alpha}{2}\sup_{t\in[0,T]}d_{bC^{\alpha}}(\mu_t,\nu_t).
$$
which proves contraction for small valus of $T$ for $\mathcal{L}$, now by a classical sequencial argument we may conclude that there exists a unique flow of marginals $(\mu_t)_{t\in[0,T]}$ such that $[X^\mu_t]=\mu_t,\ \forall t\in[0,T]$ and thus $X^\mu$ is the unique weak solution to the MKV SDE.
\end{proof}
\section{Appendix}
\begin{theorem}\label{ap1}
The bounded $\alpha$-Holder distance metrizes weak convergence of measures. More precisely given $(\mu_n)_{n\in\N}$ and $\mu$ probability measures
$$
d_{bC^{\alpha}}(\mu_n,\mu)\rightarrow 0 \Leftrightarrow \mu_n\stackrel{d}{\rightarrow}\mu.
$$
\end{theorem}
\begin{proof}
The proof will be divided in two steps and is mostly taken from \url{https://sites.stat.washington.edu/jaw/COURSES/520s/522/HO.522.20/ch11c.pdf}

1) First we prove that
$$
\mu_n\stackrel{d}{\rightarrow}\mu \Leftrightarrow \int f d\mu_n\rightarrow\int f d\mu,\ \forall f\in bC^{\alpha}.
$$
If $\mu_n\stackrel{d}{\rightarrow}\mu$ then equivalently $\int f\mu_n\rightarrow\int f\mu$ for any function $f\in bC$ which in particular means that it is true for any $f\in bC^{\alpha}$. The converse is true because if $\int f\mu_n\rightarrow\int f\mu$ for any function $f\in bC^{\alpha}$ then in particular it is true for any $f\in bLip$ which by Portmanteau's theorem implies weak convergence.

2) We will now prove that
$$
\int f d\mu_n\rightarrow\int fd\mu \ \forall f\in bC^{\alpha} \Leftrightarrow d_{bC^{\alpha}}(\mu_n,\mu)\rightarrow0.
$$
The easy implication is the right-to-left one: indeed by comparison theorem
$$
\lim_n \int f(x) \left( \mu_n(dx)-\mu(dx) \right)\leq ||f||_{bC^\alpha}\lim_n\sup_{||g||_{bC^{\alpha}}\leq 1}\left|\int g(x) \left( \mu_n(dx)-\mu(dx) \right)\right| = ||f||_{bC^\alpha}\lim_n d_{bC^{\alpha}}(\mu_n,\mu)\rightarrow0.
$$
The other way is more challenging, first by continuity from below of probability measures for any fixed $\epsilon>0$ there exists $K$ a compact set such that $\mu(K)>1-\epsilon$. Let $\mathcal{H}=\left\lbrace f\in bC^{\alpha}\ |\ ||f||_{bC^{\alpha}}\leq 1 \right\rbrace$, if we restrict each of these functions on $K$ we have that $\mathcal{H}\vert_K$ is totally bounded with respect to the $||\cdot||_\infty$ norm by Ascoli-Arzelà's theorem, in particular $\exists k$ finite and $f_1,\cdots f_k\in \mathcal{H}\vert_K$ such that for any $f\in\mathcal{H}$ $\exists j$ such that $\sup_K|f-f_j|\leq\epsilon$.

Now if we consider $d(x,y)=|x-y|$ and $K^\epsilon=\left\lbrace x\in\R^N\ |\ d(x,K)\leq\epsilon \right\rbrace$ and $f,f_j$ as before we have
$$
\sup_{x\in K^\epsilon}|f(x)-f_j(x)|\leq \sup_{x\in K^\epsilon}\left( |f(x)-f(y_x)| + |f(y_x)-f_j(y_x)| + |f_j(y_x)-f_j(x)|\right)\leq \sup_{x\in K^\epsilon}\left( 2\epsilon^\alpha + \epsilon \right)\leq C_\alpha\epsilon^\alpha.
$$
where $y_x$ is a point in $K$ such that $d(x,y_x)<\epsilon$. $C_\alpha$ may be taken uniformly of $\epsilon$ as long as $\epsilon\leq 1$.

Let $g(x)=\max\left( 0, 1-\frac{d(x,K)}{\epsilon} \right)$, evidently $g\in bLip\subseteq bC^{\alpha}$ and $\mathds{1}_K\leq g\leq \mathds{1}_{K^\epsilon}$. Thus by taking $n$ big enough we have by convergence against $bC^{\alpha}$ functions that
$$
\mu_n(K^\epsilon)\geq\int g(x) \mu_n(dx) > 1 - 2\epsilon.
$$
Thus by taking $f\in\mathcal{H}$ and the associated $f_j$ we have
\begin{align*}
\left| \int f(x)\left(\mu_n(dx)-\mu(dx)\right) \right|&= \left| \int (f(x)-f_j(x))\left(\mu_n(dx)-\mu(dx)\right) \right| + \left| \int f_j(x)\left(\mu_n(dx)-\mu(dx)\right) \right|\\
&\leq \left| \int (f(x)-f_j(x))\mu_n(dx) \right| + \left| \int (f(x)-f_j(x))\mu(dx) \right| + \left| \int f_j(x)\left(\mu_n(dx)-\mu(dx)\right) \right|\\
&\leq \left| \int_{K^\epsilon} (f(x)-f_j(x))\mu_n(dx) \right| + \left| \int_{(K^\epsilon)^c} (f(x)-f_j(x))\mu_n(dx) \right| + \\ &\qquad +\left| \int_{K^\epsilon} (f(x)-f_j(x))\mu(dx) \right| +  \left| \int_{(K^\epsilon)^c} (f(x)-f_j(x))\mu(dx) \right| + \\ &\qquad + \left| \int f_j(x)\left(\mu_n(dx)-\mu(dx)\right) \right|\\
&\leq C_\alpha\epsilon^\alpha + 4\epsilon + C_\alpha\epsilon^\alpha+2\epsilon+\epsilon\leq C_\alpha\epsilon^\alpha,
\end{align*}
where the last term gets bounded by taking $n$ big enough and by using convergence against $bC^{\alpha}$ functions, this gives us the final result.
\end{proof}
\begin{proposition}\label{ap2}
If $f$ is a $bC^\alpha$ function with $||f||_{bC^\alpha}\leq 1$ then
\begin{enumerate}
\item $||\bigtriangledown\ V^{s,t}_\nu f||_{L^\infty}\leq\frac{C}{\sqrt{t-s}},$
\item $||Hess\ V^{s,t}_\nu f||_{L^\infty}\leq \frac{C}{|t-s|^{1-\frac{\alpha}{2}}}$.
\end{enumerate}
\end{proposition}
\begin{proof}
\begin{enumerate}
\item Simply by gaussian estimates (Theorem 20.2.5 of \cite{Pascu_PB2})
\begin{align*}
\left| \partial_{x_i}V^{s,t}_\nu f (x) \right|&=\left| \int\partial_{x_i}p^\nu(s,x;t,y)f(y)dy \right| \\
&\leq \frac{C}{\sqrt{t-s}}\int \Gamma^+(t-s,x-y)|f(y)|dy\leq \frac{C}{\sqrt{t-s}}
\end{align*}
\item For the second derivates it is imperative the use of the $\alpha$-Holderian structure of $f$
\begin{align*}
\left| \partial_{x_ix_j}V^{s,t}_\nu f (x) \right|&=\left| \int\partial_{x_i x_j}p^\nu(s,x;t,y)f(y)dy \right| \\
&\leq \left| \int\partial_{x_i x_j}p^\nu(s,x;t,y)(f(y)-f(x))dy \right|+ \left| \int\partial_{x_i x_j}p^\nu(s,x;t,y)dyf(x) \right|\\
&\leq \int |\partial_{x_ix_j}p^\nu(s,x;t,y)||x-y|^\alpha dy + \left| \partial_{x_ix_j}\underbrace{\int p^\nu(s,x;t,y) dy}_{=1} \right||f(x)|\\
&\leq \frac{C}{|t-s|}\int \Gamma^+(t-s,x-y)|x-y|^\alpha dy + 0\leq \frac{C_{B,\alpha}}{|t-s|^{1-\frac{\alpha}{2}}}.
\end{align*}
by Theorem 20.2.5 and Lemma 20.3.4 of \cite{Pascu_PB2}.
\end{enumerate}
\end{proof}

\begin{thebibliography}{90}
%\bibitem{Bogachev} Vladimir I. Bogachev - Measure Theory (2007)
\bibitem{YAOZHONG} Yaozhong Hu, Michael A. Kouritzin, Jiayu Zheng -  Nonlinear McKean-Vlasov diffusions under the weak Hormander condition with quantile-dependent coefficients (2021), https://arxiv.org/abs/2101.04080
\bibitem{Kolokoltsov} Vassili N. Kolokoltsov -  Nonlinear Diffusions and Stable-Like Processes
with Coefficients Depending on the Median or VaR (2013)
\bibitem{LucePagliaPascu} G. Lucertini, A. Pagliarani, A. Pascucci - Optimal regularity for degenerate Kolmogorov equations in non-divergence form with rough-in-time coefficients (2024), https://doi.org/10.1007/s00028-023-00916-9
\bibitem{Pascu_PB2} A. Pascucci - Probability Theory II (2024), https://doi.org/10.1007/978-3-031-63193-1
%\bibitem{Villani} Cédric Villani - Optimal Transport Old and New (2009)
\end{thebibliography}
\end{document}