RESONANCE-FREE REGIONS FOR DIFFRACTIVE TRAPPING BY CONORMAL POTENTIALS

Size: px
Start display at page:

Download "RESONANCE-FREE REGIONS FOR DIFFRACTIVE TRAPPING BY CONORMAL POTENTIALS"

Transcription

1 RESONANCE-FREE REGIONS FOR DIFFRACTIVE TRAPPING BY CONORMAL POTENTIALS ORAN GANNOT AND JARED WUNSCH Abstract. We consider te Scrödinger operator P = 2 g + V on R n equipped wit a metric g tat is Euclidean outside a compact set. Te real-valued potential V is assumed be compactly supported and smoot except at conormal singularities of order 1 α along a compact ypersurface Y. For α > 2 (or even α > 1 if te classical flow is unique), we sow tat if E 0 is a non-trapping energy for te classical flow, ten te operator P as no resonances in a region [E 0 δ, E 0 + δ] i[0, ν 0 log(1/)]. Te constant ν 0 is explicit in terms of α and dynamical quantities. We also sow tat te size of tis resonance-free region is optimal for te class of piecewise-smoot potentials on te line. 1. Introduction 1.1. Main results. Let X = R n, equipped wit a smoot Riemannian metric g suc tat g ij = δ ij outside a compact set. Wit g denoting te nonnegative Laplacian on (X, g), consider te semiclassical Scrödinger operator wit compactly supported potential, P = 2 g + V. We assume tat V I [ 1 α] (Y ) is conormal to a compact ypersurface Y X wit α > 1. Tis notation means te following: V is C away from Y. In local coordinates (x 1, x ) near Y, wit Y given by {x 1 = 0}, V (x 1, x ) = e ix1 ξ 1 v(x, ξ 1 ) dξ 1, v S 1 α (R n x; R ξ 1). Suc a potential V is at least C 1,γ for some γ > 0. Let p = σ (P ) denote te semiclassical principal symbol of P, p = ξ 2 g + V. Te Hamilton vector field H p is continuous, ence always as global solutions. We furter assume tat H p as unique integral curves; tis is always true if α > 2 (were 1

2 2 ORAN GANNOT AND JARED WUNSCH H p is Lipscitz) but in general fails in te range α (1, 2]. In particular, tere is a well-defined flow ρ exp thp (ρ) on T X, wic is tangent to eac energy surface {p = E}. Let E 0 > 0 be a non-trapping energy level for te H p flow, i.e., assume tat x in bot directions along all integral curves of H p in {p = E 0 }. We sow tat for a suitable ν 0 > 0 and, δ > 0 sufficiently small, tere are no resonances of P in te spectral window [E 0 δ, E 0 + δ] i[0, ν 0 log(1/)]. Te quantity ν 0 > 0 as a dynamical caracterization wic we discuss next. Let H E T Y denote te set of yperbolic points for p E; tese are te points in pase space were te flow is transverse to Y. More precisely, if we introduce normal coordinates (x 1, x ) for g wit respect to Y, so tat Y = {x 1 = 0} locally, ten we can write p(x, ξ) E = (ξ 1 ) 2 + K(x)ξ, ξ + V (x) E (ξ 1 ) 2 r(x, ξ, E) for a positive definite matrix K(x). In tese coordinates, (x, ξ ) H E r(0, x, ξ, E) > 0. We also remark for later use tat te glancing set G E T Y is defined in coordinates by te equation r(0, x, ξ, E) = 0. For x Y, let π : T x X T x Y denote te canonical projection, wic in local coordinates is just te map (0, x, ξ 1, ξ ) (x, ξ ). Note tat π is two-to-one over H E and one-to-one over G E. Given E R, we introduce te affine lengt of te longest H p trajectory connecting two yperbolic points: diam E (Y ) = sup{ t : tere exists ρ π 1 (H E ) wit exp thp (ρ) π 1 (H E )}. (1.1) Since Y is compact, diam E (Y ) is finite for a non-trapping energy E. For I R we also define diam I (Y ) = sup{diam E (Y ) : E I}. Tis is again finite if I = [E 0 δ, E 0 + δ] for E 0 non-trapping and δ > 0 sufficiently small. Teorem 1. Let (X, g) and V I [ α 1] (Y ) be as above. If E 0 > 0 is non-trapping, ten tere exists δ 0 > 0 wit te following property. Given δ (0, δ 0 ) and α 0 < ν 0 < diam [E0 δ,e 0 +δ](y ),

3 tere exists 0 > 0 suc tat P as no resonance z wit DIFFRACTIVE TRAPPING 3 z [E 0 δ, E 0 + δ] i[0, ν 0 log(1/)] for (0, 0 ). If diam [E0 δ,e 0 +δ](y ) = 0, ten te conclusion is valid for any ν 0 (0, ). Wen V is smoot (so tat α can be taken arbitrarily large), te fact tat Teorem 1 olds for arbitrary ν 0 (0, ) is originally due to Martinez [11]. Te size of te resonance-free region in Teorem 1 is already optimal for te class of compactly supported piecewise-smoot potentials on R. Consider te operator P = (D x ) 2 + V, were V satisfies te following properties: (1) Tere exists L > 0 suc tat supp V [0, L]. (2) Te restriction of V to [0, L] is smoot. If V vanises to order k at x = 0 and to order l at x = L, ten V I [ 1 min(k,l)] ({0, L}). For use in Teorem 2 below define te quantity φ = (1/2π) arg(v (k) (0 + ) V (l) (L )), wic arises as te pase sift in a Bor Sommerfeld type formula. Consider a spectral interval [a, b], were a > sup V. Certainly any energy E [a, b] is nontrapping for te H p flow. Define (alf) te action and period by S(E) = L 0 (E V (s)) 1/2 ds, T (E) = L 0 1 ds. (1.2) 2(E V (s)) 1/2 Of course tis is a sligt abuse of terminology, since te endpoints are not turning points for te classical dynamics. Observe tat in te notation of (1.1). T (E) = diam E ({0, L}) Noting tat S(E) is increasing on any interval [a, b] as above, define [α, β] = S([a, b]) and ten set N() = {n Z : π(n + (l k)/4 + φ) 1 [α, β]}. Note tat N() 1 β α /π. We ten ave te following semiclassical analogue of [14, Teorem 6] on te existence and asymptotics of resonances. Teorem 2. Let k and l denote te orders of vanising of V at x = 0 and x = L, respectively. For eac n N() let E n = S 1 (π(n + (l k)/4 + φ)).

4 4 ORAN GANNOT AND JARED WUNSCH Tere exists 0 > 0 suc tat for eac (0, 0 ) and n N() tere is a unique resonance z n satisfying i(l + k) z n = E n 2T (E n ) log(1/) + i 2T (E n ) + O( 2 log(1/) 2 ). ( log V (k) (0 + ) V (l) (L ) (1/2)(l + k + 4) log(4e n ) ) (1.3) Furtermore, if M > 0 is sufficiently large, ten {z n : n N()} are all te resonances in [a, b] + i[ M log(1/), 0]. Because E S(E) = T (E) and T (E) is uniformly positive for E [a, b] it follows tat z n z m C n m for some C > 0 and eac n, m N(). Tus te z n are all distinct, wic gives 1 β α /π as an asymptotic formula for te number of resonances in [a, b]+[ im log(1/), 0] Context and previous work. Wen diffraction of singularities is te only trapping, an increasing body of work suggests tat resonances may occur at Im z C log(1/) for various values of C > 0, but no closer to te real axis. Tis is farter into te lower alf-plane tan te resonances occurring in cases of elliptic trapping (were tey rapidly approac te real axis as 0) or even for yperbolic trapped sets, were tere is an O() resonance-free region (see [13] for references on te subject of classical dynamical trapping and resonances). In te diffractive case, te regularization of te trapped wave wit eac successive diffraction is by contrast responsible for te faster rate of decay as measured by resonance widt. Most of te literature substantiating tis euristic is in te omogeneous rater tan te semiclassical setting; tere te analog of resonances at E 0 iν 0 log(1/) are resonances close to a curve Im λ C log Re λ (1.4) as λ. Tat diffraction of singularities can in fact create strings of resonances along suc log curves was demonstrated in te omogeneous setting by Zworski [14] (see also earlier work by Regge [12]). In te setting of diffraction by analytic corners in te plane, Burq [3] likewise sowed tat resonances lie asymptotically on families of curves (1.4) for various values of C > 0. In te setting of manifolds wit conic singularities, were similar diffractive propagation occurs, a number of recent teorems ave explored te same teme. Baskin Wunsc [1] sowed on te one and tat some nontrivial region of te form Im λ > ν 0 log Re λ, λ > R

