Numerical Analysis for Nonlinear Eigenvalue Problems
|
|
- Frank Miller
- 5 years ago
- Views:
Transcription
1 Numerical Analysis for Nonlinear Eigenvalue Problems D. Bindel Department of Computer Science Cornell University 14 Sep 29
2 Outline Nonlinear eigenvalue problems Resonances via nonlinear eigenproblems Sensitivity and backward error analysis Resonances via Perfectly Matched Layers Perturbation bounds
3 Outline Nonlinear eigenvalue problems Resonances via nonlinear eigenproblems Sensitivity and backward error analysis Resonances via Perfectly Matched Layers Perturbation bounds
4 Linear eigenvalue problems Consider the linear system of differential equations y Ay = f, y() =. Taking the Laplace transform gives (s A)Y (s) = F(s). Singularities of s A correspond to solutions y(t) = e λt v to the unforced equation, where (A λi)v =.
5 Damping and delay Free vibrations correspond to nonlinear eigenvalues in more complicated settings. With simple damping: For systems with delay Mu + Bu + Ku = f. (s 2 M + sb + K )v =. u (t) Au(t) Bu(t τ) =. (s A e τs B)v =. General picture: study asymptotic free vibrations in terms of a nonlinear eigenvalue problem A(s)v = where A : C C n n is meromorphic.
6 Spectral Schur complements Nonlinear eigenvalue problems even arise from linear problems: [ ] A11 A A = 12. A 21 A 22 The spectral Schur complement is the inverse of a piece of the resolvent R(z) = (A zi) 1 : S(z) = (R 11 (z)) 1 = A 11 zi A 12 (A 22 zi) 1 A 21. Can use to reduce a large linear eigenvalue problem to a smaller nonlinear eigenvalue problem. Bounds on (A 22 zi) 1 lead to information about spectrum of A. Think of this as looking at a (block) transfer function.
7 Resonances and nonlinear eigenvalues For problems on unbounded domains: Can still have ordinary eigenvalues (bound states) Also have continuous spectrum (scattering states) Explore the continuous spectrum with a scattering matrix S(ω) relating incoming to outgoing waves S(ω) admits analytic continuation from R into C Places where (continued) S(ω) is singular are resonances Read off information about scattering from resonances
8 Example: Resonances and transmission ( d 2 ) dx 2 + V (x) ψ =, V (x) = 5 χ [,L] (x). T (E) imag(e) real(e)
9 Outline Nonlinear eigenvalue problems Resonances via nonlinear eigenproblems Sensitivity and backward error analysis Resonances via Perfectly Matched Layers Perturbation bounds
10 Simple 1D Problem Consider 1D Schrödinger (V nice, supp(v ) [a, b]): Hψ = ( d 2 ) dx 2 + V (x) ψ = Eψ. H self-adjoint with discrete spectrum on E <, continuous spectrum on E. Continuous spectrum is a branch cut for the resolvent. (Think χ(h E) 1 χ, χ a smooth cutoff, χ([a, b]) = 1. One moral analogue of a spectral Schur complement.) Second-sheet poles of the resolvent are resonances. Correspond to trapping, quasi-stable states.
11 Simple 1D Problem Consider 1D Schrödinger: ( d 2 ) dx 2 + V (x) ψ = Eψ. How do we: 1. Quickly compute resonances (nice enough V )? 2. Make sure the computations are correct?
12 Simple 1D Problem Consider 1D Schrödinger: If supp(v ) [a, b], write ( d 2 ) dx 2 + V (x) ψ = Eψ. ( d 2 ) dx 2 + V (x) k 2 ψ =, x (a, b) ( ) d dx ik ψ =, x = b ( ) d dx + ik ψ =, x = a E = k 2, Im k for eigenvalues, Im k < for resonances.
13 Pseudospectral Discretization Sample ψ at Chebyshev nodes and approximate dψ/dx by differentiating the interpolant: ( D 2 + V (x) k 2) ψ =, x (a, b) (D ik) ψ =, x = b (D + ik) ψ =, x = a Now linearize (introduce auxiliary variable φ = kψ) to get an ordinary generalized eigenvalue problem.
14 Is it that easy? 2 1 All eigenvalues Checked eigenvalues
15 Outline Nonlinear eigenvalue problems Resonances via nonlinear eigenproblems Sensitivity and backward error analysis Resonances via Perfectly Matched Layers Perturbation bounds
16 Backward Error Analysis 1. If ( ˆψ, Ê) is a numerical solution with above scheme, then there is some ˆV s.t. for x (a, b), (H ˆV Ê) ˆψ = ( d 2 ) dx 2 + ˆV (x) Ê ˆψ = together with corresponding radiation conditions. 2. Estimate ˆV explicitly by remapping residual to finer mesh 3. Original problem is a perturbation of computed problem. 4. Use first-order perturbation theory to correct Ê. Useful to take a variational approach.
17 More General Picture Consider Schrödinger with compactly supported V in R d. For E in resolvent set on appropriate Riemann surface, seek (H V E)ψ = f on Ω ψ B(E)ψ n = on Γ where B(E) is the Dirichlet-to-Neumann map on Ω. Solutions are stationary points for I(ψ) = 1 ( ) ( ψ) T ( ψ) + ψ(v E)ψ dω + 2 Ω 1 ψb(e)ψ dγ ψf dω. 2 Γ Ω
18 Variational Formulation Check variational formulation: I(ψ) = 1 ( ) ( ψ) T ( ψ) + ψ(v E)ψ dω 2 Ω 1 ψb(e)ψ dγ ψf dω. 2 Γ Use symmetry of form (note Γ φb(e)ψ = Γ ψb(e)φ) + integration by parts: δi(ψ) = δψ ( ψ + (V E)ψ f ) dω + Ω ( ) ψ δψ n B(E) ψ dγ. Γ Ω
19 Rayleigh Quotient Analogue Now define a residual for an approximate eigenpair: ( ) r(ψ, E) = ( ψ) T ( ψ) + ψ(v E)ψ ψb(e)ψ. Ω Take variations and use symmetry of B: δr(ψ, E) = 2 δψ [( + V E)ψ] + Ω [ ] ψ 2 δψ Γ n B(E)ψ + [ ] δe ψ 2 ψb (E)ψ Ω Γ For an eigenpair or resonance, r(ψ, E) = and δr(ψ, E) =. Γ
20 Rayleigh Quotient Analogue We now implicitly define a differentiable function Ẽ(φ) in the neighborhood of an eigenpair (ψ, E), with r(φ, E(φ)) = and Ẽ(ψ) = E. Such a function should exist if ψ 2 ψb (E)ψ Ω Stationary precisely when (ψ, E) an eigenpair. Γ
21 Sensitivity Now assume δv a compactly-supported perturbation, and look at effect of δv on Rayleigh quotient analogue. Gives that isolated eigenvalues change like Ω δv ψ2 δe = Ω ψ2 Γ ψb (E)ψ Can also write in terms of a residual for ψ as a solution for the potential V + δv : Ω ψ( + (V + δv ) E)ψ δe = Ω ψ2. Γ ψb (E)ψ
22 Backward Error Analysis Revisited 1. Compute approximate solution ( ˆψ, Ê). 2. Map ˆψ to high-resolution quadrature grid to evaluate Ω δe = ˆψ( + V Ê) ˆψ Ω ˆψ 2 ˆψB Γ (Ê) ˆψ. 3. If δe large, discard Ê as spurious; otherwise, accept E Ê + δe.
23 Some Computational Issues In general, using the domain equation + DtN map to find resonances is problematic because: 1. The DtN map is nonlocal, expensive to work with computationally. 2. The Green s function (and hence the DtN map) are hard to compute for some problems I care about (e.g. elastic half space problems). 3. Nonlinear eigenvalue problems are trickier than linear problems to solve.
24 Outline Nonlinear eigenvalue problems Resonances via nonlinear eigenproblems Sensitivity and backward error analysis Resonances via Perfectly Matched Layers Perturbation bounds
25 Perfectly Matched Layers For scattering computations / resonance computations, need an outgoing BC. We use perfectly matched layers: Complex coordinate transformation Generates a perfectly matched absorbing layer Rotates the continuous spectrum to reveal resonances Idea works with general linear wave equations Electromagnetics (Berengér, 1994) Quantum mechanics exterior complex scaling (Simon, 1979 originally used to define resonances) Elasticity in standard finite element framework (Basu and Chopra, 23)
26 Model Problem Domain: x [, ) Frequency-domain equation: d 2 û dx 2 + k 2 û = Solution: û = c out e ikx + c in e ikx
27 Model with Perfectly Matched Layer σ Regular domain x Transformed domain d x dx = λ(x) where λ(s) = 1 iσ(s) d 2 û d x 2 + k 2 û = û = c out e ik x ik x + c in e
28 Model with Perfectly Matched Layer σ Regular domain x Transformed domain d x dx = λ(x) where λ(s) = 1 iσ(s), 1 d λ dx ( 1 λ ) dû + k 2 û = dx û = c out e ikx kσ(x) + c in e ikx+kσ(x) Σ(x) = x σ(s) ds
29 Model with Perfectly Matched Layer σ Regular domain x Transformed domain If solution clamped at x = L then c in c out = O(e kγ ) where γ = Σ(L) = L σ(s) ds
30 Model Problem Illustrated 1 Outgoing exp( i x) 1 Incoming exp(i x) Transformed coordinate Re( x)
31 Model Problem Illustrated 1 Outgoing exp( i x) 3 Incoming exp(i x) Transformed coordinate Re( x)
32 Model Problem Illustrated 1 Outgoing exp( i x) 6 Incoming exp(i x) Transformed coordinate Re( x)
33 Model Problem Illustrated 1 Outgoing exp( i x) 15 Incoming exp(i x) Transformed coordinate Re( x)
34 Model Problem Illustrated 1 Outgoing exp( i x) 4 Incoming exp(i x) Transformed coordinate Re( x)
35 Model Problem Illustrated 1 Outgoing exp( i x) 1 Incoming exp(i x) Transformed coordinate Re( x)
36 Finite/Spectral Element Implementation ξ 2 x 2 x 2 x(ξ) x(x) Ω ξ 1 Ω e Ωe x 1 x 1 Combine PML and isoparametric mappings k e ij = ( N i ) T D( N j ) J dω Ω m e ij = ρn i N j J dω Ω Matrices are complex symmetric
37 DtN Approximation in PML Earlier (DtN) form - Neumann form = DtN term Γ ψbψ. Eliminate PML dofs - Neumann form DtN term Approximation is rational in k, good locally Would like to understand how good approximation to DtN + some form of stability leads to error bounds.
38 12 DtN Approximation in PML log 1 B(k) PML B(k) on circle of radius 3 in R 2. Order 3 spectral elements, PML goes [3, 4]
39 12 DtN Approximation in PML Axisymmetric resonances for ring barrier for r [1, 2]. Stars for small residual; circles for spurious resonance
40 Relation of PML to DtN Approach Numerically eliminate variables in PML = local (in k or E) rational approximation to DtN-like condition. Could also form explicit Padé approximation. Approximation is local in k where can we guarantee (for example) that we have approximated all resonances?
41 Outline Nonlinear eigenvalue problems Resonances via nonlinear eigenproblems Sensitivity and backward error analysis Resonances via Perfectly Matched Layers Perturbation bounds
42 General NEP Picture A : C C n n analytic in Ω Λ(A) := {z C : A(z) singular} Λ ɛ (A) := {z C : A(z) 1 ɛ 1 } Resonance calculation with DtN map or with eliminated PML as an example. Λ(A) and Λ ɛ (A) describe asymptotics, transients of some linear differential or difference equation. Lots of function theoretic proofs from analyzing ordinary eigenvalue problems carry over without change.
43 Counting Eigenvalues If A nonsingular on Γ, analytic inside, count eigs inside by W Γ (det(a)) = 1 d ln det(a(z)) dz 2πi Γ dz ( ) 1 = tr A(z) 1 A (z) dz 2πi Suppose E also analytic inside Γ. By continuity, W Γ (det(a)) = W Γ (det(a + se)) for s in neighborhood of such that A + se remains nonsingular on Γ. Γ
44 Function Theoretic Perturbation Recipe Winding number counts give continuity of eigenvalues = Should consider eigenvalues of A + se for s 1: Analyticity of A and E + Matrix nonsingularity test for A + se = Inclusion region for Λ(A + E) + Eigenvalue counts for connected components of region
45 Example: Matrix Rouché A 1 (z)e(z) < 1 on Γ = same eigenvalue count in Γ Proof: A 1 (z)e(z) < 1 = A(z) + se(z) invertible for s 1. (Gohberg and Sigal proved a more general version in 1971.)
46 Example: Nonlinear Gershgorin Define Then G i = z : a ii(z) < j i a ij (z) 1. Λ(A) i G i 2. Connected component m i=1 G i contains m eigs (if bounded and disjoint from Ω) Proof: Write A = D + F where D = diag(a). D + sf is diagonally dominant (so invertible) off i G i.
47 Example: Pseudospectral containment Define D = {z : E(z) < ɛ}. Then 1. Λ(A + E) Λ ɛ (A) D C 2. A bounded component of Λ ɛ (A) strictly inside D contains the same number of eigs of A and A + E.
48 Application: Lattice Schrödinger Consider the discrete analogue to Schrödinger s equation: Hψ = ( T + V )ψ = Eψ where (Hψ) k = ψ k 1 + 2ψ k ψ k+1 + V k ψ k. Assume V k = for k and k L. May be complex. Want to relate the spectrum for two variants: 1. Non-negative integers: ψ = and ψ l 2 2. Bounded: ψ k = for k = and k L + N.
49 Application: Lattice Schrödinger For V 1 =.1i and V 2 = 4, N = 2.
50 Application: Lattice Schrödinger For V 1 =.1i and V 2 = 4, N = 2.
51 Spectral Schur Complement Write H in either case as [ T11 + V H = 11 e L e1 T ] e L e1 T T 22 Then Λ(H) Λ( T 22 ) c = Λ(S), where ) S(z) = ( T 11 + V 11 ) zi (e 1 T ( T 22 zi) 1 e 1 e L el T Write S (N) (z) and S ( ) (z) for bounded and unbounded cases.
52 Spectral Schur Complement For z [, 4], choose ξ 2 (2 z)ξ + 1 =, ξ < 1. Then S ( ) (z) = ( T 11 + V 11 ) zi ξe L el T ( ) 1 ξ S (N) 2N (z) = ( T 11 + V 11 ) zi ξ 1 ξ 2(N+1) e L el T Convenient to write z = 2 ξ ξ 1, use ξ as primary variable.
53 Error Bounds Find S ( ) S (N) ɛ if ξ < ( ) log(3ɛ 1 ) = 1 O 2N + 1 ( log(ɛ 1 ) N ). Therefore, eigenvalues in bounded case (in ξ plane) either 1. Are within O(log(ɛ 1 )/N) of circle (continuous spectrum) 2. Are in Λ ɛ (S ( ) ). Get exponential convergence to discrete spectrum, linear convergence to continuous spectrum.
54 Error Estimate If S ( ) has an isolated eigenvalue at γ, then S (N) asymptotically has eigenvalues γ (N) γ with γ (N) γ = γ 2N w e L el T v L (1 γ 2 )w v w e L el T v + O(γ2N+1 ) where S ( ) (γ)v = and w S ( ) (γ) =.
55 Similar Applications Linear stability analysis for traveling waves. Bounds on distance to instability via subspace projections. Estimates of damping in MEMS resonators.
56 Conclusions Nonlinear eigenproblems tell us interesting information about dynamics. Analytic structure of the eigenproblem is key to error analysis. Variational characterization gives easy first-order perturbation theory Also get analogues to standard perturbation bounds (Rouché, Gerschgorin, pseudospectral) Get interesting estimates via approximation of spectral Schur complements
Bounds and Error Estimates for Nonlinear Eigenvalue Problems
Bounds and Error Estimates for Nonlinear Eigenvalue Problems D. Bindel Courant Institute for Mathematical Sciences New York Univerity 8 Oct 28 Outline Resonances via nonlinear eigenproblems Sensitivity
More informationNonlinear Eigenvalue Problems: Theory and Applications
Nonlinear Eigenvalue Problems: Theory and Applications David Bindel 1 Department of Computer Science Cornell University 12 January 217 1 Joint work with Amanda Hood Berkeley 1 / 37 The Computational Science
More informationComputer Aided Design of Micro-Electro-Mechanical Systems
Computer Aided Design of Micro-Electro-Mechanical Systems From Resonance Poles to Dick Tracy Watches D. Bindel Courant Institute for Mathematical Sciences New York University Columbia University, 20 March
More informationEigenvalues, Resonance Poles, and Damping in MEMS
Eigenvalues,, and Damping in D. Bindel Computer Science Division Department of EECS University of California, Berkeley Matrix computations seminar, 19 Apr 26 Outline 1 2 3 4 Outline 1 2 3 4 Resonating
More informationComputer Aided Design of Micro-Electro-Mechanical Systems
Computer Aided Design of Micro-Electro-Mechanical Systems From Energy Losses to Dick Tracy Watches D. Bindel Courant Institute for Mathematical Sciences New York University Stanford University, 12 Nov
More informationComputer Aided Design of Micro-Electro-Mechanical Systems
Computer Aided Design of Micro-Electro-Mechanical Systems From Energy Losses to Dick Tracy Watches D. Bindel Courant Institute for Mathematical Sciences New York University Temple University, 7 Nov 2007
More informationA new method for the solution of scattering problems
A new method for the solution of scattering problems Thorsten Hohage, Frank Schmidt and Lin Zschiedrich Konrad-Zuse-Zentrum Berlin, hohage@zibde * after February 22: University of Göttingen Abstract We
More informationComputer Aided Design of Micro-Electro-Mechanical Systems
Computer Aided Design of Micro-Electro-Mechanical Systems From Eigenvalues to Devices David Bindel Department of Computer Science Cornell University Fudan University, 13 Dec 2012 (Cornell University) Fudan
More informationThe Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator
The Pole Condition: A Padé Approximation of the Dirichlet to Neumann Operator Martin J. Gander and Achim Schädle Mathematics Section, University of Geneva, CH-, Geneva, Switzerland, Martin.gander@unige.ch
More informationComputer-Aided Design for Micro-Electro-Mechanical Systems
Computer-Aided Design for Micro-Electro-Mechanical Systems Eigenvalues, Energy, and Dick Tracy Watches D. Bindel Computer Science Division Department of EECS University of California, Berkeley Caltech,
More informationRiesz bases of Floquet modes in semi-infinite periodic waveguides and implications
Riesz bases of Floquet modes in semi-infinite periodic waveguides and implications Thorsten Hohage joint work with Sofiane Soussi Institut für Numerische und Angewandte Mathematik Georg-August-Universität
More informationStructure-preserving model reduction for MEMS
Structure-preserving model reduction for MEMS David Bindel Department of Computer Science Cornell University SIAM CSE Meeting, 1 Mar 2011 Collaborators Sunil Bhave Emmanuel Quévy Zhaojun Bai Tsuyoshi Koyama
More informationNonlinear Eigenvalue Problems: An Introduction
Nonlinear Eigenvalue Problems: An Introduction Cedric Effenberger Seminar for Applied Mathematics ETH Zurich Pro*Doc Workshop Disentis, August 18 21, 2010 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction
More informationADI iterations for. general elliptic problems. John Strain Mathematics Department UC Berkeley July 2013
ADI iterations for general elliptic problems John Strain Mathematics Department UC Berkeley July 2013 1 OVERVIEW Classical alternating direction implicit (ADI) iteration Essentially optimal in simple domains
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationLecture 2: Computing functions of dense matrices
Lecture 2: Computing functions of dense matrices Paola Boito and Federico Poloni Università di Pisa Pisa - Hokkaido - Roma2 Summer School Pisa, August 27 - September 8, 2018 Introduction In this lecture
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More information1 Sylvester equations
1 Sylvester equations Notes for 2016-11-02 The Sylvester equation (or the special case of the Lyapunov equation) is a matrix equation of the form AX + XB = C where A R m m, B R n n, B R m n, are known,
More informationLecture 4: Resonant Scattering
Lecture 4: Resonant Scattering Sep 16, 2008 Fall 2008 8.513 Quantum Transport Analyticity properties of S-matrix Poles and zeros in a complex plane Isolated resonances; Breit-Wigner theory Quasi-stationary
More informationBindel, Fall 2016 Matrix Computations (CS 6210) Notes for
1 Iteration basics Notes for 2016-11-07 An iterative solver for Ax = b is produces a sequence of approximations x (k) x. We always stop after finitely many steps, based on some convergence criterion, e.g.
More informationNumerical and semi-analytical structure-preserving model reduction for MEMS
Numerical and semi-analytical structure-preserving model reduction for MEMS David Bindel ENUMATH 07, 10 Sep 2007 Collaborators Tsuyoshi Koyama Sanjay Govindjee Sunil Bhave Emmanuel Quévy Zhaojun Bai Resonant
More information1.1 A Scattering Experiment
1 Transfer Matrix In this chapter we introduce and discuss a mathematical method for the analysis of the wave propagation in one-dimensional systems. The method uses the transfer matrix and is commonly
More informationEIGENVALUE PROBLEMS. Background on eigenvalues/ eigenvectors / decompositions. Perturbation analysis, condition numbers..
EIGENVALUE PROBLEMS Background on eigenvalues/ eigenvectors / decompositions Perturbation analysis, condition numbers.. Power method The QR algorithm Practical QR algorithms: use of Hessenberg form and
More informationSolving Time-Harmonic Scattering Problems by the Ultra Weak Variational Formulation
Introduction Solving Time-Harmonic Scattering Problems by the Ultra Weak Variational Formulation Plane waves as basis functions Peter Monk 1 Tomi Huttunen 2 1 Department of Mathematical Sciences University
More informationMatrix Algorithms. Volume II: Eigensystems. G. W. Stewart H1HJ1L. University of Maryland College Park, Maryland
Matrix Algorithms Volume II: Eigensystems G. W. Stewart University of Maryland College Park, Maryland H1HJ1L Society for Industrial and Applied Mathematics Philadelphia CONTENTS Algorithms Preface xv xvii
More informationApplied Nuclear Physics (Fall 2006) Lecture 2 (9/11/06) Schrödinger Wave Equation
22.101 Applied Nuclear Physics (Fall 2006) Lecture 2 (9/11/06) Schrödinger Wave Equation References -- R. M. Eisberg, Fundamentals of Modern Physics (Wiley & Sons, New York, 1961). R. L. Liboff, Introductory
More informationDelocalization for Schrödinger operators with random Dirac masses
Delocalization for Schrödinger operators with random Dirac masses CEMPI Scientific Day Lille, 10 February 2017 Disordered Systems E.g. box containing N nuclei Motion of electron in this system High temperature
More informationEigenvalue Problems in Surface Acoustic Wave Filter Simulations
Eigenvalue Problems in Surface Acoustic Wave Filter Simulations S. Zaglmayr, J. Schöberl, U. Langer FWF Start Project Y-192 3D hp-finite Elements: Fast Solvers and Adaptivity JKU + RICAM, Linz in collaboration
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
More informationImplicitly Defined High-Order Operator Splittings for Parabolic and Hyperbolic Variable-Coefficient PDE Using Modified Moments
Implicitly Defined High-Order Operator Splittings for Parabolic and Hyperbolic Variable-Coefficient PDE Using Modified Moments James V. Lambers September 24, 2008 Abstract This paper presents a reformulation
More informationThe Canonical Sheaf. Stefano Filipazzi. September 14, 2015
The Canonical Sheaf Stefano Filipazzi September 14, 015 These notes are supposed to be a handout for the student seminar in algebraic geometry at the University of Utah. In this seminar, we will go over
More informationAn optimized perfectly matched layer for the Schrödinger equation
An optimized perfectly matched layer for the Schrödinger equation Anna Nissen Gunilla Kreiss June 6, 009 Abstract A perfectly matched layer (PML) for the Schrödinger equation using a modal ansatz is presented.
More informationQMI PRELIM Problem 1. All problems have the same point value. If a problem is divided in parts, each part has equal value. Show all your work.
QMI PRELIM 013 All problems have the same point value. If a problem is divided in parts, each part has equal value. Show all your work. Problem 1 L = r p, p = i h ( ) (a) Show that L z = i h y x ; (cyclic
More informationEstimates for the resolvent and spectral gaps for non self-adjoint operators. University Bordeaux
for the resolvent and spectral gaps for non self-adjoint operators 1 / 29 Estimates for the resolvent and spectral gaps for non self-adjoint operators Vesselin Petkov University Bordeaux Mathematics Days
More informationABSOLUTELY CONTINUOUS SPECTRUM OF A POLYHARMONIC OPERATOR WITH A LIMIT PERIODIC POTENTIAL IN DIMENSION TWO. V (x) = V r (x); (2)
ABSOLUTELY CONTINUOUS SPECTRUM OF A POLYHARMONIC OPERATOR WITH A LIMIT PERIODIC POTENTIAL IN DIMENSION TWO. YULIA KARPESHINA AND YOUNG-RAN LEE Abstract. We consider a polyharmonic operator H = ) l + V
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More information1 Infinite-Dimensional Vector Spaces
Theoretical Physics Notes 4: Linear Operators In this installment of the notes, we move from linear operators in a finitedimensional vector space (which can be represented as matrices) to linear operators
More informationPart II. Riemann Surfaces. Year
Part II Year 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2018 96 Paper 2, Section II 23F State the uniformisation theorem. List without proof the Riemann surfaces which are uniformised
More informationGeneralized Burgers equations and Miura Map in nonabelian ring. nonabelian rings as integrable systems.
Generalized Burgers equations and Miura Map in nonabelian rings as integrable systems. Sergey Leble Gdansk University of Technology 05.07.2015 Table of contents 1 Introduction: general remarks 2 Remainders
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationThe PML Method: Continuous and Semidiscrete Waves.
Intro Continuous Model. Finite difference. Remedies. The PML Method: Continuous and Semidiscrete Waves. 1 Enrique Zuazua 2 1 Laboratoire de Mathématiques de Versailles. 2 Universidad Autónoma, Madrid.
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More informationKrylov methods for the computation of matrix functions
Krylov methods for the computation of matrix functions Jitse Niesen (University of Leeds) in collaboration with Will Wright (Melbourne University) Heriot-Watt University, March 2010 Outline Definition
More informationDiagonalization of the Coupled-Mode System.
Diagonalization of the Coupled-Mode System. Marina Chugunova joint work with Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaborators: Mason A. Porter, California Institute
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationLOCALIZATION THEOREMS FOR NONLINEAR EIGENVALUE PROBLEMS
LOCALIZATION THEOREMS FOR NONLINEAR EIGENVALUE PROBLEMS DAVID BINDEL AND AMANDA HOOD Abstract. Let T : Ω C n n be a matrix-valued function that is analytic on some simplyconnected domain Ω C. A point λ
More informationQuantum Physics III (8.06) Spring 2016 Assignment 3
Quantum Physics III (8.6) Spring 6 Assignment 3 Readings Griffiths Chapter 9 on time-dependent perturbation theory Shankar Chapter 8 Cohen-Tannoudji, Chapter XIII. Problem Set 3. Semi-classical approximation
More informationLecture 10: Finite Differences for ODEs & Nonlinear Equations
Lecture 10: Finite Differences for ODEs & Nonlinear Equations J.K. Ryan@tudelft.nl WI3097TU Delft Institute of Applied Mathematics Delft University of Technology 21 November 2012 () Finite Differences
More informationx n+1 = x n f(x n) f (x n ), n 0.
1. Nonlinear Equations Given scalar equation, f(x) = 0, (a) Describe I) Newtons Method, II) Secant Method for approximating the solution. (b) State sufficient conditions for Newton and Secant to converge.
More informationA diamagnetic inequality for semigroup differences
A diamagnetic inequality for semigroup differences (Irvine, November 10 11, 2001) Barry Simon and 100DM 1 The integrated density of states (IDS) Schrödinger operator: H := H(V ) := 1 2 + V ω =: H(0, V
More informationSpectral approach to scattering from spheres (scalar waves)
Spectral approach to scattering from spheres (scalar waves) November 1, 2005 1 Introduction Consider the scalar wave equation [ 2 + ɛ(r)k 2 0]ψ(r) = S(r). (1) Here S(r) is the source which is usually zero
More informationCHAPTER 8 The Quantum Theory of Motion
I. Translational motion. CHAPTER 8 The Quantum Theory of Motion A. Single particle in free space, 1-D. 1. Schrodinger eqn H ψ = Eψ! 2 2m d 2 dx 2 ψ = Eψ ; no boundary conditions 2. General solution: ψ
More information16. Local theory of regular singular points and applications
16. Local theory of regular singular points and applications 265 16. Local theory of regular singular points and applications In this section we consider linear systems defined by the germs of meromorphic
More informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
More informationTime-dependent variational forms
Time-dependent variational forms Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Oct 30, 2015 PRELIMINARY VERSION
More informationFirst-order overdetermined systems. for elliptic problems. John Strain Mathematics Department UC Berkeley July 2012
First-order overdetermined systems for elliptic problems John Strain Mathematics Department UC Berkeley July 2012 1 OVERVIEW Convert elliptic problems to first-order overdetermined form Control error via
More informationAppendix B: The Transfer Matrix Method
Y D Chong (218) PH441: Quantum Mechanics III Appendix B: The Transfer Matrix Method The transfer matrix method is a numerical method for solving the 1D Schrödinger equation, and other similar equations
More informationLOCALIZATION THEOREMS FOR NONLINEAR EIGENVALUE PROBLEMS
LOCALIZATION THEOREMS FOR NONLINEAR EIGENVALUE PROBLEMS DAVID BINDEL AND AMANDA HOOD Abstract. Let T : Ω C n n be a matrix-valued function that is analytic on some simplyconnected domain Ω C. A point λ
More informationThe boundary state from open string fields. Yuji Okawa University of Tokyo, Komaba. March 9, 2009 at Nagoya
The boundary state from open string fields Yuji Okawa University of Tokyo, Komaba March 9, 2009 at Nagoya Based on arxiv:0810.1737 in collaboration with Kiermaier and Zwiebach (MIT) 1 1. Introduction Quantum
More informationOn the singularities of non-linear ODEs
On the singularities of non-linear ODEs Galina Filipuk Institute of Mathematics University of Warsaw G.Filipuk@mimuw.edu.pl Collaborators: R. Halburd (London), R. Vidunas (Tokyo), R. Kycia (Kraków) 1 Plan
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationChapter 18. Remarks on partial differential equations
Chapter 8. Remarks on partial differential equations If we try to analyze heat flow or vibration in a continuous system such as a building or an airplane, we arrive at a kind of infinite system of ordinary
More information1 Equal-time and Time-ordered Green Functions
1 Equal-time and Time-ordered Green Functions Predictions for observables in quantum field theories are made by computing expectation values of products of field operators, which are called Green functions
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationDmitry Pelinovsky Department of Mathematics, McMaster University, Canada
Spectrum of the linearized NLS problem Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaboration: Marina Chugunova (McMaster, Canada) Scipio Cuccagna (Modena and Reggio Emilia,
More informationGENERALIZED STABILITY OF THE TWO-LAYER MODEL
GENERALIZED STABILITY OF THE TWO-LAYER MODEL The simplest mid-latitude jet model supporting the baroclinic growth mechanism is the two-layer model The equations for the barotropic and baroclinic geostrophic
More informationRiemann Surfaces and Algebraic Curves
Riemann Surfaces and Algebraic Curves JWR Tuesday December 11, 2001, 9:03 AM We describe the relation between algebraic curves and Riemann surfaces. An elementary reference for this material is [1]. 1
More informationName: INSERT YOUR NAME HERE. Due to dropbox by 6pm PDT, Wednesday, December 14, 2011
AMath 584 Name: INSERT YOUR NAME HERE Take-home Final UWNetID: INSERT YOUR NETID Due to dropbox by 6pm PDT, Wednesday, December 14, 2011 The main part of the assignment (Problems 1 3) is worth 80 points.
More informationNotes by Zvi Rosen. Thanks to Alyssa Palfreyman for supplements.
Lecture: Hélène Barcelo Analytic Combinatorics ECCO 202, Bogotá Notes by Zvi Rosen. Thanks to Alyssa Palfreyman for supplements.. Tuesday, June 2, 202 Combinatorics is the study of finite structures that
More informationSolving PDEs Numerically on Manifolds with Arbitrary Spatial Topologies
Solving PDEs Numerically on Manifolds with Arbitrary Spatial Topologies Lee Lindblom Theoretical Astrophysics, Caltech Center for Astrophysics and Space Sciences, UC San Diego. Collaborators: Béla Szilágyi,
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More informationSynopsis of Numerical Linear Algebra
Synopsis of Numerical Linear Algebra Eric de Sturler Department of Mathematics, Virginia Tech sturler@vt.edu http://www.math.vt.edu/people/sturler Iterative Methods for Linear Systems: Basics to Research
More informationNOTES FOR CARDIFF LECTURES ON MICROLOCAL ANALYSIS
NOTES FOR CARDIFF LECTURES ON MICROLOCAL ANALYSIS JARED WUNSCH Note that these lectures overlap with Alex s to a degree, to ensure a smooth handoff between lecturers! Our notation is mostly, but not completely,
More informationComputer-Aided Design for Micro-Electro-Mechanical Systems
Computer-Aided Design for Micro-Electro-Mechanical Systems Eigenvalues, Energy, and Dick Tracy Watches D. Bindel Computer Science Division Department of EECS University of California, Berkeley Purdue,
More informationDeterminant lines and determinant line bundles
CHAPTER Determinant lines and determinant line bundles This appendix is an exposition of G. Segal s work sketched in [?] on determinant line bundles over the moduli spaces of Riemann surfaces with parametrized
More informationCollocation approximation of the monodromy operator of periodic, linear DDEs
p. Collocation approximation of the monodromy operator of periodic, linear DDEs Ed Bueler 1, Victoria Averina 2, and Eric Butcher 3 13 July, 2004 SIAM Annual Meeting 2004, Portland 1=Dept. of Math. Sci.,
More informationNumerical solution of the Schrödinger equation on unbounded domains
Numerical solution of the Schrödinger equation on unbounded domains Maike Schulte Institute for Numerical and Applied Mathematics Westfälische Wilhelms-Universität Münster Cetraro, September 6 OUTLINE
More informationWeighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,
More informationDispersive Equations and Hyperbolic Orbits
Dispersive Equations and Hyperbolic Orbits H. Christianson Department of Mathematics University of California, Berkeley 4/16/07 The Johns Hopkins University Outline 1 Introduction 3 Applications 2 Main
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationCoercivity of high-frequency scattering problems
Coercivity of high-frequency scattering problems Valery Smyshlyaev Department of Mathematics, University College London Joint work with: Euan Spence (Bath), Ilia Kamotski (UCL); Comm Pure Appl Math 2015.
More informationNumerical Simulation of Spin Dynamics
Numerical Simulation of Spin Dynamics Marie Kubinova MATH 789R: Advanced Numerical Linear Algebra Methods with Applications November 18, 2014 Introduction Discretization in time Computing the subpropagators
More informationAnalysis of the anomalous localized resonance
Hyundae Lee(Inha University, Korea) Joint work with Habib Ammari, Giulio Ciraolo, Hyeonbae Kang, Graeme Milton. UCI, June 22, 2012 A conference on inverse problems in honor of Gunther Uhlmann Outline Introduction
More informationQuantum Mechanics (Draft 2010 Nov.)
Quantum Mechanics (Draft 00 Nov) For a -dimensional simple harmonic quantum oscillator, V (x) = mω x, it is more convenient to describe the dynamics by dimensionless position parameter ρ = x/a (a = h )
More informationPhysics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension
Physics 221A Fall 1996 Notes 16 Bloch s Theorem and Band Structure in One Dimension In these notes we examine Bloch s theorem and band structure in problems with periodic potentials, as a part of our survey
More informationMath Introduction to Numerical Analysis - Class Notes. Fernando Guevara Vasquez. Version Date: January 17, 2012.
Math 5620 - Introduction to Numerical Analysis - Class Notes Fernando Guevara Vasquez Version 1990. Date: January 17, 2012. 3 Contents 1. Disclaimer 4 Chapter 1. Iterative methods for solving linear systems
More informationDefect-based a-posteriori error estimation for implicit ODEs and DAEs
1 / 24 Defect-based a-posteriori error estimation for implicit ODEs and DAEs W. Auzinger Institute for Analysis and Scientific Computing Vienna University of Technology Workshop on Innovative Integrators
More informationThe Eigenvalue Problem: Perturbation Theory
Jim Lambers MAT 610 Summer Session 2009-10 Lecture 13 Notes These notes correspond to Sections 7.2 and 8.1 in the text. The Eigenvalue Problem: Perturbation Theory The Unsymmetric Eigenvalue Problem Just
More informationScattering in one dimension
Scattering in one dimension Oleg Tchernyshyov Department of Physics and Astronomy, Johns Hopkins University I INTRODUCTION This writeup accompanies a numerical simulation of particle scattering in one
More informationINTEGRATION WORKSHOP 2003 COMPLEX ANALYSIS EXERCISES
INTEGRATION WORKSHOP 23 COMPLEX ANALYSIS EXERCISES DOUGLAS ULMER 1. Meromorphic functions on the Riemann sphere It s often useful to allow functions to take the value. This exercise outlines one way to
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationSummary of free theory: one particle state: vacuum state is annihilated by all a s: then, one particle state has normalization:
The LSZ reduction formula based on S-5 In order to describe scattering experiments we need to construct appropriate initial and final states and calculate scattering amplitude. Summary of free theory:
More informationEvans 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 informationPART 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 informationExplicit Examples of Strebel Differentials
Explicit Examples of Strebel Differentials arxiv:0910.475v [math.dg] 30 Oct 009 1 Introduction Philip Tynan November 14, 018 In this paper, we investigate Strebel differentials, which are a special class
More informationHerbert Stahl s proof of the BMV conjecture
Herbert Stahl s proof of the BMV conjecture A. Eremenko August 25, 203 Theorem. Let A and B be two n n Hermitian matrices, where B is positive semi-definite. Then the function is absolutely monotone, that
More information