Uncertainty quantification and systemic risk
|
|
- Stella Elliott
- 5 years ago
- Views:
Transcription
1 Uncertainty quantification and systemic risk Josselin Garnier (Université Paris 7) with George Papanicolaou and Tzu-Wei Yang (Stanford University) April 9, 2013
2 Modeling Systemic Risk We consider evolving systems with a large number of inter-connected components, each of which can be in a normal state or in a failed state. We want to study the probability of overall failure of the system, that is, its systemic risk. There are three effects that we want to model and that contribute to the behavior of systemic risk: The intrinsic stability of each component The strength of external random perturbations to the system The degree of inter-connectedness or cooperation between components
3 Possible Applications Engineering systems with a large number of interacting parts. Components can fail but the system fails only when a large number of components fail simultaneously. Power distribution systems. Individual components of the system are calibrated to withstand fluctuations in demand by sharing loads, but sharing also increases the probability of an overall failure. Banking systems. Banks cooperate and by spreading the risk of credit shocks between them can operate with less restrictive individual risk policies. However, this increases the risk that they may all fail, that is, the systemic risk. We want to propose a simple model to explain that individual risk does not affect the systemic risk in an obvious way. In fact, it is possible to simultaneously reduce individual risk and increase the systemic risk.
4 A Basic, Bistable Mean-Field Model The (banking) system has N components (banks). The risk variable X j (t) of each component j = 1,...,N, satisfies the SDE dx j (t) = hv (X j (t))dt +σdw j (t)+θ ( X(t) Xj (t) ) dt. V(x) = x 4 /4 x 2 /2 is a potential with two stable states ±1. With θ = σ = 0, X j (t) converges to one of the states. We define 1 the normal state and +1 the failed state. h 0 is the intrinsic stability. X(t) := 1 N N i=1 X i(t) is the risk variable of the system. θ 0 is the attractive interaction. {W j (t),j = 1,...,N} are independent Brownian motions. σ 0 is the noise strength.
5 Why this Model? h, σ and θ control the three effects we want to study: intrinsic stability, random perturbations, and degree of cooperation. Why mean field interaction? Because it is the simplest interaction that models cooperative behavior. It can be generalized to include diversity as well as other more complex interactions such as hierarchical ones. Connection with UQ (Uncertainty Quantification): UQ quantifies the variance of the quantity of interest. The variance quantifies the normal fluctuations. It may not characterize a rare event. Connection with QMU (Quantification of Margins of Uncertainty): In an inter-connected system, individual components can be operating closer to their margin of failure, as they can benefit from the stability of the rest of the system. This, however, reduces the overall margin of uncertainty.
6 Schematic for the Individual Risk and the Systemic Risk Status
7 The probability-measure-valued process µ N (t) First idea: to analyze X(t) = 1 N N i=1 X i(t). Unfortunately, no closed equation for X(t), so we need to generalize this problem into a larger space. We define a probability-measure-valued process µ N (t), µ N (t,dx) := 1 N N j=1 δ X j (t)(dx). Then X(t) = yµ N (t,dy). (Dawson, 1983) µ N (t) converges in law as N to a deterministic process u(t). For t > 0, u(t) has a pdf u(t,x) solution of the nonlinear Fokker-Planck equation: t u = h [ V (x)u ] θ {[ ] 12 2 yu(t,y)dy x u }+ σ2 x x x 2u, recalling that dx j (t) = hv (X j (t))dt +θ( X(t) X j (t))dt +σdw j (t).
8 Existence of two stable equilibria u(t) converges to an equilibrium uξ e with a pdf of the form: { } uξ e (x) = 1 (x ξ)2 exp 2h Z ξ σ 2V(x), 2 σ2 2θ with the compatibility condition: ξ = m(ξ) := xue ξ (x)dx. Given θ and h, there exists a critical value σ c such that ξ = m(ξ) has two stable solutions ±ξ b if and only if σ < σ c. Simplification: for small h, we have σc 2 = 2θ/3+O(h) and ξ b = 1 3σ2 2θ +O(h). Let us say that u ξ e b is the normal state of the system, u+ξ e b is the failed state. We have µ N (t) N u(t). If u(0) = u ξ e b, then u(t) = u ξ e b for all time.
9 Simulation of X: impact of N 1.5 N=120, h=0.1, σ=1, θ=10 dm(0)/dξ=1.0086, 3σ 2 /2θ= Systemic Risk N=100, h=0.1, σ=1, θ=10 dm(0)/dξ=1.0086, 3σ 2 /2θ= t 1.5 N=80, h=0.1, σ=1, θ=10 dm(0)/dξ=1.0086, 3σ 2 /2θ= Systemic Risk Systemic Risk t t X n+1 j = Xj n hv (Xj n ) t +σ Wj n+1 +θ( 1 N N Xk n Xn j ) t. k=1
10 Simulation of X: impact of θ 1.5 N=100, h=0.1, σ=1, θ=10 dm(0)/dξ=1.0086, 3σ 2 /2θ= Systemic Risk N=100, h=0.1, σ=1, θ=6 dm(0)/dξ=1.0126, 3σ 2 /2θ= t 1.5 N=100, h=0.1, σ=1, θ=5 dm(0)/dξ=1.0141, 3σ 2 /2θ= Systemic Risk Systemic Risk t t X n+1 j = Xj n hv (Xj n ) t +σ Wj n+1 +θ( 1 N N Xk n Xn j ) t. k=1
11 Probability of system failures Assume that σ < σ c so the limit of µ N (t) has two stable equilibria u e ξ b (resp. u e ξ b ), the system s normal (resp. failed) state. For large N, let the empirical density µ N (0) u e ξ b. Then we expect that µ N (t) u e ξ b for all t > 0. However, as long as N is finite, the system collapse: µ N (0) u e ξ b, µ N (T) u e +ξ b happens in the time interval [0,T] with small but nonzero probability. We use large deviations to compute this small probability.
12 Large deviation principle (Dawson & Gärtner, 1987) Given an event A in the suitable space (smaller than C([0,T],M 1 (R))) Detail, then µ N (t,dx) = 1 N N i=1 δ X j (t)(dx) satisfies the large deviation principle with the rate function I h : Detail ( ) P(µ N A) N 1 exp N inf I h(φ). φ A The rate function I h has the following variational form: Detail I h (φ) = 1 T 2σ 2 0 sup f: φ,f 2 L h φ = 1 2 {[ σ2 2 x 2φ+θ x x φ,f (t) = φ t L h φ,f 2 x (t) 0 φ,fx 2 dt, ] } yφ(t,dy) φ +h [ V (x)φ ], x φ(t, dx)f(x).
13 How to interpret the rate function? I h (φ) = 1 T 2σ 2 0 sup f: φ,fx 2 0 φ t L h φ,f 2 / φ,fx 2 dt. Recall that the limit u(t) µ e ξ b of µ N (t) satisfies u t L h u = 0, and I h(u) = 0. So P(µ N A) N 1 exp( Ninf φ A I h (u)) = 1 if u A. Given a smooth probability-measure valued process φ u with φ(0) = µ e ξ b, let g satisfy φ t L h φ = (φg) x. (1) g is the external force/randomness making φ deviate from u. Note that by (1), sup f: φ,f 2 y 0 φ t L h φ,f 2 φ,f 2 x = sup f: φ,f 2 x 0 φ,f x g 2 φ,f 2 x = φ,g 2. Then I h (φ) = 1 0 φ,g2 dt, which is the L 2 -norm of g (the randomness) and indicates the price to pay to have µ N φ. This expression (1) exists only if (very unfortunately) φ is a density function. 2σ 2 T
14 Probability of system failures and small h analysis P(µ N A) N 1 exp( Ninf φ A I h (φ)). The rare event A of system failures is the set of all possible paths starting from u e ξ b (the normal state) to u e +ξ b (the failed state): Detail A = { φ : φ(0) = u e ξ b, φ(t) = u e +ξ b }. The major work is to compute inf φ A I h (φ), which is a nonlinear and infinite-dimensional problem. Here we assume that the intrinsic stability, h, is small so that we can make this problem tractable. Two good reasons to consider a small intrinsic stability h: 1. P(µ N A) is extremely small for large h, which is not the regime we have in mind. 2. The case h = 0 is analytically solvable and the cases h small are perturbations of it.
15 The case h = 0 In this case, inf φ A I 0 (φ) can be solved without approximation. Result: For h = 0, inf φ A I 0 (φ) has the unique minimizer p e (t,x)dx, a path of Gaussian measures: { } p e 1 (t,x) = exp (x ae (t)) 2, 2π σ2 2 σ2 2θ 2θ where a e (t) = 2ξ b t/t ξ b and ξ b = 1 3σ2 2θ. Therefore, Recall that A = inf I 0(φ) = I 0 (p e ) = 2ξ2 b φ A σ 2 T = 2 σ 2 T {φ : φ(0) = u e ξb, φ(t) = u e +ξb }. ) (1 3σ2. 2θ
16 The case h small When h is small, good candidates for the transition path of empirical densities are Gaussian with small perturbations: { [ 1 (x a(t)) 2 φ = p +hq : p(t,x) = 2πb2 (t) exp 2b 2 (t) ],a(0) = ξ b,a(t) = ξ b } For h small, the large deviation problem is solvable approximately, and the transition probability is P(µ N A) ( [ exp N ) 2 (1 3σ2 + 6h σ 2 T 2θ σ 2 Comments: ( σ 2 θ ) 2 ) ]) (1 σ2 +O(h 2 ). θ A large system is more stable than a small system. In the long run (T large), a transition will happen. Increase of the intrinsic stabilization parameter h reduces systemic risk. Mean transition times are simply related to transition probabilities in this approximation (Williams 82)
17 What about the individual risk? The risk variable X j (t) of component j satisfies dx j = h(x 3 j X j )dt +θ( X X j )dt +σdw j. Assume X j (0) = 1 (the normal state) and linearize X j around 1: X j (t) = 1+δX j (t), X(t) = 1+δ X(t) and δ X(t) = 1 N N j=1 δx j(t). δx j (t) and δ X(t) satisfy linear SDEs: dδx j = (θ+2h)δx j dt+θδ Xdt+σdW j, dδ X = 2hδ Xdt+ σ N δx j (t) is a Gaussian process with the stationary distribution σ N(0, 2 2(2h+θ) ) as N. Qualitatively speaking, the individual risk is external risk(σ 2 ) intrinsic stability(h)+risk diversification(θ), and note that σ 2 2(2h+θ) σ2 2θ for small h. N dw j. j=1
18 Why and when is risk diversification undesirable? For h small: Systemic Risk exp ( 2N σ 2 T )) (1 3σ2, Individual Risk σ2 2θ 2θ. If a banker wants to take more external risk (to earn more profits) or the economy is more uncertain (increase σ 2 ), in order to keep the individual risk low, he can diversify the risk to the other banks (increase θ). Let us assume that everyone simultaneously increase σ 2 and θ but carefully keep the ratio σ2 2θ at a low level. - The individual risk is kept low. - The systemic risk increases, although this cannot be observed by the banks (until the catastrophic transition happens).
19 Extensions Diversity : the values of the parameters θ,σ,h are component-dependent. Hierarchical model : the system is stabilized by a central component.
20 Modeling of Diversity in Cooperative Behavior t u1 t u K The cooperative behavior of components can be different across groups: dx j (t) = hv (X j (t))dt +σdw j (t)+θ j ( X (t) X j (t) ) dt. The components are partitioned into K groups. In group k, the components have cooperative parameter Θ k. In the limit N the empirical densities of each group converge to the solution of the joint Fokker-Planck equations: = h [ V ] 1 2 (x)u 1 + x 2 σ2 x 2u1 Θ1 x. = h [ V ] 1 2 (x)u K + x 2 σ2 x 2u K Θ K x where ρ k N is the size of group k. {[ y {[ y } K ρ k u k (t,y)dy x ]u 1 k=1 } K ρ k u k (t,y)dy x ]u K k=1
21 Impact of Component Diversity on the Systemic Risk Why is the diversity interesting? The model is more realistic and more widely applicable. Diversity significantly affects the system stability by reducing it. Impact from the diversity: Analytical and numerical studies show that even with the same parameters and with {θ j } whose average equals θ the system still changes significantly.
22 Simulation 2 - Impact of Diversity, Change of Θ k and ρ k N=100, h=0.1, σ=1, Θ=[10;10;10] ρ=[0.33;0.34;0.33], stdev(θ)/mean(θ)=0 dm(0)/dξ=1.0086, Σ (ρ /Θ )(3σ 2 /2Θ 1)= i i i i N=100, h=0.1, σ=1, Θ=[6;10;14] ρ=[0.33;0.34;0.33], stdev(θ)/mean(θ)= dm(0)/dξ=1.0092, Σ (ρ /Θ )(3σ 2 /2Θ 1)= i i i i N=100, h=0.1, σ=1, Θ=[2;10;18] ρ=[0.33;0.34;0.33], stdev(θ)/mean(θ)= dm(0)/dξ=1.0089, Σ (ρ /Θ )(3σ 2 /2Θ 1)= i i i i Systemic Risk Systemic Risk Systemic Risk t t t N=100, h=0.1, σ=1, Θ=[5;10;15] ρ=[0.1;0.8;0.1], stdev(θ)/mean(θ)= dm(0)/dξ=1.0089, Σ i (ρ i /Θ i )(3σ 2 /2Θ i 1)= N=100, h=0.1, σ=1, Θ=[5;10;15] ρ=[0.3;0.4;0.3], stdev(θ)/mean(θ)= dm(0)/dξ=1.0095, Σ i (ρ i /Θ i )(3σ 2 /2Θ i 1)= N=100, h=0.1, σ=1, Θ=[5;10;15] ρ=[0.45;0.1;0.45], stdev(θ)/mean(θ)= dm(0)/dξ=1.0099, Σ i (ρ i /Θ i )(3σ 2 /2Θ i 1)= Systemic Risk Systemic Risk Systemic Risk t t t
23 Analysis in the Diversity Case System with diversity have larger transition probabilities: When h is zero, the average of θ j is θ, and the diversity of θ j is small, then the system has a higher transition probability. Mathematically, if h = 0, Θ k = θ(1+δα k ) with δ 1, and K k=1 ρ kα k = 0, then ( ) P µ div N A { [ exp N ) ( K ) 2 (1 3σ2 (3σ 2δ 2 ρ σ 2 k α 2 2 k T 2θ 2θ + 1 T k=1 T 0 ( 1 e θs) 2 ds ) ]}
24 A Hierarchical Model of Systemic Risk Here we consider a hierarchical model with a central component: dx 0 = h 0 V 0 (X 0)dt θ 0 ( X0 1 N dx j = σdw j hv (X j )dt θ ( X j X 0 ) dt, N ) X j dt j=1 j = 1,...,N X 0 models the central stable component. It is intrinsically stable (h 0 > 0), and not subjected to external fluctuations. It interacts with the other components through a mean field interaction. X j, j = 1,...,N model individual components that are subjected to external fluctuations. They are (h > 0) or are not (h = 0) intrinsically stable. They interact with the central component X 0.
25 A Hierarchical Model of Systemic Risk - Analysis Nonlinear Fokker-Planck equations: In the limit N the pair ( X 0 (t), 1 N N j=1 δ X j (t)(dx) ) converges to (x 0 (t),u(t,x)dx) solution of the nonlinear Fokker-Planck equation t u = σ2 2 2 xx u + [ x hv (x)+θ(x x 0 (t))u ], with dx 0 dt = h 0V 0(x ( 0 ) θ 0 x0 xu(t,x)dx ). Existence of two equilibrium states (x 0 (t),u(t,x)) (x e,u e (x)) when σ is below a critical level: u e (x) = 1 exp ( 2hV(x)+θ(x x e) 2 ) Z e σ 2 with the compatibility equation xu e (x)dx = x e + h 0 θ 0 V 0 (x e). Large deviations principle to compute the probability of transition.
26 A Hierarchical Model of Systemic Risk - Results Exact results for h = 0 and expansions for small h (for the optimal paths and for the probability of transition). For h = 0: Resolution of a fourth-order ODE for the optimal x 0 (t) with boundary conditions. For the optimal path the mean of the individual components x(t) is ahead of x 0 (t): the individual components drive the transition. Stability increases with θ and decreases with θ 0.
27 On-going and future work Study models with diversities everywhere: The most general case (when the limit exists) is that V, σ and θ in each group can be different, i.e. if X j is in the group k, then the SDE is dx j (t) = hv k (X j (t))dt+σ k dw j (t)+θ k ( X (t) X j (t) ) dt. Using the analysis as a guide, we design importance sampling algorithms for computing efficiently (very) small systemic failure probabilities.
28 Topological spaces for the mean field model M 1 (R) is the space of probability measures on R with the Prohorov metric ρ, associated with the weak convergence. C([0,T],M 1 (R)) is the space of continuous functions from [0,T] to M 1 (R) with the metric sup 0 t T ρ(µ 1 (t),µ 2 (t)). M (R) = {µ M 1 (R), ϕ(y)µ(dy) < }, where ϕ(y) = y 4, serves as a Lyapunov function. M (R) is endowed with the inductive topology: µ n µ in M (R) if and only if µ n µ in M 1 (R) and sup n ϕ(y)µn (dy) <. C([0,T],M (R)) is the space of continuous functions from [0,T] to M (R) endowed with the topology: φ n ( ) φ( ) in C([0,T],M (R)) if and only if φ n ( ) φ( ) in C([0,T],M 1 (R)) and sup 0 t T sup n ϕ(y)φn (t,dy) <. Go Back
29 Large Deviations Principle Exact statement: inf I h (φ) liminf φ A N limsup N 1 N logp(µ N A) 1 N logp(µ N A) inf I h (φ) φ A Go Back
30 The rare event A = { φ : φ(0) = u e ξ b, φ(t) = u e +ξ b }. Definition of an enlarged rare event (with non-empty interior): A δ = { φ : φ(0) = u e ξ b, ρ(φ(t),u e +ξ b ) δ }. We have lim inf I h (φ) = inf I h(φ) δ 0φ A δ φ A Go Back
31 The classical Freidlin-Wentzell formula Let dx N = b(x N )dt + 1 N σdw t (X N (t)) t [0,T] satisfies a large deviation principle in C([0,T],R) with the rate function I((x(t)) t [0,T] ) = { 1 T 2σ 2 0 +, ( t x(t) b(x(t)) ) 2 dt, (x(t))t [0,T] H 1 otherwise. Go Back
Systemic risk and uncertainty quantification
Systemic risk and uncertainty quantification Josselin Garnier (Université Paris 7) with George Papanicolaou and Tzu-Wei Yang (Stanford University) October 6, 22 Modeling Systemic Risk We consider evolving
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 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 informationImplicit sampling for particle filters. Alexandre Chorin, Mathias Morzfeld, Xuemin Tu, Ethan Atkins
0/20 Implicit sampling for particle filters Alexandre Chorin, Mathias Morzfeld, Xuemin Tu, Ethan Atkins University of California at Berkeley 2/20 Example: Try to find people in a boat in the middle of
More informationLarge deviations and averaging for systems of slow fast stochastic reaction diffusion equations.
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
More informationStochastic differential equation models in biology Susanne Ditlevsen
Stochastic differential equation models in biology Susanne Ditlevsen Introduction This chapter is concerned with continuous time processes, which are often modeled as a system of ordinary differential
More informationEndogenous Information Choice
Endogenous Information Choice Lecture 7 February 11, 2015 An optimizing trader will process those prices of most importance to his decision problem most frequently and carefully, those of less importance
More informationHypoelliptic 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 informationStochastic Volatility and Correction to the Heat Equation
Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century
More informationUNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE
Surveys in Mathematics and its Applications ISSN 1842-6298 (electronic), 1843-7265 (print) Volume 5 (2010), 275 284 UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE Iuliana Carmen Bărbăcioru Abstract.
More informationLAN property for sde s with additive fractional noise and continuous time observation
LAN property for sde s with additive fractional noise and continuous time observation Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Samy Tindel (Purdue University) Vlad s 6th birthday,
More informationImaging with Ambient Noise II
Imaging with Ambient Noise II George Papanicolaou Stanford University Michigan State University, Department of Mathematics Richard E. Phillips Lecture Series April 21, 2009 With J. Garnier, University
More informationAn Uncertain Control Model with Application to. Production-Inventory System
An Uncertain Control Model with Application to Production-Inventory System Kai Yao 1, Zhongfeng Qin 2 1 Department of Mathematical Sciences, Tsinghua University, Beijing 100084, China 2 School of Economics
More informationNumerical Solutions of ODEs by Gaussian (Kalman) Filtering
Numerical Solutions of ODEs by Gaussian (Kalman) Filtering Hans Kersting joint work with Michael Schober, Philipp Hennig, Tim Sullivan and Han C. Lie SIAM CSE, Atlanta March 1, 2017 Emmy Noether Group
More informationStochastic Modelling in Climate Science
Stochastic Modelling in Climate Science David Kelly Mathematics Department UNC Chapel Hill dtbkelly@gmail.com November 16, 2013 David Kelly (UNC) Stochastic Climate November 16, 2013 1 / 36 Why use stochastic
More information5 Applying the Fokker-Planck equation
5 Applying the Fokker-Planck equation We begin with one-dimensional examples, keeping g = constant. Recall: the FPE for the Langevin equation with η(t 1 )η(t ) = κδ(t 1 t ) is = f(x) + g(x)η(t) t = x [f(x)p
More informationApplications of controlled paths
Applications of controlled paths Massimiliano Gubinelli CEREMADE Université Paris Dauphine OxPDE conference. Oxford. September 1th 212 ( 1 / 16 ) Outline I will exhibith various applications of the idea
More informationPoisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation
Poisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation Jingyi Zhu Department of Mathematics University of Utah zhu@math.utah.edu Collaborator: Marco Avellaneda (Courant
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationWeak Convergence of Numerical Methods for Dynamical Systems and Optimal Control, and a relation with Large Deviations for Stochastic Equations
Weak Convergence of Numerical Methods for Dynamical Systems and, and a relation with Large Deviations for Stochastic Equations Mattias Sandberg KTH CSC 2010-10-21 Outline The error representation for weak
More informationLecture 4: Ito s Stochastic Calculus and SDE. Seung Yeal Ha Dept of Mathematical Sciences Seoul National University
Lecture 4: Ito s Stochastic Calculus and SDE Seung Yeal Ha Dept of Mathematical Sciences Seoul National University 1 Preliminaries What is Calculus? Integral, Differentiation. Differentiation 2 Integral
More informationStochastic contraction BACS Workshop Chamonix, January 14, 2008
Stochastic contraction BACS Workshop Chamonix, January 14, 2008 Q.-C. Pham N. Tabareau J.-J. Slotine Q.-C. Pham, N. Tabareau, J.-J. Slotine () Stochastic contraction 1 / 19 Why stochastic contraction?
More informationWeakly interacting particle systems on graphs: from dense to sparse
Weakly interacting particle systems on graphs: from dense to sparse Ruoyu Wu University of Michigan (Based on joint works with Daniel Lacker and Kavita Ramanan) USC Math Finance Colloquium October 29,
More informationA MODEL FOR THE LONG-TERM OPTIMAL CAPACITY LEVEL OF AN INVESTMENT PROJECT
A MODEL FOR HE LONG-ERM OPIMAL CAPACIY LEVEL OF AN INVESMEN PROJEC ARNE LØKKA AND MIHAIL ZERVOS Abstract. We consider an investment project that produces a single commodity. he project s operation yields
More informationStochastic differential equations in neuroscience
Stochastic differential equations in neuroscience Nils Berglund MAPMO, Orléans (CNRS, UMR 6628) http://www.univ-orleans.fr/mapmo/membres/berglund/ Barbara Gentz, Universität Bielefeld Damien Landon, MAPMO-Orléans
More information2012 NCTS Workshop on Dynamical Systems
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz 2012 NCTS Workshop on Dynamical Systems National Center for Theoretical Sciences, National Tsing-Hua University Hsinchu,
More informationDensities for the Navier Stokes equations with noise
Densities for the Navier Stokes equations with noise Marco Romito Università di Pisa Universitat de Barcelona March 25, 2015 Summary 1 Introduction & motivations 2 Malliavin calculus 3 Besov bounds 4 Other
More informationGillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde
Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk The Three Lectures 1 Gillespie s algorithm and its relation
More informationIntroduction to asymptotic techniques for stochastic systems with multiple time-scales
Introduction to asymptotic techniques for stochastic systems with multiple time-scales Eric Vanden-Eijnden Courant Institute Motivating examples Consider the ODE {Ẋ = Y 3 + sin(πt) + cos( 2πt) X() = x
More informationGhost Imaging. Josselin Garnier (Université Paris Diderot)
Grenoble December, 014 Ghost Imaging Josselin Garnier Université Paris Diderot http:/www.josselin-garnier.org General topic: correlation-based imaging with noise sources. Particular application: Ghost
More informationChapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems
Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.
More informationThe concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More informationPath integrals for classical Markov processes
Path integrals for classical Markov processes Hugo Touchette National Institute for Theoretical Physics (NITheP) Stellenbosch, South Africa Chris Engelbrecht Summer School on Non-Linear Phenomena in Field
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 informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More information16. Working with the Langevin and Fokker-Planck equations
16. Working with the Langevin and Fokker-Planck equations In the preceding Lecture, we have shown that given a Langevin equation (LE), it is possible to write down an equivalent Fokker-Planck equation
More informationMean Field Games on networks
Mean Field Games on networks Claudio Marchi Università di Padova joint works with: S. Cacace (Rome) and F. Camilli (Rome) C. Marchi (Univ. of Padova) Mean Field Games on networks Roma, June 14 th, 2017
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
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 informationFor a stochastic process {Y t : t = 0, ±1, ±2, ±3, }, the mean function is defined by (2.2.1) ± 2..., γ t,
CHAPTER 2 FUNDAMENTAL CONCEPTS This chapter describes the fundamental concepts in the theory of time series models. In particular, we introduce the concepts of stochastic processes, mean and covariance
More informationMemory and hypoellipticity in neuronal models
Memory and hypoellipticity in neuronal models S. Ditlevsen R. Höpfner E. Löcherbach M. Thieullen Banff, 2017 What this talk is about : What is the effect of memory in probabilistic models for neurons?
More informationRegularization by noise in infinite dimensions
Regularization by noise in infinite dimensions Franco Flandoli, University of Pisa King s College 2017 Franco Flandoli, University of Pisa () Regularization by noise King s College 2017 1 / 33 Plan of
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationProduction and Relative Consumption
Mean Field Growth Modeling with Cobb-Douglas Production and Relative Consumption School of Mathematics and Statistics Carleton University Ottawa, Canada Mean Field Games and Related Topics III Henri Poincaré
More informationStrong Solutions and a Bismut-Elworthy-Li Formula for Mean-F. Formula for Mean-Field SDE s with Irregular Drift
Strong Solutions and a Bismut-Elworthy-Li Formula for Mean-Field SDE s with Irregular Drift Thilo Meyer-Brandis University of Munich joint with M. Bauer, University of Munich Conference for the 1th Anniversary
More informationRunge-Kutta Method for Solving Uncertain Differential Equations
Yang and Shen Journal of Uncertainty Analysis and Applications 215) 3:17 DOI 1.1186/s4467-15-38-4 RESEARCH Runge-Kutta Method for Solving Uncertain Differential Equations Xiangfeng Yang * and Yuanyuan
More informationSPDEs, criticality, and renormalisation
SPDEs, criticality, and renormalisation Hendrik Weber Mathematics Institute University of Warwick Potsdam, 06.11.2013 An interesting model from Physics I Ising model Spin configurations: Energy: Inverse
More informationFE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology
FE610 Stochastic Calculus for Financial Engineers Lecture 3. Calculaus in Deterministic and Stochastic Environments Steve Yang Stevens Institute of Technology 01/31/2012 Outline 1 Modeling Random Behavior
More informationIn terms of measures: Exercise 1. Existence of a Gaussian process: Theorem 2. Remark 3.
1. GAUSSIAN PROCESSES A Gaussian process on a set T is a collection of random variables X =(X t ) t T on a common probability space such that for any n 1 and any t 1,...,t n T, the vector (X(t 1 ),...,X(t
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
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 informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationRegularity of competitive equilibria in a production economy with externalities
Regularity of competitive equilibria in a production economy with externalities Elena del Mercato Vincenzo Platino Paris School of Economics - Université Paris 1 Panthéon Sorbonne QED-Jamboree Copenhagen,
More informationHJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011
Department of Probability and Mathematical Statistics Faculty of Mathematics and Physics, Charles University in Prague petrasek@karlin.mff.cuni.cz Seminar in Stochastic Modelling in Economics and Finance
More informationSub-Riemannian geometry in groups of diffeomorphisms and shape spaces
Sub-Riemannian geometry in groups of diffeomorphisms and shape spaces Sylvain Arguillère, Emmanuel Trélat (Paris 6), Alain Trouvé (ENS Cachan), Laurent May 2013 Plan 1 Sub-Riemannian geometry 2 Right-invariant
More informationImproved diffusion Monte Carlo
Improved diffusion Monte Carlo Jonathan Weare University of Chicago with Martin Hairer (U of Warwick) October 4, 2012 Diffusion Monte Carlo The original motivation for DMC was to compute averages 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 informationCourse on Inverse Problems
Stanford University School of Earth Sciences Course on Inverse Problems Albert Tarantola Third Lesson: Probability (Elementary Notions) Let u and v be two Cartesian parameters (then, volumetric probabilities
More informationMultiple Random Variables
Multiple Random Variables Joint Probability Density Let X and Y be two random variables. Their joint distribution function is F ( XY x, y) P X x Y y. F XY ( ) 1, < x
More informationChapter 2 - Survival Models
2-1 Chapter 2 - Survival Models Section 2.2 - Future Lifetime Random Variable and the Survival Function Let T x = ( Future lifelength beyond age x of an individual who has survived to age x [measured in
More informationTraveling waves dans le modèle de Kuramoto quenched
Traveling waves dans le modèle de Kuramoto quenched Eric Luçon MAP5 - Université Paris Descartes Grenoble - Journées MAS 2016 31 août 2016 Travail en commun avec Christophe Poquet (Lyon 1) [Giacomin, L.,
More informationMetastability for interacting particles in double-well potentials and Allen Cahn SPDEs
Nils Berglund nils.berglund@univ-orleans.fr http://www.univ-orleans.fr/mapmo/membres/berglund/ SPA 2017 Contributed Session: Interacting particle systems Metastability for interacting particles in double-well
More informationPath Decomposition of Markov Processes. Götz Kersting. University of Frankfurt/Main
Path Decomposition of Markov Processes Götz Kersting University of Frankfurt/Main joint work with Kaya Memisoglu, Jim Pitman 1 A Brownian path with positive drift 50 40 30 20 10 0 0 200 400 600 800 1000-10
More informationIntroduction to standard and non-standard Numerical Methods
Introduction to standard and non-standard Numerical Methods Dr. Mountaga LAM AMS : African Mathematic School 2018 May 23, 2018 One-step methods Runge-Kutta Methods Nonstandard Finite Difference Scheme
More informationLecture 4: Introduction to stochastic processes and stochastic calculus
Lecture 4: Introduction to stochastic processes and stochastic calculus Cédric Archambeau Centre for Computational Statistics and Machine Learning Department of Computer Science University College London
More informationLarge Deviations Principles for McKean-Vlasov SDEs, Skeletons, Supports and the law of iterated logarithm
Large Deviations Principles for McKean-Vlasov SDEs, Skeletons, Supports and the law of iterated logarithm Gonçalo dos Reis University of Edinburgh (UK) & CMA/FCT/UNL (PT) jointly with: W. Salkeld, U. of
More informationIn essence, Dynamical Systems is a science which studies differential equations. A differential equation here is the equation
Lecture I In essence, Dynamical Systems is a science which studies differential equations. A differential equation here is the equation ẋ(t) = f(x(t), t) where f is a given function, and x(t) is an unknown
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 informationA Class of Fractional Stochastic Differential Equations
Vietnam Journal of Mathematics 36:38) 71 79 Vietnam Journal of MATHEMATICS VAST 8 A Class of Fractional Stochastic Differential Equations Nguyen Tien Dung Department of Mathematics, Vietnam National University,
More informationThe Chaotic Character of the Stochastic Heat Equation
The Chaotic Character of the Stochastic Heat Equation March 11, 2011 Intermittency The Stochastic Heat Equation Blowup of the solution Intermittency-Example ξ j, j = 1, 2,, 10 i.i.d. random variables Taking
More informationA Lévy-Fokker-Planck equation: entropies and convergence to equilibrium
1/ 22 A Lévy-Fokker-Planck equation: entropies and convergence to equilibrium I. Gentil CEREMADE, Université Paris-Dauphine International Conference on stochastic Analysis and Applications Hammamet, Tunisia,
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 informationProperties of the Autocorrelation Function
Properties of the Autocorrelation Function I The autocorrelation function of a (real-valued) random process satisfies the following properties: 1. R X (t, t) 0 2. R X (t, u) =R X (u, t) (symmetry) 3. R
More informationStep Bunching in Epitaxial Growth with Elasticity Effects
Step Bunching in Epitaxial Growth with Elasticity Effects Tao Luo Department of Mathematics The Hong Kong University of Science and Technology joint work with Yang Xiang, Aaron Yip 05 Jan 2017 Tao Luo
More informationProblem Set 5: Solutions Math 201A: Fall 2016
Problem Set 5: s Math 21A: Fall 216 Problem 1. Define f : [1, ) [1, ) by f(x) = x + 1/x. Show that f(x) f(y) < x y for all x, y [1, ) with x y, but f has no fixed point. Why doesn t this example contradict
More informationA Parcimonious Long Term Mean Field Model for Mining Industries
A Parcimonious Long Term Mean Field Model for Mining Industries Yves Achdou Laboratoire J-L. Lions, Université Paris Diderot (joint work with P-N. Giraud, J-M. Lasry and P-L. Lions ) Mining Industries
More informationBMIR Lecture Series on Probability and Statistics Fall 2015 Discrete RVs
Lecture #7 BMIR Lecture Series on Probability and Statistics Fall 2015 Department of Biomedical Engineering and Environmental Sciences National Tsing Hua University 7.1 Function of Single Variable Theorem
More informationMean-field SDE driven by a fractional BM. A related stochastic control problem
Mean-field SDE driven by a fractional BM. A related stochastic control problem Rainer Buckdahn, Université de Bretagne Occidentale, Brest Durham Symposium on Stochastic Analysis, July 1th to July 2th,
More informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
More informationRice method for the maximum of Gaussian fields
Rice method for the maximum of Gaussian fields IWAP, July 8, 2010 Jean-Marc AZAÏS Institut de Mathématiques, Université de Toulouse Jean-Marc AZAÏS ( Institut de Mathématiques, Université de Toulouse )
More informationElectric Power Network Response to Random Demand Variations: A large deviations view on voltage collapse phenomena
Electric Power Network Response to Random Demand Variations: A large deviations view on voltage collapse phenomena Christopher DeMarco Electrical and Computer Engineering University of Wisconsin-Madison
More informationDynamic Consistency for Stochastic Optimal Control Problems
Dynamic Consistency for Stochastic Optimal Control Problems Cadarache Summer School CEA/EDF/INRIA 2012 Pierre Carpentier Jean-Philippe Chancelier Michel De Lara SOWG June 2012 Lecture outline Introduction
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 informationStochastic Integration and Stochastic Differential Equations: a gentle introduction
Stochastic Integration and Stochastic Differential Equations: a gentle introduction Oleg Makhnin New Mexico Tech Dept. of Mathematics October 26, 27 Intro: why Stochastic? Brownian Motion/ Wiener process
More informationExample 4.1 Let X be a random variable and f(t) a given function of time. Then. Y (t) = f(t)x. Y (t) = X sin(ωt + δ)
Chapter 4 Stochastic Processes 4. Definition In the previous chapter we studied random variables as functions on a sample space X(ω), ω Ω, without regard to how these might depend on parameters. We now
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 informationWeak convergence and large deviation theory
First Prev Next Go To Go Back Full Screen Close Quit 1 Weak convergence and large deviation theory Large deviation principle Convergence in distribution The Bryc-Varadhan theorem Tightness and Prohorov
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 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 informationNotes on Large Deviations in Economics and Finance. Noah Williams
Notes on Large Deviations in Economics and Finance Noah Williams Princeton University and NBER http://www.princeton.edu/ noahw Notes on Large Deviations 1 Introduction What is large deviation theory? Loosely:
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More informationErgodicity and Non-Ergodicity in Economics
Abstract An stochastic system is called ergodic if it tends in probability to a limiting form that is independent of the initial conditions. Breakdown of ergodicity gives rise to path dependence. We illustrate
More informationVariational approach to mean field games with density constraints
1 / 18 Variational approach to mean field games with density constraints Alpár Richárd Mészáros LMO, Université Paris-Sud (based on ongoing joint works with F. Santambrogio, P. Cardaliaguet and F. J. Silva)
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 informationChapter 2 Event-Triggered Sampling
Chapter Event-Triggered Sampling In this chapter, some general ideas and basic results on event-triggered sampling are introduced. The process considered is described by a first-order stochastic differential
More informationElliptic Operators with Unbounded Coefficients
Elliptic Operators with Unbounded Coefficients Federica Gregorio Universitá degli Studi di Salerno 8th June 2018 joint work with S.E. Boutiah, A. Rhandi, C. Tacelli Motivation Consider the Stochastic Differential
More informationLarge Deviations in Epidemiology : Malaria modelling
Kouegou Kamen Boris Supervised by : Etienne Pardoux Ecole de printemps de la Chaire MMB - Aussois May 13, 218 Figure: Internet... - Mosquitoes are the deadliest animals for Humans ( 745 pers./an) - Guess
More informationSUMMARY OF RESULTS ON PATH SPACES AND CONVERGENCE IN DISTRIBUTION FOR STOCHASTIC PROCESSES
SUMMARY OF RESULTS ON PATH SPACES AND CONVERGENCE IN DISTRIBUTION FOR STOCHASTIC PROCESSES RUTH J. WILLIAMS October 2, 2017 Department of Mathematics, University of California, San Diego, 9500 Gilman Drive,
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More information