5 DIFFRACTIVE TRAPPING 5 contains no resonances (subject to some genericity conditions on te relationsip among te conic singularities) tis is analogous to te gap teorem obtained ere. Galkowski [7] ten found te largest ν 0 wic could be obtained by te Vainberg parametrix metod employed in [1]: ν 0 = (n 1)/(2L 0 ), wit n te dimension and L 0 te maximal distance between cone points. Work of Hillairet Wunsc using a trace formula of Ford Wunsc [6] sowed tat tis constant was in general optimal by proving existence of resonances wit Im λ ν 0 log Re λ, wile [9] refined te description of te resonances on and near tis curve. Closely related work of Datcev Kang Kessler [4] studied te distribution of resonances on surfaces of revolution wit a cone point and a funnel (as well as oter types of infinite ends). Te autors find examples were te Laplacian admits resonances wit Im λ C log Re λ as Re λ, altoug te classical flow is non-trapping; in suc cases te metric is continuous wit a conormal singularity. Te results ere are, to te best of our knowledge, te first results on resonances generated by diffractive trapping for semiclassical Scrödinger operators. Te crucial new tecnical ingredient is te propagation of singularities results recently obtained by te autors in [8], wic include estimates on te size of te wave of diffractively reflected singularities. Tese singularities bounce back off of even a mild singularity of V in violation of te most naive application of te principle of geometric optics, wic would say tat wavefront set travels along classical trajectories. Te autors intend in future work to complement te results ere, wic sow a resonance-free region in all dimensions and existence in one-dimension, by sowing te existence of resonances just below te gap obtained in Teorem 1 in all dimensions, at least in settings were te dynamics is tractable. We tus conjecture, based on te evidence of Teorem 2, tat te ν 0 obtained in Teorem 1 is optimal in general (at least generically). Acknowledgements. Te autors are grateful to Nicolas Burq and Jeff Galkowski for elpful discussions. OG was partially supported by NSF grant DMS ; JW was partially supported by NSF grant DMS Resonances 2.1. Complex scaling. We define resonances of P by te metod of complex scaling following [5, Sections 4.5, 6.2.1]. Fix R 0 > 0 suc tat outside of B(0, R 0 ), V = 0, g ij = δ ij,

