degenerate Fokker-Planck equations with linear drift
|
|
- Paula Fisher
- 5 years ago
- Views:
Transcription
1 TECHNISCHE UNIVERSITÄT WIEN Entropy method for hypocoercive & non-symmetric Fokker-Planck equations with linear drift Anton ARNOLD with Jan Erb; Franz Achleitner Paris, April 2017 Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 1 / 32
2 degenerate Fokker-Planck equations with linear drift evolution of probability density f(x,t), x R n, t > 0: ( ) f t = div D f +Cx f f(x,0) = f 0 (x) (1) D R n n... symmetric, const in x, degenerate C R n n... const in x Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 2 / 32
3 degenerate Fokker-Planck equations with linear drift evolution of probability density f(x,t), x R n, t > 0: ( ) f t = div D f +Cx f f(x,0) = f 0 (x) (1) D R n n... symmetric, const in x, degenerate C R n n... const in x goals: existence & uniqueness of steady state (x); convergence f(t) t with sharp rates; complete theory for the equation class (1) Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 2 / 32
4 hypocoercive example from plasma physics kinetic Fokker-Planck equation for f(x,v,t), x,v R n : f t + v x f }{{} free transport x V v f }{{} influence of potential V(x) steady state: (x,v) = c e ν σ V(x)... given confinement potential [ v 2 2 +V(x) ] = σ v f }{{} +ν div v (vf) }{{} diffusion, σ>0 friction, ν>0 Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 3 / 32
5 hypocoercive example from plasma physics kinetic Fokker-Planck equation for f(x,v,t), x,v R n : f t + v x f }{{} free transport x V v f }{{} influence of potential V(x) steady state: (x,v) = c e ν σ V(x)... given confinement potential [ v 2 2 +V(x) ] rewritten (with x, v variables): [ ( ) 0 0 f t = div x,v 0 σi }{{} =:D...diffusion ( x,v f + explicit Green s function for V not quadratic! = σ v f }{{} +ν div v (vf) }{{} diffusion, σ>0 friction, ν>0 ) v x V +νv }{{} drift ] f Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 3 / 32
6 Outline: 1 hypocoercivity, prototypic examples 2 review of standard entropy method for non-degenerate Fokker-Planck equations 3 decay of modified entropy dissipation functional 4 mechanism of new method Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 4 / 32
7 (hypo)coercivity 1 example 1: standard Fokker-Planck equation on R n : ( ) f t = div f +x f =: Lf... symmetric on H := L 2 (f 1 ) (x) = ce x 2 2, kerl = span( ) L is dissipative, i.e. Lf,f H 0 f D(L) L is coercive (has a spectral gap), in the sense: Lf,f H f 2 L 2 (f 1 ) f { } Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 5 / 32
8 (hypo)coercivity 2 example 2: ( ) f t = div D f +Cx f =: Lf (2) with degenerate D is degenerate parabolic; (symmetric part of) L is not coercive. Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 6 / 32
9 (hypo)coercivity 2 example 2: ( ) f t = div D f +Cx f =: Lf (2) with degenerate D is degenerate parabolic; (symmetric part of) L is not coercive. Definition 1 (Villani 2009) Consider L on Hilbert space H with K = kerl; let H K (densely) (e.g. H... weighted L 2, H... weighted H 1 ). L is called hypocoercive on H if λ > 0, c 1: e Lt f H ce λt f H f H typically c > 1 Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 6 / 32
10 Steady state of (non)degenerate FP equations: standard Fokker-Planck equation f t = div( f +x f) : unique steady state (x) = c e x 2 /2 as a balance of drift & diffusion; sharp decay rate = 1 n = 2: x 2 drift diffusion x 1 Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 7 / 32
11 degenerate prototype: degenerate diffusion (1D Fokker-Planck) + rotation [ ( ) ( ) 1 0 x1 ωx f t = div f ωx 1 }{{}}{{} =D =Cx (x) = c e x 2 /2 ω R (unique for ω 0); sharp decay rate = 1 2 (= minrλ C) for fast enough rotation ( ω > 1 2 ) x 2 ] f x 1 equilibration by drift/diffusion Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 8 / 32
12 degenerate prototype: degenerate diffusion (1D Fokker-Planck) + rotation [ ( ) ( ) 1 0 x1 ωx f t = div f ωx 1 }{{}}{{} =D =Cx (x) = c e x 2 /2 ω R (unique for ω 0); sharp decay rate = 1 2 (= minrλ C) for fast enough rotation ( ω > 1 2 ) x 2 ] f x 1 equilibration by drift/diffusion ω = 1 2 : C has a Jordan block (sharp) decay rate = 1 2 ε Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 8 / 32
13 y degenerate prototype with ω = 1: [ ( 1 0 f t = div 0 0 ) f + ( ) x f x 2 -axis: drift characteristics of ẋ = Cx tangent to level curve of x : ] x = x + y y = x Drift characteristic Level curve of P norm x level curve of distorted vector norm x P x; P = Ref: [Dolbeault-Mouhot-Schmeiser] 2015 ( Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 9 / 32 )
14 coefficients C, D in Fokker-Planck equation ( ) f t = div D f +Cx f =: Lf Condition A: No (nontrivial) subspace of kerd is invariant under C. (equivalent: L is hypoelliptic.) Proposition 1 Let Condition A hold. a) Let f 0 L 1 (R d ) f C (R n R + ). [Hörmander 1969] b) Let f 0 L 1 +(R d ) f(x,t) > 0, t > 0. (Green s fct > 0) Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 10 / 32
15 coefficients C, D in Fokker-Planck equation ( ) f t = div D f +Cx f =: Lf Condition A: No (nontrivial) subspace of kerd is invariant under C. (equivalent: L is hypoelliptic.) Proposition 1 Let Condition A hold. a) Let f 0 L 1 (R d ) f C (R n R + ). [Hörmander 1969] b) Let f 0 L 1 +(R d ) f(x,t) > 0, t > 0. (Green s fct > 0) Condition B: Condition A + let C be positively stable (i.e. Rλ C > 0) confinement potential; drift towards x = 0. hypoelliptic + confinement = hypocoercive (for FP eq.) Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 10 / 32
16 steady state ( ) f t = div D f +Cx f (3) Theorem 2 (3) has a unique (normalized) steady state L 1 (R n ) iff Condition B holds. Then: (x) = c K e x K 1 x 2... non-isotropic Gaussian 0 < K R n n... unique solution of 2D = CK+KC (continuous Lyapunov equation) Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 11 / 32
17 normalization of Fokker-Planck equations f t = div ( D f +Cx f ) with (x) = c K e x K 1 x 2 transformations: 1 y := K 1 x g t = div y ( D y g + Cy g ) with g (x) = ce y 2 2, D = C S Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 12 / 32
18 normalization of Fokker-Planck equations f t = div ( D f +Cx f ) with (x) = c K e x K 1 x 2 transformations: 1 y := K 1 x g t = div y ( D y g + Cy g ) with g (x) = ce y 2 2, D = C S 2 rotation of y Ď = diag(d 1,..., d k, 0,..., 0) }{{} n k [normalization from now on assumed] Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 12 / 32
19 review of entropy method: linear symmetric Fokker-Planck equations evolution of probability density f(x,t), x R n, t > 0: ( ) f t = div D [ f +f A(x)] =: Lf f(x,0) = f 0 (x); f 0 L 1 +(R n ), f 0 dx = 1 f(x,t) 0 R n Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 13 / 32
20 review of entropy method: linear symmetric Fokker-Planck equations evolution of probability density f(x,t), x R n, t > 0: ( ) f t = div D [ f +f A(x)] =: Lf f(x,0) = f 0 (x); f 0 L 1 +(R n ), f 0 dx = 1 f(x,t) 0 R n (x) = e A(x)... (unique) normalized steady state ( Lf = div D f )... symmetric in L 2 (R n,f 1 ) D > 0... positive definite matrix A(x)... scalar confinement potential, i.e. A(x) as x ; idea : A(x) c x 2 Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 13 / 32
21 admissible relative entropies (for entropy method) for probability densities f 1,2 : e ψ (f 1 f 2 ) := ( f1 ψ R n ) f 2 dx 0 f 2 ψ : R + 0 R entropy generators ψ 0, ψ(1) = 0, ψ > 0, (ψ ) ψ ψ IV... relative entropy Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 14 / 32
22 admissible relative entropies (for entropy method) for probability densities f 1,2 : e ψ (f 1 f 2 ) := ( f1 ψ R n ) f 2 dx 0 f 2 ψ : R + 0 R entropy generators ψ 0, ψ(1) = 0, ψ > 0, (ψ ) ψ ψ IV... relative entropy examples: 1) ψ 1 (σ) = σlnσ σ +1 2) ψ p (σ) = σ p 1 p(σ 1), 1 < p 2 ψ(σ) ϕ...quadr. e ψ (f 1 f 2 ) = 0 f 1 = f 2 σ 1 Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 14 / 32
23 entropy dissipation Lemma 1 ( ) Let f(t) solve Fokker-Planck equation f t = div D [ f +f A(x)] d dt e ψ(f(t) ) = ψ R n ( f(t) ) f(t) D f(t) dx =: I ψ (f(t) ) 0... (negative) Fisher information Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 15 / 32
24 Step 1: exponent. decay of entropy dissipation for D const D = const. in x, (x) = e A(x) Theorem 3 Let I ψ (f 0 ) <. Let D,A satisfy a Bakry - Emery condition 2 A(x) x 2 λ 1 D 1 x R n (4) >0 I ψ (f(t) ) e 2λ 1t I ψ (f 0 ), t 0 A...uniformly convex if D = I Ref s: [Bakry-Emery] 1984/85; [Arnold-Markowich-Toscani-Unterreiter] Comm. PDE 2001 robust w.r.t. many nonlinear perturbations Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 16 / 32
25 Step 2: exponential decay of relative entropy for D const Theorem 4 Let D, A satisfy BEC 2 A(x) x 2 λ 1 D 1 e ψ (f(t) ) e 2λ 1t e ψ (f 0 ), t 0 Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 17 / 32
26 Step 2: exponential decay of relative entropy for D const Theorem 4 Let D, A satisfy BEC 2 A(x) x 2 λ 1 D 1 e ψ (f(t) ) e 2λ 1t e ψ (f 0 ), t 0 Proof: from proof of Theorem 3 : d dt I(t) 2λ 1 I(t) }{{} = e (t) t... dt Since I(t), e(t) t 0: d dt e(t) 2λ 1e(t) (5) (+ density argument) Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 17 / 32
27 problem of entropy decay (cp. standard entropy method) decay of quadratic entropy e 2 (t) = f(t) 2 L 2 (f 1 ) : standard Fokker-Planck equation: non-degenerate e(t) is convex; entropy dissip. e (t) < 0 f ; e µe possible (with µ > 0) Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 18 / 32
28 problem of entropy decay (cp. standard entropy method) decay of quadratic entropy e 2 (t) = f(t) 2 L 2 (f 1 ) : standard Fokker-Planck equation: non-degenerate e(t) is convex; entropy dissip. e (t) < 0 f ; e µe possible (with µ > 0) degenerate prototype ex.: e(t) is not convex; e (t) = 0 for some f ; e µe wrong (in general) Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 18 / 32
29 entropy decay in inhomogeneous Boltzmann equation simulation of 1+2 D Boltzmann equation: wavy entropy decay H g (t): relative entropy w.r.t. the global Maxwellian Ref: [Filbet-Mouhot-Pareschi] 2006 Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 19 / 32
30 new entropy method for degen. FP: f t = div(d f +Cx f) e (t) = 0 for some f entropy dissipation: ( ) d f dt e ψ = ψ f D R n }{{} f dx =: I ψ (f) 0 is useless as Lyapunov functional. 0 Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 20 / 32
31 new entropy method for degen. FP: f t = div(d f +Cx f) e (t) = 0 for some f entropy dissipation: ( ) d f dt e ψ = ψ f D R n }{{} f dx =: I ψ (f) 0 is useless as Lyapunov functional. define modified entropy dissipation as auxiliary functional: ( ) f S ψ (f) := f P }{{} f dx 0 R n ψ d goal: estimate between S(f(t)), dts(f(t)) for good choice of P > 0. Then: P c P D S ψ (f) c P I ψ (f) ց 0 0 >0 Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 20 / 32
32 modified entropy dissipation S ψ (f) : choice of P Lemma 2 Let Q be positively stable, i.e. µ := min{rλ Q } > 0. 1 If all λ min Q {λ σ(q) Rλ = µ} are non-defective (i.e. geometric = algebraic multiplicity) P R n n, P > 0 : PQ+Q P 2µP. Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 21 / 32
33 modified entropy dissipation S ψ (f) : choice of P Lemma 2 Let Q be positively stable, i.e. µ := min{rλ Q } > 0. 1 If all λ min Q {λ σ(q) Rλ = µ} are non-defective (i.e. geometric = algebraic multiplicity) P R n n, P > 0 : PQ+Q P 2µP. 2 If (at least) one λ min Q is defective ε > 0 P = P(ε) > 0 : PQ+Q P 2(µ ε)p. Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 21 / 32
34 modified entropy dissipation S ψ (f) : choice of P Lemma 2 Let Q be positively stable, i.e. µ := min{rλ Q } > 0. 1 If all λ min Q {λ σ(q) Rλ = µ} are non-defective (i.e. geometric = algebraic multiplicity) P R n n, P > 0 : PQ+Q P 2µP. 2 If (at least) one λ min Q is defective ε > 0 P = P(ε) > 0 : PQ+Q P 2(µ ε)p. Proof: P can be constructed explicitly; e.g. for Q non-defective / diagonalizable: n P := z j z j ; z j... eigenvectors of Q j=1 P not unique; but the decay rates µ (or µ ε) are independent of P. application with Q := C. Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 21 / 32
35 Step 1: exponential decay of auxiliary functional S ψ (f) S ψ (f) := R n ψ ( f ) f P f dx 0 Proposition 2 µ := min{rλ C }. Let f 0 satisfy: ( ) ψ f0 f 0 2 dx < ( weighted H 1 seminorm) 1 If all λ min C are non-defective S(f(t)) e 2µt S(f 0 ), t 0; Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 22 / 32
36 Step 1: exponential decay of auxiliary functional S ψ (f) S ψ (f) := R n ψ ( f ) f P f dx 0 Proposition 2 µ := min{rλ C }. Let f 0 satisfy: ( ) ψ f0 f 0 2 dx < ( weighted H 1 seminorm) 1 If all λ min C are non-defective S(f(t)) e 2µt S(f 0 ), t 0; 2 If one λ min C is defective S(f(t),ε) e 2(µ ε)t S(f 0,ε), t 0. Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 22 / 32
37 Proof of Proposition 2 modified entropy method d dt S(f(t)) = 2 ψ ( f ) [ ] u PC+C P }{{} 2µP... replaces BEC u dx Tr(XY) dx 2µS(f(t)); u = f }{{} 0 Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 23 / 32
38 Proof of Proposition 2 modified entropy method d dt S(f(t)) = 2 ψ ( f ) [ ] u PC+C P }{{} 2µP... replaces BEC u dx Tr(XY) dx 2µS(f(t)); u = f }{{} 0 use X := ( ψ ψ ψ 1 2 ψiv ) ( f ) 0, since det X = 1 2 ψ ψ IV (ψ ) 2 0 (for admissible rel. entropies) ; Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 23 / 32
39 Proof of Proposition 2 modified entropy method d dt S(f(t)) = 2 ψ ( f ) [ ] u PC+C P }{{} 2µP... replaces BEC u dx Tr(XY) dx 2µS(f(t)); u = f }{{} 0 use X := ( ψ ψ ψ 1 2 ψiv ) ( f ) 0, since det X = 1 2 ψ ψ IV (ψ ) 2 0 (for admissible rel. entropies) ; ( Tr(D u x P u x ) u D u x Pu ) Y := with Cauchy-Schwarz u D u x Pu (u Pu)(u Du) 0, Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 23 / 32
40 Step 2: exponential decay of relative entropy Theorem 5 Let f 0 satisfy: ( ) ψ f0 u 0 2 dx <. e(f(t) ) c S(f(t)) ce 2µt S(f 0 ), t 0 ( reduced rate for a defective λ min C : 2(µ ε) ) Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 24 / 32
41 Step 2: exponential decay of relative entropy Theorem 5 Let f 0 satisfy: ( ) ψ f0 u 0 2 dx <. e(f(t) ) c S(f(t)) ce 2µt S(f 0 ), t 0 ( reduced rate for a defective λ min C : 2(µ ε) ) Proof: Consider non-degenerate (auxiliary) symmetric FP equation: ( g t = div }{{} P ( g ) ) ; g = = c e x 2 /2 = c e A(x) (6) >0 It satisfies the Bakry-Emery condition 2 A x 2 = I λ P P 1. convex Sobolev inequality: e ψ (g ) 1 2λ P S ψ (g) g Remark: S ψ (g) is the true entropy dissipation for (6)! Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 24 / 32
42 (parabolic) regularization of semigroup e Lt Proposition 3 The Hörmander order m [1,n k] (k =rank D) is the minimum such that m C j D(C ) j κi for some κ > 0. j=0 (Existence of m is equivalent to Condition A, i.e. hypoellipticity of L.) S ψ (f(t)) c t (2m+1) e ψ (f 0 ), 0 < t 1. Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 25 / 32
43 (parabolic) regularization of semigroup e Lt Proposition 3 The Hörmander order m [1,n k] (k =rank D) is the minimum such that m C j D(C ) j κi for some κ > 0. j=0 (Existence of m is equivalent to Condition A, i.e. hypoellipticity of L.) S ψ (f(t)) c t (2m+1) e ψ (f 0 ), 0 < t 1. Ref s: Prop. 3 is generalization to all admissible relative entropies of: [Hérau] JFA 2007; [Villani] book 2009 (only for quadratic & logarithmic entropies) Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 25 / 32
44 exp. decay of rel. entropy for f t = div(d f +Cx f) =: Lf combination of regularization for initial time with Th.5 (entropy decay) Theorem 6 (Arnold-Erb 2014) Let L satisfy Condition B; µ := min{rλ C }. c > 0: e ψ (f(t) ) ce 2µt e ψ (f 0 ), t 0 ( reduced rate for a defective λ min C : 2(µ ε) ) Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 26 / 32
45 exp. decay of rel. entropy for f t = div(d f +Cx f) =: Lf combination of regularization for initial time with Th.5 (entropy decay) Theorem 6 (Arnold-Erb 2014) Let L satisfy Condition B; µ := min{rλ C }. c > 0: e ψ (f(t) ) ce 2µt e ψ (f 0 ), t 0 ( reduced rate for a defective λ min C : 2(µ ε) ) Proof: e(t) CSI 1 S(f(t)) decay 2λ P 1 e 2µ(t δ) S(f(δ)) regularization c(δ)e 2µt e(0) 2λ P Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 26 / 32
46 exp. decay of rel. entropy for f t = div(d f +Cx f) =: Lf combination of regularization for initial time with Th.5 (entropy decay) Theorem 6 (Arnold-Erb 2014) Let L satisfy Condition B; µ := min{rλ C }. c > 0: e ψ (f(t) ) ce 2µt e ψ (f 0 ), t 0 ( reduced rate for a defective λ min C : 2(µ ε) ) Proof: e(t) CSI 1 S(f(t)) decay 2λ P 1 e 2µ(t δ) S(f(δ)) regularization c(δ)e 2µt e(0) 2λ P Remark: Rate µ is sharp, but constant c is not. Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 26 / 32
47 kinetic Fokker-Planck eq. with non-quadratic potential f t +v x f x V(x) v f = σ v f +ν div v (vf); x,v R n steady state factors in x, v: (x,v) = c e ν σ [ v 2 2 +V(x) ] Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 27 / 32
48 kinetic Fokker-Planck eq. with non-quadratic potential f t +v x f x V(x) v f = σ v f +ν div v (vf); x,v R n steady state factors in x, v: (x,v) = c e ν σ f t = div x,v [ ( 0 σ ν I σ ν I σi [ v 2 2 +V(x) ] ) x,v ( f ) f ] Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 27 / 32
49 kinetic Fokker-Planck eq. with non-quadratic potential f t +v x f x V(x) v f = σ v f +ν div v (vf); x,v R n steady state factors in x, v: (x,v) = c e ν σ f t = div x,v [ ( 0 σ ν I σ ν I σi [ v 2 2 +V(x) ] ) x,v ( f ) f ] Theorem 7 (Arnold-Erb, Achleitner-Arnold 2015) Let n = 1; V(x) = ω2 0 2 x 2 +Ṽ(x)... given confinement potential with maxv (x) minv (x) ν e ψ (f(t) ) c S ψ (f 0 )e 2κt, t 0 Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 27 / 32
50 kinetic Fokker-Planck eq. with non-quadratic potential f t +v x f x V(x) v f = σ v f +ν div v (vf); x,v R n steady state factors in x, v: (x,v) = c e ν σ f t = div x,v [ ( 0 σ ν I σ ν I σi [ v 2 2 +V(x) ] ) x,v ( f ) f ] Theorem 7 (Arnold-Erb, Achleitner-Arnold 2015) Let n = 1; V(x) = ω2 0 2 x 2 +Ṽ(x)... given confinement potential with maxv (x) minv (x) ν e ψ (f(t) ) c S ψ (f 0 )e 2κt, t 0 Ṽ... O(1) perturbation Greens function not explicit Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 27 / 32
51 local vs. global decay rate in non-symmetric FP equations decay of logarithmic entropy for the non-symmetric FP equation [ ( ) ( ) 1/4 0 1/4x1 4x ] f t = div f + 2 f 0 1 4x 1 +x 2 with (x) = ce x 2 /2 = e A(x) : find e(f(t) ) c e(0)e λt Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 28 / 32
52 local vs. global decay rate in non-symmetric FP equations decay of logarithmic entropy for the non-symmetric FP equation [ ( ) ( ) 1/4 0 1/4x1 4x ] f t = div f + 2 f 0 1 4x 1 +x 2 with (x) = ce x 2 /2 = e A(x) : find e(f(t) ) c e(0)e λt standard entropy method: BEC 2 A = I λd 1 yields sharp x 2 local decay rate: λ l = 1 4 ; c = 1 Hypocoercive method yields sharp global decay rate: λ g = 5 8 ; c > 1 Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 28 / 32
53 local vs. global decay rate in non-symmetric FP equations decay of logarithmic entropy for the non-symmetric FP equation [ ( ) ( ) 1/4 0 1/4x1 4x ] f t = div f + 2 f 0 1 4x 1 +x 2 with (x) = ce x 2 /2 = e A(x) : find e(f(t) ) c e(0)e λt standard entropy method: BEC 2 A = I λd 1 yields sharp x 2 local decay rate: λ l = 1 4 ; c = 1 Hypocoercive method yields sharp global decay rate: λ g = 5 8 ; c > 1 Lemma: 2D, admissible P: multiplicative constant S(0) 2λ P is sharp. Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 28 / 32
54 Why does the hypocoercive method work? algebraic essence: comparison of the spectral gaps of a non-symmetric matrix Q and its symmetric part Q s := 1 2 (Q+QT ). motivation: entropy method operates mostly with quadratic functionals (e.g. e 2 (f ); I ψ = ψ ( f ) T f D f dx ) Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 29 / 32
55 Why does the hypocoercive method work? algebraic essence: comparison of the spectral gaps of a non-symmetric matrix Q and its symmetric part Q s := 1 2 (Q+QT ). motivation: entropy method operates mostly with quadratic functionals (e.g. e 2 (f ); I ψ = ψ ( f ) T f D f dx ) Lemma 3 Let Q be positively stable, i.e. µ := min{rλ Q } > 0. λ min (Q s ) µ ( typically strict inequality) Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 29 / 32
56 Why does the hypocoercive method work? algebraic essence: comparison of the spectral gaps of a non-symmetric matrix Q and its symmetric part Q s := 1 2 (Q+QT ). motivation: entropy method operates mostly with quadratic functionals (e.g. e 2 (f ); I ψ = ψ ( f ) T f D f dx ) Lemma 3 Let Q be positively stable, i.e. µ := min{rλ Q } > 0. λ min (Q s ) µ ( typically strict inequality) Lemma 4 ( = Lemma 2 for choice of P) P > 0 such that the similar matrix Q := PQ P 1 satisfies: λ min ( Q s ) = µ ( λ min ( Q s ) = µ ε in the defective case ) So: a similarity transformation such that Q, Q, Q s have the same spectral gap. Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 29 / 32
57 Application to f t = div(d f +Cxf) = Lf Decay estimate of the drift characteristics ẋ = Cx: Let x 2 P := x,px. d dt x 2 P = 2x T PCx = ( Px) T ( PC P 1 + P 1 C T ) P ( Px) }{{} :=2 C s 2µI 2µ x 2 P carries over to all invariant eigenspaces of L Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 30 / 32
58 Conclusion new entropy method for degenerate Fokker-Planck eq. (+ linear drift) key tool: modified entropy dissipation functional diff. inequal. Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 31 / 32
59 Conclusion new entropy method for degenerate Fokker-Planck eq. (+ linear drift) key tool: modified entropy dissipation functional diff. inequal. extension to kinetic Fokker-Planck eq. (with non-quadratic potential) extension to non-degenerate, non-symmetric Fokker-Planck equations sharp envelope for global decay 2017: (1+t 2 )e 2κt decay for defective case [AA-Einav-Wöhrer] Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 31 / 32
60 Conclusion new entropy method for degenerate Fokker-Planck eq. (+ linear drift) key tool: modified entropy dissipation functional diff. inequal. extension to kinetic Fokker-Planck eq. (with non-quadratic potential) extension to non-degenerate, non-symmetric Fokker-Planck equations sharp envelope for global decay 2017: (1+t 2 )e 2κt decay for defective case [AA-Einav-Wöhrer] References A. Arnold, J. Erb: Sharp entropy decay for hypocoercive and non-symmetric Fokker-Planck equations with linear drift, arxiv 2014 F. Achleitner, A. Arnold, D. Stürzer: Large-time behavior in non-symmetric Fokker-Planck equations, Rivista di Matematica della Univ. di Parma, 2015 F. Achleitner, A. Arnold, E. Carlen: On linear hypocoercive BGK models, Springer Proc. in Math. & Stat., 2016 Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 31 / 32
61 References: [Villani] 2009: exponential decay in weighted H 1, but no sharp rates [Dolbeault-Mouhot-Schmeiser] Trans. AMS 2015: kinetic models, exponential decay in modified L 2 norm; m = 1 [Gadat-Miclo] KRM 2013: sharp rates for 2 Fokker-Planck toy models [Baudoin] 2014: Γ 2 formalism, includes auxiliary gradient functional Anton ARNOLD (TU Vienna) hypocoercive Fokker-Planck/entropy meth. 32 / 32
Hypocoercivity for kinetic equations with linear relaxation terms
Hypocoercivity for kinetic equations with linear relaxation terms Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (A JOINT
More informationUne approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck
Une approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (EN
More informationEntropy-dissipation methods I: Fokker-Planck equations
1 Entropy-dissipation methods I: Fokker-Planck equations Ansgar Jüngel Vienna University of Technology, Austria www.jungel.at.vu Introduction Boltzmann equation Fokker-Planck equations Degenerate parabolic
More informationON THE RATES OF DECAY TO EQUILIBRIUM IN DEGENERATE AND DEFECTIVE FOKKER-PLANCK EQUATIONS 1. INTRODUCTION
ON THE RATES OF DECAY TO EQUILIBRIUM IN DEGENERATE AND DEFECTIVE FOKKER-PLANCK EQUATIONS ANTON ARNOLD, AMIT EINAV & TOBIAS WÖHRER ABSTRACT. We establish sharp long time asymptotic behaviour for a family
More informationHypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th
Hypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th Department of Mathematics, University of Wisconsin Madison Venue: van Vleck Hall 911 Monday,
More informationSemigroup factorization and relaxation rates of kinetic equations
Semigroup factorization and relaxation rates of kinetic equations Clément Mouhot, University of Cambridge Analysis and Partial Differential Equations seminar University of Sussex 24th of february 2014
More informationEntropic structure of the Landau equation. Coulomb interaction
with Coulomb interaction Laurent Desvillettes IMJ-PRG, Université Paris Diderot May 15, 2017 Use of the entropy principle for specific equations Spatially Homogeneous Kinetic equations: 1 Fokker-Planck:
More informationDecay rates for partially dissipative hyperbolic systems
Outline Decay rates for partially dissipative hyperbolic systems Basque Center for Applied Mathematics Bilbao, Basque Country, Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ Numerical Methods
More informationEntropy Methods for Reaction-Diffusion Equations with Degenerate Diffusion Arising in Reversible Chemistry
1 Entropy Methods for Reaction-Diffusion Equations with Degenerate Diffusion Arising in Reversible Chemistry L. Desvillettes CMLA, ENS Cachan, IUF & CNRS, PRES UniverSud, 61, avenue du Président Wilson,
More informationINTERPOLATION BETWEEN LOGARITHMIC SOBOLEV AND POINCARÉ INEQUALITIES
INTERPOLATION BETWEEN LOGARITHMIC SOBOLEV AND POINCARÉ INEQUALITIES ANTON ARNOLD, JEAN-PHILIPPE BARTIER, AND JEAN DOLBEAULT Abstract. This paper is concerned with intermediate inequalities which interpolate
More informationQuantum Fokker-Planck models: kinetic & operator theory approaches
TECHNISCHE UNIVERSITÄT WIEN Quantum Fokker-Planck models: kinetic & operator theory approaches Anton ARNOLD with F. Fagnola (Milan), L. Neumann (Vienna) TU Vienna Institute for Analysis and Scientific
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationWorkshop on CENTRAL Trends in PDEs Vienna, November 12-13, 2015
Workshop on CENTRAL Trends in PDEs Vienna, November 12-13, 2015 Venue: Sky Lounge, 12th floor, Faculty of Mathematics, Oskar-Morgenstern-Platz 1, Vienna Thursday, November 12, 2015 14:00-14:30 Registration
More informationA REVIEW ON RESULTS FOR THE DERRIDA-LEBOWITZ-SPEER-SPOHN EQUATION
1 A REVIEW ON RESULTS FOR THE DERRIDA-LEBOWITZ-SPEER-SPOHN EQUATION Ansgar Jüngel Institut für Analysis und Scientific Computing, Technische Universität Wien, Wiedner Hauptstr. 8-10, 1040 Wien, Austria;
More informationDiscrete entropy methods for nonlinear diffusive evolution equations
Discrete entropy methods for nonlinear diffusive evolution equations 1 Ansgar Jüngel Vienna University of Technology, Austria www.jungel.at.vu Joint work with E. Emmrich (TU Berlin), M. Bukal (Zagreb),
More informationHypocoercivity and exponential time decay for the linear inhomogeneous relaxation Boltzmann equation
Hypocoercivity and exponential time decay for the linear inhomogeneous relaxation Boltzmann equation Frédéric Hérau Université de Reims March 16, 2005 Abstract We consider an inhomogeneous linear Boltzmann
More informationUne approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck
Une approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (EN
More informationExplosive Solution of the Nonlinear Equation of a Parabolic Type
Int. J. Contemp. Math. Sciences, Vol. 4, 2009, no. 5, 233-239 Explosive Solution of the Nonlinear Equation of a Parabolic Type T. S. Hajiev Institute of Mathematics and Mechanics, Acad. of Sciences Baku,
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 informationSpeeding up Convergence to Equilibrium for Diffusion Processes
Speeding up Convergence to Equilibrium for Diffusion Processes G.A. Pavliotis Department of Mathematics Imperial College London Joint Work with T. Lelievre(CERMICS), F. Nier (CERMICS), M. Ottobre (Warwick),
More informationThe goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T
1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.
More informationQUANTUM MODELS FOR SEMICONDUCTORS AND NONLINEAR DIFFUSION EQUATIONS OF FOURTH ORDER
QUANTUM MODELS FOR SEMICONDUCTORS AND NONLINEAR DIFFUSION EQUATIONS OF FOURTH ORDER MARIA PIA GUALDANI The modern computer and telecommunication industry relies heavily on the use of semiconductor devices.
More informationA quantum heat equation 5th Spring School on Evolution Equations, TU Berlin
A quantum heat equation 5th Spring School on Evolution Equations, TU Berlin Mario Bukal A. Jüngel and D. Matthes ACROSS - Centre for Advanced Cooperative Systems Faculty of Electrical Engineering and Computing
More informationConvergence to equilibrium of Markov processes (eventually piecewise deterministic)
Convergence to equilibrium of Markov processes (eventually piecewise deterministic) A. Guillin Université Blaise Pascal and IUF Rennes joint works with D. Bakry, F. Barthe, F. Bolley, P. Cattiaux, R. Douc,
More informationAnomalous transport of particles in Plasma physics
Anomalous transport of particles in Plasma physics L. Cesbron a, A. Mellet b,1, K. Trivisa b, a École Normale Supérieure de Cachan Campus de Ker Lann 35170 Bruz rance. b Department of Mathematics, University
More informationMildly degenerate Kirchhoff equations with weak dissipation: global existence and time decay
arxiv:93.273v [math.ap] 6 Mar 29 Mildly degenerate Kirchhoff equations with weak dissipation: global existence and time decay Marina Ghisi Università degli Studi di Pisa Dipartimento di Matematica Leonida
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationALMOST EXPONENTIAL DECAY NEAR MAXWELLIAN
ALMOST EXPONENTIAL DECAY NEAR MAXWELLIAN ROBERT M STRAIN AND YAN GUO Abstract By direct interpolation of a family of smooth energy estimates for solutions near Maxwellian equilibrium and in a periodic
More informationCOMMENT ON : HYPERCONTRACTIVITY OF HAMILTON-JACOBI EQUATIONS, BY S. BOBKOV, I. GENTIL AND M. LEDOUX
COMMENT ON : HYPERCONTRACTIVITY OF HAMILTON-JACOBI EQUATIONS, BY S. BOBKOV, I. GENTIL AND M. LEDOUX F. OTTO AND C. VILLANI In their remarkable work [], Bobkov, Gentil and Ledoux improve, generalize and
More informationhal , version 1-22 Nov 2009
Author manuscript, published in "Kinet. Relat. Models 1, 3 8) 355-368" PROPAGATION OF GEVREY REGULARITY FOR SOLUTIONS OF LANDAU EQUATIONS HUA CHEN, WEI-XI LI AND CHAO-JIANG XU Abstract. By using the energy-type
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationLogarithmic Sobolev Inequalities
Logarithmic Sobolev Inequalities M. Ledoux Institut de Mathématiques de Toulouse, France logarithmic Sobolev inequalities what they are, some history analytic, geometric, optimal transportation proofs
More informationIntertwinings for Markov processes
Intertwinings for Markov processes Aldéric Joulin - University of Toulouse Joint work with : Michel Bonnefont - Univ. Bordeaux Workshop 2 Piecewise Deterministic Markov Processes ennes - May 15-17, 2013
More informationOn the bang-bang property of time optimal controls for infinite dimensional linear systems
On the bang-bang property of time optimal controls for infinite dimensional linear systems Marius Tucsnak Université de Lorraine Paris, 6 janvier 2012 Notation and problem statement (I) Notation: X (the
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 informationHolder regularity for hypoelliptic kinetic equations
Holder regularity for hypoelliptic kinetic equations Alexis F. Vasseur Joint work with François Golse, Cyril Imbert, and Clément Mouhot The University of Texas at Austin Kinetic Equations: Modeling, Analysis
More informationOrdinary Differential Equations II
Ordinary Differential Equations II February 23 2017 Separation of variables Wave eq. (PDE) 2 u t (t, x) = 2 u 2 c2 (t, x), x2 c > 0 constant. Describes small vibrations in a homogeneous string. u(t, x)
More informationExponential approach to equilibrium for a stochastic NLS
Exponential approach to equilibrium for a stochastic NLS CNRS and Cergy Bonn, Oct, 6, 2016 I. Stochastic NLS We start with the dispersive nonlinear Schrödinger equation (NLS): i u = u ± u p 2 u, on the
More informationConservative Control Systems Described by the Schrödinger Equation
Conservative Control Systems Described by the Schrödinger Equation Salah E. Rebiai Abstract An important subclass of well-posed linear systems is formed by the conservative systems. A conservative system
More informationanalysis for transport equations and applications
Multi-scale analysis for transport equations and applications Mihaï BOSTAN, Aurélie FINOT University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Numerical methods for kinetic equations Strasbourg
More informationSectorial Forms and m-sectorial Operators
Technische Universität Berlin SEMINARARBEIT ZUM FACH FUNKTIONALANALYSIS Sectorial Forms and m-sectorial Operators Misagheh Khanalizadeh, Berlin, den 21.10.2013 Contents 1 Bounded Coercive Sesquilinear
More informationVelocity averaging a general framework
Outline Velocity averaging a general framework Martin Lazar BCAM ERC-NUMERIWAVES Seminar May 15, 2013 Joint work with D. Mitrović, University of Montenegro, Montenegro Outline Outline 1 2 L p, p >= 2 setting
More informationWEAK POINCARÉ INEQUALITIES FOR CONVERGENCE RATE OF DEGENERATE DIFFUSION PROCESSES
Submitted to the Annals of Probability arxiv: 1703.04821 WEAK POINCARÉ INEQUALITIES FOR CONVERGENCE RATE OF DEGENERATE DIFFUSION PROCESSES By Martin Grothaus and Feng Yu Wang, Mathematics Department, TU
More informationThe Kato square root problem on vector bundles with generalised bounded geometry
The Kato square root problem on vector bundles with generalised bounded geometry Lashi Bandara (Joint work with Alan McIntosh, ANU) Centre for Mathematics and its Applications Australian National University
More informationThe Kato square root problem on vector bundles with generalised bounded geometry
The Kato square root problem on vector bundles with generalised bounded geometry Lashi Bandara (Joint work with Alan McIntosh, ANU) Centre for Mathematics and its Applications Australian National University
More informationLinear Algebra Practice Problems
Linear Algebra Practice Problems Math 24 Calculus III Summer 25, Session II. Determine whether the given set is a vector space. If not, give at least one axiom that is not satisfied. Unless otherwise stated,
More informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationDiscrete entropy methods for nonlinear diffusive equations
Discrete entropy methods for nonlinear diffusive equations 1 Ansgar Jüngel Vienna University of Technology, Austria www.jungel.at.vu Joint with C. Chainais-Hillairet (Lille), J. A. Carrillo (Imperial),
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
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 informationcurvature, mixing, and entropic interpolation Simons Feb-2016 and CSE 599s Lecture 13
curvature, mixing, and entropic interpolation Simons Feb-2016 and CSE 599s Lecture 13 James R. Lee University of Washington Joint with Ronen Eldan (Weizmann) and Joseph Lehec (Paris-Dauphine) Markov chain
More informationACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017
ACM/CMS 17 Linear Analysis & Applications Fall 217 Assignment 2: PDEs and Finite Element Methods Due: 7th November 217 For this assignment the following MATLAB code will be required: Introduction http://wwwmdunloporg/cms17/assignment2zip
More informationKinetic theory of gases
Kinetic theory of gases Toan T. Nguyen Penn State University http://toannguyen.org http://blog.toannguyen.org Graduate Student seminar, PSU Jan 19th, 2017 Fall 2017, I teach a graduate topics course: same
More informationLARGE-TIME BEHAVIOR OF NON-SYMMETRIC FOKKER-PLANCK TYPE EQUATIONS
Serials Publications www.serialspublications.com LARGE-TIME BEHAVIOR OF NON-SYMMETRIC FOKKER-PLANCK TYPE EQUATIONS ANTON ARNOLD, ERIC CARLEN, AND QIANGCHANG JU Abstract. Large time asymptotics of the solutions
More information2 Formal derivation of the Shockley-Read-Hall model
We consider a semiconductor crystal represented by the bounded domain R 3 (all our results are easily extended to the one and two- dimensional situations) with a constant (in space) number density of traps
More informationSome New Results on Global Nonexistence for Abstract Evolution Equations with Positive Initial Energy
Some New Results on Global Nonexistence for Abstract Evolution Equations with Positive Initial Energy Patrizia Pucci 1 & James Serrin Dedicated to Olga Ladyzhenskaya with admiration and esteem 1. Introduction.
More informationInégalités spectrales pour le contrôle des EDP linéaires : groupe de Schrödinger contre semigroupe de la chaleur.
Inégalités spectrales pour le contrôle des EDP linéaires : groupe de Schrödinger contre semigroupe de la chaleur. Luc Miller Université Paris Ouest Nanterre La Défense, France Pde s, Dispersion, Scattering
More informationLecture 1: Review of linear algebra
Lecture 1: Review of linear algebra Linear functions and linearization Inverse matrix, least-squares and least-norm solutions Subspaces, basis, and dimension Change of basis and similarity transformations
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More informationA particle-in-cell method with adaptive phase-space remapping for kinetic plasmas
A particle-in-cell method with adaptive phase-space remapping for kinetic plasmas Bei Wang 1 Greg Miller 2 Phil Colella 3 1 Princeton Institute of Computational Science and Engineering Princeton University
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 informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
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 informationLinear-quadratic control problem with a linear term on semiinfinite interval: theory and applications
Linear-quadratic control problem with a linear term on semiinfinite interval: theory and applications L. Faybusovich T. Mouktonglang Department of Mathematics, University of Notre Dame, Notre Dame, IN
More informationIr O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v )
Section 3.2 Theorem 3.6. Let A be an m n matrix of rank r. Then r m, r n, and, by means of a finite number of elementary row and column operations, A can be transformed into the matrix ( ) Ir O D = 1 O
More informationEigenvectors. Prop-Defn
Eigenvectors Aim lecture: The simplest T -invariant subspaces are 1-dim & these give rise to the theory of eigenvectors. To compute these we introduce the similarity invariant, the characteristic polynomial.
More informationAn Observation on the Positive Real Lemma
Journal of Mathematical Analysis and Applications 255, 48 49 (21) doi:1.16/jmaa.2.7241, available online at http://www.idealibrary.com on An Observation on the Positive Real Lemma Luciano Pandolfi Dipartimento
More informationStabilization of second order evolution equations with unbounded feedback with delay
Stabilization of second order evolution equations with unbounded feedback with delay S. Nicaise and J. Valein snicaise,julie.valein@univ-valenciennes.fr Laboratoire LAMAV, Université de Valenciennes et
More information6 Inner Product Spaces
Lectures 16,17,18 6 Inner Product Spaces 6.1 Basic Definition Parallelogram law, the ability to measure angle between two vectors and in particular, the concept of perpendicularity make the euclidean space
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationDiscrete Ricci curvature: Open problems
Discrete Ricci curvature: Open problems Yann Ollivier, May 2008 Abstract This document lists some open problems related to the notion of discrete Ricci curvature defined in [Oll09, Oll07]. Do not hesitate
More informationPerturbation theory, KAM theory and Celestial Mechanics 7. KAM theory
Perturbation theory, KAM theory and Celestial Mechanics 7. KAM theory Alessandra Celletti Department of Mathematics University of Roma Tor Vergata Sevilla, 25-27 January 2016 Outline 1. Introduction 2.
More informationNonlinear Diffusion. Journal Club Presentation. Xiaowei Zhou
1 / 41 Journal Club Presentation Xiaowei Zhou Department of Electronic and Computer Engineering The Hong Kong University of Science and Technology 2009-12-11 2 / 41 Outline 1 Motivation Diffusion process
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling
1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples
More informationORNSTEIN-UHLENBECK PROCESSES ON LIE GROUPS
ORNSTEIN-UHLENBECK PROCESSES ON LIE ROUPS FABRICE BAUDOIN, MARTIN HAIRER, JOSEF TEICHMANN Abstract. We consider Ornstein-Uhlenbeck processes (OU-processes related to hypoelliptic diffusion on finite-dimensional
More informationInvariances in spectral estimates. Paris-Est Marne-la-Vallée, January 2011
Invariances in spectral estimates Franck Barthe Dario Cordero-Erausquin Paris-Est Marne-la-Vallée, January 2011 Notation Notation Given a probability measure ν on some Euclidean space, the Poincaré constant
More informationMathematical modelling of collective behavior
Mathematical modelling of collective behavior Young-Pil Choi Fakultät für Mathematik Technische Universität München This talk is based on joint works with José A. Carrillo, Maxime Hauray, and Samir Salem
More informationEigenvalues and Eigenvectors
/88 Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Eigenvalue Problem /88 Eigenvalue Equation By definition, the eigenvalue equation for matrix
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationHeat Flows, Geometric and Functional Inequalities
Heat Flows, Geometric and Functional Inequalities M. Ledoux Institut de Mathématiques de Toulouse, France heat flow and semigroup interpolations Duhamel formula (19th century) pde, probability, dynamics
More informationClassical inequalities for the Boltzmann collision operator.
Classical inequalities for the Boltzmann collision operator. Ricardo Alonso The University of Texas at Austin IPAM, April 2009 In collaboration with E. Carneiro and I. Gamba Outline Boltzmann equation:
More informationMath 350 Fall 2011 Notes about inner product spaces. In this notes we state and prove some important properties of inner product spaces.
Math 350 Fall 2011 Notes about inner product spaces In this notes we state and prove some important properties of inner product spaces. First, recall the dot product on R n : if x, y R n, say x = (x 1,...,
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 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 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 informationA note on W 1,p estimates for quasilinear parabolic equations
200-Luminy conference on Quasilinear Elliptic and Parabolic Equations and Systems, Electronic Journal of Differential Equations, Conference 08, 2002, pp 2 3. http://ejde.math.swt.edu or http://ejde.math.unt.edu
More informationConvex Optimization M2
Convex Optimization M2 Lecture 3 A. d Aspremont. Convex Optimization M2. 1/49 Duality A. d Aspremont. Convex Optimization M2. 2/49 DMs DM par email: dm.daspremont@gmail.com A. d Aspremont. Convex Optimization
More informationTraveling waves of a kinetic transport model for the KPP-Fisher equation
Traveling waves of a kinetic transport model for the KPP-Fisher equation Christian Schmeiser Universität Wien and RICAM homepage.univie.ac.at/christian.schmeiser/ Joint work with C. Cuesta (Bilbao), S.
More informationHeat equation and the sharp Young s inequality
Noname manuscript No. will be inserted by the editor) Heat equation and the sharp Young s inequality Giuseppe Toscani the date of receipt and acceptance should be inserted later Abstract We show that the
More informationDISSIPATIVE TRANSPORT APERIODIC SOLIDS KUBO S FORMULA. and. Jean BELLISSARD 1 2. Collaborations:
Vienna February 1st 2005 1 DISSIPATIVE TRANSPORT and KUBO S FORMULA in APERIODIC SOLIDS Jean BELLISSARD 1 2 Georgia Institute of Technology, Atlanta, & Institut Universitaire de France Collaborations:
More informationFollow links Class Use and other Permissions. For more information, send to:
COPYRIGHT NOTICE: Kari Astala, Tadeusz Iwaniec & Gaven Martin: Elliptic Partial Differential Equations and Quasiconformal Mappings in the Plane is published by Princeton University Press and copyrighted,
More informationHypocoercivity. Cédric Villani
Hypocoercivity Cédric Villani Unité de Mathématiques Pures et Appliquées, UMR CNRS 5669, Ecole Normale Supérieure de Lyon, 46 allée d Italie, F-69364 Lyon Cedex 07, FRANCE E-mail address: cvillani@umpa.ens-lyon.fr
More informationEntropy method and the large-time behavior of parabolic equations
Lecture notes for XXVII Summer school in mathematical physics Ravello, Italy, September 9-1, 00 Entropy method and the large-time behavior of parabolic equations Anton Arnold Institut für Numerische und
More informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
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 information16. Local theory of regular singular points and applications
16. Local theory of regular singular points and applications 265 16. Local theory of regular singular points and applications In this section we consider linear systems defined by the germs of meromorphic
More informationj=1 [We will show that the triangle inequality holds for each p-norm in Chapter 3 Section 6.] The 1-norm is A F = tr(a H A).
Math 344 Lecture #19 3.5 Normed Linear Spaces Definition 3.5.1. A seminorm on a vector space V over F is a map : V R that for all x, y V and for all α F satisfies (i) x 0 (positivity), (ii) αx = α x (scale
More information