Large deviations and averaging for systems of slow fast stochastic reaction diffusion equations.
|
|
- Loren Rogers
- 5 years ago
- Views:
Transcription
1 Large deviations and averaging for systems of slow fast stochastic reaction diffusion equations. Wenqing Hu. 1 (Joint work with Michael Salins 2, Konstantinos Spiliopoulos 3.) 1. Department of Mathematics and Statistics, Missouri S&T. 2. Department of Mathematics and Statistics, Boston University. 3. Department of Mathematics and Statistics, Boston University.
2 Systems of slow fast stochastic reaction diffusion equations : Equation. X ε,δ t Y ε,δ t (t, x) = A 1 X ε,δ (t, x) + b 1 (x, X ε,δ (t, x), Y ε,δ (t, x)) + εσ 1 (x, X ε,δ (t, x), Y ε,δ (t, x)) W Q 1 (t, x), t (t, x) = 1 ] [A δ 2 2 Y ε,δ (t, x) + b 2 (x, X ε,δ (t, x), Y ε,δ (t, x)) + 1 δ σ 2(x, X ε,δ (t, x), Y ε,δ (t, x)) W Q 2 (t, x), t X ε,δ (0, x) = X 0 (x), Y ε,δ (0, x) = Y 0 (x), x D, N 1 X ε,δ (t, x) = N 2 Y ε,δ (t, x) = 0, t 0, x D. (1)
3 Systems of slow fast stochastic reaction diffusion equations : Equation. Here ε > 0 is a small parameter and δ = δ(ε) > 0 is such that δ 0 as ε 0. The operators A 1 and A 2 are two strictly elliptic operators and b i (x, X, Y ), i = 1, 2 are the nonlinear terms. The noise processes W Q 1 and W Q 2 are two cylindrical Wiener processes with covariance matrices Q 1 and Q 2, and σ i (x, X, Y ), i = 1, 2 give the corresponding multiplicative noises. The initial values X 0 and Y 0 are assumed to be in L 2 (D). The boundary conditions are given by operators N i, i = 1, 2 which may correspond to either Dirichlet or Neumann conditions.
4 Systems of slow fast stochastic reaction diffusion equations : Averaging and Large Deviations. Fast slow system characterized by the small parameter δ > 0 ; The noise term εσ 1 (x, X ε,δ (t, x), Y ε,δ (t, x)) W Q 1 (t, x) in t the equation for the slow process X ε,δ (t, x) has a small parameter ε, and both the deterministic as well as the noise term in the equation for the fast process Y ε,δ (t, x) contain large parameters 1 δ 2 and 1 δ. In the limit, we expect an interplay between an averaging effect in the fast process Y ε,δ (t, x) and a large deviation effect in the slow process X ε,δ (t, x).
5 Systems of slow fast stochastic reaction diffusion equations : Previous work. Large deviations of stochastic partial differential equations of reaction diffusion type has been considered in previous works, but without the effect of multiple scales. Results in the case of slow fast systems of stochastic reaction diffusion equations have been considered in dimension one, with additive noise in the fast motion and no noise component in the slow motion. In finite dimensions the large deviations problem for multiscale diffusions has been well studied. To the best of our knowledge, the problem of large deviations for multiscale stochastic reaction diffusion equations in multiple dimensions, with multiplicative noise is being considered for the first time in the present work.
6 Systems of slow fast stochastic reaction diffusion equations : Application background. Fast slow stochastic partial differential equations of reaction diffusion type have many applications in chemistry and biology. In the classical chemical kinetics, the evolution of concentrations of various components in a reaction is described by ordinary differential equations. Such a description turns out to be unsatisfactory in a number of applications, especially in biology.
7 Systems of slow fast stochastic reaction diffusion equations : Application background. There are several ways to construct a more adequate mathematical model. If the reaction is fast enough, one should take into account that the concentration is not constant in space in the volume where the reaction takes place. Then the change of concentration due to the spatial transport, as a role the diffusion, should be taken into consideration and the system of ordinary differential equations should be replaced by a system of partial differential equations of reaction diffusion type.
8 Systems of slow fast stochastic reaction diffusion equations : Application background. In some cases, one should also take into account random change in time of the rates of reaction. Then the ordinary differential equation is replaced by a stochastic differential equation. If the rates change randomly not just in time but also in space, the evolution of concentrations can be described by a system of stochastic partial differential equations. On the other hand, the rates of chemical reactions in the system and the diffusion coefficients may have, and as a rule have, different orders. Some of them are much smaller than others and this leads to mathematical models based on slow fast stochastic reaction diffusion equations.
9 Systems of slow fast stochastic reaction diffusion equations : Problem Formulation. We want to have large deviation results for the slow motion X ε,δ in terms of Laplace principle. For every bounded and continuous function h : C([0, T ]; H) R we have [ ( exp 1 )] ε h(x ε,δ ) lim ε 0 ε ln E X0,Y 0 (Laplace principle.) = inf φ C([0,T ];H) [S X 0 (φ) + h(φ)]. Calculate the action functional (rate function) S X0 (φ).
10 Large deviations in terms of Laplace principle : Weak convergence method. One of the most effective methods in analyzing large deviation effects is the weak convergence method. Roughly speaking, by a variational representation of exponential functionals of Wiener processes, one can represent the exponential functional of the slow process X ε,δ (t, x) that appears in the Laplace principle (which is equivalent to large deviations principle) as a variational infimum over a family of controlled slow processes X ε,δ,u (t, x).
11 Large deviations in terms of Laplace principle : Weak convergence method. In particular, we have, for any bounded continuous function h : C([0, T ]; L 2 (D)) R, that [ ( ε ln E exp 1 )] ε h(x ε,δ ) [ 1 T ] = inf E u(s) 2 u L 2 U ([0,T ];U) 2 ds + h(x ε,δ,u ). 0
12 Large deviations in terms of Laplace principle : Weak convergence method. Here U is the control space, and the infimum is over all controls u U with finite L 2 ([0, T ]; ) norm. The controlled slow motion X ε,δ,u that appears on the right hand side of the Laplace principle comes from a controlled slow fast system (X ε,δ,u, Y ε,δ,u ) of reaction diffusion equations corresponding to (1).
13 Large deviations in terms of Laplace principle : Associated Control System. X ε,δ,u t Y ε,δ,u t (t, x) = A 1 X ε,δ,u (t, x) + b 1 (x, X ε,δ,u (t, x), Y ε,δ,u (t, x)) +σ 1 (x, X ε,δ,u (t, x), Y ε,δ,u (t, x))(q 1 u(t))(x) + εσ 1 (x, X ε,δ,u (t, x), Y ε,δ,u (t, x)) W Q 1 (t, x), t (t, x) = 1 ] [A δ 2 2 Y ε,δ,u (t, x) + b 2 (x, X ε,δ,u (t, x), Y ε,δ,u (t, x)) + 1 δ ε σ 2(x, X ε,δ,u (t, x), Y ε,δ,u (t, x))(q 2 u(t))(x) + 1 δ σ 2(x, X ε,δ,u (t, x), Y ε,δ,u (t, x)) W Q 2 (t, x), t X ε,δ,u (0, x) = X 0 (x), Y ε,δ,u (0, x) = Y 0 (x), x D, N 1 X ε,δ,u (t, x) = N 2 Y ε,δ,u (t, x) = 0, t 0, x D. (2)
14 Large deviations in terms of Laplace principle : Associated Control System. Need to analyze the limit as ε 0 (and thus δ 0) of the controlled fast slow system (X ε,δ,u, Y ε,δ,u ) in (2). Recall that for any bounded continuous function h : C([0, T ]; L 2 (D)) R, that [ ( ε ln E exp 1 )] ε h(x ε,δ ) = inf E u L 2 ([0,T ];U) u is in terms of a feedback control! [ 1 T ] u(s) 2 U 2 ds + h(x ε,δ,u ). 0
15 Set up and assmptions. H = L 2 (D) ; W (t) = k=1 β k(t) ; Q i : R H where i = 1, 2 add color to the noise and also decide if the noise are independent ; W Q i (t, x) = Q i W (t)(x). U = {x R, < } with norm x 2 U = xi 2 i=1 space in which the control lives in. xi 2 i=1 is the
16 Limit of the controlled fast slow process. We only know that the control u L 2 ([0, T ]; U) ; This makes the passage to the limit in the controlled fast slow system (X ε,δ,u, Y ε,δ,u ) as ε 0 (and therefore δ 0) difficult. Need more Hypotheses.
17 Hypothesis 1. Hypothesis 1. For i = 1, 2, there exist complete orthonormal systems {e i,k } k N in H, and sequences of non negative real numbers {α i,k } k N, such that A i e i,k = α i,k e i,k, k 1. The covariance operators Q i : U H, i = 1, 2 are diagonalized by the same orthonormal basis {e i,k } k N in the following sense. For i = 1, 2, there exists an orthonormal set {f i,k } k N. The set of {f i,k } k N is not necessarily complete. There exist sequences of non negative real numbers {λ i,k } i=1,2 satisfying k N Q i f i,k = λ i,k e i,k. Notice that if span{f 1,k } k N span{f 2,k } k N, then the driving noises of the fast and the slow motion are independent.
18 Hypothesis 1. Hypothesis 1 (cont d). If d = 1, then we have κ i := sup λ i,k e i,k 0 <, ζ i := k N k=1 α β i i,k e i,k 2 0 < for some constant β i (0, 1), and if d 2, we have κ i := k=1 λ ρ i i,k e i,k 2 0 <, ζ i := k=1 α β i i,k e i,k 2 0 < for some constants β i (0, + ) and ρ i (2, + ) such that β i (ρ i 2) ρ i < 1. Moreover inf α i,k =: λ > 0. k N i=1,2
19 Hypothesis 2. Hypothesis The mappings b i : D R 2 R and σ i : D R 2 R are measurable, both for i = 1 and for i = 2, and (L X b i + L X σ i + L Y b i + L Y σ i ) M for some M > 0. i=1,2 Moreover, sup b 2 (x, 0, 0) <, sup σ 2 (x, 0, 0) <. x D x D 2. Recalling λ, the constant introduced in (18), we have that L Y b 2 < λ.
20 Hypothesis 2. Hypothesis 2 (cont d). 3. σ 2 grows linearly in X, but is bounded in Y. There exists c > 0, sup sup σ 2 (x, X, Y ) c(1 + X ). x D Y R 4. The Lipschitz constants L Y b 2 and L Y σ 2 are so chosen that where L Y b 2 λ + K 2 (L Y σ 2 ) 2 K 2 = 0 ( β2 s β 2 ρ 2 2 ρ 2 e λ ρ 2 +2 ρ 2 s ds =: L Y b2,σ 2 < 1 e ) β2 ρ 2 2 ρ 2 ρ ρ ζ 2 2 ρ κ 2 2 and λ, β 2, ρ 2, ζ 2, κ 2 are all from Hypothesis 1.
21 Hypothesis 3. Hypothesis 3. b 1 and σ 1 grow at most linearly in X and sublinearly in Y. To be precise, there exists 0 ζ < 1 β 1(ρ 1 2) ρ 1 and a constant C > 0 such that sup ( b 1 (x, X, Y ) + σ 1 (x, X, Y ) ) C(1 + X + Y ζ ). x D
22 Limit of the controlled fast slow process. Controlled pair of fast slow process (X ε,δ,u, Y ε,δ,u ). The right way to describe the limit is a pair (ψ, P), which we call a viable pair. ψ C([0, T ]; H) is the limit dynamics, and P P(U Y [0, T ]) is a certain invariant measure. This is an averaging procedure. Constructions of the same type have been carried out thoroughly in the case of finite dimensions, but for the first time in infinite dimensions in our work.
23 Limit of the controlled fast slow process : Tightness. We introduce the family of random occupation measures P ε, (dudydt) = 1 t+ t ) 1 du (u(s))1 dy (Y ε,δ,u (s) dsdt. After very technical proof we can show that the pair (X ε,δ,u, P ε, ) is tight in the space C([0, T ]; H) P(U Y [0, T ]). We take weak topology on U and norm topology on Y, and we make use of a priori bounds and tightness functions.
24 Limit of the controlled fast slow process : Limiting slow dynamics. Tightness guarantees that (X ε,δ,u, P ε, ) ( X, P) as ε 0 (and therefore δ 0) for some limiting ( X, P). Recall that in mild form we can write t X ε,δ,u (t) = S 1 (t)x 0 + S 1 (t s)b 1 (X ε,δ,u (s), Y ε,δ,u (s))ds 0 t + S 1 (t s)σ 1 (X ε,δ,u (s), Y ε,δ,u (s))q 1 u(s)ds 0 + ε t Limiting slow dynamics ψ(t) = S 1 (t)x 0 + where 0 S 1 (t s)σ 1 (X ε,δ,u (s), Y ε,δ,u (s))dw Q 1 (s). U Y [0,t] S 1 (t s)ξ(ψ(s), Y, u)p(dudyds), ξ(x, Y, u) = Σ 1 (X, Y )Q 1 u + B 1 (X, Y ).
25 Limit of the controlled fast slow process : Limiting occupation measure. In regards to the limit for the controlled fast process Y ε,δ,u, the limiting measure P characterizes simultaneously the structure of the invariant measure of Y ε,δ,u and the control function u. In general, these two objects are intertwined and coupled together into the measure P, so that the averaging with respect to the measure P is different and hard to perform as in the classical averaging principle.
26 Limit of the controlled fast slow process : Limiting occupation measure. P(dudYdt) = η t (du Y )µ ψt (dy u)dt, and we cannot guarantee that µ is invariant measure to some process, as this process is a controlled process and we have no regularity properties of the optimal controls known other than being square integrable. In finite dimensional case with periodic coefficients, the problem can be resolved using the characterization of optimal controls through solutions to Hamilton Jacobi Bellman equations. Such a characterization is not rigorously known in infinite dimensions and even if that becomes the case, one would need to establish sufficient properties of such equations that would then imply that the resulting controlled process has a well defined invariant measure that is regular enough.
27 Limit of the controlled fast slow process : Limiting occupation measure. Is it possible to have P(dudYdt) = η t (du Y )µ ψt (dy )dt? Recall that the fast motion can be written formally as dy ε,δ,u (t) = 1 ] [A δ 2 2 Y ε,δ,u (t) + B 2 (X ε,δ,u (t), Y ε,δ,u (t)) dt+ + 1 δ 2 δ ε Σ 2 (X ε,δ,u (t), Y ε,δ,u (t))q 2 u(t)dt+ Y ε,δ,u (0) = Y 0 H. + 1 δ Σ 2(X ε,δ,u (t), Y ε,δ,u (t))dw Q 2, Send δ ε 0 as ε 0 (and therefore δ 0).
28 Hypothesis 5. Hypothesis 5. We assume that ε 0, δ = δ(ε) 0 and = (δ, ε) 0, such that lim ε 0 δ δ ε = 0, and lim ε 0 ε = 0.
29 Limit of the controlled fast slow process : Limiting occupation measure. From Hypothesis 5 we know that lim ε 0 δ δ ε = 0, and lim ε 0 ε = 0. Fast motion without control but driven by the controlled slow process dy ε,δ (t) = 1 δ 2 [ A2 Y ε,δ (t) + B 2 (X ε,δ,u (t), Y ε,δ (t)) ] dt+ + 1 δ Σ 2(X ε,δ,u (t), Y ε,δ (t))dw Q 2, Y ε,δ (0) = Y 0 H. One can show that E 1 T as ε 0 (and therefore δ 0). 0 Y ε,δ,u (t) Y ε,δ (t) 2 H dt 0,
30 Limit of the controlled fast slow process : Limiting occupation measure. Slogan : Ergodic properties are stable with respect to mildly regular small perturbations. Knowing closeness of two processes in mean square sense is already enough to draw conclusions about the equivalency of their invariant measures. P(dudYdt) = η t (du Y )µ ψt (dy )dt.
31 Averaging of the controlled fast slow process. Definition A pair (ψ, P) C([0, T ]; L 2 (D)) P(U Y [0, T ]) will be called viable with respect to (ξ, L), or simply viable if there is no confusion, if the following are satisfied. The function ψ(t) is absolutely continuous as a function in the space C([0, T ]; H), P is square integrable in the sense that ( u 2 U + Y 2 ) H P(dudYds) < U Y [0,T ] and the following hold for all t [0, T ] : ψ(t) = S 1 (t)x 0 + S 1 (t s)ξ(ψ(s), Y, u)p(dudyds), U Y [0,t]
32 Averaging of the controlled fast slow process. Definition the measure P is such that P P = { } P P(U Y [0, T ]) : P(dudYdt) = η(du Y, t)µ(dy t)dt, µ(dy t) = µ ψ(t) (dy ) for t [0, T ] where µ X is the invariant measure for the uncontrolled fast process with frozen slow component X, and P(U Y [0, t]) = t.
33 Averaging of the controlled fast slow process. Theorem Let u P N 2 (U), (X ε,δ,u, Y ε,δ,u ) be the mild solution to controlled fast process and T <. Let also P ε, (dudydt) be the occupation measure. Assume Hypotheses 1, 2, 3 and 5, X 0 H θ 1 with θ > 0 sufficiently small and Y 0 H. Then, the family of processes X ε,δ,u is tight in C([0, T ]; H) and the family of measures P ε, is tight in P(U Y [0, T ]), where U Y [0, T ] is endowed with the weak topology on U, the norm topology on Y and the standard topology on [0, T ]. Hence, given any subsequence of {(X ε,δ,u, P ε, ), ε, δ, > 0}, there exists a subsubsequence that converges in distribution with limit ( X, P). With probability 1, the accumulation point ( X, P) is a viable pair with respect to (ξ, L).
34 Large Deviations. Back to our large deviations problem : For every bounded and continuous function h : C([0, T ]; H) R we want [ ( exp 1 )] ε h(x ε,δ ) lim ε 0 ε ln E X0,Y 0 Representation formula [ ( ε ln E exp 1 )] ε h(x ε,δ ) = inf E u L 2 ([0,T ];U) = inf φ C([0,T ];H) [S X 0 (φ)+h(φ)]. [ 1 T ] u(s) 2 U 2 ds + h(x ε,δ,u ). 0
35 Large Deviations. Theorem Let (X ε,δ, Y ε,δ ) be the mild solution to fast slow SRDE and let T <. Assume Hypothesis 1, 2, 4, and 5 and let Y 0 H and X 0 H θ 1. Define S(φ) = S X0 (φ) = inf (φ,p) V (ξ,l) [ ] 1 u 2 U 2 P(dudYdt) U Y [0,T ] with the convention that the infimum over the empty set is.,
36 Large Deviations. Theorem Then, there exists θ > 0 such that for all θ (0, θ] and X 0 H1 θ, and for every bounded and continuous function h : C([0, T ]; H) R we have lim ε 0 ε ln E X0,Y 0 [ exp ( 1 ε h(x ε,δ ) )] = inf φ C([0,T ];H) [S X 0 (φ)+h(φ)]. In particular, {X ε,δ } satisfies the large deviations principle in C([0, T ]; H) with action functional S X0 ( ), uniformly for X 0 in compact subsets of H θ 1.
37 Large Deviations : Lower bound. By Fatou s Lemma we have [ ( lim inf ( ε ln E X0 exp h(x ε,δ )]) ) ε 0 ( [ ε 1 T ] ) lim inf E X0 u ε (t) 2 U ε 0 2 dt + h(x ε,δ ) ε ( [ 0 1 T 1 t+ ]) lim inf E X0 u ε (s) 2 U ε dsdt + h(x ε,δ ) [ t ]) 1 = lim inf (E X0 u 2 U ε 0 2 Pε, (dudydt) + h(x ε,δ ) U Y [0,T ] [ 1 u 2 U ], 2 P(dudYdt) + h(φ) inf (φ,p) V (ξ,l) U Y [0,T ] which concludes the proof of the Laplace principle lower bound.
38 Large Deviations : Upper bound. In order to prove the Laplace principle upper bound, we have to construct a nearly optimal control that achieves the upper bound. The nearly optimal control will be in feedback form with respect to the Y variable, Due to infinite dimensionality we can only work it out in some special cases.
39 Hypothesis 4. Hypothesis 4. b 1 is as in Hypothesis 3 and if d = 1, then there are positive constants 0 < c 0 c 1 < such that 0 < c 0 σ 2 1 (x, X, Y ) c 1, or if d > 1 then σ 1 (x, X, Y ) = σ 1 (x, X ) is independent of Y and can grow at most linearly in X uniformly in x D.
40 Thank you for your attention!
Hypoelliptic multiscale Langevin diffusions and Slow fast stochastic reaction diffusion equations.
Hypoelliptic multiscale Langevin diffusions and Slow fast stochastic reaction diffusion equations. Wenqing Hu. 1 (Joint works with Michael Salins 2 and Konstantinos Spiliopoulos 3.) 1. Department of Mathematics
More informationAsymptotic Analysis 88 (2014) DOI /ASY IOS Press
Asymptotic Analysis 88 4) 5 DOI.333/ASY-4 IOS Press Smoluchowski Kramers approximation and large deviations for infinite dimensional gradient systems Sandra Cerrai and Michael Salins Department of Mathematics,
More informationA random perturbation approach to some stochastic approximation algorithms in optimization.
A random perturbation approach to some stochastic approximation algorithms in optimization. Wenqing Hu. 1 (Presentation based on joint works with Chris Junchi Li 2, Weijie Su 3, Haoyi Xiong 4.) 1. Department
More informationFast-slow systems with chaotic noise
Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 12, 215 Averaging and homogenization workshop, Luminy. Fast-slow systems
More informationSingular Perturbations of Stochastic Control Problems with Unbounded Fast Variables
Singular Perturbations of Stochastic Control Problems with Unbounded Fast Variables Joao Meireles joint work with Martino Bardi and Guy Barles University of Padua, Italy Workshop "New Perspectives in Optimal
More informationSome Properties of NSFDEs
Chenggui Yuan (Swansea University) Some Properties of NSFDEs 1 / 41 Some Properties of NSFDEs Chenggui Yuan Swansea University Chenggui Yuan (Swansea University) Some Properties of NSFDEs 2 / 41 Outline
More informationOn the fast convergence of random perturbations of the gradient flow.
On the fast convergence of random perturbations of the gradient flow. Wenqing Hu. 1 (Joint work with Chris Junchi Li 2.) 1. Department of Mathematics and Statistics, Missouri S&T. 2. Department of Operations
More informationRegularity of the density for the stochastic heat equation
Regularity of the density for the stochastic heat equation Carl Mueller 1 Department of Mathematics University of Rochester Rochester, NY 15627 USA email: cmlr@math.rochester.edu David Nualart 2 Department
More informationv( x) u( y) dy for any r > 0, B r ( x) Ω, or equivalently u( w) ds for any r > 0, B r ( x) Ω, or ( not really) equivalently if v exists, v 0.
Sep. 26 The Perron Method In this lecture we show that one can show existence of solutions using maximum principle alone.. The Perron method. Recall in the last lecture we have shown the existence of solutions
More informationUniformly Uniformly-ergodic Markov chains and BSDEs
Uniformly Uniformly-ergodic Markov chains and BSDEs Samuel N. Cohen Mathematical Institute, University of Oxford (Based on joint work with Ying Hu, Robert Elliott, Lukas Szpruch) Centre Henri Lebesgue,
More informationOn some relations between Optimal Transport and Stochastic Geometric Mechanics
Title On some relations between Optimal Transport and Stochastic Geometric Mechanics Banff, December 218 Ana Bela Cruzeiro Dep. Mathematics IST and Grupo de Física-Matemática Univ. Lisboa 1 / 2 Title Based
More informationOn the Goodness-of-Fit Tests for Some Continuous Time Processes
On the Goodness-of-Fit Tests for Some Continuous Time Processes Sergueï Dachian and Yury A. Kutoyants Laboratoire de Mathématiques, Université Blaise Pascal Laboratoire de Statistique et Processus, Université
More informationStrong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term
1 Strong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term Enrico Priola Torino (Italy) Joint work with G. Da Prato, F. Flandoli and M. Röckner Stochastic Processes
More informationLipschitz continuity for solutions of Hamilton-Jacobi equation with Ornstein-Uhlenbeck operator
Lipschitz continuity for solutions of Hamilton-Jacobi equation with Ornstein-Uhlenbeck operator Thi Tuyen Nguyen Ph.D student of University of Rennes 1 Joint work with: Prof. E. Chasseigne(University of
More informationParameter Estimation for Stochastic Evolution Equations with Non-commuting Operators
Parameter Estimation for Stochastic Evolution Equations with Non-commuting Operators Sergey V. Lototsky and Boris L. Rosovskii In: V. Korolyuk, N. Portenko, and H. Syta (editors), Skorokhod s Ideas in
More informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationResolving the White Noise Paradox in the Regularisation of Inverse Problems
1 / 32 Resolving the White Noise Paradox in the Regularisation of Inverse Problems Hanne Kekkonen joint work with Matti Lassas and Samuli Siltanen Department of Mathematics and Statistics University of
More informationNeighboring feasible trajectories in infinite dimension
Neighboring feasible trajectories in infinite dimension Marco Mazzola Université Pierre et Marie Curie (Paris 6) H. Frankowska and E. M. Marchini Control of State Constrained Dynamical Systems Padova,
More informationSUPPLEMENT TO CONTROLLED EQUILIBRIUM SELECTION IN STOCHASTICALLY PERTURBED DYNAMICS
SUPPLEMENT TO CONTROLLED EQUILIBRIUM SELECTION IN STOCHASTICALLY PERTURBED DYNAMICS By Ari Arapostathis, Anup Biswas, and Vivek S. Borkar The University of Texas at Austin, Indian Institute of Science
More informationStochastic Gradient Descent in Continuous Time
Stochastic Gradient Descent in Continuous Time Justin Sirignano University of Illinois at Urbana Champaign with Konstantinos Spiliopoulos (Boston University) 1 / 27 We consider a diffusion X t X = R m
More informationHölder continuity of the solution to the 3-dimensional stochastic wave equation
Hölder continuity of the solution to the 3-dimensional stochastic wave equation (joint work with Yaozhong Hu and Jingyu Huang) Department of Mathematics Kansas University CBMS Conference: Analysis of Stochastic
More informationBackward Stochastic Differential Equations with Infinite Time Horizon
Backward Stochastic Differential Equations with Infinite Time Horizon Holger Metzler PhD advisor: Prof. G. Tessitore Università di Milano-Bicocca Spring School Stochastic Control in Finance Roscoff, March
More information1 Directional Derivatives and Differentiability
Wednesday, January 18, 2012 1 Directional Derivatives and Differentiability Let E R N, let f : E R and let x 0 E. Given a direction v R N, let L be the line through x 0 in the direction v, that is, L :=
More informationPairwise Comparison Dynamics for Games with Continuous Strategy Space
Pairwise Comparison Dynamics for Games with Continuous Strategy Space Man-Wah Cheung https://sites.google.com/site/jennymwcheung University of Wisconsin Madison Department of Economics Nov 5, 2013 Evolutionary
More informationKramers formula for chemical reactions in the context of Wasserstein gradient flows. Michael Herrmann. Mathematical Institute, University of Oxford
eport no. OxPDE-/8 Kramers formula for chemical reactions in the context of Wasserstein gradient flows by Michael Herrmann Mathematical Institute, University of Oxford & Barbara Niethammer Mathematical
More informationDynamical systems with Gaussian and Levy noise: analytical and stochastic approaches
Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Noise is often considered as some disturbing component of the system. In particular physical situations, noise becomes
More informationStarting from Heat Equation
Department of Applied Mathematics National Chiao Tung University Hsin-Chu 30010, TAIWAN 20th August 2009 Analytical Theory of Heat The differential equations of the propagation of heat express the most
More informationPropagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R
Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455 Abstract We consider semilinear
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationFinding Limits of Multiscale SPDEs
Finding Limits of Multiscale SPDEs David Kelly Martin Hairer Mathematics Institute University of Warwick Coventry UK CV4 7AL dtbkelly@gmail.com November 3, 213 Multiscale Systems, Warwick David Kelly (Warwick)
More informationStochastic nonlinear Schrödinger equations and modulation of solitary waves
Stochastic nonlinear Schrödinger equations and modulation of solitary waves A. de Bouard CMAP, Ecole Polytechnique, France joint work with R. Fukuizumi (Sendai, Japan) Deterministic and stochastic front
More informationExercises from other sources REAL NUMBERS 2,...,
Exercises from other sources REAL NUMBERS 1. Find the supremum and infimum of the following sets: a) {1, b) c) 12, 13, 14, }, { 1 3, 4 9, 13 27, 40 } 81,, { 2, 2 + 2, 2 + 2 + } 2,..., d) {n N : n 2 < 10},
More informationProceedings of the 2014 Winter Simulation Conference A. Tolk, S. D. Diallo, I. O. Ryzhov, L. Yilmaz, S. Buckley, and J. A. Miller, eds.
Proceedings of the 014 Winter Simulation Conference A. Tolk, S. D. Diallo, I. O. Ryzhov, L. Yilmaz, S. Buckley, and J. A. Miller, eds. Rare event simulation in the neighborhood of a rest point Paul Dupuis
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
More informationSingularities and Laplacians in Boundary Value Problems for Nonlinear Ordinary Differential Equations
Singularities and Laplacians in Boundary Value Problems for Nonlinear Ordinary Differential Equations Irena Rachůnková, Svatoslav Staněk, Department of Mathematics, Palacký University, 779 OLOMOUC, Tomkova
More informationExam February h
Master 2 Mathématiques et Applications PUF Ho Chi Minh Ville 2009/10 Viscosity solutions, HJ Equations and Control O.Ley (INSA de Rennes) Exam February 2010 3h Written-by-hands documents are allowed. Printed
More informationScaling Limits of Waves in Convex Scalar Conservation Laws under Random Initial Perturbations
Scaling Limits of Waves in Convex Scalar Conservation Laws under Random Initial Perturbations Jan Wehr and Jack Xin Abstract We study waves in convex scalar conservation laws under noisy initial perturbations.
More informationFast-slow systems with chaotic noise
Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 1, 216 Statistical properties of dynamical systems, ESI Vienna. David
More informationAn iterative procedure for constructing subsolutions of discrete-time optimal control problems
An iterative procedure for constructing subsolutions of discrete-time optimal control problems Markus Fischer version of November, 2011 Abstract An iterative procedure for constructing subsolutions of
More informationOn semilinear elliptic equations with measure data
On semilinear elliptic equations with measure data Andrzej Rozkosz (joint work with T. Klimsiak) Nicolaus Copernicus University (Toruń, Poland) Controlled Deterministic and Stochastic Systems Iasi, July
More informationAsymptotic stability of an evolutionary nonlinear Boltzmann-type equation
Acta Polytechnica Hungarica Vol. 14, No. 5, 217 Asymptotic stability of an evolutionary nonlinear Boltzmann-type equation Roksana Brodnicka, Henryk Gacki Institute of Mathematics, University of Silesia
More informationChapter 2 Random Ordinary Differential Equations
Chapter 2 Random Ordinary Differential Equations Let (Ω,F, P) be a probability space, where F is a σ -algebra on Ω and P is a probability measure, and let η :[0, T ] Ω R m be an R m -valued stochastic
More informationOn a class of stochastic differential equations in a financial network model
1 On a class of stochastic differential equations in a financial network model Tomoyuki Ichiba Department of Statistics & Applied Probability, Center for Financial Mathematics and Actuarial Research, University
More informationHarnack Inequalities and Applications for Stochastic Equations
p. 1/32 Harnack Inequalities and Applications for Stochastic Equations PhD Thesis Defense Shun-Xiang Ouyang Under the Supervision of Prof. Michael Röckner & Prof. Feng-Yu Wang March 6, 29 p. 2/32 Outline
More informationBranching Processes II: Convergence of critical branching to Feller s CSB
Chapter 4 Branching Processes II: Convergence of critical branching to Feller s CSB Figure 4.1: Feller 4.1 Birth and Death Processes 4.1.1 Linear birth and death processes Branching processes can be studied
More informationWeak convergence and Brownian Motion. (telegram style notes) P.J.C. Spreij
Weak convergence and Brownian Motion (telegram style notes) P.J.C. Spreij this version: December 8, 2006 1 The space C[0, ) In this section we summarize some facts concerning the space C[0, ) of real
More informationScaling Limits of Waves in Convex Scalar Conservation Laws Under Random Initial Perturbations
Journal of Statistical Physics, Vol. 122, No. 2, January 2006 ( C 2006 ) DOI: 10.1007/s10955-005-8006-x Scaling Limits of Waves in Convex Scalar Conservation Laws Under Random Initial Perturbations Jan
More informationUncertainty quantification and systemic risk
Uncertainty quantification and systemic risk Josselin Garnier (Université Paris Diderot) with George Papanicolaou and Tzu-Wei Yang (Stanford University) February 3, 2016 Modeling systemic risk We consider
More informationBridging the Gap between Center and Tail for Multiscale Processes
Bridging the Gap between Center and Tail for Multiscale Processes Matthew R. Morse Department of Mathematics and Statistics Boston University BU-Keio 2016, August 16 Matthew R. Morse (BU) Moderate Deviations
More informationContinuous dependence estimates for the ergodic problem with an application to homogenization
Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,
More informationOn LDP and inviscid hydrodynamical equations
On LDP and inviscid hydrodynamical equations Annie Millet (amillet@univ-paris1.fr ) collaboration with H. Bessaih (Bessaih@uwyo.edu) SAMM, Université Paris 1 and PMA University of Wyoming Workshop Stochastic
More informationConsistency of the maximum likelihood estimator for general hidden Markov models
Consistency of the maximum likelihood estimator for general hidden Markov models Jimmy Olsson Centre for Mathematical Sciences Lund University Nordstat 2012 Umeå, Sweden Collaborators Hidden Markov models
More informationTheory and applications of random Poincaré maps
Nils Berglund nils.berglund@univ-orleans.fr http://www.univ-orleans.fr/mapmo/membres/berglund/ Random Dynamical Systems Theory and applications of random Poincaré maps Nils Berglund MAPMO, Université d'orléans,
More informationBERNOULLI ACTIONS AND INFINITE ENTROPY
BERNOULLI ACTIONS AND INFINITE ENTROPY DAVID KERR AND HANFENG LI Abstract. We show that, for countable sofic groups, a Bernoulli action with infinite entropy base has infinite entropy with respect to every
More informationWELL-POSEDNESS FOR HYPERBOLIC PROBLEMS (0.2)
WELL-POSEDNESS FOR HYPERBOLIC PROBLEMS We will use the familiar Hilbert spaces H = L 2 (Ω) and V = H 1 (Ω). We consider the Cauchy problem (.1) c u = ( 2 t c )u = f L 2 ((, T ) Ω) on [, T ] Ω u() = u H
More information1. Stochastic Processes and filtrations
1. Stochastic Processes and 1. Stoch. pr., A stochastic process (X t ) t T is a collection of random variables on (Ω, F) with values in a measurable space (S, S), i.e., for all t, In our case X t : Ω S
More informationOn the stochastic nonlinear Schrödinger equation
On the stochastic nonlinear Schrödinger equation Annie Millet collaboration with Z. Brzezniak SAMM, Paris 1 and PMA Workshop Women in Applied Mathematics, Heraklion - May 3 211 Outline 1 The NL Shrödinger
More informationEffective dynamics for the (overdamped) Langevin equation
Effective dynamics for the (overdamped) Langevin equation Frédéric Legoll ENPC and INRIA joint work with T. Lelièvre (ENPC and INRIA) Enumath conference, MS Numerical methods for molecular dynamics EnuMath
More informationHypothesis testing for Stochastic PDEs. Igor Cialenco
Hypothesis testing for Stochastic PDEs Igor Cialenco Department of Applied Mathematics Illinois Institute of Technology igor@math.iit.edu Joint work with Liaosha Xu Research partially funded by NSF grants
More informationLecture Notes in Advanced Calculus 1 (80315) Raz Kupferman Institute of Mathematics The Hebrew University
Lecture Notes in Advanced Calculus 1 (80315) Raz Kupferman Institute of Mathematics The Hebrew University February 7, 2007 2 Contents 1 Metric Spaces 1 1.1 Basic definitions...........................
More informationExercises in stochastic analysis
Exercises in stochastic analysis Franco Flandoli, Mario Maurelli, Dario Trevisan The exercises with a P are those which have been done totally or partially) in the previous lectures; the exercises with
More informationSDE Coefficients. March 4, 2008
SDE Coefficients March 4, 2008 The following is a summary of GARD sections 3.3 and 6., mainly as an overview of the two main approaches to creating a SDE model. Stochastic Differential Equations (SDE)
More informationON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS
Bendikov, A. and Saloff-Coste, L. Osaka J. Math. 4 (5), 677 7 ON THE REGULARITY OF SAMPLE PATHS OF SUB-ELLIPTIC DIFFUSIONS ON MANIFOLDS ALEXANDER BENDIKOV and LAURENT SALOFF-COSTE (Received March 4, 4)
More informationEcon 204 Differential Equations. 1 Existence and Uniqueness of Solutions
Econ 4 Differential Equations In this supplement, we use the methods we have developed so far to study differential equations. 1 Existence and Uniqueness of Solutions Definition 1 A differential equation
More informationu( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3)
M ath 5 2 7 Fall 2 0 0 9 L ecture 4 ( S ep. 6, 2 0 0 9 ) Properties and Estimates of Laplace s and Poisson s Equations In our last lecture we derived the formulas for the solutions of Poisson s equation
More informationNonlinear error dynamics for cycled data assimilation methods
Nonlinear error dynamics for cycled data assimilation methods A J F Moodey 1, A S Lawless 1,2, P J van Leeuwen 2, R W E Potthast 1,3 1 Department of Mathematics and Statistics, University of Reading, UK.
More informationGeneral Theory of Large Deviations
Chapter 30 General Theory of Large Deviations A family of random variables follows the large deviations principle if the probability of the variables falling into bad sets, representing large deviations
More informationPhase-field systems with nonlinear coupling and dynamic boundary conditions
1 / 46 Phase-field systems with nonlinear coupling and dynamic boundary conditions Cecilia Cavaterra Dipartimento di Matematica F. Enriques Università degli Studi di Milano cecilia.cavaterra@unimi.it VIII
More informationStochastic Homogenization for Reaction-Diffusion Equations
Stochastic Homogenization for Reaction-Diffusion Equations Jessica Lin McGill University Joint Work with Andrej Zlatoš June 18, 2018 Motivation: Forest Fires ç ç ç ç ç ç ç ç ç ç Motivation: Forest Fires
More informationPatrizio Frosini. Department of Mathematics and ARCES, University of Bologna
An observer-oriented approach to topological data analysis Part 2: The algebra of group invariant non-expansive operators and its application in the project GIPHOD - Group Invariant Persistent Homology
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationStochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno
Stochastic Processes M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno 1 Outline Stochastic (random) processes. Autocorrelation. Crosscorrelation. Spectral density function.
More informationb i (µ, x, s) ei ϕ(x) µ s (dx) ds (2) i=1
NONLINEAR EVOLTION EQATIONS FOR MEASRES ON INFINITE DIMENSIONAL SPACES V.I. Bogachev 1, G. Da Prato 2, M. Röckner 3, S.V. Shaposhnikov 1 The goal of this work is to prove the existence of a solution to
More information3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due ). Show that the open disk x 2 + y 2 < 1 is a countable union of planar elementary sets. Show that the closed disk x 2 + y 2 1 is a countable
More informationSolution of Stochastic Optimal Control Problems and Financial Applications
Journal of Mathematical Extension Vol. 11, No. 4, (2017), 27-44 ISSN: 1735-8299 URL: http://www.ijmex.com Solution of Stochastic Optimal Control Problems and Financial Applications 2 Mat B. Kafash 1 Faculty
More informationLarge Deviations from the Hydrodynamic Limit for a System with Nearest Neighbor Interactions
Large Deviations from the Hydrodynamic Limit for a System with Nearest Neighbor Interactions Amarjit Budhiraja Department of Statistics & Operations Research University of North Carolina at Chapel Hill
More informationWeb Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D.
Web Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D. Ruppert A. EMPIRICAL ESTIMATE OF THE KERNEL MIXTURE Here we
More informationProbability and Measure
Part II Year 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2018 84 Paper 4, Section II 26J Let (X, A) be a measurable space. Let T : X X be a measurable map, and µ a probability
More informationAutomorphic Equivalence Within Gapped Phases
1 Harvard University May 18, 2011 Automorphic Equivalence Within Gapped Phases Robert Sims University of Arizona based on joint work with Sven Bachmann, Spyridon Michalakis, and Bruno Nachtergaele 2 Outline:
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationBrownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539
Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory
More informationControllability of linear PDEs (I): The wave equation
Controllability of linear PDEs (I): The wave equation M. González-Burgos IMUS, Universidad de Sevilla Doc Course, Course 2, Sevilla, 2018 Contents 1 Introduction. Statement of the problem 2 Distributed
More informationZ. Zhou On the classical limit of a time-dependent self-consistent field system: analysis. computation
On the classical limit of a time-dependent self-consistent field system: analysis and computation Zhennan Zhou 1 joint work with Prof. Shi Jin and Prof. Christof Sparber. 1 Department of Mathematics Duke
More informationGeometry of log-concave Ensembles of random matrices
Geometry of log-concave Ensembles of random matrices Nicole Tomczak-Jaegermann Joint work with Radosław Adamczak, Rafał Latała, Alexander Litvak, Alain Pajor Cortona, June 2011 Nicole Tomczak-Jaegermann
More informationLecture 7 Introduction to Statistical Decision Theory
Lecture 7 Introduction to Statistical Decision Theory I-Hsiang Wang Department of Electrical Engineering National Taiwan University ihwang@ntu.edu.tw December 20, 2016 1 / 55 I-Hsiang Wang IT Lecture 7
More informationConvergence of Particle Filtering Method for Nonlinear Estimation of Vortex Dynamics
Convergence of Particle Filtering Method for Nonlinear Estimation of Vortex Dynamics Meng Xu Department of Mathematics University of Wyoming February 20, 2010 Outline 1 Nonlinear Filtering Stochastic Vortex
More informationFunctional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals
Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Noèlia Viles Cuadros BCAM- Basque Center of Applied Mathematics with Prof. Enrico
More informationUncertainty quantification and systemic risk
Uncertainty quantification and systemic risk Josselin Garnier (Université Paris 7) with George Papanicolaou and Tzu-Wei Yang (Stanford University) April 9, 2013 Modeling Systemic Risk We consider evolving
More informationMyopic Models of Population Dynamics on Infinite Networks
Myopic Models of Population Dynamics on Infinite Networks Robert Carlson Department of Mathematics University of Colorado at Colorado Springs rcarlson@uccs.edu June 30, 2014 Outline Reaction-diffusion
More informationErgodicity in data assimilation methods
Ergodicity in data assimilation methods David Kelly Andy Majda Xin Tong Courant Institute New York University New York NY www.dtbkelly.com April 15, 2016 ETH Zurich David Kelly (CIMS) Data assimilation
More informationPontryagin Approximations for Optimal Design
Pontryagin Approximations for Optimal Design Jesper Carlsson CSC, KTH BIT Conference, 2006 Optimal shape problem Find domain D R d such that inf F(ϕ) dx D D ff, G(ϕ) = 0 in D ad D Parameter design problem
More informationLandesman-Lazer type results for second order Hamilton-Jacobi-Bellman equations
Author manuscript, published in "Journal of Functional Analysis 258, 12 (2010) 4154-4182" Landesman-Lazer type results for second order Hamilton-Jacobi-Bellman equations Patricio FELMER, Alexander QUAAS,
More informationEcon Lecture 14. Outline
Econ 204 2010 Lecture 14 Outline 1. Differential Equations and Solutions 2. Existence and Uniqueness of Solutions 3. Autonomous Differential Equations 4. Complex Exponentials 5. Linear Differential Equations
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationON THE POLICY IMPROVEMENT ALGORITHM IN CONTINUOUS TIME
ON THE POLICY IMPROVEMENT ALGORITHM IN CONTINUOUS TIME SAUL D. JACKA AND ALEKSANDAR MIJATOVIĆ Abstract. We develop a general approach to the Policy Improvement Algorithm (PIA) for stochastic control problems
More informationStability of Stochastic Differential Equations
Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010
More informationVariational and Topological methods : Theory, Applications, Numerical Simulations, and Open Problems 6-9 June 2012, Northern Arizona University
Variational and Topological methods : Theory, Applications, Numerical Simulations, and Open Problems 6-9 June 22, Northern Arizona University Some methods using monotonicity for solving quasilinear parabolic
More information