Key words. Optimal mass transport, Wasserstein distance, gradient flow, Schrödinger bridge.
|
|
- Jane Craig
- 5 years ago
- Views:
Transcription
1 ENTROPIC INTERPOLATION AND GRADIENT FLOWS ON WASSERSTEIN PRODUCT SPACES YONGXIN CHEN, TRYPHON GEORGIOU AND MICHELE PAVON Abstract. We show that the entropic interpolation between two given marginals provie by the Schröinger brige may be characterize as the curve in the Wasserstein space W which minimizes a suitable action. We also stuy the relative entropy evolution between an uncontrolle an a controlle ranom evolution as a graient flow on W W. Key wors. Optimal mass transport, Wasserstein istance, graient flow, Schröinger brige. AMS subject classifications. 47H07, 47H09, 60J5, 34A34, 49J0. Introuction. In the Schröinger brige problem (SBP) 3, one seeks the ranom evolution (a probability measure on path-space) which is closest in the relative entropy sense to a prior Markov iffusion evolution an has certain prescribe initial an final marginals µ an ν. As alreay observe by Schröinger 37, 38, the problem may be reuce to a static problem which, except for the cost, resembles the Kantorovich relaxe formulation of the optimal mass transport problem (OMT). Consiering that since (OMT) also has a ynamic formulation, we have two problems which amit equivalent static an ynamic versions 4. Moreover, in both cases, the solution entails a flow of one-time marginals joining µ an ν. The OMT yiels a isplacement interpolation flow whereas the SBP provies an entropic interpolation flow. Trough the work of Mikami, Mikami-Thieullen an Leonar 6, 7, 8, 3, 4, we know that the OMT may be viewe as a zero-noise limit of SBP when the prior is a sort of uniform measure on path space with vanishing variance. This connection has been extene to more general prior evolutions in 8, 9. Moreover, we know that, thanks to a very useful intuition by Otto 3, the isplacement interpolation flow {µ t ; 0 t } may be viewe as a constant-spee geoesic joining µ an ν in Wasserstein space 40. What can be sai from this geometric viewpoint of the entropic flow? It cannot be a geoesic, but can it be characterize as a curve minimizing a suitable action? In this paper, we show that this is inee the case resorting to a timesymmetric flui ynamic formulation of SBP. This characterization of the Schröinger brige answers as a byprouct a question pose by Carlen 4, pp Department of Mechanical Engineering, University of Minnesota, Minneapolis, Minnesota MN 55455, USA; chen468@umn.eu Department of Electrical an Computer Engineering, University of Minnesota, Minneapolis, Minnesota MN 55455, USA; tryphon@umn.eu Dipartimento i Matematica, Università i Paova, via Trieste 63, 35 Paova, Italy; pavon@math.unip.it Research partially supporte by the NSF, the AFOSR an by the University of Paova Research Project CPDA
2 SBP may be also formulate as a stochastic control problem with atypical bounary constraints. It is therefore interesting to compare the flow associate to the uncontrolle evolution (prior) to the optimal one. In particular, it is interesting to stuy the evolution of the relative entropy on the prouct Wasserstein space. The paper is outline as follows.... Elements of optimal mass transport theory. The literature on this problem is by now so vast an our egree of competence is such that we shall not even attempt here to give a reasonable an/or balance introuction to the various fascinating aspects of this theory. Fortunately, there exist excellent monographs an survey papers on this topic, see 35,, 40,, 4, 33, to which we refer the reaer. We shall only briefly review some concepts an results which are relevant for the topics of this paper... The static problem. Let µ an ν be probability measures on the measurable spaces X an Y, respectively. Let c : X Y 0, + ) be a measurable map with c(x, y) representing the cost of transporting a unit of mass from location x to location y. Let T µν be the family of measurable maps T : X Y such that T #µ = ν, namely such that ν is the push-forwar of µ uner T. Then Monge s optimal mass transport problem (OMT) is (.) inf c(x, T (x))µ(x). T T µν X Y As is well known, this problem may be unfeasible, namely the family T µν may be empty. This is never the case for the relaxe version of the problem stuie by Kantorovich in the 940 s (.) inf π Π(µ,ν) c(x, y)π(x, y) X Y where Π(µ, ν) are couplings of µ an ν, namely probability istributions on X Y with marginals µ an ν. Inee, Π(µ, ν) always contains the prouct measure µ ν. Let us specialize the Monge-Kantorovich problem (.) to the case X = Y = an c(x, y) = x y. Then, if µ oes not give mass to sets of imension n, by Brenier s theorem 40, p.66, there exists a unique optimal transport plan π (Kantorovich) inuce by a µ a.e. unique map T (Monge), T = ϕ, ϕ convex, an we have (.3) π = (I ϕ)#µ, ϕ#µ = ν. Among the extensions of this result, we mention that to strictly convex, superlinear costs c by Gangbo an McCann 5. The optimal transport problem may be use to introuce a useful istance between probability measures. Inee, let P ( ) be the set of probability measures µ on with finite secon moment. For µ, ν P ( ), the Kantorovich-Rubinstein istance, usually calle Wasserstein (Vasershtein) quaratic istance, is efine by ( (.4) W (µ, ν) = inf π Π(µ,ν) ) / x y π(x, y).
3 As is well known 40, Theorem 7.3, W is a bona fie istance. Moreover, it provies a most natural way to metrize weak convergence in P ( ) 40, Theorem 7.,, Proposition 7..5 (the same applies to the case p replacing with p everywhere). The Wasserstein space W is efine as the metric space ( P ( ), W ). It is a Polish space, namely a separable, complete metric space... The ynamic problem. So far, we have ealt with the static optimal transport problem. Nevertheless, in, p.378 it is observe that...a continuum mechanics formulation was alreay implicitly containe in the original problem aresse by Monge... Eliminating the time variable was just a clever way of reucing the imension of the problem. Thus, a ynamic version of the OMT problem was alreay in fieri in Gaspar Monge s 78 mémoire sur la théorie es éblais et es remblais! It was elegantly accomplishe by Benamou an Brenier in by showing that (.5a) W (µ, ν) = inf v(x, t) µ t (x)t, (µ,v) 0 (.5b) µ + (vµ) = 0, t (.5c) µ 0 = µ, µ = ν. Here the flow {µ t ; 0 t } varies over continuous maps from 0, to P ( ) an v over smooth fiels. In 4, Villani states at the beginning of Chapter 7 that two main motivations for the time-epenent version of OMT are a time-epenent moel gives a more complete escription of the transport; the richer mathematical structure will be useful later on. We can a three further reasons: it opens the way to establish a connection with the Schröinger brige problem (see Section 9 below), where the latter appears as a regularization of the former 6, 7, 8, 3, 4, 8, 9; it allows to view the optimal transport problem as an (atypical) optimal control problem 6-9. In some applications such as interpolation of images 0 or spectral morphing, the interpolating flow is essential! Let {µ t ; 0 t } an {v (x, t); (x, t) 0, } be optimal for (.5). Then µ t = ( t)i + t ϕ #µ, with T = ϕ solving Monge s problem, provies, in McCann s language, the isplacement interpolation between µ an ν. Then {µ t ; 0 t } may be viwe as a constant-spee geoesic joining µ an ν in Wasserstein space (Otto). This formally enows W with a pseuo Riemannian structure. McCann iscovere 5 that certain functionals are isplacement convex, namely convex along Wasserstein geoesics. µ k converges weakly to µ if fµ k fµ for every continuous, boune function f. 3
4 This has le to a variety of applications. Following one of Otto s main iscoveries, 3, it turns out that a large class of PDE s may be viewe as graient flows on the Wasserstein space W. This interpretation, because of the isplacement convexity of the functionals, is well suite to establish uniqueness an to stuy energy issipation an convergence to equilibrium. A rigorous setting in which to make sense of the Otto calculus has been evelope by Ambrosio, Gigli an Savaré for a suitable class of functionals. Convexity along geoesics in W also leas to new proofs of various geometric an functional inequalities 5, 40, Chapter 9. Finally, we mention that, when the space is not flat, qualitative properties of optimal transport can be quantifie in terms of how bouns on the Ricci-Curbastro curvature affect the isplacement convexity of certain specific functionals 4, Part II. The tangent space of P ( ) at a probability measure µ, enote by T µ P ( ) may be ientifie with the closure in L µ of the span of { ϕ : ϕ Cc }, where Cc is the family of smooth functions with compact support. It is naturally equippe with the scalar prouct of L µ. 3. The Fokker-Planck equation as a graient flow on Wasserstein space. Let us review the variational formulation of the Fokker-Planck equation as a graient flow on Wasserstein space, 40. Consier a physical system with phase space an with Hamiltonian H : x H(x) = E x. The thermoynamic states of the system are given by the family P( ) of probability istributions P on amitting ensity ρ. On P( ), we efine the internal energy as the expecte value of the Energy observable in state P (3.) U(H, ρ) = E P {H} = H(x) ρ(x)x = H, ρ. Let us also introuce the (ifferential) Gibbs entropy (3.) S(p) = k ln ρ(x)ρ(x)x, where k is Boltzmann s constant. S is strictly concave on P( ). Accoring to the Gibbsian postulate of classical statistical mechanics, the equilibrium state of a microscopic system at constant absolute temperature T an with Hamiltonian function H is necessarily given by the Boltzmann istribution law with ensity (3.3) ρ(x) = Z exp H(x) kt where Z is the partition function. Let us introuce the Free Energy functional F efine by (3.4) F (H, ρ, T ) := U(H, ρ) T S(ρ). Since S is strictly concave on S an U(E, ) is linear, it follows that F is strictly convex on the state space P( ). By Gibbs variational principle, the Boltzmann The letter Z was chosen by Boltzmann to inicate zustänige Summe (pertinent sum- here integral). 4
5 istribution ρ is a minimum point of the free energy F on P( ). Also notice that D(ρ ρ) = log ρ(x) R ρ(x) ρ(x)x N = S(ρ) + log Z + H(x)ρ(x)x = F (H, ρ, T ) + log Z. k kt R kt N Since Z oes not epen on ρ, we conclue that Gibb s principle is a trivial consequence of the fact that ρ minimizes D(ρ ρ) on D( ). Consier now an absolutely continuous curve µ t :, t W. Then, Chapter 8, there exist velocity fiel v t L µ t such that the following continuity equation hols on (0, T ) t µ t + (v t µ t ) = 0. Suppose µ t = x, so that the continuity equation (3.5) ρ t + (vρ) = 0 hols. We want to stuy the free energy functional F (H,, T ) or, equivalently, D( ρ), along the flow { ; t t }. Using (3.5), we get t D( ρ) = + log + R kt H(x) ρt t x N = + log + (3.6) kt H(x) (v )x. Integrating by parts, if the bounary terms at infinity vanish, we get (3.7) t D( ρ) = log + kt H(x) v x = log + kt H(x), v L. Thus, the Wasserstein graient of D( ρ) is The corresponing graient flow is (3.8) t = W D( ρ) = log + kt H(x). ( log + kt H(x) ) = kt H(x) +. But this is precisely the Fokker-Planck equation corresponing to the iffusion process (3.9) X t = kt H(X t)t + W t where W is a stanar n-imensional Wiener process. The process (3.9) has the Boltzmann istribution (3.3) as invariant ensity. Recall that, p.0 F (H,, T ) or, equivalently, D( ρ) are isplacement convex an have therefore a unique minimizer. 5
6 Remark. It seems worthwhile investigating to what extent the funamental assumption of statistical mechanics that the variables with longer relaxation time form a vector Markov process having (3.3) as invariant ensity is equivalent to the requirement that the flow of one-time ensities be a graient flow in Wasserstein space for the free energy. Let us finally plug the steepest escent (3.8) into (3.6). We get, after integrating by parts, the well known formula 6 (3.0) t D( ρ) = = + log + + log + = log ( ρt ρ kt H(x) ρt t x kt H(x) ) x. kt H(x) + x The last integral in (3.0) is sometimes calle the relative Fisher information of with respect to ρ 40, p Relative entropy as a functional on Wasserstein prouct spaces. Consier now two absolutely continuous curves µ t :, t W an µ t :, t W an their velocity fiels v t L µ t an ṽ t L µ t. Then, on (0, T ) (4.) (4.) t µ t + (v t µ t ) = 0, t µ t + (ṽ t µ t ) = 0. Let us suppose that µ t = (x)x an µ t = (x)x, for all t, t. Then (4.)-(4.) become (4.3) ρ + (vρ) = 0, t (4.4) ρ + (ṽ ρ) = 0, t where the fiels v an ṽ satisfy v(x, t) (x)x <, ṽ(x, t) (x)x <. The ifferentiability of the Wasserstein istance W (, ) has been stuie 4, Theorem 3.9. Consier instea the relative entropy functional on W W ( ) ( ) ρt ρ D( ) = h(, )x = log x, h( ρ, ρ) = log ρ. R ρ N t ρ Relative entropy functionals D( γ), where γ is a fixe probability measure (ensity), have been stuie as geoesically convex functionals on P ( ), see, Section 9.4. Our stuy of the evolution of D( ) is motivate by problems on a finite time interval such as the Schröinger brige problem (Section 6) an stochastic control 6
7 problems (Section 8) where it is important to evaluate relative entropy on two flows of marginals. We get t D( ) = (4.5) = h ρ ρ t + h ρ x ρ t ( ( + log log ) ( (ṽ ) + ρ ) t ( (v ) x After an integration by parts, assuming that the bounary terms at infinity vanish, we get ( ) t D( ρ ρt t ) = log ṽ ρ t v x R ρ N t ( ( )) log ρt (ṽ ) ρ ρt (4.6) = t x. ρt v Notice that the last expression looks like ( ( )) log ρt (ṽ ), ρt v L L. Thus, we ientify the graient of the functional D( ρ ρ) on W W as ( ) ( ( )) (4.7) W D( ρ ρ) log ρ ρ =. W D( ρ ρ) ρ Let us now compute the graient flow on W W corresponing to graient (4.7). We get ( ) ) ( ρt (4.8) log ρt ( ) t ρt = 0. Since ( ) ( ) ρt ρt J = log = = J, we observe the remarkable property that in the steepest escent (4.8) on the prouct Wasserstein space the fluxes are opposite an, therefore, ρ t = ρ t. If we plug the steepest escent (4.8) into (4.5), we get what appears to be a new formula ( t D( ) = + log log + ρ ) t ρ x R ρ N t t ( ) ( ) ρt = log ρt + x R ρ N t ( = + ρ ) ( ) t ρt (4.9) log x, 7
8 which shoul be compare to (3.0). Let us return to equation (4.6). By multiplying an iviing by in the last term of the mile expression, we get ( ) (4.0) t D( ρ ρt t ) = log (ṽ v) x which is precisely the expression obtaine in 34, Theorem III.. 5. Elements of Nelson-Föllmer kinematics of finite-energy iffusion processes. Let (Ω, F, P) be a complete probability space. A stochastic process {ξ(t); t t } is calle a finite-energy iffusion with constant iffusion coefficient σ I N if the paths ξ(ω) belong to C (, t ; R ) N (N-imensional continuous functions) an (5.) ξ(t) ξ(s) = t s β(τ)τ + σw + (t) W + (s), s < t t, where β(t) is at each time t a measurable function of the past {ξ(τ); τ t} an W is a stanar N-imensional Wiener process. Moreover, the rift β satisfies the finite energy conition { t } E β τ <. In, Föllmer has shown that a finite-energy iffusion also amits a reverse-time Ito ifferential. Namely, there exists a measurable function γ(t) of the future {ξ(τ); t τ t } calle backwar rift an another Wiener process W such that (5.) ξ(t) ξ(s) = Moreover, γ satisfies t s γ(τ)τ + σw (t) W (s), s < t t. { t } E γ τ <. Let us agree that t always inicate a strictly positive variable. For any function f :, t R let + f(t) = f(t + t) f(t) be the forwar increment at time t an let f(t) = f(t) f(t t) be the backwar increment at time t. For a finite-energy iffusion, Föllmer has also shown in that forwar an backwar rifts may be obtaine as Nelson s conitional erivatives 9 { } + ξ(t) β(t) = lim E ξ(τ), τ t, t 0 t { } ξ(t) γ(t) = lim E ξ(τ), t τ t, t 0 t the limits being taken in L N (Ω, F, P). It was finally shown in that the one-time probability ensity ( ) of ξ(t) (which exists for every t. ) is absolutely continuous on an the following uality relation hols t > 0 (5.3) E {β(t) γ(t) ξ(t)} = σ log ρ(ξ(t), t), a.s.. 8
9 Let us introuce the fiels b + (x, t) = E {β(t) ξ(t) = x}, b (x, t) = E {γ(t) ξ(t) = x}. Then, Ito s rule for the forwar an backwar ifferential of ξ imply that satisfies the two Fokker-Planck equations (5.4) ρ t + (b +ρ) σ ρ = 0, (5.5) ρ t + (b ρ) + σ ρ = 0. Following Nelson, let us introuce the current an osmotic rift of ξ by (5.6) v(t) = β(t) + γ(t), u(t) = β(t) γ(t), respectively. Clearly v is similar to the classical velocity, whereas u is the velocity ue to the noise which tens to zero when σ tens to zero. Let us also introuce v(x, t) = E {v(t) ξ(t) = x} = b +(x, t) + b (x, t). Then, combining (5.4) an (5.5), we get (5.7) ρ + (vρ) = 0, t which has the form of a continuity equation expressing conservation of mass. When ξ is Markovian with β(t) = b + (ξ(t), t) an γ(t) = b (ξ(t), t), (5.3) reuces to Nelson s relation (5.8) b + (x, t) b (x, t) = σ log (x). Then (5.7) hols with (5.9) v(x, t) = b + (x, t) σ log (x). 6. Schröinger briges an entropic interpolation. Let Ω = C(, t ; ) be the space of value continuous functions. Let Wx σ enote Wiener measure on Ω with variance σ I N starting at the point x at time. If, instea of a Dirac measure concentrate at x, we give the volume measure as initial conition, we get the unboune measure on path space W σ = W σ x x. It is a useful tool to introuce the family of istributions P on Ω which are equivalent to it. Let P P represent an a priori ranom evolution an let Q P. Then, we have D(Q P ) = E Q log Q t (6.a) = D(q 0 p 0 ) + E Q P σ βq β P t (6.b) = D(q p ) + E Q t σ γq γ P t Here q 0, q are the marginal ensities of Q at an t, respectively. Similarly, p 0, p are the marginal ensities of P. Let ρ 0 an ρ be two probability ensities. Let 9.
10 P(ρ 0, ρ ) enote the set of istributions in P having the prescribe marginal ensities at an t. Given P P, we consier the following problem: (6.) Minimize D(Q P ) over Q P(ρ 0, ρ ). Conitions for existence an uniqueness for this problem an properties of the minimizing measure have been stuie by many authors, most noticeably by Fortet, Beurlin, Jamison an Föllmer 4, 3, 0, 3. If there is at least one Q in D(ρ 0 ρ ) such that D(Q P ) <, there exists a unique minimizer Q in P(ρ 0 ρ ) calle the Schröinger brige from ρ 0 to ρ over P 3. Existence is guarantee by uner conitions on P, ρ 0 an ρ, see 4, Proposition.5. We shall tacitly assume henceforth that they are satisfie so that Q is well efine. In view of Sanov s theorem 36, solving the maximum entropy problem (6.) is equivalent to a problem of large eviations of the empirical istribution as showe by Föllmer 3 recovering Schröinger s original motivation. Using (5.6) an observing that D(q 0 p 0 ) an D(q p ) are constant for Q varying in P(ρ 0 ρ ), we get that problem (6.) is equivalent to (6.3) t Minimize E Q σ vq v P + σ uq u P t over Q P(ρ 0, ρ ). Suppose now that the prior measure P is Markovian. In this case, the classical results of Jamison 0 imply that the solution of (6.) or, equivalently, of (6.3) is also Markovian. Thus, we can restrict our search to P M (ρ 0, ρ ), namely Markovian measures in P(ρ 0, ρ ). Taking (5.8) into account, we can rewrite (6.3) in the form (6.4) Minimize + σ σ log ρq t ρ P t E Q t (ξ(t), t) t σ vq (ξ(t), t) v P (ξ(t), t) over Q P M (ρ 0 ρ ). Finally, in view of (5.7), we get the following equivalent flui ynamic formulation t (6.5a) Minimize (ρ,v) R σ v(x, t) vp (x, t) N + σ ρ log (x, t) ρ(x, t)xt, 8 ρp (6.5b) ρ + (vρ) = 0, t (6.5c) ρ = 0 0, ρ =. Comparing (6.5) to the OMT with prior formulate an stuie in 8, see also 9 for the Gauss-Markov case, we see that the essential ifference is that there is here an extra term in the action functional which has the form of a relative Fisher information of with respect to the prior one-time ensity ρ P t (Dirichlet form) 40, p.78 integrate over time. Also notice that, in view of (4.7), we have (6.6) log ρ P t = W D( ρ P t ). 0
11 To fin the connection to the classical OMT, let us specialize to the situation where the prior P = W σ. In that case, v P = u P = 0 an, multiplying the criterion by σ, we get the problem t (6.7a) Minimize (ρ,v) R v(x, t) + σ4 8 log (x) (x)xt, N (6.7b) ρ + (vρ) = 0, t (6.7c) ρ = 0 0, ρ =. Again, specializing (6.6) to the case ρ P t (x), we get (6.8) log = W k S(). If σ 0 it appears that in the limit we get the Benamou-Brenier formulation of OMT (.5). This is inee the case, see 6, 7, 8, 3, 4 an 8, 9 for the case with prior. 7. Optimal transport an Nelson s stochastic mechanics. There has been some interest in connecting optimal transport with Nelson s stochastic mechanics 4, 4, p.707 or irectly with the Schröinger equation 4. Consier the case of a free, non-relativistic particle of mass m. Then, a variational principle leaing to the Schröinger equation, can be base on the Guerra-Morato action functional 5 which, in flui ynamic form, is (7.) where A GM (, t ) = = t t m v(x, t) u(x, t) (x)xt m v(x, t) (x)x 8m I() RN ρ (7.) I(ρ) = x ρ t is the Fisher information of ρ since, for the Nelson process, σ = /m. Instea, the Yasue action 44 in flui-ynamic form is (7.3) A Y (, t ) = = t t m v(x, t) + u(x, t) (x)xt m v(x, t) (x)x + 8m I() t. In 4, p.3, Eric Carlen poses the question of minimizing the Yasue action subject to the continuity equation (6.7b) for given initial an final marginals (6.7c) ) stating that...the Euler-Lagrange equations for it are not easy to unerstan. In view of Section 6, we alreay know the solution to this problem: It is provie by the current velocity an the flow of one-time ensities of the Schröinger brige with (6.7c) an stationary Wiener measure as a prior.
12 8. Relative entropy prouction for controlle evolution. Consier on, t a finite-energy Markov process taking values in with forwar Ito ifferential (8.) ξ = b + (ξ(t), t)t + σw +. Let (x) be the probability ensity of ξ(t). Consier also the feeback controlle process ξ u with forwar ifferential (8.) ξ u = b + (ξ u (t), t)t + u(ξ u (t), t)t + σw +. Here the control u is aapte to the past an is such that ξ u is a finite-energy iffusion. Let ρ u t (x) be the probability ensity of ξ u (t). We are intereste in the evolution of D(ρ u t ). By (5.9)-(5.7), the ensities satisfy (8.3) (8.4) ρ u t + (vu ρ u ) = 0, By (4.0), we now get t D(ρu t ) = (8.5) = ρ + (vρ) = 0, t v(x, t) = b +(x, t) σ log (x) v u (x, t) = b + (x, t) + u(x, t) σ log ρu t (x). ( ρ u log t ( ρ u log t ) (v u v) ρ u t x ) ( )) (u σ ρ u log t ρ u t x. Suppose now ρ u t = ρ 0 t is also uncontrolle an iffers from only because of the initial conition at t =. Then (8.5) gives the well known formula generalizing (3.0) ( ) ( ) (8.6) t D(ρ0 t ) = σ ρ 0 log t ρ 0 log t ρ u t x R ρ N t which shows that two solutions of the same Fokker-Plank equation ten to get closer. 9. Schröinger briges as controlle evolution. Consier now the situation where ξ(t) represents a prior evolution on, t an the controlle evolution ξ u = ξ is the solution of the Schröinger brige problem for a pair of initial an final marginals ρ 0 an ρ 3, 43. Then the ifferential of ξ is given by (9.) ξ = b + (ξ (t), t)t + σ log ϕ(ξ (t), t)t + σw + where ϕ is space-time harmonic for the prior evolution, namely it satisfies (9.) ϕ t + b + ϕ + σ ϕ = 0. Let ρ ϕ be the ensity of ξ. Let us single out a special case of the Schröinger brige problem where relative entropy on path space is minimise uner the only constraint that the final marginal ensity be ρ. In such case, we have ρ ϕ t (x) = (x) ϕ(x, t).
13 Then (8.5) gives (9.3) t D(ρϕ t ) = σ log ϕ log ϕ ρ ϕ t x. This shows that D(ρ ϕ t ) increases up to time t = t. It represents the intuitive fact that the brige evolution has to be as close as possible to the prior but the final value of the relative entropy must be the positive quantity D(ρ ). Thus, D(ρ ϕ t ) approaches this positive quantity from below. Result (9.3) may be viewe as a reverse-time H-theorem, as the brige an the reference evolution have the same backwar rift 3. REFERENCES L. Ambrosio, N. Gigli an G. Savaré, Graient Flows in Metric Spaces an in the Space of Probability Measures, Lectures in Mathematics ETH Zürich, Birkhäuser Verlag, Basel, n e J. Benamou an Y. Brenier, A computational flui mechanics solution to the Monge- Kantorovich mass transfer problem, Numerische Mathematik, 84, no. 3, 000, A. Beurling, An automorphism of prouct measures, Ann. Math. 7 (960), E. Carlen, Stochastic mechanics: A look back an a look ahea, in Diffusion, quantum theory an raically elementary mathematics, W. G. Faris e., Mathematical Notes 47, Princeton University Press, 006, Y. Chen an T.T. Georgiou, Stochastic briges of linear systems, preprint, abs/407.34, IEEE Trans. Aut. Control, to appear. 6 Y. Chen, T.T. Georgiou an M. Pavon, Optimal steering of a linear stochastic system to a final probability istribution,, Part I, Aug. 04, IEEE Trans. Aut. Control, to appear. 7 Y. Chen, T.T. Georgiou an M. Pavon, Optimal steering of a linear stochastic system to a final probability istribution, part II, Oct. 04, IEEE Trans. Aut. Control, to appear. 8 Y. Chen, T.T. Georgiou an M. Pavon, On the relation between optimal transport an Schröinger briges: A stochastic control viewpoint Dec. 04, arxiv :4.4430v, submitte for publication. 9 Y. Chen, T.T. Georgiou an M. Pavon, Optimal transport over a linear ynamical system, Feb. 05, arxiv:50.065v, submitte for publication. 0 Y. Chen, T.T. Georgiou an M. Pavon, Entropic an isplacement interpolation: a computational approach using the Hilbert metric, June 05, arxiv: v, submitte for publication. L. C. Evans, Partial ifferential equations an Monge-Kantorovich mass transfer, in Current evelopments in mathematics, 977, Int. Press, Boston, 999, H. Föllmer, in Stochastic Processes - Mathematics an Physics, Lecture Notes in Mathematics (Springer-Verlag, New York,986), Vol. 58, pp H. Föllmer, Ranom fiels an iffusion processes, in: Ècole Ètè e Probabilitès e Saint- Flour XV-XVII, eite by P. L. Hennequin, Lecture Notes in Mathematics, Springer- Verlag, New York, 988, vol.36, R. Fortet, Résolution un système equations e M. Schröinger, J. Math. Pure Appl. IX (940), W. Gangbo an R. McCann, The geometry of optimal transportation. Acta Mathematica, 77 (996), no., 3?6. 6 R. Graham, Path integral methos in nonequilibrium thermoynamics an statistics, in Stochastic Processes in Nonequilibrium Systems, L. Garrio, P. Seglar an P.J.Shepher Es., Lecture Notes in Physics 84, Springer-Verlag, New York, 978, F. Guerra an L. Morato, Quantization of ynamical systems an stochastic control theory, Physical Review D, 7 (8), U.G.Haussmann an E.Paroux, Time reversal of iffusions, The Annals of Probability 4, 986, D.B.Hernanez an M. Pavon, Equilibrium escription of a particle system in a heat bath, Acta Applicanae Mathematicae 4 (989),
14 0 B. Jamison, The Markov processes of Schröinger, Z. Wahrscheinlichkeitstheorie verw. Gebiete 3 (975), X. Jiang, Z. Luo, an T. Georgiou, Geometric methos for spectral analysis, IEEE Transactions on Signal Processing, 60, no. 3, 0, R. Joran, D. Kinerlehrer, an F. Otto, The variational formulation of the Fokker-Planck equation, SIAM J. Math. Anal. 9: -7 (998). 3 C. Léonar, A survey of the Schroeinger problem an some of its connections with optimal transport, Discrete Contin. Dyn. Syst. A, 04, 34 (4): C. Léonar, From the Schröinger problem to the Monge-Kantorovich problem, J. Funct. Anal., 0, 6, R. McCann, A convexity principle for interacting gases, Av. Math., 8 (997),, T. Mikami, Monge s problem with a quaratic cost by the zero-noise limit of h-path processes, Probab. Theory Relat. Fiels, 9, (004), T. Mikami an M. Thieullen, Duality theorem for the stochastic optimal control problem., Stoch. Proc. Appl., 6, 85?835 (006). 8 T. Mikami an M. Thieullen, Optimal Transportation Problem by Stochastic Optimal Control, SIAM Journal of Control an Optimization, 47, N. 3, 7-39 (008). 9 E. Nelson, Dynamical Theories of Brownian Motion, Princeton University Press, Princeton, L. Ning, X. Jiang an T. Georgiou, On the geometry of covariance matrices, IEEE Signal Processing Letters, vol. 0, Augus3, L. Ning, T.T. Georgiou, an A. Tannenbaum, Matrix-value Monge-Kantorovich Optimal Mass Transport, IEEE Trans. Aut. Contr., to appear. 3 F. Otto, The geometry of issipative evolution equations: the porous meium equation, Comm. Partial Differential Equations 6,? (00), 0? Y. Ollivier, H. Pajot an C. Villani, Optimal Transportation. Theory an Applications, Lonon Mathematical Society Lecture Notes Series 43, Cambrige University Press, M. Pavon an F. Ticozzi, On entropy prouction for controlle Markovian evolution, J. Math. Phys., 47, 06330, oi:0.063/.0776 (006). 35 T. Rachev an L. Rüschenorf, Mass Transportation Problems, Vol. I: Theory, Vol. II: Applications. Probability an its applications. Springer- Verlag, New York, I. S. Sanov, On the probability of large eviations of ranom magnitues (in Russian), Mat. Sb. N. S., 4 (84) (957) 44. Select. Transl. Math. Statist. Probab.,, 3-44 (96). 37 E. Schröinger, Über ie Umkehrung er Naturgesetze, Sitzungsberichte er Preuss Aka. Wissen. Berlin, Phys. Math. Klasse (93), E. Schröinger, Sur la théorie relativiste e l electron et l interpretation e la mécanique quantique, Ann. Inst. H. Poincaré, 69 (93). 39 E. Tannenbaum, T. Georgiou an A. Tannenbaum, Signals an Control Aspects of Optimal Mass Transport an the Boltzmann Entropy, Proc. 49th IEEE-CDC conference, Atlanta, GA, 00, C. Villani, Topics in optimal transportation, AMS, 003, vol C. Villani, Optimal transport. Ol an new. Grunlehren er Mathematischen Wissenschaften 338. Springer-Verlag, M. K. von Renesse, An optimal transport view on Schröinger s equation, Cana. Math. Bull, 55 (4), 0, A. Wakolbinger, Schroeinger briges from 93 to 99, in Proc. of the 4th Latin American Congress in Probability an Mathematical Statistics, Mexico City 990, Contribuciones en probabilia y estaistica matematica, 3 (99), K. Yasue, Stochastic calculus of variations, J. Funct. Analysis, 4, (98). 4
WUCHEN LI AND STANLEY OSHER
CONSTRAINED DYNAMICAL OPTIMAL TRANSPORT AND ITS LAGRANGIAN FORMULATION WUCHEN LI AND STANLEY OSHER Abstract. We propose ynamical optimal transport (OT) problems constraine in a parameterize probability
More informationSecond order differentiation formula on RCD(K, N) spaces
Secon orer ifferentiation formula on RCD(K, N) spaces Nicola Gigli Luca Tamanini February 8, 018 Abstract We prove the secon orer ifferentiation formula along geoesics in finite-imensional RCD(K, N) spaces.
More informationThe dynamics of Schrödinger bridges
Stochastic processes and statistical machine learning February, 15, 2018 Plan of the talk The Schrödinger problem and relations with Monge-Kantorovich problem Newton s law for entropic interpolation The
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationA new proof of the sharpness of the phase transition for Bernoulli percolation on Z d
A new proof of the sharpness of the phase transition for Bernoulli percolation on Z Hugo Duminil-Copin an Vincent Tassion October 8, 205 Abstract We provie a new proof of the sharpness of the phase transition
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationIntroduction to Markov Processes
Introuction to Markov Processes Connexions moule m44014 Zzis law Gustav) Meglicki, Jr Office of the VP for Information Technology Iniana University RCS: Section-2.tex,v 1.24 2012/12/21 18:03:08 gustav
More informationSome geometric consequences of the Schrödinger problem
Some geometric consequences of the Schrödinger problem Christian Léonard To cite this version: Christian Léonard. Some geometric consequences of the Schrödinger problem. Second International Conference
More informationEntropic and Displacement Interpolation: Tryphon Georgiou
Entropic and Displacement Interpolation: Schrödinger bridges, Monge-Kantorovich optimal transport, and a Computational Approach Based on the Hilbert Metric Tryphon Georgiou University of Minnesota IMA
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationThe continuity equation
Chapter 6 The continuity equation 61 The equation of continuity It is evient that in a certain region of space the matter entering it must be equal to the matter leaving it Let us consier an infinitesimal
More informationProblem set 2: Solutions Math 207B, Winter 2016
Problem set : Solutions Math 07B, Winter 016 1. A particle of mass m with position x(t) at time t has potential energy V ( x) an kinetic energy T = 1 m x t. The action of the particle over times t t 1
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationMean field games via probability manifold II
Mean field games via probability manifold II Wuchen Li IPAM summer school, 2018. Introduction Consider a nonlinear Schrödinger equation hi t Ψ = h2 1 2 Ψ + ΨV (x) + Ψ R d W (x, y) Ψ(y) 2 dy. The unknown
More informationSINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
Communications on Stochastic Analysis Vol. 2, No. 2 (28) 289-36 Serials Publications www.serialspublications.com SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationOn the relation between optimal transport and Schrödinger bridges: A stochastic control viewpoint
DECEMBER 4, 4 On the relation between optimal transport and Schrödinger bridges: A stochastic control viewpoint arxiv:4.443v [cs.sy 5 Dec 4 Yongxin Chen, Tryphon Georgiou and Michele Pavon Abstract We
More informationWELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL PRESSURE IN SOBOLEV SPACES
Electronic Journal of Differential Equations, Vol. 017 (017), No. 38, pp. 1 7. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationConvergence of Langevin MCMC in KL-divergence
Convergence of Langevin MCMC in KL-ivergence Xiang Cheng x.cheng@berkeley.eu an Peter Bartlett peter@berkeley.eu Eitor: Abstract Langevin iffusion is a commonly use tool for sampling from a given istribution.
More informationOn the number of isolated eigenvalues of a pair of particles in a quantum wire
On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in
More informationLocal semiconvexity of Kantorovich potentials on non-compact manifolds
Local semiconvexity of Kantorovich potentials on non-compact manifolds Alessio Figalli, Nicola Gigli Abstract We prove that any Kantorovich potential for the cost function c = d / on a Riemannian manifold
More informationFrom slow diffusion to a hard height constraint: characterizing congested aggregation
3 2 1 Time Position -0.5 0.0 0.5 0 From slow iffusion to a har height constraint: characterizing congeste aggregation Katy Craig University of California, Santa Barbara BIRS April 12, 2017 2 collective
More informationTOWARDS THERMOELASTICITY OF FRACTAL MEDIA
ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More informationConvergence of Random Walks
Chapter 16 Convergence of Ranom Walks This lecture examines the convergence of ranom walks to the Wiener process. This is very important both physically an statistically, an illustrates the utility of
More informationGradient flow of the Chapman-Rubinstein-Schatzman model for signed vortices
Graient flow of the Chapman-Rubinstein-Schatzman moel for signe vortices Luigi Ambrosio, Eoaro Mainini an Sylvia Serfaty Deicate to the memory of Michelle Schatzman (1949-2010) Abstract We continue the
More informationarxiv:nlin/ v1 [nlin.cd] 21 Mar 2002
Entropy prouction of iffusion in spatially perioic eterministic systems arxiv:nlin/0203046v [nlin.cd] 2 Mar 2002 J. R. Dorfman, P. Gaspar, 2 an T. Gilbert 3 Department of Physics an Institute for Physical
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationA family of functional inequalities: lojasiewicz inequalities and displacement convex functions
17-66 A family of functional inequalities: lojasiewicz inequalities an isplacement convex functions Arien Blanchet an Jérôme Bolte A FAMILY OF FUNCTIONAL INEQUALITIES: LOJASIEWICZ INEQUALITIES AND DISPLACEMENT
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationarxiv: v1 [math-ph] 5 May 2014
DIFFERENTIAL-ALGEBRAIC SOLUTIONS OF THE HEAT EQUATION VICTOR M. BUCHSTABER, ELENA YU. NETAY arxiv:1405.0926v1 [math-ph] 5 May 2014 Abstract. In this work we introuce the notion of ifferential-algebraic
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationDiscrete Ricci curvature via convexity of the entropy
Discrete Ricci curvature via convexity of the entropy Jan Maas University of Bonn Joint work with Matthias Erbar Simons Institute for the Theory of Computing UC Berkeley 2 October 2013 Starting point McCann
More informationA nonlinear inverse problem of the Korteweg-de Vries equation
Bull. Math. Sci. https://oi.org/0.007/s3373-08-025- A nonlinear inverse problem of the Korteweg-e Vries equation Shengqi Lu Miaochao Chen 2 Qilin Liu 3 Receive: 0 March 207 / Revise: 30 April 208 / Accepte:
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationAbstract A nonlinear partial differential equation of the following form is considered:
M P E J Mathematical Physics Electronic Journal ISSN 86-6655 Volume 2, 26 Paper 5 Receive: May 3, 25, Revise: Sep, 26, Accepte: Oct 6, 26 Eitor: C.E. Wayne A Nonlinear Heat Equation with Temperature-Depenent
More informationSurvival exponents for fractional Brownian motion with multivariate time
Survival exponents for fractional Brownian motion with multivariate time G Molchan Institute of Earthquae Preiction heory an Mathematical Geophysics Russian Acaemy of Science 84/3 Profsoyuznaya st 7997
More informationAlmost everywhere well-posedness of continuity equations with measure initial data
Almost everywhere well-poseness of continuity equations with measure initial ata Luigi Ambrosio Alessio Figalli Abstract The aim of this note is to present some new results concerning almost everywhere
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationREVERSIBILITY FOR DIFFUSIONS VIA QUASI-INVARIANCE. 1. Introduction We look at the problem of reversibility for operators of the form
REVERSIBILITY FOR DIFFUSIONS VIA QUASI-INVARIANCE OMAR RIVASPLATA, JAN RYCHTÁŘ, AND BYRON SCHMULAND Abstract. Why is the rift coefficient b associate with a reversible iffusion on R given by a graient?
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationLower bounds on Locality Sensitive Hashing
Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,
More informationThe effect of dissipation on solutions of the complex KdV equation
Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,
More informationDiscrete Operators in Canonical Domains
Discrete Operators in Canonical Domains VLADIMIR VASILYEV Belgoro National Research University Chair of Differential Equations Stuencheskaya 14/1, 308007 Belgoro RUSSIA vlaimir.b.vasilyev@gmail.com Abstract:
More informationSOME RESULTS ON THE GEOMETRY OF MINKOWSKI PLANE. Bing Ye Wu
ARCHIVUM MATHEMATICUM (BRNO Tomus 46 (21, 177 184 SOME RESULTS ON THE GEOMETRY OF MINKOWSKI PLANE Bing Ye Wu Abstract. In this paper we stuy the geometry of Minkowski plane an obtain some results. We focus
More informationON ISENTROPIC APPROXIMATIONS FOR COMPRESSIBLE EULER EQUATIONS
ON ISENTROPIC APPROXIMATIONS FOR COMPRESSILE EULER EQUATIONS JUNXIONG JIA AND RONGHUA PAN Abstract. In this paper, we first generalize the classical results on Cauchy problem for positive symmetric quasilinear
More informationAsymptotic estimates on the time derivative of entropy on a Riemannian manifold
Asymptotic estimates on the time erivative of entropy on a Riemannian manifol Arian P. C. Lim a, Dejun Luo b a Nanyang Technological University, athematics an athematics Eucation, Singapore b Key Lab of
More informationA LIMIT THEOREM FOR RANDOM FIELDS WITH A SINGULARITY IN THE SPECTRUM
Teor Imov r. ta Matem. Statist. Theor. Probability an Math. Statist. Vip. 81, 1 No. 81, 1, Pages 147 158 S 94-911)816- Article electronically publishe on January, 11 UDC 519.1 A LIMIT THEOREM FOR RANDOM
More informationarxiv: v4 [math.pr] 27 Jul 2016
The Asymptotic Distribution of the Determinant of a Ranom Correlation Matrix arxiv:309768v4 mathpr] 7 Jul 06 AM Hanea a, & GF Nane b a Centre of xcellence for Biosecurity Risk Analysis, University of Melbourne,
More informationA transmission problem for the Timoshenko system
Volume 6, N., pp. 5 34, 7 Copyright 7 SBMAC ISSN -85 www.scielo.br/cam A transmission problem for the Timoshenko system C.A. RAPOSO, W.D. BASTOS an M.L. SANTOS 3 Department of Mathematics, UFSJ, Praça
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationarxiv: v1 [physics.class-ph] 20 Dec 2017
arxiv:1712.07328v1 [physics.class-ph] 20 Dec 2017 Demystifying the constancy of the Ermakov-Lewis invariant for a time epenent oscillator T. Pamanabhan IUCAA, Post Bag 4, Ganeshkhin, Pune - 411 007, Inia.
More informationORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS. Gianluca Crippa
Manuscript submitte to AIMS Journals Volume X, Number 0X, XX 200X Website: http://aimsciences.org pp. X XX ORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS Gianluca Crippa Departement Mathematik
More informationApproximate reduction of dynamic systems
Systems & Control Letters 57 2008 538 545 www.elsevier.com/locate/sysconle Approximate reuction of ynamic systems Paulo Tabuaa a,, Aaron D. Ames b, Agung Julius c, George J. Pappas c a Department of Electrical
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationarxiv: v1 [math-ph] 2 May 2016
NONLINEAR HEAT CONDUCTION EQUATIONS WITH MEMORY: PHYSICAL MEANING AND ANALYTICAL RESULTS PIETRO ARTALE HARRIS 1 AND ROBERTO GARRA arxiv:165.576v1 math-ph] May 16 Abstract. We stuy nonlinear heat conuction
More informationSYMPLECTIC GEOMETRY: LECTURE 3
SYMPLECTIC GEOMETRY: LECTURE 3 LIAT KESSLER 1. Local forms Vector fiels an the Lie erivative. A vector fiel on a manifol M is a smooth assignment of a vector tangent to M at each point. We think of M as
More informationMonotonicity for excited random walk in high dimensions
Monotonicity for excite ranom walk in high imensions Remco van er Hofsta Mark Holmes March, 2009 Abstract We prove that the rift θ, β) for excite ranom walk in imension is monotone in the excitement parameter
More informationConservation Laws. Chapter Conservation of Energy
20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action
More informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More information1 Heisenberg Representation
1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.
More informationANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS
ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics
More informationarxiv: v1 [cond-mat.stat-mech] 9 Jan 2012
arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,
More informationA REPRESENTATION FOR THE KANTOROVICH RUBINSTEIN DISTANCE DEFINED BY THE CAMERON MARTIN NORM OF A GAUSSIAN MEASURE ON A BANACH SPACE
Theory of Stochastic Processes Vol. 21 (37), no. 2, 2016, pp. 84 90 G. V. RIABOV A REPRESENTATION FOR THE KANTOROVICH RUBINSTEIN DISTANCE DEFINED BY THE CAMERON MARTIN NORM OF A GAUSSIAN MEASURE ON A BANACH
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationARBITRARY NUMBER OF LIMIT CYCLES FOR PLANAR DISCONTINUOUS PIECEWISE LINEAR DIFFERENTIAL SYSTEMS WITH TWO ZONES
Electronic Journal of Differential Equations, Vol. 2015 (2015), No. 228, pp. 1 12. ISSN: 1072-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu ftp eje.math.txstate.eu ARBITRARY NUMBER OF
More informationAn Optimal Algorithm for Bandit and Zero-Order Convex Optimization with Two-Point Feedback
Journal of Machine Learning Research 8 07) - Submitte /6; Publishe 5/7 An Optimal Algorithm for Banit an Zero-Orer Convex Optimization with wo-point Feeback Oha Shamir Department of Computer Science an
More informationII. First variation of functionals
II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent
More information12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes
Differentiation of vectors 12.5 Introuction The area known as vector calculus is use to moel mathematically a vast range of engineering phenomena incluing electrostatics, electromagnetic fiels, air flow
More informationApproximate Reduction of Dynamical Systems
Proceeings of the 4th IEEE Conference on Decision & Control Manchester Gran Hyatt Hotel San Diego, CA, USA, December 3-, 6 FrIP.7 Approximate Reuction of Dynamical Systems Paulo Tabuaa, Aaron D. Ames,
More informationChaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena
Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More information1. Aufgabenblatt zur Vorlesung Probability Theory
24.10.17 1. Aufgabenblatt zur Vorlesung By (Ω, A, P ) we always enote the unerlying probability space, unless state otherwise. 1. Let r > 0, an efine f(x) = 1 [0, [ (x) exp( r x), x R. a) Show that p f
More informationApproximations of displacement interpolations by entropic interpolations
Approximations of displacement interpolations by entropic interpolations Christian Léonard Université Paris Ouest Mokaplan 10 décembre 2015 Interpolations in P(X ) X : Riemannian manifold (state space)
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson
JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises
More informationOn Characterizing the Delay-Performance of Wireless Scheduling Algorithms
On Characterizing the Delay-Performance of Wireless Scheuling Algorithms Xiaojun Lin Center for Wireless Systems an Applications School of Electrical an Computer Engineering, Purue University West Lafayette,
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationCalculus of Variations
Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,
More informationQuantum mechanical approaches to the virial
Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from
More informationConvergence rates of moment-sum-of-squares hierarchies for optimal control problems
Convergence rates of moment-sum-of-squares hierarchies for optimal control problems Milan Kora 1, Diier Henrion 2,3,4, Colin N. Jones 1 Draft of September 8, 2016 Abstract We stuy the convergence rate
More informationBEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi
BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica
More informationThe Generalized Incompressible Navier-Stokes Equations in Besov Spaces
Dynamics of PDE, Vol1, No4, 381-400, 2004 The Generalize Incompressible Navier-Stokes Equations in Besov Spaces Jiahong Wu Communicate by Charles Li, receive July 21, 2004 Abstract This paper is concerne
More informationCHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold
CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the
More informationIntroduction. A Dirichlet Form approach to MCMC Optimal Scaling. MCMC idea
Introuction A Dirichlet Form approach to MCMC Optimal Scaling Markov chain Monte Carlo (MCMC quotes: Metropolis et al. (1953, running coe on the Los Alamos MANIAC: a feasible approach to statistical mechanics
More informationCurvature, Conformal Mapping, and 2D Stationary Fluid Flows. Michael Taylor
Curvature, Conformal Mapping, an 2D Stationary Flui Flows Michael Taylor 1. Introuction Let Ω be a 2D Riemannian manifol possibly with bounary). Assume Ω is oriente, with J enoting counterclockwise rotation
More informationDissipative numerical methods for the Hunter-Saxton equation
Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a
More information3 The variational formulation of elliptic PDEs
Chapter 3 The variational formulation of elliptic PDEs We now begin the theoretical stuy of elliptic partial ifferential equations an bounary value problems. We will focus on one approach, which is calle
More informationON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract
More informationFLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction
FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS ALINA BUCUR, CHANTAL DAVID, BROOKE FEIGON, MATILDE LALÍN 1 Introuction In this note, we stuy the fluctuations in the number
More information