The mean field Schr\"odinger problem: ergodic behavior, entropy estimates and functional inequalities

We study the mean field Schr\"odinger problem (MFSP), that is the problem of finding the most likely evolution of a cloud of interacting Brownian particles conditionally on the observation of their initial and final configuration. Its rigorous formulation is in terms of an optimization problem with marginal constraints whose objective function is the large deviation rate function associated with a system of weakly dependent Brownian particles. We undertake a fine study of the dynamics of its solutions, including quantitative energy dissipation estimates yielding the exponential convergence to equilibrium as the the time between observations grows larger and larger, as well as a novel class of functional inequalities involving the mean field entropic cost (i.e. the optimal value in (MFSP)). Our strategy unveils an interesting connection between forward backward stochastic differential equations and the Riemannian calculus on the space of probability measures introduced by Otto, which is of independent interest.


Introduction and statement of the main results
In the seminal works [44,45] E. Schrödinger addressed the problem of finding the most likely evolution of a cloud of independent Brownian particles conditionally on the observation of their initial and final configuration. In modern language this is an entropy minimization problem with marginal constraints. The aim of this work is to take the first steps in the understanding of the Mean Field Schrödinger Problem, obtained by replacing in the above description the independent particles by interacting ones.
To obtain an informal description of the problem, consider N Brownian particles (X i,N t ) t∈[0,T ],1≤i≤N interacting through a pair potential W Their evolution is encoded in the random empirical path measure At a given time T , the configuration of the particle system is visible to an external observer that finds it close to an "unexpected" (écart spontané et considerable in [45]) probability measure µ fin , It is a classical result [47,20,4] that the sequence of empirical path measures (2) obeys the large deviations principle (LDP). Thus, the problem of finding the most likely evolution conditionally on the observations is recast as the problem of minimizing the LDP rate function among all path measures whose marginal at time 0 is µ in and whose marginal at time T is µ fin . This is the mean field Schrödinger problem (MFSP). Extending naturally the classical terminology we say that an optimal path measure is a mean field Schrödinger bridge (henceforth MFSB) and the optimal value is the mean field entropic cost. The latter generalizes both the Wasserstein distance and the entropic cost.
The classical Schrödinger problem has been the object of recent intense research activity (see [36]). This is due to the computational advantages deriving from introducing an entropic penalization in the Monge-Kantrovich problem [19] or to its relations with functional inequalities, entropy estimates and the geometrical aspects of optimal transport. Our article contributes to this second line of research, recently explored by the papers [17,28,34,42,32]. Leaving all precise statements to the main body of the introduction, let us give a concise summary of our contributions.
Dynamics of mean field Schrödinger bridges Our mean field version of the Schrödinger problem stems from fundamental results in large deviations for weakly interacting particle systems such as [47,20] and shares some analogies with the control problems considered in [16] and with the article [2] in which an entropic formulation of second order variational mean field games is studied. Among the more fundamental results we establish for the mean field Schrödinger problem, we highlight • the existence of MFSBs and, starting from the original large deviations formulation, the derivation of both an equivalent reformulation in terms of a McKean-Vlasov control problem as well as a Benamou-Brenier formula, • establishing that MFSBs solve forward backward stochastic differential equations (FBSDE) of McKean-Vlasov type (cf. [10,11]).
The proof strategy we adopt in this article combines ideas coming from large deviations and stochastic calculus of variations, see [18,49,23]. Another interesting consequence of having a large deviations viewpoint is that we can also exhibit some regularity properties of MFSBs, taking advantage of Föllmer's results [26] on time reversal. Building on [17,28] we establish a link between FBSDEs and the Riemannian calculus on probability measures introduced by Otto [40] that is of independent interest and underlies our proof strategies. In a nutshell, the seminal article [31] established that the heat equation is the gradient flow of the relative entropy w.r.t. the squared Wasserstein distance. Thus, classical first order SDEs yield probabilistic representations for first order ODEs in the Riemannian manifold of optimal transport. Our observation may be seen as the second order counterpart to the results of [31]: indeed we will present an heuristic strongly supporting the fact that Markovian solutions of "second order" pathwise equations (FBSDEs) yield probabilistic representations for second order ODEs in the Riemannian manifold of optimal transport.
Ergodicity of Schrödinger bridges and functional inequalities Consider again (1) and assume that W is convex so that the particle system is rapidly mixing and there is a well defined notion of equilibrium configuration µ ∞ . If N and T are large, one expects that (i) The configurations 1 N N i=1 δ X i t at times t = 0, T /2, T are almost independent.
(ii) The configuration at T /2 is with high probability very similar to µ ∞ .
Because of (i), even when the external observer acquires the information (3), he/she still expects (ii) to hold. Thus mean field Schrödinger bridges are to spend most of their time around the equilibrium configuration. All our quantitative results originate in an attempt to justify rigorously this claim.
In this work we obtain a number of precise quantitative energy dissipation estimates. These lead us to the main quantitative results of the article: • we characterize the long time behavior of MFSBs, proving exponential convergence to equilibrium with sharp exponential rates, • we derive a novel class of functional inequalities involving the mean field entropic cost. Precisely, we obtain a Talagrand inequality and an HWI inequality that generalize those previously obtained in [12] by Carrillo, McCann and Villani.
Regarding the second point above, we can in fact retrieve (formally) the inequalities in [12] by looking at asymptotic regimes for the mean field Schrödinger problem. Besides the intrinsic interest and their usefulness in establishing some of our main results, our functional inequalities may have consequences in terms of concentration of measure and hypercontractivity of non linear semigroups, but this is left to future work. The fact that optimal curves of a given optimal control problem spend most of their time around an equilibrium is known in the literature as the turnpike property. The first turnpike theorems have been established in the 60's for problems arising in econometry [38]; general results for deterministic finite dimensional problems are by now available, see [48]. In view of the McKean-Vlasov formulation of the mean field Schrödinger problem, some of our results may be viewed as turnpike theorems as well, but for a class of infinite dimensional and stochastic problems. An interesting feature is that, by exploiting the specific structure of our setting, we are able to establish the turnpike property in a quantitative, rather than qualitative form. The McKean-Vlasov formulation also connects our findings with the study of the long time behavior of mean field games [8,7,5,9].
Concerning the proof methods, our starting point is Otto calculus and the recent rigorous results of [17] together with the heuristics put forward in [28]. The first new ingredient of our proof strategy is the above mentioned connection between FBSDEs and Otto calculus that plays a key role in turning the heuristics into rigorous statements. It is worth remarking that using a pathwise approach does not just provide with a way of making some heuristics rigorous, but it also permits to obtain a stronger form of some of the results conjectured in [28] which then simply follow by averaging pathwise estimates. The second new ingredient in our proofs involves a conserved quantity that plays an analogous role to the total energy of a phyisical system. For such quantity we derive a further functional inequality which seems to be novel already in the classical Schrödinger problem (i.e. for independent particles) and allows to establish the turnpike property.
Structure of the article In the remainder introductory section we state and comment our main results. In Section 2 we provide a geometrical interpretation sketching some interesting heuristic connections between optimal transport and stochastic calculus. The material of this section is not used later on; therefore the reader who is not interested in optimal transport may avoid it. Sections 3 and 4 contain the proofs of our main results, the former being devoted to the results concerning the dynamics of MFSBs and the latter one dealing with the ergodic results. Finally an appendix section contains some technical results.

Frequently used notation
• (Ω, F t , F T ) is the canonical space of R d -valued continuous paths on [0, T ], so {F t } t≤T is the coordinate filtration. Ω is endowed with the uniform topology.
• P(Ω) and P(R d ) denote the set of Borel probability measures on Ω and R d respectively.
• R µ is the Wiener measure with starting distribution µ.
• H(P|Q) denotes the relative entropy of P with respect to Q, defined as E P log dP dQ if P ≪ Q and +∞ otherwise.
• P t denotes the marginal distribution of a measure P ∈ P(Ω) at time t.
• P β (Ω) is the set of measures on Ω for which sup t≤T |X t | β is integrable. P β (R d ) is the set of measures on R d for which the function | · | β is integrable.
• The β-Wasserstein distance on P β (Ω) is defined by With a slight abuse of notation we also denote by W β the β-Wasserstein distance on P β (R d ) defined analogously.
• For a given measurable marginal flow the space of square integrable functions from [0, T ] × R d to R d associated to the reference measure µ t (dx)dt and the corresponding almost-sure identification. We consider likewise the Hilbert space defined as the closure in L 2 ((µ t ) t∈[0,T ] ) of the smooth subspace • γ and λ are respectively the standard Gaussian and Lebesgue measure in R d .
is the space of real-valued smooth functions on [0, T ] × R d with compact support. The gradient ∇ and Laplacian ∆ act only in the space variable.
• If f is a function and µ a measure, its convolution is x → f * µ(x) := f (x − y)µ(dy).