6 6 ORAN GANNOT AND JARED WUNSCH and let R 1 > R 0. We define te complex scaled operator P θ using a contour Γ θ = f θ (R n ), were f θ (x) = x + i x F (x), and F is a smoot convex function satisfying F = 0 near B(0, R 1 ), F = tan θ 2 x 2 for x 2R 1. Here we take any fixed θ (0, π/2). In particular, te semiclassical principal symbol of p θ is given by { (1 + i 2 F (x)) 2 ξ, ξ, x > R 0, p θ (x, ξ) = ξ 2 g + V (x) x < R 1. Referring to te proof of [5, Proposition 6.10], we record te following important observation: for an interval I, exp thre pθ (ρ) {Im p θ = 0} for all t I = exp thre pθ (ρ) = exp thp (ρ) for all t I Te time-dependent Scrödinger equation. Define M = R t X, Y M = R t Y. (2.1) As a preliminary step, we pass to te time-dependent semiclassical Scrödinger operator on M given by Q = D t + P. If τ denotes te momentum dual to t, ten te principal symbol of Q (wic does not depend on t) is q(x, τ, ξ) = τ + p(x, ξ). Note tat τ is conserved under te H q flow, and t evolves at unit speed. We will also work wit te complex-scaled operator Q θ = D t + P θ, writing q θ = τ + p θ for its principal symbol. Since Im q θ (τ, x, ξ) = Im p θ (x, ξ) for all τ, observe tat (2.1) also olds wit q replacing p. Given E R, consider te joint energy surface {q = 0, τ = E}. In general, tere may be non-trivial beavior of solutions to Qw = 0 in te caracteristic set of (τ, ξ) 2 q at fiber-infinity in T M, since Q is not elliptic in te non-semiclassical sense. However, we explicitly avoid suc issues by restricting to τ = E. Tus (t, x, E, ξ) {q = 0, τ = E} (x, ξ) {p = E}, wit similar observations for p θ and q θ. We now consider propagation of singularities for te operator Q near T Y M M. Define te lift of H E by Ĥ E = {(t, x, E, ξ ) : (x, ξ ) H E } T Y M,

7 DIFFRACTIVE TRAPPING 7 wit te analogous definition for ĜE. We also write ˆπ : TmM TmY M for te projection wenever m Y M. Adapting te results of [8], we ave te following teorem on propagation of singularities for Q, stated locally: Teorem 3. Let w = w() be -tempered in H 1,loc (M) suc tat Qw = O( ) L 2 loc. (1) If q ĤE, let µ ± {q = 0, τ = E} be te preimages of q under ˆπ wit opposite normal momenta. If µ + WF s (w) for some s R, ten tere exists ε > 0 suc tat exp thq (µ + ) WF s (w), or exp thq (µ ) WF s α (w), or bot, for all t (0, ε). (2) If q ĜE, let µ {q = 0, τ = E} be te unique preimage of q under ˆπ wit vanising normal momentum. If µ WF s (w) for some s R, ten tere exists ε > 0 suc tat exp thq (µ) WF s (w) for all t (0, ε). Te (minor) modifications to [8] needed to prove Teorem 3 are outlined in Appendix A. Note tat singularities propagate straigt troug ˆπ 1 (ĜE). Combined wit ordinary propagation of singularities for Q θ away from T Y M M, we obtain te following: Corollary 4. Let w = w() be -tempered in H 1,loc (M) suc tat Q θw = O( ) L 2 loc. Let s R and T 0. If µ {q θ = 0, τ = E} is suc tat exp thre qθ (µ) is disjoint from ˆπ 1 (ĤE) for eac t [0, T ], ten exp T HRe qθ (µ) / WF s (v) = µ / WF s (v). Proof. Suppose conversely tat µ WF s (v), and define T 0 = sup{t 0 : exp t H Re qθ (µ) WF s (v) for all t [0, t]}. Note tat µ = exp T0 H Re qθ (µ) WF s (v) since wavefront set is closed. We claim tat T 0 T, wic tus completes te proof. Indeed, if T 0 < T, ten by ypotesis µ / T Y M M or µ ˆπ 1 (ĜE). In te first case Q θ is smoot in a neigborood of µ, and since Im q θ 0 we can apply propagation of singularities forward along te H Re qθ flow to deduce a contradiction. In te second case q = q θ in a neigborood of µ, so we can apply te second part of Teorem 3 to deduce a contradiction. We make te following dynamical definitions, recalling te definition of R 1 from Section 2.1.

8 8 ORAN GANNOT AND JARED WUNSCH Definition 2.1. Let E R and ρ {p = E} \ π 1 (H E ). We say tat ρ IC E if tere exists T 0 0 suc tat exp thp (µ) is disjoint from π 1 (H E ) for eac t [0, T 0 ] and x(exp T0 H p (ρ)) 2R 1. If ρ {p = E} \ π 1 (H E ), ten we say tat ρ Γ ± E exp thp (ρ) π 1 (H E ). if tere exists t > 0 suc tat Tus, IC E consists of points tat are incoming from infinity witout interaction wit te yperbolic region, wile Γ ± E represent points itting te yperbolic region in backward-/forward-time (using notation in loose analogy wit tat of unstable/stable manifolds). Note tat { x 2R 1 } {Im p θ 0}. (2.2) Te next observation follows immediately from Corollary 4. Lemma 5. Let w = w() be -tempered in H 1,loc (M) suc tat Q θw = O( ) L 2 loc, and let E R. If (x, ξ) IC E and we set µ = (t, x, E, ξ) {q θ = 0, τ = E}, ten µ / WF s (w) for all s R. Proof. Let T 0 > 0 be as in Definition 2.1. It cannot be tat exp thre qθ (µ) {Im q θ = 0} for all t [0, T 0 ], since using Im q θ (τ, x, ξ) = Im p θ (x, ξ) and (2.1), we would conclude tat exp thre qθ (µ) = exp thq (µ) for all t [0, T 0 ], (2.3) and in particular exp T0 H q (µ) {Im q θ = 0}; tis contradicts te definition of T 0 according to (2.2). Tus we can define T = inf{t [0, T 0 ] : exp T HRe qθ (µ) {Im q θ 0}}. Note tat µ = exp T HRe qθ (µ) / T Y M M since Im q θ = 0 in a neigborood of T Y M Y. By semiclassical ellipticity in te smoot setting and te definition of T, tere exists δ > 0 suc tat µ = exp (T +δ)hre qθ (µ) / WF (w), and ence by forward propagation in te smoot setting we ave µ / WF (w) as well. On te oter and, by te definition of T we must ave tat (2.3) olds for all t [0, T ]. Since (x, ξ) IC E, it follows tat te backward H Re qθ flow from µ to µ is disjoint from ˆπ 1 (ĤE), and ence µ / WF s (w) for all s R by Corollary 4. Now we make te assumption tat E is a non-trapping energy level. Lemma 6. Let w = w() be -tempered in H 1,loc (M) suc tat Q θw = O( ) L 2 loc, and let E > 0 be non-trapping. If tere exists µ WF s (w) {q θ = 0, τ = E}

9 DIFFRACTIVE TRAPPING 9 for some s R, ten tere exists r R and µ = (t, x, E, ξ) WF r (w) {q θ = 0} suc tat (x, ξ) Γ E. Proof. First assume µ / ˆπ 1 (ĤE). Arguing precisely as in Lemma 5, te backwards flow exp thre qθ (µ) must encounter a yperbolic point µ + for some t > 0. Oterwise, if µ ˆπ 1 (ĤE) to begin wit, simply set µ + = µ. Now let µ project to te same yperbolic point as µ +, but wit opposite normal momentum. We know from te first part of Teorem 3 tat tere exists µ wit te requisite properties, obtained by flowing backwards along H q from eiter µ + (taking r = s) or µ (taking r = s α) for a sort time ε > 0, noting tat q = Re q θ = q θ near µ ±. Note in Lemma 6 we assume tat µ is in bot WF s (w) and te caracteristic set {q θ = 0}, since te inclusion of former in te latter is only guaranteed for a certain range of s owing to te singularity of V (see [8, Proposition 7.5]). Lemma 7. Let E > 0 be non-trapping and ρ Γ E. If diam E(Y ) > 0 and T > diam E (Y ), ten If diam E (Y ) = 0, ten ρ IC E. exp T Hp (ρ) IC E. Proof. Since ρ Γ E, tere exists t 0 > 0 suc tat ρ 0 = exp t0 H p (ρ) π 1 (H E ). Clearly if diam E (Y ) = 0, ten ρ IC E. Oterwise, if 0 < diam E (Y ) < T, set ρ = exp T Hp (ρ 0 ) = exp (t0 T )H p (ρ). Ten exp thp (ρ ) is disjoint from π 1 (H E ) for all t 0, and since E is non-trapping te proof is finised. We record a key lemma, wic is a refinement of Lemma 6. Lemma 8. Let w = w() be -tempered in H 1,loc (M) suc tat Q θw = O( ) L 2 loc, and let E > 0 be nontrapping. If (x, ξ) Γ E and µ = (0, x, E, ξ) WF s (w) {q θ = 0} for some s R, ten for all δ > 0 tere exists T 0 diam E (Y ) + δ and suc tat (x 0, ξ 0 ) Γ E. µ = ( T 0, x 0, E, ξ 0 ) WF s α (w) {q θ = 0}

10 10 ORAN GANNOT AND JARED WUNSCH Proof. Define t 0 sup { t 0: exp t H Re qθ (µ) WF s (w) for all t [0, t] }. We claim tat t 0 always exists and t 0 diam E (Y ). To see tis, first note tat for t 0, exp t H Re qθ (µ) WF s (w) for all t [0, t] = exp t H Re qθ (µ) = exp t H q (µ) for all t [0, t]. Indeed, if te conclusion fails, ten by (2.3) tere must be t [0, t] suc tat exp t H Re qθ (µ) Im{q θ 0}. But in te neigborood of suc a point q θ is smoot, and ence by semiclassical ellipticity exp t H Re qθ (µ) / WF s (w), wic is a contradiction. Furtermore, by Lemmas 5 and 7, if T > diam E (Y ), ten wic sows tat t < T. exp T Hq (µ) / WF s (w), Now set µ + = exp t0 H Re qθ (µ); we ave µ + WF s (w) since wavefront set is closed. Also by definition of t 0, tere is a sequence ε n > 0 wit ε n 0 suc tat exp εnh Re qθ (µ + ) / WF s (w), By ordinary propagation of singularities of Q θ away from T Y M Y we must ave µ + T Y M M, and by Teorem 3, ˆπ(µ + ) must be a yperbolic point. Let µ denote te point in ˆπ 1 (ˆπ(µ + )) wit opposite normal momentum. Ten Teorem 3 sows tat µ exp δhre qθ (µ ) WF s α (w) for δ > 0 arbitrarily small. Now T 0 = t(µ ) = t 0 + δ diam E (Y ) + δ as desired. We are now ready to prove Teorem Resonance widts. We begin by working quite generally, witout prejudice as to te operator P. Suppose tat u = u() is -tempered in H,loc 1 (X). Given a family z = z() C, set Im z() ν() = log(). Now form te functions w(t, x) = e izt/ u(x), and note tat a sufficient condition to guarantee tat w = w() is -tempered in H 1,loc (M) is ν() = O(1). Furtermore, if (P z)u = O( ) L 2 loc, ten Qw = O( ) L 2 loc. We need a simple lemma comparing te wavefront of w wit its restriction to a fixed time slice {t = t 0 }. For our purposes it will suffice to consider wavefront set away

11 DIFFRACTIVE TRAPPING 11 from fiber-infinity. In addition, to te assumption ν() = O(1), we also assume tat z() = E + o(1) for some E R. Note tat WF s (w) T M {τ = E} for all s, since (D t + z)w = 0 and D t + z is elliptic wen τ E. Of course tis leaves open te possibility tat w as wavefront set at fiber-infinity in a direction wit τ = 0. We need to be sligtly more precise: let χ 0 Cc (R) ave support near E and χ 1 Cc (R) ave support near t 0. Ten for eac K X compact, after integrating by parts. For eac fixed t 0 R, introduce te set (1 χ 0 (D t ))χ 1 (t)w L 2 (R K) = O( ) (2.4) W s (t 0 ) = {(x, ξ) T X : (t 0, x, E, ξ) WF s (w)}. We ave te following lifting lemma: Lemma 9. If z = z() C satisfies z() = E + o(1) and ν() ν 0 for some E R and ν 0 > 0, ten for eac s, t 0 R and δ > 0 W s δ (t 0 ) WF s ( ν()t 0 u) T X W s (t 0 ). Proof. First observe tat WF s (w(t 0, )) = WF s ( ν()t 0 u) for eac fixed t 0 R. Now suppose tat (x 0, ξ 0 ) T X satisfies (t 0, x 0, E, ξ 0 ) / WF s (w). For appropriate cutoffs φ 0, φ 1 C c (X), and χ 1 C c (R) supported near ξ 0, x 0, and t 0, respectively, (2.4) implies tat φ 0 (D x )φ 1 (x)χ 1 (t)w = O( s ) L 2 (M), wic immediately implies tat (x 0, ξ 0 ) / WF s (w(t 0, )). (x 0, ξ 0 ) / WF s (w(t 0, )) T X, so in te notation above, e izt 0/ φ 0 (D x )φ 1 (x)u = O( s ) L 2 (X). Conversely, suppose tat Now if χ 1 as above as sufficiently small support depending on δ > 0, ten we can arrange tat φ 0 (D x )φ 1 (x)χ 1 (t)w = O( s δ ) L 2 (M). Again using (2.4), tis implies tat (x 0, ξ 0 ) / W s δ (t 0 ) for eac δ > 0. Now we place ourselves in te setting of Teorem 1. If tere is a resonance in te set [E 0 δ, E 0 + δ] i[0, ν 0 log(1/)] for some ν 0 > 0, ten we can find a sequence k 0, complex numbers z( k ) satisfying Re z( k ) = E + o(1), Im z( k ) ν 0 k log(1/ k ),

12 12 ORAN GANNOT AND JARED WUNSCH for some E E 0 < δ, and eigenfunctions u( k ) L 2 (X) of te complex-scaled operator, (P θ z( k ))u( k ) = 0. Trougout, we suppress te index k. Normalize u by u L 2 = 1. In order to eventually apply Teorem 3, we verify tat te H 1 (in fact, te H2 ) norms of u are also uniformly bounded in. Tis follows from standard semiclassical analysis by observing tat P θ V is smoot, tat uniformly on X, and tat V L (X). σ (P θ V )(x, ξ) δ ξ 2 C 0 Observe tat WF s (u) for eac s > 0, since u L 2 = 1. Taking s = α, we conclude tat WF α (u) { ξ 2 (p θ E) = 0} Indeed, away from Y tis is ordinary semiclassical elliptic regularity (ence applies wit any s), wereas near T Y X, were p θ = p, we can apply [8, Proposition 7.5] (wose proof applies nearly verbatim even wen z is not real valued). Of course since p θ is elliptic at fiber-infinity, tis can be rewritten as WF α (u) {p θ = E}. Next, we let w = e izt/ u wit u = u() our family of eigenfunctions as above, and observe tat Lemma 9 applies to te family w = w(). Proof of Teorem 1. As noted above, tere exists (x 0, ξ 0 ) WF α (u) {p θ = E}. Let µ 0 = (0, x 0, E, ξ 0 ) {q θ = 0, τ = E} and note tat µ 0 WF α (w) by Lemma 9. Applying Lemmas 6 and 9 we conclude tere exists ρ WF r (u) Γ E for some r R. Tis is already a contradiction if diam E (Y ) = 0, since by Lemmas 5, 7, and 9 we ave ρ IC E and ence ρ / WF r (u). Oterwise, we derive a lower bound for ν 0 as follows. Since WF 0 (u) =, we may define s 0 = inf{s : WF s (u) Γ E }. Pick any s > s 0, so tere is (x, ξ) WF s (u) Γ E, and apply Lemma 8 to µ = (0, x, E, ξ), wic by Lemma 9 indeed satisfies µ WF s (w) {q θ = 0}. Tus for eac δ > 0 we can find T < diam E (Y ) + δ suc tat µ = ( T, x 1, E, ξ) WF s α (w) {q θ = 0}

13 DIFFRACTIVE TRAPPING 13 and (x 1, ξ 1 ) Γ E. Furtermore, by Lemma 9, (x 1, ξ 1 ) WF s α+ν 0T +δ (u), and ence (x 1, ξ 1 ) WF s α+ν 0T +δ (u) Γ E. By definition of s 0, we must ave s α + ν 0 T + δ > s 0. But s > s 0 and δ > 0 are bot arbitrary, wic implies tat ν 0 α/ diam Y (E). Tis completes te proof of te Teorem Existence of resonances in one dimension 3.1. WKB solutions. Trougout tis section we adopt te notation of Teorem 2 and te paragrap preceding it. Since supp V [0, L], a complex number z C is a resonance of P = (D x ) 2 + V precisely if it satisfies (P z)u = 0, D x u(0) + z 1/2 u(0) = 0, (3.1) D x u(l) z 1/2 u(l) = 0. Fix an interval [a, b] wit a > sup V. Given M > 0, set Ω M () = [a, b] + i[ M log(1/), 0]. Te fact tat I = [0, L] is a classically allowed region implies te following: Lemma 10. If u solves te equation (P z)u = 0 wit initial conditions u(0) = O(1), D x u(0) = O(1), ten u and x u are polynomially bounded on I = [0, L] uniformly in z Ω M () Proof. Tis follows from a semiclassical version of te energy estimates in [10, Section 23.2]; cf. [5, Lemma E.60]. At most polynomial growt in arises from te fact tat te imaginary part of z can be of size O( log(1/)). Next we consider approximate WKB solutions to (P z)u = 0. Following [2], for tis problem it is convenient to consider approximate solutions in exponential form. Define ψ 0 = (z V ) 1/2, and ten set ψ +,0 = ψ,0 = ψ 0. For i 1, define ψ ±,i recursively by ψ ±,k (x) = ± i 2ψ 0 (x) xψ ±,k 1 (x) 1 k 1 ψ ±,j (x)ψ ±,k j (x). (3.2) 2ψ 0 (x) j=1

14 14 ORAN GANNOT AND JARED WUNSCH Eac ψ ±,k is smoot on I and depends olomorpically on z [a, b] + i[ C 0, 0]. Let ψ ± be a function admitting an asymptotic expansion ψ ± j ψ ±,j j=0 and depending olomorpically on z. We ten set ϕ ± (x) = x 0 ψ ± (s) ds, u ± = e ±iϕ ±/. If z Ω M (), ten u ± are polynomially bounded on I, and ence te recursion relation (3.2) guarantees tat (P z)u ± = O( ) C (I), z Ω M () for any fixed M > 0. Te usefulness of tis exponential form comes from te following observation (cf. [2, Appendix 2]). Lemma 11. For eac j, ψ ±,j (x) = ( ) j ±i ψ (j) 0 (x) + F j (ψ 0 (x), ψ 2ψ 0 (x) 0(x),..., ψ (j 1) 0 (x)) for a smoot function F j (t 0, t 1,..., t j 1 ) suc tat F j (t 0, 0,..., 0) = 0 for all t 0 R. Proof. Tis follows by induction from te recursion relations (3.2). We ten obtain te following corollary: if V vanises to order k 1 at a point x 0 I, ten ψ ±,j (x 0 ) = 0 for j = 1,..., k 1, and ψ ±,k (x 0 ) = i ±k (2z 1/2 ) k 1 V (k) (x 0 ). (3.3) If x 0 is an endpoint of I = [0, L], ten te same formulas old in te sense of one-sided limits. Indeed, V (x) = c 0 (x x 0 ) k + O( x x 0 k+1 ), c 0 = V (k) 0 (x 0 )/k!, wic sows tat ψ 0 (x) = z 1/2 (c 0 /2)z 1/2 (x x 0 ) k + O( x x 0 k+1 ). Te formula (3.3) follows immediately from tis expression.

15 DIFFRACTIVE TRAPPING Outgoing condition. Observe tat u ± (0) = 1 and D x u ± (0) = ±ψ ± (0), so if we form te function v = u z1/2 ψ (0) z 1/2 + ψ + (0) u +, ten D x v(0) = z 1/2 v(0), and v is polynomially bounded on I for z Ω M (). Lemma 12. Let u solve te equation (P z)u = 0 wit initial data u(0) = v(0) and D x u(0) = z 1/2 u(0). Ten uniformly for z Ω M (). u = v + O( ) C (I) Proof. Tis follows from te identity u = W (u, u + )u W (u, u )u +, W (u, u + ) were W (f, g) = f ( x g) ( x f) g is te semiclassical Wronskian. Indeed, te fact tat u and u ± are polynomially bounded on I for z Ω M () implies tat all Wronskians appearing in te formula above are constant modulo O( ). It is ten a straigtfoward computation of te Wronskians at x = 0 using te specified initial conditions. Fix M > 0. In view of (3.1) and Lemma 12, resonances z Ω M () are caracterized as solutions to an equation of te form ( ) ( ) ψ+ (L) z 1/2 ψ (0) z 1/2 e iϕ +(L)/ e iϕ (L)/ = O( ). (3.4) ψ + (L) + z 1/2 ψ (0) + z 1/2 Bot sides of tis equation are olomorpic in z Ω M (). We replace tis wit a simpler expression by inserting te asymptotics of ψ ±, making sure to only incur errors tat are olomorpic in z Ω M (); we will continue to use ordinary Landau notation to denote tese errors. First, note tat according to (3.3), ψ + (L) z 1/2 = i l l (2z 1/2 ) l 1 V (l) (L )(1 + O()) ψ (0) z 1/2 = i k k (2z 1/2 ) k 1 V (k) (0 + )(1 + O()), ψ + (L) + z 1/2 = 2z 1/2 + O( l ), ψ (0) + z 1/2 = 2z 1/2 + O( k ). We also multiply (3.4) troug by e iϕ (L)/ = O( N ) (for some N N). Define te pase function ϕ(x) = x 0 (z V (s)) 1/2 ds, (3.5) and observe tat ϕ + (L) + ϕ (L) = 2ϕ(L) + O( 2 ) since ψ +,1 = ψ,1. In particular, e i(ϕ +(L)+ϕ (L))/ = e 2iϕ(L)/ (1 + O()).

16 16 ORAN GANNOT AND JARED WUNSCH Inserting tis information into (3.4), we obtain an equivalent equation i l k l+k (2z 1/2 ) l k 4 (V (k) (0 + ) V (l) (L ))e 2iϕ(L)/ 1 = O(). Tis equation already implies te necessity of Teorem 2 (i.e., te fact tat resonances in Ω M () may only be of te form (1.3)) by considering te modulus and argument of tis equation. To sow existence of tese resonances, consider te function F (w, E, ) = i l k l+k (2E 1/2 ) l k 4 (V (k) (0 + ) V (l) (L ))e 2i(S(E)+wT (E))/ 1. Here E [a, b] is treated as a parameter, and F (w, E, ) is olomorpic in w C. We ave used te notation S(E) and T (E) from Teorem 2. Note tat if z = E + w Ω M (), ten ϕ(l) = S(E) + wt (E) + O( w 2 ). If w and w 2 / are bot small, ten z is resonance if and only if w satisfies an equation of te form F (w, E, ) = O( + w + w 2 /) (3.6) wit bot sides olomorpic in w. Now if n N() and E n are as in Teorem 2, ten F (w n, E n, ) = 0 for te coice w n = i (l + k) log(1/) 2T (E n ) + i ( log V (k) (0 + )V (l) (L ) (1/2)(l + k + 4) log(4e n ) ). 2T (E n ) Furtermore, w F (w n, E n, ) = (2i/)T (E n ), and given ε > 0 tere exists C > 0 suc tat n N() and w w n ε = 2 wf (w, E n, ) C 2. Tus by Taylor s teorem, for any A > 0 we can find 0, C 0 > 0 suc tat F (w, E n, ) A log(1/) 2 for all (0, 0 ) and n N() wenever w w n = C 0 2 log(1/) 2. On te oter and, tere exists B > 0 independent of C 0 suc tat rigt and side of (3.6) is bounded by B log(1/) 2 wenever w w n = C 0 2 log(1/) 2 and (0, 0 ) uniformly in n N(), srinking 0 > 0 depending on C 0 if necessary. Tus we first fix A > B, coose 0, C 0 > 0 as above, and ten apply Roucé s teorem. It follows tat for eac (0, 0 ) and n N() tere exists a unique resonance z n Ω M () satisfying z n = E n + w n + O( 2 log(1/) 2 ), tus completing te proof of Teorem 2.

17 DIFFRACTIVE TRAPPING 17 Appendix A. Proof of Teorem 3 In [8], propagation of singularities for operators of te form P = 2 + V wit V real-valued was discussed. Altoug Q is not elliptic at fiber-infinity (unlike P ), since we are only considering propagation of singularities in compact subsets of T M viewed as te interior of T M, tere is little difference in te proofs. Tere are two ingredients: (1) Propagation of singularities along generalized broken bicaracteristics (GBBs). (2) Diffractive improvements at yperbolic and glancing points. Te main difference is in te propagation arguments along GBBs. Te preliminary material in [8, Section 4] goes troug uncanged provided p is replaced wit q = τ + ξ 2 g + V. In te neigborood of a point m Y M we can use coordinates (t, x 1, x ), were (x 1, x ) are normal coordinates on Y wit respect to g. We ten replace P wit Q = Q (D x1 ) (D x1 ) in tese local coordinates. Te elliptic, yperbolic, and glancing sets are defined in te obvious way; we denote tese Ê, Ĥ, Ĝ. In te notation of Section 2.2, Ĥ E = Ĥ {τ = E}, Ĝ E = Ĝ {τ = E}. Note tat Ĥ Ĝ may intersect fiber-infinity of T Y M in directions wit τ 0, but of course ĤE ĜE are compact subsets of T Y M. Te study of te elliptic set Ê in [8, Section 5.1] needs only minor modifications provided we stay away from fiber infinity; tis amounts to working wit compactly microlocalized b-pseudodifferential operators only. Te most important point is to replace te Diriclet form associated to P in [8, Lemma 5.3] wit te expression 2 d X Aw 2 g + V Aw 2 + D t Aw Aw d g, M were d X is te differential on X lifted to M = R X by te product structure, and g is te product metric dt 2 + g. Te proof of [8, Lemma 5.3] applies witout cange. Modifying te Diriclet form affects [8, Eq. 5.3]; in te notation tere, te relevant replacement is an estimate of te form Aw H 1 C GQw H 1 + C Gw H 1 + C 0 ( Aw L 2 + D t Aw L 2) + O( ) w H 1,

18 18 ORAN GANNOT AND JARED WUNSCH wit C 0 > 0 independent of A, namely one must add D t Aw L 2 to te rigt and side. On te oter and, if A as compact b-microsupport, ten we can estimate D t Aw L 2 C 1 Aw L 2 + O( ) w L 2, were C 1 > 0 depends only on te size of te b-microsupport of A; tus te important estimate [8, Eq. 5.3] is still valid. Te proofs of [8, Lemmas 5.5, 5.6] go troug, provided V is replaced wit V + τ. Wit tese modifications, b-elliptic regularity ([8, Proposition 5.2]) continues to old, at least away from fiber-infinity. Te analysis at Ĥ is essentially uncanged, provided one only consider propagation of singularities away from fiber-infinity; tis is certainly te case near ĤE. One needs only to account for te additional localization in (t, τ), and in te tird step of te proof of [8, Lemma 5.10] we make te replacement 2 k = 2 k + D t. In tat case, [8, Proposition 5.8] still olds. (A.1) Finally, we come to te analogues of [8, Teorems 2 and 3], wic are proved using ordinary pseudodifferential operators. Again, since in Teorem 3 only points µ in te finite parts of te fiber of T M (rater tan fiber-infinity in te compactification) are considered, te proofs in [8, Section 7] are still valid provided we make te replacement (A.1). References [1] Dean Baskin and Jared Wunsc. Resolvent estimates and local decay of waves on conic manifolds. J. Differential Geom., 95(2): , [2] Micael V Berry. Semiclassically weak reflections above analytic and non-analytic potential barriers. Journal of Pysics A: Matematical and General, 15(12):3693, [3] Nicolas Burq. Pôles de diffusion engendrés par un coin. Astérisque, (242):ii+122, [4] Kiril R Datcev, Daniel D Kang, and Andre P Kessler. Nontrapping surfaces of revolution wit long living resonances. arxiv: , [5] Semyon Dyatlov and Maciej Zworski. Matematical teory of scattering resonances. mat.mit. edu/~dyatlov/res/. [6] G Austin Ford and Jared Wunsc. Te diffractive wave trace on manifolds wit conic singularities. Advances in Matematics, 304: , [7] Jeffrey Galkowski. A quantitative Vainberg metod for black box scattering. Communications in Matematical Pysics, 349(2): , [8] Oran Gannot and Jared Wunsc. Semiclassical diffraction by conormal potential singularities. arxiv: , [9] Luc Hillairet and Jared Wunsc. On resonances generated by conic diffraction. arxiv: , [10] Lars Hörmander. Te Analysis of Linear Partial Differential Operators, volume 3. Springer- Verlag, Berlin, Heidelberg, New York, Tokyo, [11] André Martinez. Resonance free domains for non globally analytic potentials. Annales Henri Poincaré, 3(4): , 2002.

19 DIFFRACTIVE TRAPPING 19 [12] Tullio Regge. Analytic properties of te scattering matrix. Nuovo Cimento (10), 8: , [13] Jared Wunsc. Resolvent estimates wit mild trapping, Journées Équations aux Dérivées Partielles. [14] Maciej Zworski. Distribution of poles for scattering on te real line. J. Funct. Anal., 73: , Department of Matematics, Lunt Hall, Nortwestern University, Evanston, IL 60208, USA address: and

A SHORT INTRODUCTION TO BANACH LATTICES AND

A SHORT INTRODUCTION TO BANACH LATTICES AND CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,

More information

Poisson Equation in Sobolev Spaces

Poisson Equation in Sobolev Spaces Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on

More information

ERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*

ERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER* EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T

More information

SEMICLASSICAL DIFFRACTION BY CONORMAL POTENTIAL SINGULARITIES

SEMICLASSICAL DIFFRACTION BY CONORMAL POTENTIAL SINGULARITIES SEMICLASSICAL DIFFRACTION BY CONORMAL POTENTIAL SINGULARITIES ORAN GANNOT AND JARED WUNSCH Abstract. We establis propagation of singularities for te semiclassical Scrödinger equation, were te potential

More information

MA455 Manifolds Solutions 1 May 2008

MA455 Manifolds Solutions 1 May 2008 MA455 Manifolds Solutions 1 May 2008 1. (i) Given real numbers a < b, find a diffeomorpism (a, b) R. Solution: For example first map (a, b) to (0, π/2) and ten map (0, π/2) diffeomorpically to R using

More information

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx. Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions

More information

Differentiation in higher dimensions

Differentiation in higher dimensions Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends

More information

ch (for some fixed positive number c) reaching c

ch (for some fixed positive number c) reaching c GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan

More information

BOUNDARY REGULARITY FOR SOLUTIONS TO THE LINEARIZED MONGE-AMPÈRE EQUATIONS

BOUNDARY REGULARITY FOR SOLUTIONS TO THE LINEARIZED MONGE-AMPÈRE EQUATIONS BOUNDARY REGULARITY FOR SOLUTIONS TO THE LINEARIZED MONGE-AMPÈRE EQUATIONS N. Q. LE AND O. SAVIN Abstract. We obtain boundary Hölder gradient estimates and regularity for solutions to te linearized Monge-Ampère

More information

Math 161 (33) - Final exam

Math 161 (33) - Final exam Name: Id #: Mat 161 (33) - Final exam Fall Quarter 2015 Wednesday December 9, 2015-10:30am to 12:30am Instructions: Prob. Points Score possible 1 25 2 25 3 25 4 25 TOTAL 75 (BEST 3) Read eac problem carefully.

More information

Polynomial Interpolation

Polynomial Interpolation Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc

More information

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of

More information

Chapter 5 FINITE DIFFERENCE METHOD (FDM)

Chapter 5 FINITE DIFFERENCE METHOD (FDM) MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential

More information

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra

More information

Weierstrass and Hadamard products

Weierstrass and Hadamard products August 9, 3 Weierstrass and Hadamard products Paul Garrett garrett@mat.umn.edu ttp://www.mat.umn.edu/ garrett/ [Tis document is ttp://www.mat.umn.edu/ garrett/m/mfms/notes 3-4/b Hadamard products.pdf].

More information

Diffraction by Edges. András Vasy (with Richard Melrose and Jared Wunsch)

Diffraction by Edges. András Vasy (with Richard Melrose and Jared Wunsch) Diffraction by Edges András Vasy (with Richard Melrose and Jared Wunsch) Cambridge, July 2006 Consider the wave equation Pu = 0, Pu = D 2 t u gu, on manifolds with corners M; here g 0 the Laplacian, D

More information

Numerical Differentiation

Numerical Differentiation Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function

More information

A Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems

A Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems

More information

arxiv: v1 [math.dg] 4 Feb 2015

arxiv: v1 [math.dg] 4 Feb 2015 CENTROID OF TRIANGLES ASSOCIATED WITH A CURVE arxiv:1502.01205v1 [mat.dg] 4 Feb 2015 Dong-Soo Kim and Dong Seo Kim Abstract. Arcimedes sowed tat te area between a parabola and any cord AB on te parabola

More information

Fractal Weyl Laws and Wave Decay for General Trapping

Fractal Weyl Laws and Wave Decay for General Trapping Fractal Weyl Laws and Wave Decay for General Trapping Jeffrey Galkowski McGill University July 26, 2017 Joint w/ Semyon Dyatlov The Plan The setting and a brief review of scattering resonances Heuristic

More information

Finite Difference Method

Finite Difference Method Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According

More information

4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.

4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these. Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra

More information

HOMEWORK HELP 2 FOR MATH 151

HOMEWORK HELP 2 FOR MATH 151 HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,

More information

1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)

