Spatial decay of rotating waves in parabolic systems

Size: px
Start display at page:

Download "Spatial decay of rotating waves in parabolic systems"

Transcription

1 Spatial decay of rotating waves in parabolic systems Nonlinear Waves, CRC 701, Bielefeld, June 19, 2013 Denny Otten Department of Mathematics Bielefeld University Germany June 19, 2013 CRC 701 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

2 Outline 1 Introduction: Rotating pattern in R d 2 Main result: Exponential decay of v 3 Outline of proof: Exponential decay of v 4 Spectrum of rotating waves Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

3 Rotating Patterns in R d Consider a reaction diffusion system (1) u t (x,t) = A u(x,t)+f(u(x,t)), t > 0, x R d, d 2, u(x,0) = u 0 (x), t = 0, x R d. where u : R d [0, [ R N, A R N,N, f C 2 (R N,R N ). Assume a rotating wave solution u : R d [0, [ R N of (1) u (x,t) = v (e ts x) v : R d R N profile (pattern), 0 S R d,d skew-symmetric. Transformation (into a rotating frame): v(x,t) = u(e ts x,t) solves (2) v t (x,t) = A v(x,t)+ Sx, v(x,t) +f(v(x,t)), t > 0, x R d, d 2, v(x,0) = u 0 (x), t = 0, x R d. Sx, v(x) := d d S ij x j D i v(x) (drift term). i=1 j=1 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

4 Rotating Patterns in R d Consider a reaction diffusion system (1) u t (x,t) = A u(x,t)+f(u(x,t)), t > 0, x R d, d 2, u(x,0) = u 0 (x), t = 0, x R d. where u : R d [0, [ R N, A R N,N, f C 2 (R N,R N ). Assume a rotating wave solution u : R d [0, [ R N of (1) u (x,t) = v (e ts x) v : R d R N profile (pattern), 0 S R d,d skew-symmetric. Transformation (into a rotating frame): v(x,t) = u(e ts x,t) solves (2) v t (x,t) = A v(x,t)+ Sx, v(x,t) +f(v(x,t)), t > 0, x R d, d 2, v(x,0) = u 0 (x), t = 0, x R d. Sx, v(x) := d d S ij x j D i v(x) (drift term). i=1 j=1 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

5 Rotating Patterns in R d Consider a reaction diffusion system (1) u t (x,t) = A u(x,t)+f(u(x,t)), t > 0, x R d, d 2, u(x,0) = u 0 (x), t = 0, x R d. where u : R d [0, [ R N, A R N,N, f C 2 (R N,R N ). Assume a rotating wave solution u : R d [0, [ R N of (1) u (x,t) = v (e ts x) v : R d R N profile (pattern), 0 S R d,d skew-symmetric. Transformation (into a rotating frame): v(x,t) = u(e ts x,t) solves (2) v t (x,t) = A v(x,t)+ Sx, v(x,t) +f(v(x,t)), t > 0, x R d, d 2, v(x,0) = u 0 (x), t = 0, x R d. Sx, v(x) := d d S ij x j D i v(x) (drift term). i=1 j=1 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

6 Rotating Patterns in R d Consider a reaction diffusion system (1) u t (x,t) = A u(x,t)+f(u(x,t)), t > 0, x R d, d 2, u(x,0) = u 0 (x), t = 0, x R d. where u : R d [0, [ R N, A R N,N, f C 2 (R N,R N ). Assume a rotating wave solution u : R d [0, [ R N of (1) u (x,t) = v (e ts x) v : R d R N profile (pattern), 0 S R d,d skew-symmetric. Transformation (into a rotating frame): v(x,t) = u(e ts x,t) solves (2) v t (x,t) = A v(x,t)+ Sx, v(x,t) +f(v(x,t)), t > 0, x R d, d 2, v(x,0) = u 0 (x), t = 0, x R d. Sx, v(x) := d 1 d i=1 j=i+1 S ij (x j D i x i D j )v(x) (rotational term). Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

7 Rotating Patterns in R d Consider a reaction diffusion system (1) u t (x,t) = A u(x,t)+f(u(x,t)), t > 0, x R d, d 2, u(x,0) = u 0 (x), t = 0, x R d. where u : R d [0, [ R N, A R N,N, f C 2 (R N,R N ). Assume a rotating wave solution u : R d [0, [ R N of (1) u (x,t) = v (e ts x) v : R d R N profile (pattern), 0 S R d,d skew-symmetric. Transformation (into a rotating frame): v(x,t) = u(e ts x,t) solves (2) v t (x,t) = A v(x,t)+ Sx, v(x,t) +f(v(x,t)), t > 0, x R d, d 2, v(x,0) = u 0 (x), t = 0, x R d. v is a stationary solution of (2). Question: How to show exponential decay of v at x =? Consequence: Exponentially small error by restriction to bounded domain. Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

8 Rotating Patterns in R d Consider a reaction diffusion system (1) u t (x,t) = A u(x,t)+f(u(x,t)), t > 0, x R d, d 2, u(x,0) = u 0 (x), t = 0, x R d. where u : R d [0, [ R N, A R N,N, f C 2 (R N,R N ). Assume a rotating wave solution u : R d [0, [ R N of (1) u (x,t) = v (e ts x) v : R d R N profile (pattern), 0 S R d,d skew-symmetric. Transformation (into a rotating frame): v(x,t) = u(e ts x,t) solves (2) v t (x,t) = A v(x,t)+ Sx, v(x,t) +f(v(x,t)), t > 0, x R d, d 2, v(x,0) = u 0 (x), t = 0, x R d. v is a stationary solution of (2). d = 2: Spectral stability implies nonlinear stability. [BL] W.-J. Beyn, J. Lorenz Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

9 Example Consider the quintic complex Ginzburg-Landau equation (QCGL): ( u t = α u +u µ+β u 2 +γ u 4), u = u(x,t) C with u : R d [0, [ C, d {2,3}. For the parameters α = i, β = 5 2 +i, γ = i, µ = 1 2 this equation exhibits so called spinning soliton solutions. Applications: superconductivity, superfluidity, nonlinear optical systems. Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

10 Example Consider the quintic complex Ginzburg-Landau equation (QCGL): ( u t = α u +u µ+β u 2 +γ u 4), u = u(x,t) C with u : R d [0, [ C, d {2,3}. For the parameters α = i, β = 5 2 +i, γ = i, µ = 1 2 this equation exhibits so called spinning soliton solutions. [CMM] L.-C. Crasovan, B.A. Malomed, D. Mihalache Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

11 Example Consider the quintic complex Ginzburg-Landau equation (QCGL): ( u t = α u +u µ+β u 2 +γ u 4), u = u(x,t) C with u : R d [0, [ C, d {2,3}. For the parameters α = i, β = 5 2 +i, γ = i, µ = 1 2 this equation exhibits so called spinning soliton solutions. Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

12 Main result: Exponential decay of v Theorem: (Exponential Decay of v ) Let f(v ) = 0 and Reσ(Df(v )) < 0. Under further assumptions holds: For every 1 < p <, 0 < ϑ < 1 and for every radially nondecreasing weight function θ C(R d,r) of exponential growth rate η 0 with 0 η 2 ϑ 2 3 a 0 b 0 a 2 max p2 there exists K 1 = K 1 (A,f,v,d,p,θ,ϑ) > 0 with the following property: Every classical solution v of A v(x)+ Sx, v(x) +f(v(x)) = 0, x R d, such that v v L p (R d,r N ) and sup x R 0 v (x) v K 1 for some R 0 > 0 satisfies v v W 1,p θ (R d,r N ) (weighted Sobolev space). Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

13 Main result: The assumptions Reλ > 0 and Imλ p 2 2 p 1Reλ λ σ(a) Imλ 2 < p < 1 < p < 2 ϕ = arctan ( ) 2 p 1 p 2 Reλ A,Df(v ) R N,N simultaneously diagonalizable over C a 0 Reλ, λ a max λ σ(a) Reµ b 0 < 0 µ σ(df(v )) Σ p Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

14 Main result: The assumptions A positive function θ C(R d,r) is called a weight function of exponential growth rate η 0 provided that C θ > 0 : θ(x +y) C θ θ(x)e η y x,y R d. [ZM] S. Zelik, A. Mielke Examples: µ R, x R d ( θ 1 (x) = exp( µ x ), θ 3 (x) = exp µ ) x 2 +1, ( θ 2 (x) = cosh(µ x ), θ 4 (x) = cosh µ ) x Exponentially weighted Sobolev spaces: 1 p, k N 0 L p θ (Rd,R N ) := { v L 1 loc (Rd,R N ) θv L p < }, W k,p θ (R d,r N ) := { v L p θ (Rd,R N ) D β u L p θ (Rd,R N ) β k }. Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

15 Outline of proof: Exponential Decay of v Consider the nonlinear problem A v (x)+ Sx, v (x) +f(v (x)) = 0, x R d, d Far-Field Linearization: f C 1, Taylor s theorem, f(v ) = 0 a(x) = 1 0 Df(v +t(v (x) v ))dt, w(x) := v (x) v A w(x)+ Sx, w(x) +a(x)w(x) = 0, x R d. Q(x) Q c(x) K 1 Q ε(x) R 0 x = R Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

16 Outline of proof: Exponential Decay of v Consider the nonlinear problem A v (x)+ Sx, v (x) +f(v (x)) = 0, x R d, d Decomposition of a: Df(v )+Q(x) = 1 0 Df(v +t(v (x) v ))dt, w(x) := v (x) v A w(x)+ Sx, w(x) + (Df(v )+Q(x))w(x) = 0, x R d. Q(x) Q c(x) K 1 Q ε(x) R 0 x = R Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

17 Outline of proof: Exponential Decay of v Consider the nonlinear problem A v (x)+ Sx, v (x) +f(v (x)) = 0, x R d, d Decomposition of a: Df(v )+Q(x) = 1 0 Df(v +t(v (x) v ))dt, w(x) := v (x) v A w(x)+ Sx, w(x) + (Df(v )+Q ε (x)+q c (x))w(x) = 0, x R d. Q(x) 3. Decomposition of Q: Q c(x) Q(x) = Q ε (x)+q c (x), K 1 Q ε(x) R 0 x = R Q,Q ε,q c L (R d,r N,N ), Q ε small, i.e. Q ε L < K 1, Q c compactly supported. Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

18 Outline of proof: Exponential Decay of v Consider the nonlinear problem A v (x)+ Sx, v (x) +f(v (x)) = 0, x R d, d Decomposition of a: Df(v )+Q(x) = 1 0 Df(v +t(v (x) v ))dt, w(x) := v (x) v A w(x)+ Sx, w(x) + (Df(v )+Q ε (x)+q c (x))w(x) = 0, x R d. Q(x) 3. Decomposition of Q: Q c(x) Q(x) = Q ε (x)+q c (x), K 1 Q ε(x) R 0 x = R Q,Q ε,q c L (R d,r N,N ), Q ε small, i.e. Q ε L < K 1, Q c compactly supported. Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

19 Outline of proof: Exponential Decay of v Consider the nonlinear problem A v (x)+ Sx, v (x) +f(v (x)) = 0, x R d, d Decomposition of a: Df(v )+Q(x) = 1 0 Df(v +t(v (x) v ))dt, w(x) := v (x) v A w(x)+ Sx, w(x) + (Df(v )+Q ε (x)+q c (x))w(x) = 0, x R d. Q(x) 3. Decomposition of Q: Q c(x) Q(x) = Q ε (x)+q c (x), K 1 Q ε(x) R 0 x = R Q,Q ε,q c L (R d,r N,N ), Q ε small, i.e. Q ε L < K 1, Q c compactly supported. Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

20 Exponential Decay: To show exponential decay for the solution v of A v(x)+ Sx, v(x) +f(v(x)) = 0, x R d, d 2 investigate the far-field linearization (w.l.o.g. v = 0) A v(x)+ Sx, v(x) +(Df(v )+Q ε (x)+q c (x))v(x) = 0, x R d, d 2. Operators: Study the following operators L Q v :=A v + S, v +Df(v )v +Q ε v +Q c v, L Qε v :=A v + S, v +Df(v )v +Q ε v, L v :=A v + S, v +Df(v )v, L 0 v :=A v + S, v (Ornstein-Uhlenbeck operator). (max. domain) Ornstein-Uhlenbeck Operator Let P,B R d,d, P = P T, P > 0 and B 0. d d d d T P v(x)+ Bx, v(x) = D i (P ij D j v(x))+ D i v(x)b ij x j, x R d Here: P = I d and B = S. i=1 j=1 i=1 j=1 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

21 Exponential Decay: To show exponential decay for the solution v of A v(x)+ Sx, v(x) +f(v(x)) = 0, x R d, d 2 investigate the far-field linearization (w.l.o.g. v = 0) A v(x)+ Sx, v(x) +(Df(v )+Q ε (x)+q c (x))v(x) = 0, x R d, d 2. Operators: Study the following operators L Q v :=A v + S, v +Df(v )v +Q ε v +Q c v, L Qε v :=A v + S, v +Df(v )v +Q ε v, L v :=A v + S, v +Df(v )v, L 0 v :=A v + S, v (Ornstein-Uhlenbeck operator). (max. domain) Ornstein-Uhlenbeck Operator Let P,B R d,d, P = P T, P > 0 and B 0. d d d d T P v(x)+ Bx, v(x) = D i (P ij D j v(x))+ D i v(x)b ij x j, x R d Here: P = I d and B = S. i=1 j=1 i=1 j=1 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

22 Exponential Decay: To show exponential decay for the solution v of A v(x)+ Sx, v(x) +f(v(x)) = 0, x R d, d 2 investigate the far-field linearization (w.l.o.g. v = 0) A v(x)+ Sx, v(x) +(Df(v )+Q ε (x)+q c (x))v(x) = 0, x R d, d 2. Operators: Study the following operators L Q v :=A v + S, v +Df(v )v +Q ε v +Q c v, L Qε v :=A v + S, v +Df(v )v +Q ε v, L v :=A v + S, v +Df(v )v, (exp. decay) L 0 v :=A v + S, v (Ornstein-Uhlenbeck operator). (max. domain) Ornstein-Uhlenbeck Operator Let P,B R d,d, P = P T, P > 0 and B 0. d d d d T P v(x)+ Bx, v(x) = D i (P ij D j v(x))+ D i v(x)b ij x j, x R d Here: P = I d and B = S. i=1 j=1 i=1 j=1 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

23 Exponential Decay: To show exponential decay for the solution v of A v(x)+ Sx, v(x) +f(v(x)) = 0, x R d, d 2 investigate the far-field linearization (w.l.o.g. v = 0) A v(x)+ Sx, v(x) +(Df(v )+Q ε (x)+q c (x))v(x) = 0, x R d, d 2. Operators: Study the following operators L Q v :=A v + S, v +Df(v )v +Q ε v +Q c v, L Qε v :=A v + S, v +Df(v )v +Q ε v, (exp. decay) L v :=A v + S, v +Df(v )v, (exp. decay) L 0 v :=A v + S, v (Ornstein-Uhlenbeck operator). (max. domain) Ornstein-Uhlenbeck Operator Let P,B R d,d, P = P T, P > 0 and B 0. d d d d T P v(x)+ Bx, v(x) = D i (P ij D j v(x))+ D i v(x)b ij x j, x R d Here: P = I d and B = S. i=1 j=1 i=1 j=1 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

24 Exponential Decay: To show exponential decay for the solution v of A v(x)+ Sx, v(x) +f(v(x)) = 0, x R d, d 2 investigate the far-field linearization (w.l.o.g. v = 0) A v(x)+ Sx, v(x) +(Df(v )+Q ε (x)+q c (x))v(x) = 0, x R d, d 2. Operators: Study the following operators L Q v :=A v + S, v +Df(v )v +Q ε v +Q c v, (exp. decay) L Qε v :=A v + S, v +Df(v )v +Q ε v, (exp. decay) L v :=A v + S, v +Df(v )v, (exp. decay) L 0 v :=A v + S, v (Ornstein-Uhlenbeck operator). (max. domain) Ornstein-Uhlenbeck Operator Let P,B R d,d, P = P T, P > 0 and B 0. d d d d T P v(x)+ Bx, v(x) = D i (P ij D j v(x))+ D i v(x)b ij x j, x R d Here: P = I d and B = S. i=1 j=1 i=1 j=1 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

25 Exponential Decay: To show exponential decay for the solution v of A v(x)+ Sx, v(x) +f(v(x)) = 0, x R d, d 2 investigate the far-field linearization (w.l.o.g. v = 0) A v(x)+ Sx, v(x) +(Df(v )+Q ε (x)+q c (x))v(x) = 0, x R d, d 2. Operators: Study the following operators L Q v :=A v + S, v +Df(v )v +Q ε v +Q c v, (exp. decay) L Qε v :=A v + S, v +Df(v )v +Q ε v, (exp. decay) L v :=A v + S, v +Df(v )v, (exp. decay) L 0 v :=A v + S, v (Ornstein-Uhlenbeck operator). (max. domain) Ornstein-Uhlenbeck Operator Let P,B R d,d, P = P T, P > 0 and B 0. d d d d T P v(x)+ Bx, v(x) = D i (P ij D j v(x))+ D i v(x)b ij x j, x R d Here: P = I d and B = S. i=1 j=1 i=1 j=1 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

26 Exponential Decay: To show exponential decay for the solution v of A v(x)+ Sx, v(x) +f(v(x)) = 0, x R d, d 2 investigate the far-field linearization (w.l.o.g. v = 0) A v(x)+ Sx, v(x) +(Df(v )+Q ε (x)+q c (x))v(x) = 0, x R d, d 2. Operators: Study the following operators L Q v :=A v + S, v +Df(v )v +Q ε v +Q c v, (exp. decay) L Qε v :=A v + S, v +Df(v )v +Q ε v, (exp. decay) L v :=A v + S, v +Df(v )v, (exp. decay) L 0 v :=A v + S, v (Ornstein-Uhlenbeck operator). (max. domain) [MPV] G. Metafune, D. Pallara, V. Vespri [MPRS] G. Metafune Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

27 The operator L 0 Ornstein-Uhlenbeck operator [L 0v](x) = A v(x)+ Sx, v(x), x R d, d 2. ( Heat kernel H 0(x,ξ,t) = (4πtA) d 2 exp (4tA) 1 e ts x ξ 2), x,ξ R d, t > 0. semigroup theory ւ unique solv. of resolvent equ., 1 p < (λi A p)v = g L p. Semigroup in L p (R d,r N ), 1 p [T 0(t)v](x) = H 0(x,ξ,t)v(ξ)dξ, t > 0. R d strong continuity Infinitesimal generator (A p,d(a p)), 1 p <. A-priori estimates exponential decay, 1 p < v W 1,p ց L 0: L p -resolvent est. max. domain and max. realization, 1 < p < A p = L 0 on D(A p) = D p (L 0). θ. Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

28 Spectrum of localized rotating waves Linearization about the profile of the rotating wave [Lv](x) = A v(x)+ Sx, v(x) +Df(v (x))v(x), x R d, d 2. Eigenvalue problem [Lv](x) = λv(x), x R d, d 2, λ C. A rotating wave solution u (x,t) = v ( e ts x ) is spectrally stable if Decompose the spectrum σ(l) into σ(l) {λ C Reλ 0}. σ(l) = σ ess (L) σ pt (L), with essential spectrum σ ess (L) and point spectrum σ pt (L). Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

29 Illustration: Point spectrum of L λ (σ(s) {λ+µ λ,µ σ(s), λ µ}) σ pt (L) with algebraic multiplicity Imλ Imλ 1 Imλ i(σ 1 +σ 2 ) 1 Imλ i(σ 1 +σ 2 ) 1 iσ 1 2 iσ 1 1 i(σ 1 σ 2 ) 1 i(σ 1 σ 2 ) 1 iσ 1 2 iσ 1 1 iσ 2 2 iσ Reλ iσ Reλ iσ Reλ iσ Reλ iσ 2 1 i(σ 1 σ 2 ) 1 i(σ 1 σ 2 ) 1 iσ 1 2 iσ 1 1 i(σ 1 +σ 2 ) 1 i(σ 1 +σ 2 ) d = 2 dimse(2) = 3 d = 3 dimse(3) = 6 d = 4 dimse(4) = 10 d = 5 dimse(5) = 15 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

30 Point spectrum of L Theorem: (Point spectrum of L) Let A R N,N satisfy the main assumptions, f C 2 (R N,R N ), 1 < p < and let v be a classical solution of [L 0 v](x)+f(v(x)) = 0. Then v(x) = C rot x +I d C tra, v (x), x R d, C rot so(d), C tra R d solves Lv = λv, whenever (λ,(c rot,c tra )) solves In particular: λc rot = SC rot + ( SC rot) T, λc tra = SC tra. σ(s) {λ+µ λ,µ σ(s), λ µ} σ pt (L), v(x) = Sx, v (x) eigenfunction of λ = 0 for every d 2, group action σ pt (L), Theorem also valid for spiral waves, scroll waves, scroll rings. Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

31 Point spectrum of L Theorem: (Point spectrum of L) Let A R N,N satisfy the main assumptions, f C 2 (R N,R N ), 1 < p < and let v be a classical solution of [L 0 v](x)+f(v(x)) = 0. Then v(x) = C rot x +I d C tra, v (x), x R d, C rot so(d), C tra R d solves Lv = λv, whenever (λ,(c rot,c tra )) solves In particular: λc rot = SC rot + ( SC rot) T, λc tra = SC tra. σ(s) {λ+µ λ,µ σ(s), λ µ} σ pt (L), v(x) = Sx, v (x) eigenfunction of λ = 0 for every d 2, group action σ pt (L), Theorem also valid for spiral waves, scroll waves, scroll rings. Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

32 Exponential decay of v Theorem: (Exponential decay of v) Let the assumptions of the main result be satisfied. The every classical solution v L p (R d,c N ) of A v(x)+ Sx, v(x) +Df(v (x))v(x) = λv(x), x R d, with λ C and Reλ b0 3 satisfies v W 1,p θ (R d,c N ) v exp. localized v exp. localized (with same rate) Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

33 Outline of proof: Point spectrum of L (3) 0 = A v (x)+ Sx, v (x) +f(v (x)), x R d, d Group action: Apply a(r,τ) to (3) 0 = a(r,τ)(a v (x)+ Sx, v (x) +f(v (x))) 2. Derivative d d(r,τ) at (R,τ) = (I d,0) leads to d(d+1) 2 equations 0 =(x j D i x i D j )(A v (x)+ Sx, v (x) +f(v (x))) 0 =D l (A v (x)+ Sx, v (x) +f(v (x))) for i = 1,...,d 1, j = i +1,...,d, l = 1,...,d. 3. Commutator relations for differential operators yield, D (ij) := x j D i x i D j ( ) 0 =L D (ij) v (x) + 0 =L(D l v (x)) d S in D (jn) v (x) n=1 n j d S ln D n v (x) n=1 d S jn D (in) v (x) n=1 n i Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

34 Outline of proof: Point spectrum of L (3) 0 = A v (x)+ Sx, v (x) +f(v (x)), x R d, d Group action: Apply a(r,τ) to (3) 0 = a(r,τ)(a v (x)+ Sx, v (x) +f(v (x))) 2. Derivative d d(r,τ) at (R,τ) = (I d,0) leads to d(d+1) 2 equations 0 =(x j D i x i D j )(A v (x)+ Sx, v (x) +f(v (x))) 0 =D l (A v (x)+ Sx, v (x) +f(v (x))) for i = 1,...,d 1, j = i +1,...,d, l = 1,...,d. 3. Commutator relations for differential operators yield, D (ij) := x j D i x i D j ( ) 0 =L D (ij) v (x) + 0 =L(D l v (x)) d S in D (jn) v (x) n=1 n j d S ln D n v (x) n=1 d S jn D (in) v (x) n=1 n i Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

35 Outline of proof: Point spectrum of L (3) 0 = A v (x)+ Sx, v (x) +f(v (x)), x R d, d Group action: Apply a(r,τ) to (3) 0 = a(r,τ)(a v (x)+ Sx, v (x) +f(v (x))) 2. Derivative d d(r,τ) at (R,τ) = (I d,0) leads to d(d+1) 2 equations 0 =(x j D i x i D j )(A v (x)+ Sx, v (x) +f(v (x))) 0 =D l (A v (x)+ Sx, v (x) +f(v (x))) for i = 1,...,d 1, j = i +1,...,d, l = 1,...,d. 3. Commutator relations for differential operators yield, D (ij) := x j D i x i D j ( ) 0 =L D (ij) v (x) + 0 =L(D l v (x)) d S in D (jn) v (x) n=1 n j d S ln D n v (x) n=1 d S jn D (in) v (x) n=1 n i Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

36 Outline of proof: Point spectrum of L 3. Commutator relations for differential operators yield, D (ij) := x j D i x i D j ( ) 0 =L D (ij) v (x) + 0 =L(D l v (x)) d S in D (jn) v (x) n=1 n j d S ln D n v (x) n=1 d S jn D (in) v (x) n=1 n i 4. Finite-dimensional eigenvalue problem: The ansatz v(x) = d 1 d i=1 j=i+1 reduces Lv = λv to a Note: S is unitary diagonalizable. C rot ij (x j D i x i D j )v (x)+ d l=1 λc rot = SC rot + ( SC rot) T, λc tra = SC tra. Cl tra D l v (x), Cij rot,cl tra C Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

37 Outline of proof: Point spectrum of L 3. Commutator relations for differential operators yield, D (ij) := x j D i x i D j ( ) 0 =L D (ij) v (x) + 0 =L(D l v (x)) d S in D (jn) v (x) n=1 n j d S ln D n v (x) n=1 d S jn D (in) v (x) n=1 n i 4. Finite-dimensional eigenvalue problem: The ansatz v(x) = d 1 d i=1 j=i+1 reduces Lv = λv to a Note: S is unitary diagonalizable. C rot ij (x j D i x i D j )v (x)+ d l=1 λc rot = SC rot + ( SC rot) T, λc tra = SC tra. Cl tra D l v (x), Cij rot,cl tra C Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

38 Illustration: Essential spectrum of L ( λ(ω) + i k X l=1 2 nl σl λ(ω) eigenvalue of ω A Df (v ) d = 2 or 3 Denny Otten (Bielefeld University) d = 4 (not dense) Spatial decay of rotating waves in parabolic systems ) σess (L) d = 4 (dense) Bielefeld / 22

39 Essential spectrum of L Theorem: (Essential spectrum of v) Let the assumptions of the main result be satisfied. Moreover, let ±iσ 1,...,±iσ k denote the nonzero eigenvalues of S and let λ(ω) denote an eigenvalue of ω 2 A Df(v ) for some ω R. Then { in L p (R d,c N ). λ = λ(ω) i } k n l σ l C n l Z, ω R σ ess (L) l=1 Far-field linearization σ ess (L) Dispersion relation: λ σ ess (L) if ( ) k det λi N +ω 2 A+i n l σ l I N Df(v ) = 0 for some κ R. l=1 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

40 Outline of proof: Essential spectrum of L Linearization at the profile v : [Lv](x) = A v(x)+ Sx, v(x) +Df(v )v(x)+q(x)v(x) Q(x) := Df(v (x)) Df(v ), sup Q(x) 2 0 as R x R 1. Orthogonal transformation: S R d,d, S T = S, S = PΛ S block PT. T 1 (x) = Px yields with [L 1 v](x) = A v(x)+ Λ S block x, v(x) +Df(v )v(x)+q(t 1 (x))v(x) Λ S block x, v(x) = k σ l (x 2l D 2l 1 x 2l 1 D 2l )v(x). l=1 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

41 Outline of proof: Essential spectrum of L Orthogonal transformation: [L 1 v](x) = A v(x)+ Λ S block x, v(x) +Df(v )v(x)+q(t 1 (x))v(x) Λ S block x, v(x) = k σ l (x 2l D 2l 1 x 2l 1 D 2l )v(x) l=1 2. Several planar polar coordinates: Transformation ( ) ( ) x2l 1 rl cosφ = T(r x l,φ l ) := l, l = 1,...,k, φ 2l r l sinφ l ] π,π], r l > 0. l yields for ξ = (r 1,φ 1,...,r k,φ k,x 2k+1,...,x d ) with total transformation T 2 (ξ), Q(ξ) := Q(T 1 (T 2 (ξ))) [L 2 v](x) =A [ k l=1 ( 2 r l + 1 r l rl + 1 rl 2 2 φ l )+ d l=2k+1 2 x l ]v(ξ) k σ l φl v(ξ)+df(v )v(ξ)+q(ξ)v(ξ), l=1 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

42 Outline of proof: Essential spectrum of L Several planar polar coordinates: [L 2 v](ξ) =A [ k l=1 ( 2 r l + 1 r l rl + 1 rl 2 2 φ l )+ d l=2k+1 2 x l ]v(ξ) k σ l φl v(ξ)+df(v )v(ξ)+q(ξ)v(ξ), l=1 ξ = (r 1,φ 1,...,r k,φ k,x 2k+1,...,x d ), Q(ξ) := Q(T 1 (T 2 (ξ))) 3. Simplified operator (far-field linearization): Neglecting O ( 1 r) terms yields [ L sim 2 v ] [ k (x) = A r 2 l + l=1 d l=2k+1 2 x l ]v(ξ) k σ l φl v(ξ)+df(v )v(ξ). l=1 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

43 Outline of proof: Essential spectrum of L Simplified operator (far-field linearization): [ L sim 2 v ] [ k (ξ) = A r 2 l + l=1 d l=2k+1 2 x l ]v(ξ) k σ l φl v(ξ)+df(v )v(ξ) l=1 4. Angular Fourier decomposition: ( ( k k v(ξ) = exp iω r l )exp i n l φ l )ˆv,n l Z, ω R, ˆv C N, ˆv = 1 l=1 l=1 φ l ] π,π], r l > 0, l = 1,...,k, yields [( λi L sim 2 ) v ] (ξ) = (λi N +ω 2 A+i ) k n l σ l I N Df(v ) v(ξ). l=1 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

44 Outline of proof: Essential spectrum of L Angular Fourier decomposition: [( λi L sim 2 ) v ] (ξ) = (λi N +κ 2 A+i ) k n l σ l I N Df(v ) v(ξ). l=1 n l Z, κ R, ±iσ l nonzero eigenvalues of S R d,d 5. Finite-dimensional eigenvalue problem: [( ) ] λi L sim 2 v (ξ) = 0 for every ξ if λ C satisfies ( ω 2 A Df(v ) )ˆv ( k = λ+i n l σ l )ˆv, for some ω R. l=1 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

45 Example Consider the quintic complex Ginzburg-Landau equation (QCGL): ( u t = α u +u µ+β u 2 +γ u 4), u = u(x,t) C with u : R d [0, [ C, d {2,3}. For the parameters α = i, β = 5 2 +i, γ = i, µ = 1 2 this equation exhibits so called spinning soliton solutions σess approx σ approx pt σ ess(l) σ(s) {0} Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

46 Figure: Eigenfunctions of QCGL for a spinning soliton with d = 3 Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

47 Work in progress exponential decay in space of continuous functions rotating waves in bounded domains approximation theorem for rotating waves asymptotic boundary conditions numerical computations (interaction of multisolitons) Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

48 Work in progress exponential decay in space of continuous functions rotating waves in bounded domains approximation theorem for rotating waves asymptotic boundary conditions numerical computations (interaction of multisolitons) Denny Otten (Bielefeld University) Spatial decay of rotating waves in parabolic systems Bielefeld / 22

Spatial Decay of Rotating Waves and Restrictions on Finite Disks.

Spatial Decay of Rotating Waves and Restrictions on Finite Disks. University of New Mexico UNM Digital Repository Mathematics & Statistics ETDs Electronic Theses and Dissertations 9-1-2015 Spatial Decay of Rotating Waves and Restrictions on Finite Disks. Sahitya Konda

More information

Spatial Decay of Rotating Waves in Reaction Diffusion Systems

Spatial Decay of Rotating Waves in Reaction Diffusion Systems Satial Decay of Rotating Waves in Reaction Diffusion Systems Wolf-Jürgen Beyn 1 Deartment of Mathematics Bielefeld University 33501 Bielefeld Germany Denny Otten Deartment of Mathematics Bielefeld University

More information

LINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM

LINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM LINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM Unless otherwise stated, all vector spaces in this worksheet are finite dimensional and the scalar field F is R or C. Definition 1. A linear operator

More information

Elliptic Operators with Unbounded Coefficients

Elliptic 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 information

Spectral stability of periodic waves in dispersive models

Spectral stability of periodic waves in dispersive models in dispersive models Collaborators: Th. Gallay, E. Lombardi T. Kapitula, A. Scheel in dispersive models One-dimensional nonlinear waves Standing and travelling waves u(x ct) with c = 0 or c 0 periodic

More information

Nonlinear Stability of Rotating Patterns

Nonlinear Stability of Rotating Patterns Dynamics of PDE, Vol.5, No.4, 349-4, 28 Nonlinear Stability of Rotating Patterns Wolf Jürgen Beyn and Jens Lorenz Communicated by Yannis Kevrekidis, received August 21, 28. Abstract. We consider 2D localized

More information

Analyticity of semigroups generated by Fleming-Viot type operators

Analyticity of semigroups generated by Fleming-Viot type operators Analyticity of semigroups generated by Fleming-Viot type operators Elisabetta Mangino, in collaboration with A. Albanese Università del Salento, Lecce, Italy s Au(x) = x i (δ ij x j )D ij u + b i (x)d

More information

1. Find the solution of the following uncontrolled linear system. 2 α 1 1

1. Find the solution of the following uncontrolled linear system. 2 α 1 1 Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +

More information

Eigenvalues and Eigenvectors

Eigenvalues 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 information

Global 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 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 information

MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.

MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. Adjoint operator and adjoint matrix Given a linear operator L on an inner product space V, the adjoint of L is a transformation

More information

Evans function review

Evans function review Evans function review Part I: History, construction, properties and applications Strathclyde, April 18th 2005 Simon J.A. Malham http://www.ma.hw.ac.uk/ simonm/talks/ Acknowledgments T.J. Bridges, C.K.R.T.

More information

Nonlinear convective stability of travelling fronts near Turing and Hopf instabilities

Nonlinear convective stability of travelling fronts near Turing and Hopf instabilities Nonlinear convective stability of travelling fronts near Turing and Hopf instabilities Margaret Beck Joint work with Anna Ghazaryan, University of Kansas and Björn Sandstede, Brown University September

More information

Nonlinear Stability of Rotating Patterns

Nonlinear Stability of Rotating Patterns Nonlinear Stability of Rotating Patterns Wolf Jürgen Beyn Fakultät für Mathematik, Postfach 1131 Universität Bielefeld, 3351 Bielefeld Jens Lorenz Department of Mathematics and Statistics UNM, Albuquerque,

More information

Introduction to Matrix Algebra

Introduction to Matrix Algebra Introduction to Matrix Algebra August 18, 2010 1 Vectors 1.1 Notations A p-dimensional vector is p numbers put together. Written as x 1 x =. x p. When p = 1, this represents a point in the line. When p

More information

Lecture Notes on PDEs

Lecture 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 information

Partial Differential Equations

Partial Differential Equations M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1

More information

Intertwinings for Markov processes

Intertwinings 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 information

Problems in Linear Algebra and Representation Theory

Problems in Linear Algebra and Representation Theory Problems in Linear Algebra and Representation Theory (Most of these were provided by Victor Ginzburg) The problems appearing below have varying level of difficulty. They are not listed in any specific

More information

Speeding up Convergence to Equilibrium for Diffusion Processes

Speeding 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 information

Sturm-Liouville operators have form (given p(x) > 0, q(x)) + q(x), (notation means Lf = (pf ) + qf ) dx

Sturm-Liouville operators have form (given p(x) > 0, q(x)) + q(x), (notation means Lf = (pf ) + qf ) dx Sturm-Liouville operators Sturm-Liouville operators have form (given p(x) > 0, q(x)) L = d dx ( p(x) d ) + q(x), (notation means Lf = (pf ) + qf ) dx Sturm-Liouville operators Sturm-Liouville operators

More information

REPRESENTATION THEORY WEEK 7

REPRESENTATION THEORY WEEK 7 REPRESENTATION THEORY WEEK 7 1. Characters of L k and S n A character of an irreducible representation of L k is a polynomial function constant on every conjugacy class. Since the set of diagonalizable

More information

Definite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra

Definite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra Definite versus Indefinite Linear Algebra Christian Mehl Institut für Mathematik TU Berlin Germany 10th SIAM Conference on Applied Linear Algebra Monterey Bay Seaside, October 26-29, 2009 Indefinite Linear

More information

Oblique derivative problems for elliptic and parabolic equations, Lecture II

Oblique derivative problems for elliptic and parabolic equations, Lecture II of the for elliptic and parabolic equations, Lecture II Iowa State University July 22, 2011 of the 1 2 of the of the As a preliminary step in our further we now look at a special situation for elliptic.

More information

Nonlinear stability of time-periodic viscous shocks. Margaret Beck Brown University

Nonlinear stability of time-periodic viscous shocks. Margaret Beck Brown University Nonlinear stability of time-periodic viscous shocks Margaret Beck Brown University Motivation Time-periodic patterns in reaction-diffusion systems: t x Experiment: chemical reaction chlorite-iodite-malonic-acid

More information

Stability 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 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 information

Exponential Stability of the Traveling Fronts for a Pseudo-Para. Pseudo-Parabolic Fisher-KPP Equation

Exponential Stability of the Traveling Fronts for a Pseudo-Para. Pseudo-Parabolic Fisher-KPP Equation Exponential Stability of the Traveling Fronts for a Pseudo-Parabolic Fisher-KPP Equation Based on joint work with Xueli Bai (Center for PDE, East China Normal Univ.) and Yang Cao (Dalian University of

More information

Gene Regulatory Network and Functional Analysis: Hopf Bifurcation causes by molecular movement

Gene Regulatory Network and Functional Analysis: Hopf Bifurcation causes by molecular movement Gene Regulatory Network and Functional Analysis: Hopf Bifurcation causes by molecular movement Mariya Ptashnyk Division of Mathematics, University of Dundee joined work with Mark Chaplain and Marc Sturrock

More information

Math Matrix Algebra

Math Matrix Algebra Math 44 - Matrix Algebra Review notes - (Alberto Bressan, Spring 7) sec: Orthogonal diagonalization of symmetric matrices When we seek to diagonalize a general n n matrix A, two difficulties may arise:

More information

MATH 304 Linear Algebra Lecture 33: Bases of eigenvectors. Diagonalization.

MATH 304 Linear Algebra Lecture 33: Bases of eigenvectors. Diagonalization. MATH 304 Linear Algebra Lecture 33: Bases of eigenvectors. Diagonalization. Eigenvalues and eigenvectors of an operator Definition. Let V be a vector space and L : V V be a linear operator. A number λ

More information

Eddy-diffusivity fields in turbulent transport: explicit expressions

Eddy-diffusivity fields in turbulent transport: explicit expressions DICCA University of Genova Italy Eddy-diffusivity fields in turbulent transport: explicit expressions Andrea Mazzino andrea.mazzino@unige.it Lagrangian transport: from complex flows to complex fluids Lecce

More information

TRAVELLING WAVES. Morteza Fotouhi Sharif Univ. of Technology

TRAVELLING WAVES. Morteza Fotouhi Sharif Univ. of Technology TRAVELLING WAVES Morteza Fotohi Sharif Univ. of Technology Mini Math NeroScience Mini Math NeroScience Agst 28 REACTION DIFFUSION EQUATIONS U = DU + f ( U ) t xx x t > U n D d 1 = d j > d n 2 Travelling

More information

Online Exercises for Linear Algebra XM511

Online Exercises for Linear Algebra XM511 This document lists the online exercises for XM511. The section ( ) numbers refer to the textbook. TYPE I are True/False. Lecture 02 ( 1.1) Online Exercises for Linear Algebra XM511 1) The matrix [3 2

More information

g(t) = f(x 1 (t),..., x n (t)).

g(t) = f(x 1 (t),..., x n (t)). Reading: [Simon] p. 313-333, 833-836. 0.1 The Chain Rule Partial derivatives describe how a function changes in directions parallel to the coordinate axes. Now we shall demonstrate how the partial derivatives

More information

Before you begin read these instructions carefully:

Before you begin read these instructions carefully: NATURAL SCIENCES TRIPOS Part IB & II (General) Tuesday, 30 May, 2017 9:00 am to 12:00 pm MATHEMATICS (1) Before you begin read these instructions carefully: You may submit answers to no more than six questions.

More information

Physics 505 Homework No. 8 Solutions S Spinor rotations. Somewhat based on a problem in Schwabl.

Physics 505 Homework No. 8 Solutions S Spinor rotations. Somewhat based on a problem in Schwabl. Physics 505 Homework No 8 s S8- Spinor rotations Somewhat based on a problem in Schwabl a) Suppose n is a unit vector We are interested in n σ Show that n σ) = I where I is the identity matrix We will

More information

j=1 x j p, if 1 p <, x i ξ : x i < ξ} 0 as p.

j=1 x j p, if 1 p <, x i ξ : x i < ξ} 0 as p. LINEAR ALGEBRA Fall 203 The final exam Almost all of the problems solved Exercise Let (V, ) be a normed vector space. Prove x y x y for all x, y V. Everybody knows how to do this! Exercise 2 If V is a

More information

An Introduction to Stability Theory for Nonlinear PDEs

An Introduction to Stability Theory for Nonlinear PDEs An Introduction to Stability Theory for Nonlinear PDEs Mathew A. Johnson 1 Abstract These notes were prepared for the 217 Participating School in Analysis of PDE: Stability of Solitons and Periodic Waves

More information

ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations

ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.01 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential

More information

November 18, 2013 ANALYTIC FUNCTIONAL CALCULUS

November 18, 2013 ANALYTIC FUNCTIONAL CALCULUS November 8, 203 ANALYTIC FUNCTIONAL CALCULUS RODICA D. COSTIN Contents. The spectral projection theorem. Functional calculus 2.. The spectral projection theorem for self-adjoint matrices 2.2. The spectral

More information

Math Bootcamp An p-dimensional vector is p numbers put together. Written as. x 1 x =. x p

Math Bootcamp An p-dimensional vector is p numbers put together. Written as. x 1 x =. x p Math Bootcamp 2012 1 Review of matrix algebra 1.1 Vectors and rules of operations An p-dimensional vector is p numbers put together. Written as x 1 x =. x p. When p = 1, this represents a point in the

More information

MATH 583A REVIEW SESSION #1

MATH 583A REVIEW SESSION #1 MATH 583A REVIEW SESSION #1 BOJAN DURICKOVIC 1. Vector Spaces Very quick review of the basic linear algebra concepts (see any linear algebra textbook): (finite dimensional) vector space (or linear space),

More information

Matrix Theory. A.Holst, V.Ufnarovski

Matrix Theory. A.Holst, V.Ufnarovski Matrix Theory AHolst, VUfnarovski 55 HINTS AND ANSWERS 9 55 Hints and answers There are two different approaches In the first one write A as a block of rows and note that in B = E ij A all rows different

More information

Math 108b: Notes on the Spectral Theorem

Math 108b: Notes on the Spectral Theorem Math 108b: Notes on the Spectral Theorem From section 6.3, we know that every linear operator T on a finite dimensional inner product space V has an adjoint. (T is defined as the unique linear operator

More information

MTH 464: Computational Linear Algebra

MTH 464: Computational Linear Algebra MTH 464: Computational Linear Algebra Lecture Outlines Exam 4 Material Prof. M. Beauregard Department of Mathematics & Statistics Stephen F. Austin State University April 15, 2018 Linear Algebra (MTH 464)

More information

Slow Modulation & Large-Time Dynamics Near Periodic Waves

Slow Modulation & Large-Time Dynamics Near Periodic Waves Slow Modulation & Large-Time Dynamics Near Periodic Waves Miguel Rodrigues IRMAR Université Rennes 1 France SIAG-APDE Prize Lecture Jointly with Mathew Johnson (Kansas), Pascal Noble (INSA Toulouse), Kevin

More information

Criterions on periodic feedback stabilization for some evolution equations

Criterions on periodic feedback stabilization for some evolution equations Criterions on periodic feedback stabilization for some evolution equations School of Mathematics and Statistics, Wuhan University, P. R. China (Joint work with Yashan Xu, Fudan University) Toulouse, June,

More information

The following definition is fundamental.

The following definition is fundamental. 1. Some Basics from Linear Algebra With these notes, I will try and clarify certain topics that I only quickly mention in class. First and foremost, I will assume that you are familiar with many basic

More information

Math 46, Applied Math (Spring 2008): Final

Math 46, Applied Math (Spring 2008): Final Math 46, Applied Math (Spring 2008): Final 3 hours, 80 points total, 9 questions, roughly in syllabus order (apart from short answers) 1. [16 points. Note part c, worth 7 points, is independent of the

More information

4. Special Solutions and Stability

4. Special Solutions and Stability 4. Special Solutions and Stability Most PDEs in science and engineering are nonlinear. They cannot be solved either explicitly or semiexplicitly (separation of variables, transform methods, Green s functions,

More information

Compact symetric bilinear forms

Compact symetric bilinear forms Compact symetric bilinear forms Mihai Mathematics Department UC Santa Barbara IWOTA 2006 IWOTA 2006 Compact forms [1] joint work with: J. Danciger (Stanford) S. Garcia (Pomona

More information

Nash Type Inequalities for Fractional Powers of Non-Negative Self-adjoint Operators. ( Wroclaw 2006) P.Maheux (Orléans. France)

Nash Type Inequalities for Fractional Powers of Non-Negative Self-adjoint Operators. ( Wroclaw 2006) P.Maheux (Orléans. France) Nash Type Inequalities for Fractional Powers of Non-Negative Self-adjoint Operators ( Wroclaw 006) P.Maheux (Orléans. France) joint work with A.Bendikov. European Network (HARP) (to appear in T.A.M.S)

More information

1 Math 241A-B Homework Problem List for F2015 and W2016

1 Math 241A-B Homework Problem List for F2015 and W2016 1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let

More information

PART IV Spectral Methods

PART IV Spectral Methods PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,

More information

i=1 α i. Given an m-times continuously

i=1 α i. Given an m-times continuously 1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable

More information

Dynamic programming using radial basis functions and Shepard approximations

Dynamic programming using radial basis functions and Shepard approximations Dynamic programming using radial basis functions and Shepard approximations Oliver Junge, Alex Schreiber Fakultät für Mathematik Technische Universität München Workshop on algorithms for dynamical systems

More information

On Schrödinger equations with inverse-square singular potentials

On Schrödinger equations with inverse-square singular potentials On Schrödinger equations with inverse-square singular potentials Veronica Felli Dipartimento di Statistica University of Milano Bicocca veronica.felli@unimib.it joint work with Elsa M. Marchini and Susanna

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

Introduction to Semigroup Theory

Introduction to Semigroup Theory Introduction to Semigroup Theory Franz X. Gmeineder LMU München, U Firenze Bruck am Ziller / Dec 15th 2012 Franz X. Gmeineder Introduction to Semigroup Theory 1/25 The Way Up: Opening The prototype of

More information

P(E t, Ω)dt, (2) 4t has an advantage with respect. to the compactly supported mollifiers, i.e., the function W (t)f satisfies a semigroup law:

P(E t, Ω)dt, (2) 4t has an advantage with respect. to the compactly supported mollifiers, i.e., the function W (t)f satisfies a semigroup law: Introduction Functions of bounded variation, usually denoted by BV, have had and have an important role in several problems of calculus of variations. The main features that make BV functions suitable

More information

Nonlinear Control. Nonlinear Control Lecture # 6 Passivity and Input-Output Stability

Nonlinear Control. Nonlinear Control Lecture # 6 Passivity and Input-Output Stability Nonlinear Control Lecture # 6 Passivity and Input-Output Stability Passivity: Memoryless Functions y y y u u u (a) (b) (c) Passive Passive Not passive y = h(t,u), h [0, ] Vector case: y = h(t,u), h T =

More information

OR MSc Maths Revision Course

OR MSc Maths Revision Course OR MSc Maths Revision Course Tom Byrne School of Mathematics University of Edinburgh t.m.byrne@sms.ed.ac.uk 15 September 2017 General Information Today JCMB Lecture Theatre A, 09:30-12:30 Mathematics revision

More information

Differential equations

Differential equations Differential equations Math 7 Spring Practice problems for April Exam Problem Use the method of elimination to find the x-component of the general solution of x y = 6x 9x + y = x 6y 9y Soln: The system

More information

2 Lyapunov Stability. x(0) x 0 < δ x(t) x 0 < ɛ

2 Lyapunov Stability. x(0) x 0 < δ x(t) x 0 < ɛ 1 2 Lyapunov Stability Whereas I/O stability is concerned with the effect of inputs on outputs, Lyapunov stability deals with unforced systems: ẋ = f(x, t) (1) where x R n, t R +, and f : R n R + R n.

More information

Stability Analysis of Stationary Solutions for the Cahn Hilliard Equation

Stability Analysis of Stationary Solutions for the Cahn Hilliard Equation Stability Analysis of Stationary Solutions for the Cahn Hilliard Equation Peter Howard, Texas A&M University University of Louisville, Oct. 19, 2007 References d = 1: Commun. Math. Phys. 269 (2007) 765

More information

Qualification Exam: Mathematical Methods

Qualification Exam: Mathematical Methods Qualification Exam: Mathematical Methods Name:, QEID#41534189: August, 218 Qualification Exam QEID#41534189 2 1 Mathematical Methods I Problem 1. ID:MM-1-2 Solve the differential equation dy + y = sin

More information

Positive Stabilization of Infinite-Dimensional Linear Systems

Positive Stabilization of Infinite-Dimensional Linear Systems Positive Stabilization of Infinite-Dimensional Linear Systems Joseph Winkin Namur Center of Complex Systems (NaXys) and Department of Mathematics, University of Namur, Belgium Joint work with Bouchra Abouzaid

More information

1 Definition and Basic Properties of Compa Operator

1 Definition and Basic Properties of Compa Operator 1 Definition and Basic Properties of Compa Operator 1.1 Let X be a infinite dimensional Banach space. Show that if A C(X ), A does not have bounded inverse. Proof. Denote the unit ball of X by B and the

More information

Dynamic interpretation of eigenvectors

Dynamic interpretation of eigenvectors EE263 Autumn 2015 S. Boyd and S. Lall Dynamic interpretation of eigenvectors invariant sets complex eigenvectors & invariant planes left eigenvectors modal form discrete-time stability 1 Dynamic interpretation

More information

Stability of Parameter Adaptation Algorithms. Big picture

Stability of Parameter Adaptation Algorithms. Big picture ME5895, UConn, Fall 215 Prof. Xu Chen Big picture For ˆθ (k + 1) = ˆθ (k) + [correction term] we haven t talked about whether ˆθ(k) will converge to the true value θ if k. We haven t even talked about

More information

Asymptotic stability for solitons of the Gross-Pitaevskii and Landau-Lifshitz equations

Asymptotic stability for solitons of the Gross-Pitaevskii and Landau-Lifshitz equations Asymptotic stability for solitons of the Gross-Pitaevskii and Landau-Lifshitz equations Philippe Gravejat Cergy-Pontoise University Joint work with F. Béthuel (Paris 6), A. de Laire (Lille) and D. Smets

More information

Asymptotic behaviour of the heat equation in twisted waveguides

Asymptotic behaviour of the heat equation in twisted waveguides Asymptotic behaviour of the heat equation in twisted waveguides Gabriela Malenová Faculty of Nuclear Sciences and Physical Engineering, CTU, Prague Nuclear Physics Institute, AS ČR, Řež Graphs and Spectra,

More information

PH.D. PRELIMINARY EXAMINATION MATHEMATICS

PH.D. PRELIMINARY EXAMINATION MATHEMATICS UNIVERSITY OF CALIFORNIA, BERKELEY SPRING SEMESTER 207 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem

More information

Symmetry Properties and Exact Solutions of the Fokker-Planck Equation

Symmetry Properties and Exact Solutions of the Fokker-Planck Equation Nonlinear Mathematical Physics 1997, V.4, N 1, 13 136. Symmetry Properties and Exact Solutions of the Fokker-Planck Equation Valery STOHNY Kyïv Polytechnical Institute, 37 Pobedy Avenue, Kyïv, Ukraïna

More information

Linear Algebra. Workbook

Linear Algebra. Workbook Linear Algebra Workbook Paul Yiu Department of Mathematics Florida Atlantic University Last Update: November 21 Student: Fall 2011 Checklist Name: A B C D E F F G H I J 1 2 3 4 5 6 7 8 9 10 xxx xxx xxx

More information

1 Linear transformations; the basics

1 Linear transformations; the basics Linear Algebra Fall 2013 Linear Transformations 1 Linear transformations; the basics Definition 1 Let V, W be vector spaces over the same field F. A linear transformation (also known as linear map, or

More information

Partial Differential Equations

Partial 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 information

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is

More information

BIHARMONIC WAVE MAPS INTO SPHERES

BIHARMONIC WAVE MAPS INTO SPHERES BIHARMONIC WAVE MAPS INTO SPHERES SEBASTIAN HERR, TOBIAS LAMM, AND ROLAND SCHNAUBELT Abstract. A global weak solution of the biharmonic wave map equation in the energy space for spherical targets is constructed.

More information

Mathematical Analysis of a Generalized Chiral Quark Soliton Model

Mathematical Analysis of a Generalized Chiral Quark Soliton Model Symmetry, Integrability and Geometry: Methods and Applications Vol. 2 (2006), Paper 018, 12 pages Mathematical Analysis of a Generalized Chiral Quark Soliton Model Asao ARAI Department of Mathematics,

More information

Iné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. 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 information

Maximum Principles for Parabolic Equations

Maximum Principles for Parabolic Equations Maximum Principles for Parabolic Equations Kamyar Malakpoor 24 November 2004 Textbooks: Friedman, A. Partial Differential Equations of Parabolic Type; Protter, M. H, Weinberger, H. F, Maximum Principles

More information

Linear DifferentiaL Equation

Linear DifferentiaL Equation Linear DifferentiaL Equation Massoud Malek The set F of all complex-valued functions is known to be a vector space of infinite dimension. Solutions to any linear differential equations, form a subspace

More information

SPECTRAL METHODS: ORTHOGONAL POLYNOMIALS

SPECTRAL METHODS: ORTHOGONAL POLYNOMIALS SPECTRAL METHODS: ORTHOGONAL POLYNOMIALS 31 October, 2007 1 INTRODUCTION 2 ORTHOGONAL POLYNOMIALS Properties of Orthogonal Polynomials 3 GAUSS INTEGRATION Gauss- Radau Integration Gauss -Lobatto Integration

More information

Robust control and applications in economic theory

Robust control and applications in economic theory Robust control and applications in economic theory In honour of Professor Emeritus Grigoris Kalogeropoulos on the occasion of his retirement A. N. Yannacopoulos Department of Statistics AUEB 24 May 2013

More information

EE5120 Linear Algebra: Tutorial 6, July-Dec Covers sec 4.2, 5.1, 5.2 of GS

EE5120 Linear Algebra: Tutorial 6, July-Dec Covers sec 4.2, 5.1, 5.2 of GS EE0 Linear Algebra: Tutorial 6, July-Dec 07-8 Covers sec 4.,.,. of GS. State True or False with proper explanation: (a) All vectors are eigenvectors of the Identity matrix. (b) Any matrix can be diagonalized.

More information

Polarization Optics. N. Fressengeas

Polarization Optics. N. Fressengeas Polarization Optics N. Fressengeas Laboratoire Matériaux Optiques, Photonique et Systèmes Unité de Recherche commune à l Université de Lorraine et à Supélec Download this document from http://arche.univ-lorraine.fr/

More information

SPRING 2006 PRELIMINARY EXAMINATION SOLUTIONS

SPRING 2006 PRELIMINARY EXAMINATION SOLUTIONS SPRING 006 PRELIMINARY EXAMINATION SOLUTIONS 1A. Let G be the subgroup of the free abelian group Z 4 consisting of all integer vectors (x, y, z, w) such that x + 3y + 5z + 7w = 0. (a) Determine a linearly

More information

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points Nonlinear Control Lecture # 3 Stability of Equilibrium Points The Invariance Principle Definitions Let x(t) be a solution of ẋ = f(x) A point p is a positive limit point of x(t) if there is a sequence

More information

Maths for Signals and Systems Linear Algebra in Engineering

Maths for Signals and Systems Linear Algebra in Engineering Maths for Signals and Systems Linear Algebra in Engineering Lectures 13-14, Tuesday 11 th November 2014 DR TANIA STATHAKI READER (ASSOCIATE PROFFESOR) IN SIGNAL PROCESSING IMPERIAL COLLEGE LONDON Eigenvectors

More information

Pseudospectra and Nonnormal Dynamical Systems

Pseudospectra and Nonnormal Dynamical Systems Pseudospectra and Nonnormal Dynamical Systems Mark Embree and Russell Carden Computational and Applied Mathematics Rice University Houston, Texas ELGERSBURG MARCH 1 Overview of the Course These lectures

More information

Nonlinear stability of semidiscrete shocks for two-sided schemes

Nonlinear stability of semidiscrete shocks for two-sided schemes Nonlinear stability of semidiscrete shocks for two-sided schemes Margaret Beck Boston University Joint work with Hermen Jan Hupkes, Björn Sandstede, and Kevin Zumbrun Setting: semi-discrete conservation

More information

Stabilization of Distributed Parameter Systems by State Feedback with Positivity Constraints

Stabilization of Distributed Parameter Systems by State Feedback with Positivity Constraints Stabilization of Distributed Parameter Systems by State Feedback with Positivity Constraints Joseph Winkin Namur Center of Complex Systems (naxys) and Dept. of Mathematics, University of Namur, Belgium

More information

w T 1 w T 2. w T n 0 if i j 1 if i = j

w T 1 w T 2. w T n 0 if i j 1 if i = j Lyapunov Operator Let A F n n be given, and define a linear operator L A : C n n C n n as L A (X) := A X + XA Suppose A is diagonalizable (what follows can be generalized even if this is not possible -

More information

Magnetic wells in dimension three

Magnetic wells in dimension three Magnetic wells in dimension three Yuri A. Kordyukov joint with Bernard Helffer & Nicolas Raymond & San Vũ Ngọc Magnetic Fields and Semiclassical Analysis Rennes, May 21, 2015 Yuri A. Kordyukov (Ufa) Magnetic

More information

On some nonlinear parabolic equation involving variable exponents

On some nonlinear parabolic equation involving variable exponents On some nonlinear parabolic equation involving variable exponents Goro Akagi (Kobe University, Japan) Based on a joint work with Giulio Schimperna (Pavia Univ., Italy) Workshop DIMO-2013 Diffuse Interface

More information

Strong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term

Strong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term 1 Strong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term Enrico Priola Torino (Italy) Joint work with G. Da Prato, F. Flandoli and M. Röckner Stochastic Processes

More information

MATH 423 Linear Algebra II Lecture 10: Inverse matrix. Change of coordinates.

MATH 423 Linear Algebra II Lecture 10: Inverse matrix. Change of coordinates. MATH 423 Linear Algebra II Lecture 10: Inverse matrix. Change of coordinates. Let V be a vector space and α = [v 1,...,v n ] be an ordered basis for V. Theorem 1 The coordinate mapping C : V F n given

More information

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Stability 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 information

Interest Rate Models:

Interest Rate Models: 1/17 Interest Rate Models: from Parametric Statistics to Infinite Dimensional Stochastic Analysis René Carmona Bendheim Center for Finance ORFE & PACM, Princeton University email: rcarmna@princeton.edu

More information