The mean field Schrödinger problem and its equivalent formulations
We are given a so-called interaction potential W : R d → R, for which we assume W is of class C 2 (R d ; R) and symmetric, i.e. W (·) = W (−·), (H1) Besides the interaction potential, the data of the problem are a pair of probability measures µ in , µ fin on which we impose where the free energy or entropy functionalF is defined for µ ∈ P 2 (R d ) bỹ In the above, and in the rest of the article, we shall make no distinction between a measure and its density against Lebesgue measure λ, provided it exists.
We recall that the McKean Vlasov dynamics is the non linear SDE Under the hypothesis (H1), it is a classical result (see e.g. [13,Thm 2.6]) that (5) admits a unique strong solution whose law we denote P MKV . The functionalF plays a crucial role in the sequel. For the moment, let us just remark that the marginal flow of the McKean-Vlasov dynamics may be viewed as the gradient flow of 1 2F in the Wasserstein space (P 2 (R d ), W 2 (·, ·)). If P ∈ P 1 (Ω) is given, then the stochastic differential equation admits a unique strong solution (cf. Section 3.2) whose law we denote Γ(P). With this we can now introduce the main object of study of the article: Definition 1.1. The mean field Schrödinger problem is inf H(P|Γ(P)) : P ∈ P 1 (Ω), P 0 = µ in , P T = µ fin . (MFSP) Its optimal value, denoted C T (µ in , µ fin ), is called mean field entropic transportation cost. Its optimizers are called mean field Schrödinger bridges (MFSB).
Remark 1.1. Of course, the choice of the pair potential W as interaction mechanism is a particular one. Thus (MFSP) is by far not the only mean field Schrödinger problem of interest. Also, it would have been easy to include in the dynamics a confinement (single-site) potential. However, since one of the goals of this article is to understand the role of the pair potential, we preferred not to do that. Indeed, in many situations, the single site potential may be the one that determines the long time behavior of mean field Schrödinger bridges, thus obscuring the role of the pair potential.
It is not difficult to provide existence of optimizers for (MFSP). In the classical case, uniqueness is an easy consequence of the convexity of the entropy functional. However the rate function H(P|Γ(P)) is not convex in general. Therefore, it is not clear whether or not (MFSP) admits a unique solution.

Large deviations principle (LDP)
We start by deriving the LDP interpretation of (MFSP). We consider for N ∈ N the interacting particle system where {B i : i = 1, . . . , N } are independent Brownian motions and {X i,N 0 : i = 1, . . . , N } are independent to each other and to the Brownian motions. The theory of stochastic differential equations guarantees the strong existence and uniqueness for the above interacting particle system under (H1),(H2). In the next theorem we obtain a LDP for the sequence of empirical path measures; in view of the classical results of [20], it is not surprising that the LDP holds. However, even the most recent works on large deviations for weakly interacting particle systems such as [4] do not seem to cover the setting and scope of Theorem 1.1. Essentially, this is because in those references the LDPs are obtained for a topology that is weaker than the W 1 -topology, that is what we need later on. Theorem 1.1. In addition to (H1),(H2) assume that Then the sequence of empirical measures satisfies the LDP on P 1 (Ω) equipped with the W 1 -topology, with good rate function given by P 1 (Ω) ∋ P → I (P) := H(P|Γ(P)) , P ≪ Γ(P), +∞ , otherwise.
In fact we will prove in Section 3 a strengthened version of Theorem 1.1 where the drift term is much more general. For this, we will follow Tanaka's elegant reasoning [47].

Remark 1.2. Having a rate function implies heuristically
Hence Problem (MFSP) has the desired interpretation of finding the most likely evolution of the particle system conditionally on the observations (when N is very large). This is in fact a consequence of the Gibbs conditioning principle; see [21] and the references therein.

McKean-Vlasov control and Benamou-Brenier formulation
We now reinterpret the mean field Schrödinger problem (MFSP) in terms of McKean-Vlasov stochastic control (also known as mean field control). Lemma 1.1. Let P be admissible for (MFSP). There exists a predictable process (α P t ) t∈[0,T ] s.t.
and so that has law R µ in under P. The problem (MFSP) is equivalent to as well as to The formulations (11)- (12) can be seen as McKean-Vlasov stochastic control problems. In the first case one is steering through α P part of the drift of a McKean-Vlasov SDE. In the second case one is controlling the drift Φ of a standard SDE but the optimization cost depends non-linearly on the law of the controlled process. In both cases, the condition P T = µ fin is rather unconventional. By analogy with the theory of mean field games, one could refer to (11)-(12) as planning McKean-Vlasov stochastic control problems, owing to this type of terminal condition.
The third and last formulation of (MFSP) we propose relates to the well known fluid dynamics representation of the Monge Kantorovich distance due to Benamou and Brenier (cf. [50]) that has been recently extended to the standard entropic transportation cost [27,15]. The interest of this formula is twofold: on the one hand it clearly shows that (MFSP) is equivalent and gives a rigorous meaning to some of the generalized Schrödinger problems formally introduced in [28,34]. On the other hand, it allows to interpret (MFSP) as a control problem in the Riemannian manifold of optimal transport. This viewpoint, that we shall explore in more detail in Section 2, provides with a strong guideline towards the study of the long time behavior of Schrödinger bridges.
We define the set A as the collection of all absolutely continuous curves (µ t ) t∈[0,T ] ⊂ P 2 (R d ) (cf. Section 4.2) such that µ 0 = µ in , µ T = µ fin and We then define Theorem 1.2. Let (H1),(H2) hold. Then If P is optimal for (MFSP) and the latter is finite, then (P t ) t∈[0,T ] is optimal in (13) and its associated tangent vector field w is given by where Ψ is as in Theorem 1.3 below. Conversely, if (µ t ) t∈[0,T ] is optimal for C BB T (µ in , µ fin ) and the latter is finite, then there exists an optimizer of C T (µ in , µ fin ) whose marginal flow equals (µ t ) t∈[0,T ] .