1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x) Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of

More information

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note

More information

Average Rate of Change

Average Rate of Change Te Derivative Tis can be tougt of as an attempt to draw a parallel (pysically and metaporically) between a line and a curve, applying te concept of slope to someting tat isn't actually straigt. Te slope

More information

Preconditioning in H(div) and Applications

Preconditioning in H(div) and Applications 1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition

More information

THE IMPLICIT FUNCTION THEOREM

THE IMPLICIT FUNCTION THEOREM THE IMPLICIT FUNCTION THEOREM ALEXANDRU ALEMAN 1. Motivation and statement We want to understand a general situation wic occurs in almost any area wic uses matematics. Suppose we are given number of equations

More information

POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY

POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY APPLICATIONES MATHEMATICAE 36, (29), pp. 2 Zbigniew Ciesielski (Sopot) Ryszard Zieliński (Warszawa) POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY Abstract. Dvoretzky

More information

Click here to see an animation of the derivative

Click here to see an animation of the derivative Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,

More information

arxiv: v1 [math.ag] 13 Jan 2019

arxiv: v1 [math.ag] 13 Jan 2019 Adiabatic Limit and Deformations of Complex Structures Dan Popovici ariv:1901.04087v1 [mat.ag] 13 Jan 2019 Abstract. Based on our recent adaptation of te adiabatic limit construction to te case of complex