Mean field Schrödinger bridges
Leveraging the stochastic control interpretation, and building on the stochastic calculus of variations perspective, we obtain the following necessary optimality conditions for (MFSP). Theorem 1.3. Assume (H1),(H2) and let P be optimal for (MFSP). Then there exist Ψ ∈ where (α P t ) t∈[0,T ] is related to P as in Lemma 1.1. The process t → Ψ t (X t ) is continuous 1 and the process (M t ) t∈[0,T ] defined by is a continuous martingale under P on [0, T [, where (X t ) t∈[0,T ] is an independent copy of (X t ) t∈[0,T ] defined on some probability space (Ω,F,P) andẼP denotes the expectation on (Ω,F ,P).
We shall refer to Ψ as the corrector of P. Correctors will play an important role in the ergodic results. In this part, we give an interpretation of Theorem 1.1 in terms of stochastic analysis (FBSDEs) and partial differential equations.

Planning McKean-Vlasov FBSDE for MFSB
We consider the following McKean Vlasov forward-backward stochastic differential equation (FB-SDE) in the unknowns (X, Y, Z): As in the stochastic control interpretation of the mean field Schrödinger problem, here too the terminal condition X T ∼ µ fin is somewhat unconventional. We hence call this forward-backward system the planning McKean-Vlasov FBSDE.
Thanks to the results in Section 1.2.2 we can actually solve (16). If P is optimal for (MFSP) with associated Ψ as recalled in Theorem 1.3 above, all we need to do is take Y t := Ψ t (X t ) and reinterpret (10) for the dynamics of the canonical process X and (15) for the dynamics of Y (in the latter case using martingale representation).
One remarkable aspect of this connection between Schrödinger problems and FBSDEs is that one can prove existence of solutions to such FBSDEs by a purely variational method. Indeed, we remark that (16) is beyond the scope of existing FBSDE theory, such as Carmona and Delarue's [11,Theorem 5.1]. Further, we also obtained for free an extra bit of information: the constructed process Y lives in H −1 ((P t ) t∈[0,T ] ). This is in tandem with the usual heuristic relating FBSDEs and PDEs (where Y is conjectured to be an actual gradient) as explained in Carmona and Delarue's [10,Remark 3.1]. In fact, if we make the additional assumption that Y t = ∇ψ t (X t ) for some timedependent sufficiently smooth potential ψ t (x), and we set µ t = (X t ) # P, then the first equation in (16) tells that ∂ t µ t − 1 2 ∆µ t + ∇ · ((−∇W * µ t + ∇ψ t )µ t ) = 0, whereas the second equation of (16) and a standard argument with Itô formula yields that for each where L is the non-linear generator 1 2 ∆ − ∇W * µ t − ∇ψ t · ∇. All in all we arrive at the PDE system 2 : It is worth stressing that the link between (16) and (17) can be established if the FBSDE solution Y t is a gradient vector field depending only on t and X t . We have gathered preliminary evidence that (16) admits non Markovian solutions even in the simple case when W = 0. This is in contrast with what is expected for standard FBSDEs [10, Lemma 3.5] whose boundary conditions are not of planning type.

Schrödinger potentials and the mean field planning PDE system
The PDE system (17) is the literal translation of the planning McKean-Vlasov FBSDE in the case when the process Y is an actual gradient, Y = ∇ψ. In the next corollary we show that if this is the case, and if ψ is sufficiently regular, then (17) can be rewritten as a system of two coupled PDEs, the first being a Hamilton-Jacobi-Bellman equation for ψ, and the second one being a Fokker-Planck equation. This type of PDE system is the prototype of a planning mean field game [33]. Corollary 1.1. Let P be an optimizer for (MFSP), Ψ · (·) be as in Theorem 1.3 and set µ t = P t for all t ∈ [0, T ]. If µ · (·) is everywhere positive and of class C 1,2 ([0, T ] × R d ; R) and Ψ · (·) is of A fundamental result [25] concerning the structure of optimizers in the classical Schrödinger problem is that their density takes a product form, i.e.
where ϕ t (x), ψ t (x) solve respectively the forward and backward Hamilton Jacobi Bellman equation It is interesting to see that this structure is preserved in (MFSP), at least formally. The effect of having considered interacting Brownian particles instead of independent ones is reflected in the fact that the two Hamilton Jacobi Bellman PDEs are coupled not only through the boundary conditions but also through their dynamics. Corollary 1.2. Using the same notation and under the same hypotheses of Remark 1.3. In the article [34] the authors implement another approach in order to generalize the decomposition (19) to interacting systems based on the Hopf-Cole transform, whereas ours relies on time reversal. In particular, in their Example 4, they deal with a generalized Schrödinger problem that would correspond to erasing the term 1 2 ∇ log µ t in (13). It appears that the variables η, η * produced by the Hopf-Cole transformation do not coincide with the variables ψ, ϕ we use here. Therefore these two approaches may be seen as complementary.

Convergence to equilibrium and functional inequalities
Our aim is to show that MFSBs spend most of their time in a small neighborhood of the equilibrium configuration µ ∞ , to study their long time behavior, and to derive a new class of functional inequalities involving the mean field entropic cost C T (µ in , µ fin ).
Throughout this section we make the assumption that W is uniformly convex, ie. that where the inequality above has to be understood as an inequality between quadratic forms. Under (H3) the McKean Vlasov dynamics associated with the particle system (6) converges in the limit as T → +∞ to an equilibrium measure µ ∞ , that is found by minimizing the functionalF over the elements of P 2 (R d ) whose mean is the same as µ in . Existence and uniqueness of µ ∞ has been proven in [37]. We shall often assume that µ in and µ fin have the same mean: is a classical one ensuring exponential convergence rates for the McKean-Vlasov dynamics. It may be weakened in various ways, see the work [12] by Carrillo, McCann and Villani or the more recent [3] by Bolley, Gentil and Guillin, for instance. It is an interesting question to determine which of the results of this section still hold in the more general setup. Hypothesis (H4) can be easily removed using the fact that the mean evolves linearly along any Schrödinger bridge (see Lemma 4.2 below); the reason for imposing (H4) is that it allows to simplify the presentation of our results. We insist that all our quantitative results below use the technical assumptions (H1) and (H2) only in order to make sure that all objects we are dealing with are well defined. If one can ensure this by other means, the results are still valid under the key assumption (H3).
Long time behavior of mean field games: The articles [8,7,5,9] study the asymptotic behaviour of dynamic mean field games showing convergence towards an ergodic mean field game with exponential rates. Following [33], we can associate to (18) the ergodic problem with unknowns (λ, ψ, µ) This PDE system expresses optimality conditions for the ergodic control problem corresponding to (11). It is easy to see that (0, 0, µ ∞ ) is a solution of (20). Therefore, we are addressing the same questions studied in the above mentioned articles. However, the equations we are looking at are quite different. A fundamental difference is that the coupling terms in (11) are not monotone in the sense of [6, Eq.(7) pg. 8].

Exponential convergence to equilibrium and the turnpike property
A key step towards the forthcoming quantitative estimates is to consider the time-reversed version of our mean field Schrödinger problem. For Q ∈ P(Ω) the time reversalQ is the law of the time reversed process (X T −t ) t∈[0,T ] . In Lemma 4.5 we prove that if P is an optimizer for (MFSP), then P optimizes The optimality ofP implies the existence of an associated processΨ as described in Theorem 1.3. We show at Theorem 1.7 below that the function is a constant, that we denote E P (µ in , µ fin ) and call the conserved quantity. Naturally this quantity depends also on T but we omit this from the notation. Theorem 1.4 confirms the intuition that mean field Schrödinger bridges are localized around µ ∞ providing an explicit upper bound for F (P t ) along any MFSB, where We recall that µ ∞ is found by minimizingF among all elements of P 2 (R d ) whose mean is the same as µ. IfF is thought of as a free energy, then F should be thought of as a divergence (from equilibrium). 3 Theorem 1.4. Assume (H1)-(H4) and let P be an optimizer for (MFSP).
In particular, since F (P θT ) dominates W 2 (P θT , µ ∞ ) (see e.g. [12, (ii), Thm 2.2 1]), we obtain that P θT converges exponentially fast to µ ∞ and that the exponential rate is proportional to κ. The proof of (24) is done by bounding the second derivative of the function t → F (P t ) along Schrödinger bridges with the help of the logarithmic Sobolev inequality established in [12]. To obtain (25) from (24) we use a functional inequality for the conserved quantity and a Talagrand inequality for C T (µ in , µ fin ), that are the content of Theorem 1.7 and Corollary 1.3 below. It is worth mentioning that the estimates (24),(25) (as well as (35) below) appear to be new even for the classical Schrödinger bridge problem and have not been anticipated by the heuristic articles [28,34]. Conversely, the above mentioned estimates admit a geometrical interpretation in the framework of Otto calculus that allows to formally extend their validity to the whole class of problems studied in [28].
Remark 1.5. The exponential rate in (25) has a sharp dependence on κ. To see this, fix µ in and choose µ fin = P MKV T . Then it is easy to see that the restriction of P MKV to the interval [0, T ] is an optimizer for (MFSP). Setting θ = 1/2 and considering (25) for T = 2t we arrive at Thus, we obtain the same exponential rate as in [12] 4 , that is easily seen to be optimal under the assumption that W is κ-convex. A similar argument can be used to show the optimal dependence of the rate in θ.
In the previous theorem we showed that, when looking at a timescale that is proportional to T , the marginal distribution of any Schrödinger bridge is exponentially close to µ ∞ . Here we show that for a fixed value of t, we have an exponential convergence towards the law of the McKean-Vlasov dynamics P MKV , see (5).
Theorem 1.5. Assume (H1)-(H4) and let P be an optimizer for (MFSP). For all t ∈ [0, T ] we have In particular, the above theorem tells that W 2 2 (P t , P MKV t ) decays asymptotically at least as fast as exp(−2κT ) when T is large.

Functional inequalities for the mean field entropic cost
It is well known that analysing the evolution of entropy-like functionals along the so-called displacement interpolation of optimal transport has far reaching consequences in terms functional inequalities [52]. Since (MFSP) provides with an alternative way of interpolating between probability measures, it is tempting to see if it leads to new functional inequalities involving the cost C T (µ in , µ fin ). Here, we present a Talagrand and an HWI inequality that we used in order to study the long time behavior of MFSBs. They generalize their respective counterparts in [46], [41]. Both inequalities are based on another upper bound for the evolution of F along MFSBs that we present now. This bound generalizes to the mean field setup the results of [17]: The above result may be seen as a rigorous version of some of the heuristic arguments put forward in [28] and [34], upon slightly modifying the definition of C T .
The following Talagrand inequality, obtained as a corollary of Theorem 1.6, tells that the mean field entropic cost grows at most linearly with F .
In particular, choosing µ fin = µ ∞ leads to Theorem 1.6 also yields an entropic HWI inequality where the Wasserstein distance is replaced by the conserved quantity E P in the first term on the rhs and by the mean field entropic cost in the second term. An extra positive contribution 1 4 I F is present in the first term. Our interpretation is that this compensates for the fact that in the "gain" term we put the cost C T , that is larger than the squared Wasserstein distance. In order to state the HWI inequality, we introduce the non linear Fisher information functional I F defined for µ ∈ P 2 (R d ) by where by ∇ log µ ∈ L 2 µ we mean µ ≪ λ and that log µ is an absolutely continuous function on R d whose derivative is in L 2 µ . The non linear Fisher information can be seen to be equal to the derivative of the free energyF along the marginal flow of the McKean Vlasov dynamics. Corollary 1.4 (An HWI Inequality). Assume (H1)-(H3) and choose µ fin = µ ∞ . If P is an optimizer for (MFSP) and t → I F (P t ) is continuous in a right neighbourhood of 0, then It is worth noticing that by letting T → +∞ in the above HWI inequality we obtain the logarithmic Sobolev inequality [12,Thm 2.2]. Indeed, C T (µ in , µ ∞ ) is always non negative and we shall see at Thoerem 1.7 is the standard entropic cost and we have under suitable hypothesis on µ in (see [39]) The heuristic arguments put forward in [28] tell that (32) is expected to be true even when W is a general potential satisfying (H1). Following again (32), one also expects that Putting (32) and (33) together we obtain an heuristic justification of the fact that in the limit as T → 0 (31) becomes the classical HWI inequality put forward in [12], namely Remark 1.6. In Corollary 1.4 we asked t → I F (P t ) to be continuous in a right neighbourhood of 0. This is somewhat annoying. However, we were not able conclude that (H1) and (H2) imply the desired continuity, although we could establish that they grant the continuity of Our last result is a functional inequality that establishes a hierarchical relation between the conserved quantity and the mean field entropic cost: the former is exponentially small in T and κ in comparison with the latter. We may refer to this as an energy-transport inequality since the conserved quantity may be geometrically interpreted as the total energy of a physical system (cf [17, Corollary 1.1]). Theorem 1.7. Assume (H1)-(H4) and let P be an optimizer. Then the function is constant. Denoting this constant by E P (µ in , µ fin ), we have 2 Connections with optimal transport In this section we shall see how the results of this article relate to the Riemannian calculus on P 2 (R d ) introduced by Otto [40], at least formally. The link is rooted in a seemingly novel connection between (McKean-Vlasov) FBSDEs and second order ODEs in the Riemannian manifold of optimal transport that we find of independent interest. To better understand this connection, let us begin by recalling that in the seminal article [31] it is proven that the marginal flow of the pathwise SDE can be interpreted as the gradient flow of the entropy functional w.r.t. the 2-Wasserstein metric. Thus, first order Itô SDEs provide with probabilistic representations for first order ODEs in the Riemannian manifold of optimal transport. Of course, since a path measure is not fully determined by its one time marginals, the SDE (36) contains more information than the gradient flow equation. It has been shown in [17] that the marginal flow of a classical Schrödinger bridge satisfies a second order ODE, more precisely a Newton's law in which the acceleration field is the Wasserstein gradient of the Fisher information functional. The natural question is then: What pathwise (second order) SDE governs the dynamics of a Schrödinger bridge and yields a probabilistic representation for the associated Newton's law?
In order to answer this, let us first recall some notions of Otto calculus.
2.1 Second order calculus on P 2 (R d ).
In the next lines, we sketch the ideas behind the Riemannian calculus on P 2 (R d ). It would be impossible to provide a self-contained introduction in this work and we refer to [50] or [29] for detailed accounts. The main idea is to equip P 2 (R d ) with a Riemannian metric such that the associated geodesic distance is W 2 (·, ·). To do this, one begins by identifying the tangent space T µ P 2 at µ ∈ P 2 (R d ) as the space closure in L 2 µ of the subspace of gradient vector fields The velocity (first derivative) of a sufficiently regular curve [0, T ] ∋ t → µ t ∈ P 2 (R d ) is then defined by looking at the only solution v t (x) of the continuity equation such that v t ∈ T µt P 2 for all t ∈ [0, T ]. Finally, the Riemannian metric (Otto metric) ·, · TµP2 is defined by It can be seen that the constant speed geodesic curves associated to the Riemannian metric we have introduced coincide with the displacement interpolations of optimal transport and that the corresponding geodesic distance is indeed W 2 (·, ·). This makes it possible to carry out several explicit calculations. In particular, we can compute the gradient grad W2 F and the Hessian Hess W2 F of a smooth functional F : P 2 (R d ) → R. At least formally, we have , where we used the notation # for the push forward. In particular, setting W = 0 for simplicity in (30) we obtain that the classical Fisher information functional I has a gradient that can be computed with the rules above. One obtains that (cf. [51]) The Levi-Civita connection associated to the Riemannian metric (37) can also be explicitly computed with the help of the orthogonal projection operator Π µ : L 2 µ → T µ P 2 . To do this, consider a regular curve (µ t ) t∈[0,T ] with velocity (v t ) t∈[0,T ] and a tangent vector field t → u t ∈ T µt P 2 along then this covariant derivative satisfies the compatibility with the metric and the torsion-free identity, i.e. it is the Levi-Civita connection. The acceleration of the curve (µ t ) t∈[0,T ] is then the covariant derivative of the velocity along the curve, i.e.

Newton's laws and FBSDEs
According to the above discussion the Newton's law in (P 2 (R d ), ., . T·P2 ) provides with a geometrical interpretation for the PDE system (see [17] for more details) where to derive the latter equation we observe that the requirement that v t ∈ T µt P 2 for all t ∈ [0, T ] is formally equivalent to v t = ∇φ t for some time dependent potential (t, x) → φ t (x).
As we have seen in section 1.3.1, solutions of the FBSDE having the additional property that Y t = ∇ψ t (X t ) yield a probabilistic representation for Some tedious though standard calculations allow to see that the change of variable φ t = − 1 2 log µ t + ψ t transforms the PDE system (42) in (40). Summing up, we have obtained the following Informal statement We have: varying potential ψ, then the marginal flow (µ t ) t∈[0,T ] of X t is a solution for the Newton's law (39).
(ii) If P is the (classical) Schrödinger bridge between µ in and µ fin , then under P the canonical process We leave it to future work to prove a rigorous version of the informal statement above. On the formal level, there is no conceptual difficulty in extending it to include the interaction potential W . Essentially, the only difference is that one has to deal with the non linear Fisher information functional I F instead of I.
Beside its intrinsic interest, the parallelism between Newton's laws and FBSDEs is very useful when studying the long time behavior of the latter. Indeed, the Riemannian structure underlying (39) allows to find tractable expressions for the first and second derivative of entropy-like functionals along the marginal flow of the FBSDE.
Remark 2.1. Classical Schrödinger bridges are h−transforms in the sense of Doob [22]. Therefore, one can also describe their dynamics with a first order SDE and a PDE that encodes the evolution of the drift field. This is not strictly speaking a probabilistic representation of (39) since there is already a PDE involved. Our FBSDE approach may be viewed as a way to interpret in a pathwise sense the PDE governing the dirft in the h−transform representation.

The mean field Schrödinger problem and its equivalent formulations: proofs
In this part we complement the discussion undertaken in Section 1.2 and provide the proofs of the results stated therein. This section is organized into four subsections so that • Subsection 3.1 contains the proof of a more general version of Theorem 3.1, that is made with the help of several lemmas, • Subsection 3.2 is where we prove Proposition 1.1, Lemma 1.1 and Theorem 1.3.
• Theorem 1.2 is proven in subsection 3.3.
In the whole section, apart from subsection 3.1 that has its own assumptions, we always assume that (H1),(H2) are in force, even if we do not write them down explicitly in the statements of the lemmas and propositions.

A large deviations principle for particles interacting through their drifts
We consider for N ∈ N the interacting particle system where for some constant C > 0 and all (t, ω 1 , ω 2 ,ω 1 ,ω 2 ) ∈ [0, T ] × Ω 4 . Finally, regarding the measure µ in we assume that We stress that the usual theory of stochastic differential equations guarantees the strong existence and uniqueness for the above interacting particle system. Furthermore, if P ∈ P 1 (Ω) then the same arguments show that the stochastic differential equation admits a unique strong solution. We denote Γ(P) the law of X P . We can now state the main result of this part, which contains Theorem 1.1 as a very particular case.
The result is sharp, in that it fails for β = 2; see [53]. We follow Tanaka's reasoning [47] in order to establish this large deviations result. We remark that the assumption on exponential moments (46) is only used in the proof of Theorem 3.1, and not in the results preceding this proof.
For Q ∈ P β (Ω) we consider the equation t (ω) := ω 0 , Q ∈ P β (Ω), and consider the iterations Then a) The iteration is well-defined ω-by-ω (in particular, the Q-integrals are well-defined and finite) and in fact sup n E Q sup t≤T |Y n t | β is finite.
b) For each ω ∈ Ω the sequence {Y (n) (ω)} n∈N is convergent in the sup-norm to some limiting is adapted to the canonical filtration.
Proof. From the Lipschitz assumption on b we first derive Raising this to β, taking expectations and using Jensen's inequality, we derive where C ′ only depends on T and β. From this we establish for some R ≥ 0 that Observe that ∆ 1 t ≤ 2 sup s≤T |ω s − ω 0 | + C, so from the above inequality we obtain by induction that } n∈N is (for each ω) summable in n, so the same happens to {∆ n T } n∈N and therefore the uniform limit of the Y (n) exists for all ω. We denote by Y (∞) this limit. By Fatou's lemma E Q sup t≤T |Y Lemma 3.2. For any Q ∈ P β there exists a unique adapted continuous process satisfying (48) pointwise. Denoting Y Q this process, we further have Proof. If X and Y are solutions, then the Lipschitz assumption on b implies so from Grönwall we derive E Q sup s≤T |Y s − X s | = 0. With this, and using again the Lipschitz assumption on b, we find Thanks to this result we can define the operator where Y Q denotes the unique solution of (48).
Furthermore, if {X i,N : i ≤ N, N ∈ N} is the aforementioned interacting particle system, which is driven by {B i : i ∈ N} independent Brownian motions stared like µ in , then Proof. That Y R µ in is a solution to the McKean-Vlasov SDE is clear since ω is a Brownian motion under R µ in . That the solution is unique follows by observing that the drift in this SDE is Lipschitz jointly in Z and P = Law(Z), from where usual arguments apply. For the second point, consider first ω 1 , . . . , ω N continuous paths and define Q = 1 Replacing the deterministic paths ω 1 , . . . , ω N by those of B 1 , . . . , B N we conclude.
The key observation is that 1 N N i=1 δ B i satisfies a large deviations principle in P β (Ω) equipped with the W β topology, with good rate function given by the relative entropy H(·|R µ in ). This is true for β < 2 under our exponential moments assumption (46), but fails for β = 2, as follows easily from [53]. By Lemma 3. 3

we may derive, via the contraction principle ([21, Theorem 4.2.1]) a large deviations principle for
if we could only establish the continuity of Θ. This is our next step. Proof. We first prove the Lipschitz property. Let π be a coupling with first marginal Q and second marginal P. Denoting (ω,ω) the canonical process on Ω × Ω, and by the Lipschitz assumption on b, we have By Grönwall we have so taking infimum over such π we conclude that W β (Θ(Q), Θ(P)) ≤ K ′ W β (Q, P).
,ω)P(dω) ds, and the same holds forQ instead of Q. Denoting b(s, ω,ω)P(dω) ds, we therefore have We can now provide the proof of Theorem 3.1: Proof of Theorem 3.1. As we have observed, if {B i : i ∈ N} is and iid sequence of R µ in -distributed processes, then 1 N N i=1 δ B i · satisfies a large deviations principle in P β (Ω) equipped with the W β topology, with good rate function given by the relative entropy H(·|R µ in ). By (51), and since Θ : N ∈ N} satisfies a large deviations principle in P β (Ω) equipped with the W β topology. Since Θ is injective the good rate function is given bỹ In fact observe that if P ∈ range(Θ) and Θ −1 (P) ≪ R µ in then 5 P ≪ R µ in , sõ I (P) := H(Θ −1 (P)|R µ in ) if P ∈ range(Θ) and P ≪ R µ in +∞ otherwise.
Now take P ∈ range(Θ) and call Q = Θ −1 (P). It is immediate by the definition of Γ(·) that Γ(P) = R µ in • (Y Q ) −1 . On the other hand observe that the filtration generated by Y Q is equal to the canonical filtration: indeed Y Q is adapted and conversely and thereforeĨ (P) = H(P|Γ(P)) if P ∈ range(Θ) and P ≪ R µ in +∞ otherwise.
On the other hand call X P the unique adapted pointwise solution to which exists by Lemma 3.2 applied to the drift b(·, ·,ω)P(dω). As we recall in Lemma 5.4, X P and τ are P-a.s. inverses if P ≪ R µ in , since the above drift is Lipschitz. Now introduce Q := P • (τ ) −1 , so that Q • (X P ) −1 = P and in particular We have arrived atĨ To obtain the desired form (47) of the rate function it suffices to use Lemma 5.2 in the Appendix. 5 Let P = Θ(Q) for Q ≪ R µ in . The process Y Q satisfies pointwise dY Q t = dωt +b(t, Y Q )dt, whereb(t, y) = b(t, y, Y Q (ω))Q(dω). We have R µ in • (Y Q ) −1 ≪ R µ in since in fact their relative entropy is finite. Hence, if R µ in (A) = 0 then R µ in ((Y Q ) −1 (A)) = 0, and so Q ≪ R µ in implies Q((Y Q ) −1 (A)) = 0 therefore P(A) = 0 as desired.
Proof. We begin by proving continuity. Fix t, x and (t n , x n ) → (t, x). The sequence ∇W (x n − X tn ) converges pointwise to ∇W (x − X t ), since X is the (continuous) canonical process. By the fundamental theorem of calculus and (H1) we have |∇W (x n − X tn )| ≤ C 1 + C 2 sup s∈[0,T ] |X s |. Since P ∈ P 1 (Ω), we may use dominated convergence to conclude The space Lipschitzianity of −∇W * P t follows from (H1). Space differentiability follows similarly from (H1) and dominated convergence.
We will often make use of the next technical lemma, whose proof we defer to the appendix: for some C < +∞. DefineR as the law of the SDE and let P ∈ P(Ω) with X 0 ∼ µ. The following are equivalent (i) H(P|R) < +∞.
(ii) There exist a P-a.s. defined adapted process (ᾱ t ) t∈[0,T ] such that is a Brownian motion under P.