More information

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department

More information

Introduction to Derivatives

Introduction to Derivatives Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))

More information

Chapter 2 Limits and Continuity

Chapter 2 Limits and Continuity 4 Section. Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 6) Quick Review.. f () ( ) () 4 0. f () 4( ) 4. f () sin sin 0 4. f (). 4 4 4 6. c c c 7. 8. c d d c d d c d c 9. 8 ( )(

More information

LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION

LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y

More information

RIGIDITY OF TIME-FLAT SURFACES IN THE MINKOWSKI SPACETIME

RIGIDITY OF TIME-FLAT SURFACES IN THE MINKOWSKI SPACETIME RIGIDITY OF TIME-FLAT SURFACES IN THE MINKOWSKI SPACETIME PO-NING CHEN, MU-TAO WANG, AND YE-KAI WANG Abstract. A time-flat condition on spacelike 2-surfaces in spacetime is considered ere. Tis condition

More information

1. Introduction. Consider a semilinear parabolic equation in the form

1. Introduction. Consider a semilinear parabolic equation in the form A POSTERIORI ERROR ESTIMATION FOR PARABOLIC PROBLEMS USING ELLIPTIC RECONSTRUCTIONS. I: BACKWARD-EULER AND CRANK-NICOLSON METHODS NATALIA KOPTEVA AND TORSTEN LINSS Abstract. A semilinear second-order parabolic

More information

Polynomial Interpolation

Polynomial Interpolation Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x

More information

Nonlinear elliptic-parabolic problems

Nonlinear elliptic-parabolic problems Nonlinear elliptic-parabolic problems Inwon C. Kim and Norbert Požár Abstract We introduce a notion of viscosity solutions for a general class of elliptic-parabolic pase transition problems. Tese include

More information

SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY

SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY (Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative

More information

A = h w (1) Error Analysis Physics 141

A = h w (1) Error Analysis Physics 141 Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.

More information

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225 THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:

More information

Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems

Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for

More information

Convexity and Smoothness

Convexity and Smoothness Capter 4 Convexity and Smootness 4.1 Strict Convexity, Smootness, and Gateaux Differentiablity Definition 4.1.1. Let X be a Banac space wit a norm denoted by. A map f : X \{0} X \{0}, f f x is called a

More information

Gluing semiclassical resolvent estimates via propagation of singularities

Gluing semiclassical resolvent estimates via propagation of singularities Gluing semiclassical resolvent estimates via propagation of singularities Kiril Datchev (MIT) and András Vasy (Stanford) Postdoc Seminar MSRI Kiril Datchev (MIT) Semiclassical resolvent estimates 7 December

More information

The derivative function

The derivative function Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL IFFERENTIATION FIRST ERIVATIVES Te simplest difference formulas are based on using a straigt line to interpolate te given data; tey use two data pints to estimate te derivative. We assume tat

More information

Generic maximum nullity of a graph

Generic maximum nullity of a graph Generic maximum nullity of a grap Leslie Hogben Bryan Sader Marc 5, 2008 Abstract For a grap G of order n, te maximum nullity of G is defined to be te largest possible nullity over all real symmetric n

More information

Notes on wavefunctions II: momentum wavefunctions

Notes on wavefunctions II: momentum wavefunctions Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles

More information

Quantum Numbers and Rules

Quantum Numbers and Rules OpenStax-CNX module: m42614 1 Quantum Numbers and Rules OpenStax College Tis work is produced by OpenStax-CNX and licensed under te Creative Commons Attribution License 3.0 Abstract Dene quantum number.

More information

Subdifferentials of convex functions

Subdifferentials of convex functions Subdifferentials of convex functions Jordan Bell jordan.bell@gmail.com Department of Matematics, University of Toronto April 21, 2014 Wenever we speak about a vector space in tis note we mean a vector

More information

lecture 26: Richardson extrapolation

lecture 26: Richardson extrapolation 43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)

More information

On convexity of polynomial paths and generalized majorizations

On convexity of polynomial paths and generalized majorizations On convexity of polynomial pats and generalized majorizations Marija Dodig Centro de Estruturas Lineares e Combinatórias, CELC, Universidade de Lisboa, Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal

More information

NUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,

NUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example, NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing

More information

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems 5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we

More information

Math Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim

Math Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z

More information

INTRODUCTION AND MATHEMATICAL CONCEPTS

INTRODUCTION AND MATHEMATICAL CONCEPTS Capter 1 INTRODUCTION ND MTHEMTICL CONCEPTS PREVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips

More information

DYNAMICAL ZETA FUNCTIONS FOR ANOSOV FLOWS VIA MICROLOCAL ANALYSIS

DYNAMICAL ZETA FUNCTIONS FOR ANOSOV FLOWS VIA MICROLOCAL ANALYSIS DYNAMICAL ZETA FUNCTIONS FOR ANOSOV FLOWS VIA MICROLOCAL ANALYSIS SEMYON DYATLOV AND MACIEJ ZWORSKI Abstract. Te purpose of tis paper is to give a sort microlocal proof of te meromorpic continuation of

More information

arxiv: v2 [math.na] 5 Jul 2017

arxiv: v2 [math.na] 5 Jul 2017 Trace Finite Element Metods for PDEs on Surfaces Maxim A. Olsanskii and Arnold Reusken arxiv:1612.00054v2 [mat.na] 5 Jul 2017 Abstract In tis paper we consider a class of unfitted finite element metods

More information

More on generalized inverses of partitioned matrices with Banachiewicz-Schur forms

More on generalized inverses of partitioned matrices with Banachiewicz-Schur forms More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,

More information

Exam 1 Review Solutions

Exam 1 Review Solutions Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),

More information

Continuity and Differentiability of the Trigonometric Functions

Continuity and Differentiability of the Trigonometric Functions [Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te

More information

Journal of Computational and Applied Mathematics

Journal of Computational and Applied Mathematics Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing

More information

LECTURE NOTES ON QUANTUM CHAOS COURSE , MIT, APRIL 2016, VERSION 2

LECTURE NOTES ON QUANTUM CHAOS COURSE , MIT, APRIL 2016, VERSION 2 LECTUE NOTES ON QUANTUM CHAOS COUSE 18.158, MIT, APIL 2016, VESION 2 SEMYON DYATLOV Abstract. We give an overview of te quantum ergodicity result. 1. Quantum ergodicity in te pysical space 1.1. Concentration

More information

Digital Filter Structures

Digital Filter Structures Digital Filter Structures Te convolution sum description of an LTI discrete-time system can, in principle, be used to implement te system For an IIR finite-dimensional system tis approac is not practical

More information

Some Error Estimates for the Finite Volume Element Method for a Parabolic Problem

Some Error Estimates for the Finite Volume Element Method for a Parabolic Problem Computational Metods in Applied Matematics Vol. 13 (213), No. 3, pp. 251 279 c 213 Institute of Matematics, NAS of Belarus Doi: 1.1515/cmam-212-6 Some Error Estimates for te Finite Volume Element Metod

More information

DIGRAPHS FROM POWERS MODULO p

DIGRAPHS FROM POWERS MODULO p DIGRAPHS FROM POWERS MODULO p Caroline Luceta Box 111 GCC, 100 Campus Drive, Grove City PA 1617 USA Eli Miller PO Box 410, Sumneytown, PA 18084 USA Clifford Reiter Department of Matematics, Lafayette College,

More information

Recall from our discussion of continuity in lecture a function is continuous at a point x = a if and only if

Recall from our discussion of continuity in lecture a function is continuous at a point x = a if and only if Computational Aspects of its. Keeping te simple simple. Recall by elementary functions we mean :Polynomials (including linear and quadratic equations) Eponentials Logaritms Trig Functions Rational Functions

More information

Diffraction. S.M.Lea. Fall 1998

Diffraction. S.M.Lea. Fall 1998 Diffraction.M.Lea Fall 1998 Diffraction occurs wen EM waves approac an aperture (or an obstacle) wit dimension d > λ. We sall refer to te region containing te source of te waves as region I and te region

More information

Stationary Gaussian Markov Processes As Limits of Stationary Autoregressive Time Series

Stationary Gaussian Markov Processes As Limits of Stationary Autoregressive Time Series Stationary Gaussian Markov Processes As Limits of Stationary Autoregressive Time Series Lawrence D. Brown, Pilip A. Ernst, Larry Sepp, and Robert Wolpert August 27, 2015 Abstract We consider te class,

More information

POINTWISE BOUNDS ON QUASIMODES OF SEMICLASSICAL SCHRÖDINGER OPERATORS IN DIMENSION TWO

POINTWISE BOUNDS ON QUASIMODES OF SEMICLASSICAL SCHRÖDINGER OPERATORS IN DIMENSION TWO POINTWISE BOUNDS ON QUASIMODES OF SEMICLASSICAL SCHRÖDINGER OPERATORS IN DIMENSION TWO HART F. SMITH AND MACIEJ ZWORSKI Abstract. We prove optimal pointwise bounds on quasimodes of semiclassical Schrödinger

More information

Andrea Braides, Anneliese Defranceschi and Enrico Vitali. Introduction

Andrea Braides, Anneliese Defranceschi and Enrico Vitali. Introduction HOMOGENIZATION OF FREE DISCONTINUITY PROBLEMS Andrea Braides, Anneliese Defrancesci and Enrico Vitali Introduction Following Griffit s teory, yperelastic brittle media subject to fracture can be modeled

More information

Mathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative

Mathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x

More information

3 Parabolic Differential Equations

3 Parabolic Differential Equations 3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded

More information

Combining functions: algebraic methods

Combining functions: algebraic methods Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)

More information

MIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp

MIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a

More information

Approximation of the Viability Kernel

Approximation of the Viability Kernel Approximation of te Viability Kernel Patrick Saint-Pierre CEREMADE, Université Paris-Daupine Place du Marécal de Lattre de Tassigny 75775 Paris cedex 16 26 october 1990 Abstract We study recursive inclusions

More information

A h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation

A h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt

More information

Mass Lumping for Constant Density Acoustics

Mass Lumping for Constant Density Acoustics Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial

More information

Chapter Seven The Quantum Mechanical Simple Harmonic Oscillator

Chapter Seven The Quantum Mechanical Simple Harmonic Oscillator Capter Seven Te Quantum Mecanical Simple Harmonic Oscillator Introduction Te potential energy function for a classical, simple armonic oscillator is given by ZÐBÑ œ 5B were 5 is te spring constant. Suc

More information

The Schrödinger propagator for scattering metrics

The Schrödinger propagator for scattering metrics The Schrödinger propagator for scattering metrics Andrew Hassell (Australian National University) joint work with Jared Wunsch (Northwestern) MSRI, May 5-9, 2003 http://arxiv.org/math.ap/0301341 1 Schrödinger

More information

. If lim. x 2 x 1. f(x+h) f(x)

. If lim. x 2 x 1. f(x+h) f(x) Review of Differential Calculus Wen te value of one variable y is uniquely determined by te value of anoter variable x, ten te relationsip between x and y is described by a function f tat assigns a value

More information

THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein

THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein Worksop on Transforms and Filter Banks (WTFB),Brandenburg, Germany, Marc 999 THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS L. Trautmann, R. Rabenstein Lerstul

More information

Linear Algebra and its Applications

Linear Algebra and its Applications Linear Algebra and its Applications 466 (15) 1 116 Nonlinear Analysis 16 (15) 17 Contents lists available at ScienceDirect Contents lists available at ScienceDirect Linear Algebra and its Applications

More information

0.1 Differentiation Rules

0.1 Differentiation Rules 0.1 Differentiation Rules From our previous work we ve seen tat it can be quite a task to calculate te erivative of an arbitrary function. Just working wit a secon-orer polynomial tings get pretty complicate

More information

Analysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem

Analysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem Analysis of A Continuous inite Element Metod for Hcurl, div)-elliptic Interface Problem Huoyuan Duan, Ping Lin, and Roger C. E. Tan Abstract In tis paper, we develop a continuous finite element metod for

More information

An approximation method using approximate approximations

An approximation method using approximate approximations Applicable Analysis: An International Journal Vol. 00, No. 00, September 2005, 1 13 An approximation metod using approximate approximations FRANK MÜLLER and WERNER VARNHORN, University of Kassel, Germany,

More information

arxiv:math/ v1 [math.ca] 1 Oct 2003

arxiv:math/ v1 [math.ca] 1 Oct 2003 arxiv:mat/0310017v1 [mat.ca] 1 Oct 2003 Cange of Variable for Multi-dimensional Integral 4 Marc 2003 Isidore Fleiscer Abstract Te cange of variable teorem is proved under te sole ypotesis of differentiability

More information

Chapter 2. Limits and Continuity 16( ) 16( 9) = = 001. Section 2.1 Rates of Change and Limits (pp ) Quick Review 2.1

Chapter 2. Limits and Continuity 16( ) 16( 9) = = 001. Section 2.1 Rates of Change and Limits (pp ) Quick Review 2.1 Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 969) Quick Review..... f ( ) ( ) ( ) 0 ( ) f ( ) f ( ) sin π sin π 0 f ( ). < < < 6. < c c < < c 7. < < < < < 8. 9. 0. c < d d < c

More information

Minimal surfaces of revolution

Minimal surfaces of revolution 5 April 013 Minimal surfaces of revolution Maggie Miller 1 Introduction In tis paper, we will prove tat all non-planar minimal surfaces of revolution can be generated by functions of te form f = 1 C cos(cx),

More information

Math 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006

Math 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2

More information

arxiv: v1 [math.pr] 28 Dec 2018

arxiv: v1 [math.pr] 28 Dec 2018 Approximating Sepp s constants for te Slepian process Jack Noonan a, Anatoly Zigljavsky a, a Scool of Matematics, Cardiff University, Cardiff, CF4 4AG, UK arxiv:8.0v [mat.pr] 8 Dec 08 Abstract Slepian

More information

Order of Accuracy. ũ h u Ch p, (1)

Order of Accuracy. ũ h u Ch p, (1) Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical

More information

University Mathematics 2

University Mathematics 2 University Matematics 2 1 Differentiability In tis section, we discuss te differentiability of functions. Definition 1.1 Differentiable function). Let f) be a function. We say tat f is differentiable at

More information

Stability properties of a family of chock capturing methods for hyperbolic conservation laws

Stability properties of a family of chock capturing methods for hyperbolic conservation laws Proceedings of te 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August 0-, 005 (pp48-5) Stability properties of a family of cock capturing metods for yperbolic conservation

More information

Symmetry Labeling of Molecular Energies

Symmetry Labeling of Molecular Energies Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry

More information

Semigroups of Operators

Semigroups of Operators Lecture 11 Semigroups of Operators In tis Lecture we gater a few notions on one-parameter semigroups of linear operators, confining to te essential tools tat are needed in te sequel. As usual, X is a real

More information

Functions of the Complex Variable z

Functions of the Complex Variable z Capter 2 Functions of te Complex Variable z Introduction We wis to examine te notion of a function of z were z is a complex variable. To be sure, a complex variable can be viewed as noting but a pair of

More information

A NON HOMOGENEOUS EXTRA TERM FOR THE LIMIT OF DIRICHLET PROBLEMS IN PERFORATED DOMAINS

A NON HOMOGENEOUS EXTRA TERM FOR THE LIMIT OF DIRICHLET PROBLEMS IN PERFORATED DOMAINS J. Casado Diaz & A. Garroni A NON HOMOGENEOUS EXTRA TERM FOR THE LIMIT OF DIRICHLET PROBLEMS IN PERFORATED DOMAINS J. CASADO DIAZ A. GARRONI Abstract: We study te asymptotic beaviour of Diriclet problems

More information

New Streamfunction Approach for Magnetohydrodynamics

New Streamfunction Approach for Magnetohydrodynamics New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite

More information