Moreover, under condition (ii)
and In particular under either (i) or (ii) we have P ∈ P 2 (Ω).
We turn to proving Lemma 1.1 stated in the introduction: Proof of Lemma 1.1. Define the vector fieldb(t, z) := −∇W * P t (z). Lemma 3.5 grants thatb fulfills the hypotheses of Lemma 3.6, giving the desired conclusions.
We can prove Proposition 1.1 of the introduction, concerning the existence of MFSBs. Recall the definition of Γ(P) and (MFSP) from the introduction.
Proof of Proposition 1.1. Let R µ in be the law of the Brownian motion started at µ in . (H2) grants that the classical Schrödinger problem (namely wrt. Brownian motion) is admissible. To see this, it suffices to verify that the coupling µ in ⊗ µ fin is admissible for the static version of the Schrödinger problem [36, Def 2.2] and then use the equivalence between the static and dynamic versions [36,Prop 2.3]. Therefore, there exist some P ∈ P(Ω) such that P 0 = µ in and H(P|R µ in ) < +∞. Lemma 3.6 (or its specialization Lemma 5.1 in the appendix) yields that P ∈ P 1 (Ω). On the other hand Lemma 5.2 in the appendix proves that for any P ∈ P(Ω) H(P|Γ(P)) < +∞ if and only if H(P|R µ in ) < +∞. Thus (MFSP) is admissible as well. Now observe that P → H(P|Γ(P)) is lower semicontinuous in P β (Ω), since on the one hand the relative entropy is jointly lower semicontinuous in the weak topology, and on the other hand Γ is readily seen to be continuous in P 1 (Ω). Recalling the definition of the operator Θ given in (50), to finish the proof we only need to justify that is relatively compact in P 1 (Ω) for each M , since the proof of Theorem 3.1 established 6 that  2 2 µ in (dx) < ∞. As Θ is per Lemma 3.4 Lipschitz in P 1 (Ω), it remains to prove that {H(Q|R γ ) ≤M } is W 1compact. This can be easily done by hand, or by invoking Sanov Theorem in the W 1 -topology for independent particles distributed according to R γ (see e.g. [53]), finishing the proof.
Proof of Theorem 1. 3 We split the proof into two propositions, namely Propositions 3.1 and 3.2. We begin by addressing the issue of Markovianity of the minimizers. Recall the definition of H −1 ((µ t ) t∈[0,T ] ) given under 'frequently used notation.' We rely strongly on the work [14] by Cattiaux and Léonard for the proof of the following result: where (α P t ) t∈[0,T ] is given in Lemma 1.1.
Proof. If P be optimal for (MFSP), then it is also optimal for inf {H(Q|Γ(P)) : Q ∈ P 1 (Ω), Q t = P t for all t ∈ [0, T ]} , since Γ(P) only depends on the marginals of P. The above problem is an instance of [14], ie. its optimizer is a so-called critical Nelson process. However, the drift of the path-measure Γ(P) may not fulfill the hypotheses in [14]. For this reason we need to make a slight detour. Let θ n ∈ C ∞ c ([0, T ] × R d ) and R n be defined as in Lemma 5.3 in the appendix, meaning that ∇θ n · (·) converges to −∇W * P t (z) in H −1 ((P t ) t∈[0,T ] ) and that R n is the law of For any n consider the problem Using [14, Lemma 3.1,Theorem 3.6] we obtain that for all n the unique optimizerP of (60) is the same for all n, and is such that is a Brownian motion underP. Lemma 3.5 grants that if we setb(t, z) = −∇W * P t (z) then the hypotheses of Lemma 3.6 are met. Since H(P|Γ(P)) < +∞, we derive from (57) therein that Using the implication (ii) ⇒ (i) of Lemma 3.6 we finally obtain that H(P|Γ(P)) < +∞ and therefore that we can use Lemma 5.3 for the choice Q =P therein. Now consider Q admissible for (59) and such that H(Q|Γ(P)) < +∞. Using Lemma 5. ThusP is also an optimizer for (59). But thenP = P since (59) can have at most one minimizer by strict convexity of the entropy and convexity of the admissible region. Combining (61) with (10) we get that t 0 −∇W * P s (X s ) + α P s − Φ s (X s ) ds is a continuous martingale with finite variation. But then it is constant P-a.s. The conclusion follows setting Ψ t (z) := Φ t (z) + ∇W * P t (z) and observing that ∇W * P · (·) ∈ H −1 ((P t ) r∈[0,T ] ).
Notice that the above proposition proves the first half of Theorem 1.3 from the introduction. We now establish the second half of this result: Proposition 3.2. Assume that P is optimal for (MFSP). Then Ψ t (X t ) has a continuous version adapted to the P-augmented canonical filtration, and the process (M t ) t∈[0,T ] defined by is a continuous martingale under P on [0, T [ and satisfies E P T 0 |M t | 2 dt < +∞.
To carry out the proof, we will use a well-known characterization of martingales (see e.g. [23]) which is as follows: an adapted process Proof. Define (M t ) t∈[0,T ] via (63). Using (H1), (9) and (57) for some C > 0. Consider now a process (h t ) t∈[0,T ] satisfying (65) and (66) and for ε > 0 define the shift transformation Under the current assumptions, τ ε admits an adapted inverse Y ε , i.e. there exists an adapted Indeed, since H(P|Γ(P)) < +∞, Lemma 5.2 in the appendix yields that P ≪ R µ in ; this entitles us to apply Lemma 5.4 in the same section, providing the existence of the inverse Y ε . If we set P ε = P • (τ ε ) −1 we have that P ε ∈ P 1 (Ω) is admissible for (MFSP), thanks to (65). Moreover, Lemma 1.1 and (68) imply that is a Brownian motion under P ε . Combining (9), (66) and (H1) we get that Lemma 3.5 grants thatb(t, x) = −∇W * P ε t (x) fulfills the hypothesis of Lemma 3.6 and (69) allows to use the implication (ii) ⇒ (i) which yields that H(P ε |Γ(P ε )) is finite and equals the left hand side of (69). Using the definition of P ε , we can rewrite H(P ε |Γ(P ε )) as Imposing optimality of P and letting ε to zero, using Taylor's expansion In the above equation, (X t ,h t ) t∈[0,T ] is an independent copy of (X t , h t ) t∈[0,T ] defined on some probability space (Ω,F ,P) andẼP denotes the expectation on (Ω,F ,P). Moreover, the exchange of limit and expectation is justified by (52), (9) and the dominated convergence theorem. Using the symmetry of W , and taking ±h, we can rewrite the latter condition as By integration by parts and the boundary condition (65), we arrive at proving the desired martingale property. By [43,Theorem IV.36.5] we know that a martingale in an augmented Brownian filtration admits a continuous version. Using again Lemma 5.2 we have that P ≪ R µ , and we so obtain a continuous version of our martingale (63), and a fortiori of Ψ t (X t ).

Benamou-Brenier formulation
We finally turn to the Benamou-Brenier formulation. Recall that C T (µ in , µ fin ) denotes the optimal value of the mean field Schrödinger problem. We define the set A as the collection of all absolutely continuous curves Recall from the introduction the problem (70) In (70), solutions to the continuity equation ∂ t µ t + ∇ · (w t µ t ) = 0 are meant in the weak sense.
Proof of Theorem 1.2. We first show that C T (µ in , µ fin ) ≥ C BB T (µ in , µ fin ). To this end, we may assume that the l.h.s. if finite and denote P an optimizer. As established in Theorem 1.3, the drift of X under P is equal to t 0 Ψ s (X s ) − ∇W * P s (X s )ds, where Ψ ∈ H −1 ((P t ) t∈[0,T ] ) and As we will see in Lemma 4.4 and Remark 4.1, the flow of marginals (P t ) t∈[0,T ] is absolutely continuous and its tangent velocity field v is given by Hence We conclude the desired inequality by noticing that ∇ log P t ∈ L 2 (dP t dt) and ∇W * P t ∈ L 2 (dP t dt). To wit, the first statement follows from [26, Thm 3.10] combined with Lemma 5.2 in our appendix, and the second from (57) used withb = −∇W * P t (z). We now establish C T (µ in , µ fin ) ≤ C BB T (µ in , µ fin ), so we may assume that (µ t ) t∈[0,T ] is feasible for the r.h.s. and leads to a finite value. Denote byṽ its tangent velocity field. We define Φ t (z) :=ṽ t (z) + 1 2 ∇ log µ t (z), so from the continuity equation for (µ t ) t∈[0,T ] we deduce the following equation in the distributional sense Observing that Φ ∈ H −1 ((µ t ) t∈[0,T ] ), we may apply the equivalence "(a) iff (c)" in [14, Theorem 3.4] 7 . We thus obtain a measure P whose marginals are exactly (µ t ) t∈[0,T ] , and by the uniqueness statement in [14,Theorem 3.4] we also know that the drift of X under P is precisely Φ s (X s ). Hence where the inequality follows from the equivalent expression of C T (µ in , µ fin ) given in (11).
We have proven C T (µ in , µ fin ) = C BB T (µ in , µ fin ), and the other statements follow from the previous arguments.

Schrödinger potentials and mean field PDE system: proofs
Proof of Corollary 1.1. We know by Theorem 1.1 that Ψ belongs to H((µ t )) t∈[0,T ] . The regularity hypothesis imposed on Ψ t (x) and µ t (x) allow us to conclude that Ψ is a true gradient, i.e. there exist ψ such that Ψ t (x) = ∇ψ t (x) for all (t, x) ∈ [0, T ] × R d . Lemma 1.1 together with Theorem 1.3 yield that µ t is a weak solution of the Fokker Planck equation in (18). Because of the regularity assumptions we made on Ψ and µ, we can conclude that µ t is indeed a classical solution. For the same reasons, we can turn the martingale condition (1.3) into the system of PDEs where L is the generator 1 2 ∆ + (∇(−W * µ t + ψ t )) · ∇. After some tedious but standard calculations we can rewrite the above as Since ψ is defined up to the addition of a function that depends on time only, the conclusion follows. Corollary 1.2 can be proven with a direct calculation using the definition of ϕ t and (18).

Convergence to equilibrium and functional inequalities: proofs
In this part we complement the discussion undertaken in Section 1.4 and provide proofs for the results stated therein. This section is organized as follows: • Subsections 4.1,4.2,4.3 are devoted to stating and proving some preparatory results that we shall use at different times in the proofs of the main results.
• In Subsection 4.4 we prove Theorem 1.7 and Theorem 1.6 together with its corollaries, ie. the Talagrand inequality (Corollary 1.3) and the HWI inequality (Corollary 1.4).
• Finally, in Subsection 4.5 we prove Theorem 1.4 and Theorem 1.5.
In all the lemmas and theorems in this subsection we always assume (H1)-(H2) to hold, and throughout P, α P , Ψ, M are as given in Theorem 1.3. We refer to Sections 1.2 and 1.4 for any unexplained notation.

Exponential upper bound for the corrector
Recall that we called Ψ the corrector. The goal of this part is to quantify the size of the corrector, as stated in Lemma 4.3 below. Before doing this we prove two preliminary lemmas. As usual, we denote by · the quadratic variation of a semimartingale.
Proof. We have shown at Theorem 3.2 that E P T 0 |M t | 2 dt < +∞ which gives that E P |M t | 2 < +∞ for almost every t ∈ [0, T [. But since E P |M t | 2 is an increasing function of t, we get E P |M t | 2 < +∞ for all t ∈ [0, T [. To complete the proof of (71) it suffices to observe that by definition of quadratic variation and since M t is an To prove the continuity of t → E P [ M t ] we start by observing that since M t is a continuous martingale, then M t has continuous and increasing paths. Thus, we obtain by monotone convergence that E P [ M t+h ] → E P [ M t ] as h ↓ 0, which gives the desired result. The proof of (72) follows from (63), the bounded Hessian of W (see(H1)) and the first part of Theorem 1.3.
Proof. Using the symmetry of W and the martingale property (63) it is easily derived that E P [Ψ t (X t )] is constant as a function of t. Therefore we get for all t ∈ [0, T ] Using again the symmetry of W we get that t 0 E P [∇W * P s (X s )]ds = 0, from which the conclusion follows.
We can now provide some key estimates on the corrector: Lemma 4.3. Assume (H1)-(H4). If P is an optimizer for (MFSP) and Ψ the associated corrector, then for any t ∈ (0, T ) we have Proof. Consider the function t → ϕ(t) defined by Fubini's Theorem allows to interchange the time integral and the expectation to get that ϕ is an absolutely continuous function with derivative From Itô's formula and Theorem 1.
We observe that the fact that M t is a martingale together with (72) and (71) make sure that E P t 0 Ψ r (X r ) · dM r = 0. Thus, taking expectation on both sides of the above equation yields Because of (72) we can use Fubini's Theorem and write where we used the symmetry of W to obtain the last expression. Plugging it back in (76) and using that t → ϕ ′ (t) is • increasing on [0, T [ since W is convex and the quadratic variation is an increasing process, we conclude that t → ϕ ′ (t) is absolutely continuous on the same interval. Moreover, using the κ-convexity of W and again the fact that the quadratic variation is an increasing process we get where to establish the last inequality we used that the hypothesis on µ in and µ fin together with Lemma 4.2 imply E P [Ψ t (X t )] = 0. The bound (73) follows by integrating the differential inequality (78) as done for instance in Lemma 5.5 in the Appendix, and observing that 1 2 E P T 0 |Ψ r (X r )| 2 dr = H(P|Γ(P)). To prove (74), we begin by observing that (78) also yields that Next, by definition of entropic cost we get the trivial bound The desired conclusion follows by plugging (79) in the above equation and some standard calculations.

First derivative of F
We compute the first derivative of F along the marginal flow of Q, assuming that H(Q|Γ(Q)) < +∞ and that Q is Markovian. To do this, we use an approach based on optimal transport, and some results of [26]. To be self-contained, we recall the basic notions of optimal transport we need to state the results. We refer to [1] for more details.
Tangent space Let µ ∈ P 2 (R d ). The tangent space Tan µ P 2 at µ is the closure in L 2 µ of ∇ψ; ψ ∈ C ∞ c (R d ) .
Absolutely continuous curves and velocity field Following [1, Th 8.3.1], we say that a curve (µ t ) t∈[0,T ] ⊆ P 2 (R d ) is absolutely continuous if there exists a Borel measurable vector field (t, z) → w t (z) such that • (w t ) t∈[0,T ] solves (in the sense of distributions) the continuity equation • w t satisfies the integrability condition Consider an absolutely continuous curve (µ t ) t∈[0,T ] . It is a consequence of the results in Chapter 8, and in particular of Proposition 8.4.5 of [1], that there exist a unique Borel measurable vector field v t (z) solving (80) and such that z → v t (z) belongs to the tangent space Tan µt P 2 for almost every t ∈ [0, T ]. We call such v t the (tangent) velocity field of (µ t ) t∈[0,T ] .
Remark 4.1. Let (µ t ) t∈[0,T ] be an absolutely continuous curve and w t (z) be in H −1 ((µ t ) t∈[0,T ] ). It is rather easy to see that z → w t (z) belongs to Tan µt P 2 for almost every t ∈ [0, T ].
Throughout the rest of the paper, if Q ∈ P(Ω) is such that H(Q|Γ(Q)) < +∞, we say that Q is Markovian if α Q t is σ(X t )-measurable for all t ∈ [0, T ], (α Q t ) t∈[0,T ] being defined by (9). In that case we write Ξ Q t (X t ) instead of α Q t .
Lemma 4.4. Let Q be such that H(Q|Γ(Q)) < +∞ and Markovian. Then is an absolutely continuous curve. Its tangent velocity field is given by Moreover, (ii) The function t → F (P t ) is absolutely continuous and Proof.
Proof of (i) To show that (Q t ) t∈[0,T ] is absolutely continuous it suffices to show that there exists a distributional solution of the continuity equation Let now ϕ ∈ C ∞ c (]0, T [×R d ). Using Itô's formula and taking expectation we obtain Lemma 5.2 in the appendix grants that under the current assumptions H(Q|R µ in ) < +∞, where R µ in is the Wiener measure started at µ in . But then, using [26, Thm 3.10] 8 we obtain that log Q t is an absolutely continuous function for almost every t and that (t, z) → ∇ log Q t (z) belongs to H −1 ((Q t ) t∈[0,T ] ). Therefore we can use integration by parts in (86) to obtain which gives, using the definition of the projection operator Π Qt , that the rhs of (81) solves the continuity equation in the sense of distributions. Next, we observe that (9)  Proof of (ii) From point (i) we know that z → ∇ log Q t (z) belongs to L 2 Qt for almost every t; this implies that ∇ log Q t +2∇W * Q t belongs to the subdifferential of F at Q t for almost every t (see e.g. [1,Thm. 10.4.13]). The chain rule [1, sec. E, pg. 233-234] gives the desired result (83), provided its hypothesis are verified. We have to check that (a) (Q t ) t∈[0,T ] is an absolutely continuous curve and F (Q t ) < +∞ for all t ∈ [0, T ], (b) F (·) is displacement λ-convex for some λ ∈ R, and (c) that To wit, (a) follows from point (i) and the fact that H(Q|Γ(Q)) < +∞, and (b) is a consequence of displacement convexity of the entropy and (H1). Finally, (c) is granted by (82) and the fact that ∇ log Q t (z) + 2∇W * Q t (z) belongs to H −1 ((Q t ) t∈[0,T ] ) (see the proof of (i)).

Time reversal
For Q ∈ P(Ω) the time reversalQ is the law of the time reversed process (X T −t ) t∈[0,T ] . In this section we derive an expression for H(Q|Γ(Q)) and use it to derive the bound (93) below, which plays a fundamental role in the proof of Theorem 1.6.
A consequence of Proposition 4.1 is that optimality is preserved under time reversal.

Functional inequalities: proofs
The goal of this section is to prove Theorem 1.7 and Theorem 1.6 together with its corollaries.

Proof of Theorem 1.6
Using a time reversal argument, we prove the bound (93) which is the key ingredient of the proof of Theorem 1.6 together with the bound for the correction term (73). Tthe backward correctorΨ is obtained by the same argument used in Proposition 4.1, replacing (MFSP) with (91) to obtain that there exist a Borel measurable vector fieldΨ t (z) ∈ H −1 ((P t ) t∈[0,T ] ) such that X t − t 0 −∇W * P s (X s ) +Ψ s (X s ) ds is a Brownian motion underP. Moreover, the following relation holdŝ Lemma 4.6. Assume (H3),(H4) and let P be an optimizer for (MFSP). Then hold for all r ∈ (0, T ).
Proof. Using (81) we can rewrite the above equation (92) aŝ From Proposition 4.5 we also know thatP is optimal for (91) and hence that H(P|Γ(P)) = C T (µ fin , µ in ). Therefore, by inverting again the roles of µ in and µ fin , we can use Lemma 4.3 for the problem (91) setting t = T − r to derive that Thanks to (95) we can write The bound (93) follows by plugging in (97) into (96) using the above equation, (87) and recalling that H(P|Γ(P)) = C T (µ in , µ fin ). The proof of (94) goes along the same lines: SinceP is optimal for (91) we also get from Lemma 4.3, and in particular from (74) for the choice t = T − r that Now the proof of Theorem 1.6 and its corollaries in the introduction is an easy task, given all the preparatory work.
Proof of Theorem 1.6. It amounts to add (73) and (93) with the choice r = t, and use the relation Proof of Corollary 1.3. It follows from Theorem 1.6 (Eq. (27)), observing that F (P t ) ≥ 0.
Proof of Corollary 1.4. Combining (92),(95), (34) we get that Using the above relation, Cauchy Schwartz inequality and the continuity of I F (P t ) in a neighborhood of 0, (83) we get that Consider now the bound (27). Observing that F (µ ∞ ) = 0, subtracting F (µ in ) on both sides, dividing by t, letting t → 0, using (98) and finally rearranging the resulting terms we get (31).

Proof of Theorem 1.7
We prove here Theorem 1.7 of the introduction. In the proof we will write which is used in the rest of the article. This is done in order to better deal with time reversal.
Proof of Theorem 1.7. Let M t be the martingale defined at (15). SinceP is optimal for (91), from Proposition 3.2 we get that )P s (dy) ds is an L 2 -martingale on [0, T [ underP. We define the stochastic processes We have, using the Markovianity of both P andP, that Taking conditional expectation w.r.t. σ(X [0,t] ) and using that both A t andÂ T −t are X [0,t] -measurable we get that the above expression equals Using the fact that W is symmetric and the definition of A t ,Â T −t , one easily obtains that the latter expression is worth 0. Indeed it holds that The proof that the function (34) is constant on (0, T ) is now concluded. In order to establish (35) we set t = T /2 in (34) and Cauchy Schwartz to get that |E P (µ in , µ fin )| ≤ E P [|Ψ T /2 (X T /2 )| 2 ]E P [|Ψ T /2 (X T /2 )| 2 ]  = E P [(∇ log P t (X t ) + 2∇W * P t (X t )) · v t (X t )] (92)+(95) Going along the same lines as the proof of Lemma 4.3 we get that both E P |Ψ(t, X)| 2 and E P |Ψ T −t (X T −t )| 2 are differentiable as functions of t in the open interval (0, T ). Moreover = κ E P |∇ log P t (X t ) + 2∇W * P t (X t )| 2 − 2κ E P (µ in , µ f in ).
The κ-convexity of W and the fact that the center of mass E P [X t ] is constant (see Lemma 4.2) allow to use the logarithmic Sobolev inequality [12, (ii), Thm 2.2] to obtain 9 κE P |∇ log P t (X t ) + 2∇W * P t (X t )| 2 ≥ 4κ 2 F (P t ).

Proof of Theorem 1.5
Proof of Theorem 1.5. Let P be optimal for (MFSP) and Ψ be given by Proposition 3.1. Then if we define B t := X t − t 0 ∇W * P s (X s ) + Ψ s (X s )ds the process (B t ) t∈[0,T ] is a Brownian motion under P. Since the McKean Vlasov SDE admits a unique strong solution, there exists a map Y : Ω −→ Ω such that Y • B · := Y satisfies Y 0 = X 0 (P − a.s.) and P − a.s.
Integrating the differential inequality and using that δ(0) = 0: ≤ t 1/2 2 exp(2κt)−1 exp(2κT )−1 C T (µ in , µ fin ) The conclusion follows from (28) and the observation that W 2 2 (P t , P MKV t ) ≤ δ(t). 9 Again, the apparent mismatch between the constant in the Log Sobolev inequality from [12] and the one we use here is due to the fact that in our definition ofF , there is no 1/2 in front of W .

Appendix
The following lemma is well known. For simple proofs see [35] or the appendix of [24].
(ii) There exist a P-a.s. defined adapted process (α t ) t∈[0,T ] such that E P T 0 |α t | 2 dt < +∞ and X t − t 0 α s ds is a Brownian motion under P.

Moreover, if (ii) holds, then
and With the help of Lemma 5.1 we can readily prove its generalization given in Lemma 3.6 concerning the case when R µ is replaced by the law of a diffusion with Lipschitz drift.
Under the current hypotheses on b, R µ andR are mutually absolutely continuous and dR µ dR = exp − T 0b (t, X t ) · dX t + 1 2 T 0 |b(t, X t )| 2 dt .
Therefore, using some standard calculations and (ii) we get Since R µ andR are mutually absolutely continuous and H(P|R µ ) < +∞ we get that P ≪R and H(P|R) = E P log dP dR µ + log dR µ dR (105) Thus, H(P|R) < +∞ if both H(P|R µ ) and E P log dR µ dR are finite. But this is indeed the case, thanks to (104),(103). The proof that (ii) ⇒ (i) is now complete. The desired form of the relative entropy follows by plugging in (103) and (104) into (105). Finally (57) follows from (102) and (101), and the last statement from (57).
Proof. By Lemma 5.2 we have H(Q|R µ in ) < +∞, where R µ in is the law of the Brownian motion with initial law µ in . Using implication (i) ⇒ (ii) from Lemma 5.1 and then (ii) ⇒ (i) together with (56) from Lemma 3.6 for the choiceb = ∇θ n , we get for all n ∈ N: