Numerical solution of the Schrödinger equation on unbounded domains
|
|
- Lionel Boone
- 5 years ago
- Views:
Transcription
1 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
2 OUTLINE Transparent boundary conditions (TBC) OUTLINE
3 OUTLINE Transparent boundary conditions (TBC) Analytical transparent boundary conditions for the Schrödinger equation OUTLINE -A
4 OUTLINE Transparent boundary conditions (TBC) Analytical transparent boundary conditions for the Schrödinger equation Derivation of discrete transparent boundary conditions (DTBC) for the Schrödinger equation for a higher order compact 9-point scheme OUTLINE -B
5 OUTLINE Transparent boundary conditions (TBC) Analytical transparent boundary conditions for the Schrödinger equation Derivation of discrete transparent boundary conditions (DTBC) for the Schrödinger equation for a higher order compact 9-point scheme Approximation of DTBC OUTLINE -C
6 OUTLINE Transparent boundary conditions (TBC) Analytical transparent boundary conditions for the Schrödinger equation Derivation of discrete transparent boundary conditions (DTBC) for the Schrödinger equation for a higher order compact 9-point scheme Approximation of DTBC Simulation of quantum wave guides OUTLINE -D
7 OUTLINE Transparent boundary conditions (TBC) Analytical transparent boundary conditions for the Schrödinger equation Derivation of discrete transparent boundary conditions (DTBC) for the Schrödinger equation for a higher order compact 9-point scheme Approximation of DTBC Simulation of quantum wave guides Extension of the DTBC to (nearly) arbitrary potentials OUTLINE -E
8 OUTLINE Transparent boundary conditions (TBC) Analytical transparent boundary conditions for the Schrödinger equation Derivation of discrete transparent boundary conditions (DTBC) for the Schrödinger equation for a higher order compact 9-point scheme Approximation of DTBC Simulation of quantum wave guides Extension of the DTBC to (nearly) arbitrary potentials more realistic simulations OUTLINE -G
9 OUTLINE Transparent boundary conditions (TBC) Analytical transparent boundary conditions for the Schrödinger equation Derivation of discrete transparent boundary conditions (DTBC) for the Schrödinger equation for a higher order compact 9-point scheme Approximation of DTBC Simulation of quantum wave guides Extension of the DTBC to (nearly) arbitrary potentials more realistic simulations solving the Schrödinger equation on circular domains OUTLINE -H
10 TRANSPARENT BOUNDARY CONDITIONS (IN GENERAL) Γ R n ψ Ω Definition (TBC): Consider a given whole-space initial value problem (IVP) on R n and Ω R n, Γ = Ω. We are interested in the solution of the IVP on Ω. Therefore we need new artificial boundary conditions on Γ. We call these artificial BC transparent, if the solution of the IVBP on Ω corresponds to the whole-space solution of the IVP restricted on Ω. TRANSPARENT BOUNDARY CONDITIONS (IN GENERAL) 3
11 ANALYTICAL TBC FOR THE SCHRÖDINGER EQUATION IVP: time-dependent Schrödinger equation (here: D) i «t ψ(x,t) = + V (x, t) ψ(x,t), x R, t > m x ψ(x,) = ψ I (x) L (R) on a domain of interest Ω = {x R < x < X}. ANALYTICAL TBC FOR THE SCHRÖDINGER EQUATION 4
12 ANALYTICAL TBC FOR THE SCHRÖDINGER EQUATION IVP: time-dependent Schrödinger equation (here: D) i «t ψ(x,t) = + V (x, t) ψ(x,t), x R, t > m x ψ(x,) = ψ I (x) L (R) on a domain of interest Ω = {x R < x < X}. Assumptions: supp ψ I Ω potential V (., t) L (R), V (x,.) is piecewise continuous V constant on R\Ω (here: V (x, t) = for x and V (x, t) = V X for x X) Goal: Calculate the solution ψ(x,t) C on Ω with TBC at x = and x = X. ANALYTICAL TBC FOR THE SCHRÖDINGER EQUATION 4-A
13 DERIVATION OF TBC ψ V = const V = const G Ω G ψ I X x x G : x G : i ψ t = m ψ xx + V (x, t)ψ i ψ t = m ψ xx + V (x, t)ψ ψ(x, ) = ψ I (x) v(x, ) = ψ x (, t) = (T ψ)(,t) v(x, t) = Φ(t) t >, Φ() = ψ x (X, t) = (T X ψ)(x,t) lim v(x, t) = x v x (, t) = (T X Φ)(t) DERIVATION OF TBC 5
14 Laplace-transformation on the exterior domains: ˆv xx (x, s) + im s + iv X ˆv(x, s) = x > X ˆv(X, s) = ˆΦ(s) lim ˆv(x, s) = x ˆv x (X, s) = (T X Φ)(s) r solution: ˆv(x, s) = e i + im `s+ iv X (x X) ˆΦ(s) r r m (TX Φ)(s) = e iπ 4 + s + iv X ˆΦ(s) With the inverse Laplace-Transformation follows the analytical TBC ψ x (X, t) = r π V πm e i 4 e i X t d dt Z t ψ(x,τ)e i V X τ τ t dτ. [J. S. Papadakis (98)] DERIVATION OF TBC 6
15 Application of DTBC for the Schrödinger equation: Simulation of quantum transistors in quantum waveguides (with inhomogeneous DTBC for the D Schrödinger equation) Analyse steady states and transient behaviour y V control Y ψ inc Γ Ω Γ X x Y y ψ inc V control Γ Ω Γ ψ pw X x DERIVATION OF TBC 7
16 FORMER STRATEGIES: Discretization of the analytic TBCs with an numerical approximation of the convolution integral [e.g. B. Mayfield (989)] only conditionally stable, not transparent! FORMER STRATEGIES: 8
17 FORMER STRATEGIES: Discretization of the analytic TBCs with an numerical approximation of the convolution integral [e.g. B. Mayfield (989)] only conditionally stable, not transparent! Create a buffer zone Θ of the length d with a complex potential V (X) = W ia around the computational domain Ω with Dirichlet -BC at Θ and absorbing boundary conditions on Ω [e.g. L. Burgnies (997)] ψ= Θ Ω X Θ ψ= unconditionally stable, unphysical reflections at the boundary, huge numerical costs FORMER STRATEGIES: 8-A
18 SUCCESSFUL STRATEGIES Family of absorbing BCs (also for the non-linear Schrödinger equation, wave equation) [J. Szeftel (5)] discretize the whole space problem with an unconditionally stable scheme (e.g. Crank-Nicolson finite difference scheme) and calculate new discrete transparent boundary conditions for the full discretized Schrödinger equation DTBC: dx=.65, dt=e 5, t= ψ x [A. Arnold, M. Ehrhardt (since 995)] SUCCESSFUL STRATEGIES 9
19 DERIVATION OF DISCRETE TRB Discretize D Schrödinger equation: Crank-Nicolson scheme in time, t n = n t, n N, H := + V Hamilton-Operator, Ω = [,X] [, Y ] + ih t «ψ(x,y, t + t) = ih t «ψ(x,y, t) x + y «`ψn+ (x, y) + ψ n (x, y) = g n+ (x, y) `ψn+ (x, y) + ψ n (x, y) + Wψ n (x, y) compact 9-point scheme in space, x j = j x, y k = k y with j Z, k K DTBC at x = and x J = X = J x with J Z -BC at y = and y = Y = K y with K N DERIVATION OF DISCRETE TRB
20 9-point discretization scheme: (j+,k+) (j,k) (j+,k) D x + Dy + x + y Dx Dy = ψ n+ j,k I + x D x + y D y «h V n+ j,k ψ n+ j,k id + t ψ n j,k i with ψ n+ j,k = `ψn+ j,k + ψj,k n D + t ψj,k n = ψn+ j,k t ψn j,k, n D x ψ n j,k = ψn j,k ψ n j,k + ψ n j+,k x, j Z D y ψ n j,k = ψn j,k ψ n j,k + ψ n j,k+ y, k N DERIVATION OF DISCRETE TRB
21 Discrete Sine-Transformation in y-direction: ˆψ n j,m := K K X k= «πkm ψj,k n sin K m =,...,K Motivation: Solve discrete stationary Schrödinger equation in D: y χ m j,k = E m χ m j,k, k =,...,K χ m j, = χ m j,k =. The eigenfunctions χ m j,k = sin ` πkm K E m = πm cos y K provide the energies. Hence follows for m =,...,K y ψ n j,k ψ n j,k + ψ n j,k+ b m = πm cos y K ˆψn j,m. DERIVATION OF DISCRETE TRB
22 Sine-Transformation of the discrete Schrödinger equation on the exterior domains j, j J yields for the modes m =,...,K: C m j+ ˆψ n+ j+,m + Cm j ˆψ n+ j,m + Rm j ˆψ n+ j,m = (D C m j+) ˆψ n j+,m + (D C m j ) ˆψ n j,m + (B m j R m j ) ˆψ n j,m. DERIVATION OF DISCRETE TRB 3
23 Sine-Transformation of the discrete Schrödinger equation on the exterior domains j, j J yields for the modes m =,...,K: C m j+ ˆψ n+ j+,m + Cm j ˆψ n+ j,m + Rm j ˆψ n+ j,m = (D C m j+) ˆψ n j+,m + (D C m j ) ˆψ n j,m + (B m j R m j ) ˆψ n j,m. Definition [Z - Transformation]: The Z-Transformation of a sequence (ψ n ) n N is given by Z {ψ n } = Ψ(z) := ψ n z n z C, z >. n= DERIVATION OF DISCRETE TRB 3-A
24 One can show: Z = z ˆψ j,m + zψ m j (z) ˆψn+ j,m ψ J+,k = ψ J,k = ψ J,k = for k =,...,K V j constant for j, j J C m j = C m, R m j = R m, B m j = B m Ψ J+ (z) +» R z + R B Ψ J (z) + Ψ J (z) =. C z + C D DERIVATION OF DISCRETE TRB 4
25 One can show: Z = z ˆψ j,m + zψ m j (z) ˆψn+ j,m ψ J+,k = ψ J,k = ψ J,k = for k =,...,K V j constant for j, j J C m j = C m, R m j = R m, B m j = B m Ψ J+ (z) +» R z + R B Ψ J (z) + Ψ J (z) =. C z + C D This difference equation with constant coefficients is solved by Ψ j (z) = ν j (z):» R z + R B ν (z) + ν(z) + =, C z + C D Physical background forces decay of the solution for j ν(z) > und ν(z)ψ J (z) = Ψ J (z) Z-transformed DTBC DERIVATION OF DISCRETE TRB 4-A
26 Theorem [DTBC for the D Schrödinger equation]: Discretize the D Schrödinger-Equation with the compact 9-point difference scheme in space and with the Crank-Nicolson scheme in time. Then the DTBC at x J = J x and x = for n read ˆψ n,m s (),m ˆψ n,m = ˆψ n J,m s () J,m ˆψ n J,m = n ν= n ν= s (n ν),m ˆψ ν,m a m s (n ν) J,m ˆψ ν J,m a m J ˆψ n,m, ˆψ n J,m. The convolution coefficients s (n) j,m can be calculated by s (n) j,m = αm j (λ m j ) n n [ ] P n (µ m j ) P n (µ m j ) with the Legendre-Polynomials P n ( P P ). DERIVATION OF DISCRETE TRB 5
27 Advantages of the new DTBC: no numerical reflections 3-point recursion for s (n) these DTBC have exactly the same structure like the DTBC calculated with the 5-point scheme [A. Arnold, M. Ehrhardt] convergence: O( x 4 + y 4 + t ) same numerical effort like discretized analytical TRB s (n) j,m = O(n 3/ ) CN-FD scheme with DTBC is unconditionally stable, with A := I + x D x + y D y follows: D + t ψ A = with ψ A := ψ, Aψ. () DERIVATION OF DISCRETE TRB 6
28 some drawbacks: DTBC are non-local in time high memory costs: the solution ψ n j,m x and x J for all time steps n =,,... has to be saved in in each time step n =,,... you have to calculate K convolutions of the length n (FFT not useable, O(Kn )) DERIVATION OF DISCRETE TRB 7
29 some drawbacks: DTBC are non-local in time high memory costs: the solution ψ n j,m x and x J for all time steps n =,,... has to be saved in in each time step n =,,... you have to calculate K convolutions of the length n (FFT not useable, O(Kn )) These DTBC are given in the Sine-transformed form: BC of one mode is a linear combination of all other boundary points diagonal structure of the system matrix is destroyed DERIVATION OF DISCRETE TRB 7-A
30 Sparsity pattern of the system matrix for J = K = DERIVATION OF DISCRETE TRB 8
31 Example : free D Schrödinger equation on Ω = [, ] [, ] with the initial data: + y ik x+ik y 6 x ) ( ( x y ), (x, y) Ω ψ I (x, y) = e it = it = 6 it = x y T = D ERIVATION OF DISCRETE TRB y T = 6 t x x y T = 8 t 9
32 * APPROXIMATION OF DTBC Idea: Approximate the convolution coefficients s (n) j,m exponentials: by a sum of s (n) s (n) = L l= b l q n l, n N, q l >, L 4 APPROXIMATION OF DTBC
33 * APPROXIMATION OF DTBC Idea: Approximate the convolution coefficients s (n) j,m exponentials: by a sum of s (n) s (n) = L l= b l q n l, n N, q l >, L 4 b l, q l are calculated by the Padé - Approximation of f(x) = L n= s (n) x n, x [A.Arnold, M. Ehrhardt, I. Sofronov (3)] APPROXIMATION OF DTBC -A
34 * Recursion formula for the convolution coefficients: n X t= s (n t) ψ t = LX l= c (n) l with c (n) l = q l c () l = c (n ) l + b l q l ψ n, n =,...,N APPROXIMATION OF DTBC
35 * Recursion formula for the convolution coefficients: n X t= s (n t) ψ t = LX l= c (n) l with c (n) l = q l c () l = c (n ) l + b l q l ψ n, n =,...,N Advantages: If you have calculated b l, q l once for a set x, y, t, V, you ll easily derive b l, q l for any x, y, t, V by q l = q lā b a q l b b l = b l q l aā b b (a q l b)(q l ā b) Numerical effort: O(Kn ) O(KLn) Memory: O(Kn) O(KL) + q l + q l APPROXIMATION OF DTBC -A
36 * Example : free D Schrödinger-Equation in Ω = [, ] [, ] Initial function: ψ I (x, y) = sin(πy)e ik xx 6(x ), (x, y) Ω IT=5 IT=3 IT= psi(x,y), t=t.6.4 psi(x,y), t=t.6.4 psi(x,y), t=t.6.4 L = 5 : y T = 5 t x y T = 3 t x y T = 5 t x IT=5 IT=3 IT= psi(x,y), t=t.6.4 psi(x,y), t=t.6.4 psi(x,y), t=t.6.4 L = : y T = 5 t x y T = 3 t x y T = 5 t x APPROXIMATION OF DTBC
37 SIMULATION OF QUANTUM WAVEGUIDES y 3nm V control 4.5nm Y ψ inc Γ nm Ω Γ nm X x Incoming wave at x = : ψ inc (,y, t) = sin(πy)e iet, E = 9.9meV inhomogeneous DTBC at x =, DTBC at x = X SIMULATION OF QUANTUM WAVEGUIDES 3
38 Little trick to suppress oscillations in time: ψ inc oscillates like e iet in time. E big: fast oscillation of the solution in time small time step size is necessary high numerical effort for the analysis of steady state and long-time behaviour Define ϕ(x, y, t) := e iωt ψ(x,y, t) mit ω = E. ϕ solve the modified Schrödinger equation i ϕ t = m (ϕ xx + ϕ yy ) + (V ω ) {z } =:Ṽ ϕ SIMULATION OF QUANTUM WAVEGUIDES 4
39 .4 V= V= E t [ s] f (t) = ψ (x, y, t) ψ ref (x, y, t) f (t) = ψ (x, y, t) ψ ref (x, y, t) ψ is calculated with Ṽ =, ψ with Ṽ = E and ψ ref is a numerical reference solution, which has been calculated with high accuracy. SIMULATION OF QUANTUM WAVEGUIDES 5
40 Initial function: X=, Y=.86667, dx=.6667, dy=.6667, dt=., V= Quantenwellen.m, solution after T =.7998 ps, it = Quantenwellen.m, solution after T =.998 ps, it = psi(x,y), t=.5 ψ(x,y,t).5 ψ(x,y,t) x T = t y y [nm] T = t y [nm] T = 9 t Quantenwellen.m, solution after T =. ps, it = Quantenwellen.m, solution after T =.8 ps, it = 4 Quantenwellen.m, solution after T =.9 ps, it = ψ(x,y,t).5 ψ(x,y,t).5 ψ(x,y,t) y [nm] T = t y [nm] T = 4 t y [nm] T = 95 t SIMULATION OF QUANTUM WAVEGUIDES 6
41 EXTENSION OF THE DTBC TO MORE ARBITRARY POTENTIALS Drawbacks of the simulations: Hard walls (zero Dirichlet BCs) and edges are not practicable in industry! Geometry of computational domain shall be realized by potentials also in the simulations. How to choose the incoming wave and the initial function then? Potentials are NOT constant in the exterior domains! Decoupling of the modes for V (x, y) const. after Sine- Transformation is not possible EXTENSION OF THE DTBC TO MORE ARBITRARY POTENTIALS 7
42 EXTENSION OF THE DTBC TO MORE ARBITRARY POTENTIALS Drawbacks of the simulations: Hard walls (zero Dirichlet BCs) and edges are not practicable in industry! Geometry of computational domain shall be realized by potentials also in the simulations. How to choose the incoming wave and the initial function then? Potentials are NOT constant in the exterior domains! Decoupling of the modes for V (x, y) const. after Sine- Transformation is not possible Discrete Schrödinger equation: 5-point scheme in space, Crank-Nicolson in time: i D + t ψj,k n = `D m x + Dy ψn+ j,k + V n+ j,k ψ n+ j,k Calculate new Eigenfunctions, which take the potential into account! EXTENSION OF THE DTBC TO MORE ARBITRARY POTENTIALS 7-A
43 Solve the eigenvalue equation on the exterior domains (j, j J): with `χn y j,k,m χ n j,k,m + χ n n+ j,k+,m + V k χ n j,k,m = Ej,mχ n n j,k,m y KX χ n j,k,m = and χ n j,,m = χ n j,k,m = k= for k, m K, n >. Transformation w.r.t. the eigenfunctions yields ˆψ n j,m = y K X k= χ n j,k,mψ n j,k, m K. i D + t ˆψ j,m n = m D ˆψ n+ x j,m + En ˆψ n+ j,m j,m same structure as the sine-transformed Schrödinger equation! [N. Ben Abdallah, M.S. (5)] EXTENSION OF THE DTBC TO MORE ARBITRARY POTENTIALS 8
44 Example 3: D channel Potential V(,y) Eigenfunction χ for j=, χ j,k,m of the m = st mode with the Eigenvalue λ = V(,y k ) χ j,k,m y k y k potential V (y) = 4y( y) eigenfunction χ of m = initial function: incoming wave: ψ I j,k = e ik xj x χ j,k,m ψ inc j,k,n = e ik xj x χ j,k,me ie xn t with E x = cos(k x x) x EXTENSION OF THE DTBC TO MORE ARBITRARY POTENTIALS 9
45 t = 5 t =. t= ψ(x,y,t) ψ(x,y,t) ψ(x,y,t) x y.8 x T = t t = 5 t = T = 5 t y y ψ(x,y,t) ψ(x,y,t) t = 7 t =.34.5 x.6 T = 5 t.4.4 x..5 y. T = 5 t t = 5 t =.5 ψ(x,y,t) t = 5 t =.3 x.8.5 y T = 5 t E XTENSION OF THE DTBC TO MORE ARBITRARY POTENTIALS x.8 y T = 7 t 3
46 SCHRÖDINGER EQUATION ON CIRCULAR DOMAINS Solve the Schrödinger equation (in polar coordinates) iψ t = r (rψ r) r + «r ψ θθ +V (r,θ, t)ψ, r >, θ π, t > on a circular domain Ω = [,R] [,π] with TBC at x = R. Problems: solve a second order difference equation with varying coefficients: a j Ψ J+ (z) + b j (z)ψ J (z) + c j Ψ J (z) = calculation of the convolution coefficients for the DTBC "recursion from infinity" singularity at r = not equidistant offset-grid r j = r j+ approximation of the convolution coefficients and the -sum [A.Arnold, M. Ehrhardt, M. S., I. Sofronov (6)] SCHRÖDINGER EQUATION ON CIRCULAR DOMAINS 3
47 Example 4: free Schrödinger equation on unit disc Ω = [, ] [, π] ψ I (r,θ) = αx α y e ik xr cos θ+ik y r sin θ (r cos θ) (r sin θ) αx αy initial function on Ω = [,]x[,π], r = /8, θ = π /8 solution on Ω = [,]x[,π] at time t=.5, t = r = /8, θ = π /8 solution on Ω = [,]x[,π] at time t=.5, t = r = /8, θ = π / ψ(r,θ,) ψ(r,θ,t) 4 3 ψ(r,θ,t) initial function T =.5 T =.5 initial function on Ω = [,]x[,π], r = /8, θ = π /8 solution on Ω = [,]x[,π] at time t=.5, t = r = /8, θ = π /8 solution on Ω = [,]x[,π] at time t=.5, t = r = /8, θ = π / SCHRÖDINGER EQUATION ON CIRCULAR DOMAINS 3
48 Error due to the scheme/tbc: L(ψ, ϕ, t n, Ω) := ψ(r j, θ k, t n ) ϕ(r j, θ k, t n ) Ω, ϕ(r j, θ k,t n ) Ω, ψ: numerical solution ϕ: exact solution or numerical reference solution relative L error due to the scheme J=K=64, t=/64 J=K=8, t=/8 J=K=56, t=/56 relative L error due to the TBC t 5 J=K=64, t=/64 J=K=8, t=/8 J=K=56, t=/ t SCHRÖDINGER EQUATION ON CIRCULAR DOMAINS 33
Discrete Transparent Boundary Conditions for the Schrödinger Equation. Matthias Ehrhardt, Anton Arnold
Discrete Transparent Boundary Conditions for the Schrödinger Equation Matthias Ehrhardt, Anton Arnold Fakultät für Mathematik und Informatik, Universität des Saarlandes, Postfach 15 11 5, D-6641 Saarbrücken,
More informationArtificial boundary conditions for dispersive equations. Christophe Besse
Artificial boundary conditions for dispersive equations by Christophe Besse Institut Mathématique de Toulouse, Université Toulouse 3, CNRS Groupe de travail MathOcéan Bordeaux INSTITUT de MATHEMATIQUES
More informationPerfectly matched layers versus discrete transparent boundary conditions in quantum device simulations
ASC Report No. 43/213 Perfectly matched layers versus discrete transparent boundary conditions in quantum device simulations J.-F. Mennemann, A. Jüngel Institute for Analysis and Scientific Computing Vienna
More informationNumerical aspects of the nonlinear Schro dinger equation. Christophe Besse
Numerical aspects of the nonlinear Schro dinger equation by Christophe Besse Laboratoire Paul Painleve, Universite Lille 1, CNRS, Inria Simpaf Team Lille Outline of the talk 1 Motivation 2 Numerical methods
More informationAnton ARNOLD. with N. Ben Abdallah (Toulouse), J. Geier (Vienna), C. Negulescu (Marseille) TU Vienna Institute for Analysis and Scientific Computing
TECHNISCHE UNIVERSITÄT WIEN Asymptotically correct finite difference schemes for highly oscillatory ODEs Anton ARNOLD with N. Ben Abdallah (Toulouse, J. Geier (Vienna, C. Negulescu (Marseille TU Vienna
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 informationMethod of infinite system of equations for problems in unbounded domains
Method of infinite system of equations for problems in unbounded domains Dang Quang A Institute of Information Technology 18 Hoang Quoc Viet, Cau giay, Hanoi, Vietnam Innovative Time Integration, May 13-16,
More information6 Non-homogeneous Heat Problems
6 Non-homogeneous Heat Problems Up to this point all the problems we have considered for the heat or wave equation we what we call homogeneous problems. This means that for an interval < x < l the problems
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 informationLecture 19: Heat conduction with distributed sources/sinks
Introductory lecture notes on Partial Differential Equations - c Anthony Peirce. Not to be copied, used, or revised without explicit written permission from the copyright owner. 1 ecture 19: Heat conduction
More informationThe Factorization Method for Inverse Scattering Problems Part I
The Factorization Method for Inverse Scattering Problems Part I Andreas Kirsch Madrid 2011 Department of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research Center
More informationLecture 7. More dimensions
Lecture 7 More dimensions 67 68 LECTURE 7. MORE DIMENSIONS 7.1 Introduction In this lecture we generalize the concepts introduced so far to systems that evolve in more than one spatial dimension. While
More information4 Power Series Solutions: Frobenius Method
4 Power Series Solutions: Frobenius Method Now the ODE adventure takes us to series solutions for ODEs, a technique A & W, that is often viable, valuable and informative. These can be readily applied Sec.
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationf xx g dx, (5) The point (i) is straightforward. Let us check the point (ii). Doing the integral by parts, we get 1 f g xx dx = f g x 1 0 f xg x dx =
Problem 19. Consider the heat equation u t = u xx, (1) u = u(x, t), x [, 1], with the boundary conditions (note the derivative in the first one) u x (, t) =, u(1, t) = (2) and the initial condition u(x,
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 21 Quantum Mechanics in Three Dimensions Lecture 21 Physics 342 Quantum Mechanics I Monday, March 22nd, 21 We are used to the temporal separation that gives, for example, the timeindependent
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
More informationQuantum Light-Matter Interactions
Quantum Light-Matter Interactions QIC 895: Theory of Quantum Optics David Layden June 8, 2015 Outline Background Review Jaynes-Cummings Model Vacuum Rabi Oscillations, Collapse & Revival Spontaneous Emission
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationQuantum Mechanics Exercises and solutions
Quantum Mechanics Exercises and solutions P.J. Mulders Department of Physics and Astronomy, Faculty of Sciences, Vrije Universiteit Amsterdam De Boelelaan 181, 181 HV Amsterdam, the Netherlands email:
More informationApproximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation
Approximation and Fast Cacuation of Non-oca Boundary Conditions for the Time-dependent Schrödinger Equation Anton Arnod, Matthias Ehrhardt 2, and Ivan Sofronov 3 Universität Münster, Institut für Numerische
More informationVibrating-string problem
EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u
More informationUnderstand the basic principles of spectroscopy using selection rules and the energy levels. Derive Hund s Rule from the symmetrization postulate.
CHEM 5314: Advanced Physical Chemistry Overall Goals: Use quantum mechanics to understand that molecules have quantized translational, rotational, vibrational, and electronic energy levels. In a large
More informationMATH 220: THE FOURIER TRANSFORM TEMPERED DISTRIBUTIONS. Beforehand we constructed distributions by taking the set Cc
MATH 220: THE FOURIER TRANSFORM TEMPERED DISTRIBUTIONS Beforehand we constructed distributions by taking the set Cc ( ) as the set of very nice functions, and defined distributions as continuous linear
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationThe one and three-dimensional particle in a box are prototypes of bound systems. As we
6 Lecture 10 The one and three-dimensional particle in a box are prototypes of bound systems. As we move on in our study of quantum chemistry, we'll be considering bound systems that are more and more
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 informationUNIVERSITY OF SURREY FACULTY OF ENGINEERING AND PHYSICAL SCIENCES DEPARTMENT OF PHYSICS. BSc and MPhys Undergraduate Programmes in Physics LEVEL HE2
Phys/Level /1/9/Semester, 009-10 (1 handout) UNIVERSITY OF SURREY FACULTY OF ENGINEERING AND PHYSICAL SCIENCES DEPARTMENT OF PHYSICS BSc and MPhys Undergraduate Programmes in Physics LEVEL HE PAPER 1 MATHEMATICAL,
More informationStability estimates in stationary inverse transport
Stability estimates in stationary inverse transport Guillaume Bal and Alexandre Jollivet April 8, 28 Abstract We study the stability of the reconstruction of the scattering and absorption coefficients
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, x < 1, 1 + x, 1 < x < 0, ψ(x) = 1 x, 0 < x < 1, 0, x > 1, so that it verifies ψ 0, ψ(x) = 0 if x 1 and ψ(x)dx = 1. Consider (ψ j )
More informationChemistry 532 Problem Set 7 Spring 2012 Solutions
Chemistry 53 Problem Set 7 Spring 01 Solutions 1. The study of the time-independent Schrödinger equation for a one-dimensional particle subject to the potential function leads to the differential equation
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationMathematical Methods in Quantum Mechanics With Applications to Schrödinger Operators 2nd edition
Errata G. Teschl, Mathematical Methods in Quantum Mechanics With Applications to Schrödinger Operators nd edition Graduate Studies in Mathematics, Vol. 157 American Mathematical Society, Providence, Rhode
More informationSAMPLE FINAL EXAM SOLUTIONS
LAST (family) NAME: FIRST (given) NAME: ID # : MATHEMATICS 3FF3 McMaster University Final Examination Day Class Duration of Examination: 3 hours Dr. J.-P. Gabardo THIS EXAMINATION PAPER INCLUDES 22 PAGES
More information17 Source Problems for Heat and Wave IB- VPs
17 Source Problems for Heat and Wave IB- VPs We have mostly dealt with homogeneous equations, homogeneous b.c.s in this course so far. Recall that if we have non-homogeneous b.c.s, then we want to first
More information1 A complete Fourier series solution
Math 128 Notes 13 In this last set of notes I will try to tie up some loose ends. 1 A complete Fourier series solution First here is an example of the full solution of a pde by Fourier series. Consider
More informationDynamically assisted Sauter-Schwinger effect
Dynamically assisted Sauter-Schwinger effect Ralf Schützhold Fachbereich Physik Universität Duisburg-ssen Dynamically assisted Sauter-Schwinger effect p.1/16 Dirac Sea Schrödinger equation (non-relativistic)
More informationPartial Differential Equations for Engineering Math 312, Fall 2012
Partial Differential Equations for Engineering Math 312, Fall 2012 Jens Lorenz July 17, 2012 Contents Department of Mathematics and Statistics, UNM, Albuquerque, NM 87131 1 Second Order ODEs with Constant
More informationLinear Differential Equations. Problems
Chapter 1 Linear Differential Equations. Problems 1.1 Introduction 1.1.1 Show that the function ϕ : R R, given by the expression ϕ(t) = 2e 3t for all t R, is a solution of the Initial Value Problem x =
More informationLecture #5: Begin Quantum Mechanics: Free Particle and Particle in a 1D Box
561 Fall 017 Lecture #5 page 1 Last time: Lecture #5: Begin Quantum Mechanics: Free Particle and Particle in a 1D Box 1-D Wave equation u x = 1 u v t * u(x,t): displacements as function of x,t * nd -order:
More informationCollection of formulae Quantum mechanics. Basic Formulas Division of Material Science Hans Weber. Operators
Basic Formulas 17-1-1 Division of Material Science Hans Weer The de Broglie wave length λ = h p The Schrödinger equation Hψr,t = i h t ψr,t Stationary states Hψr,t = Eψr,t Collection of formulae Quantum
More informationKernel B Splines and Interpolation
Kernel B Splines and Interpolation M. Bozzini, L. Lenarduzzi and R. Schaback February 6, 5 Abstract This paper applies divided differences to conditionally positive definite kernels in order to generate
More informationLINEAR RESPONSE THEORY
MIT Department of Chemistry 5.74, Spring 5: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff p. 8 LINEAR RESPONSE THEORY We have statistically described the time-dependent behavior
More informationSeparation of Variables
Separation of Variables A typical starting point to study differential equations is to guess solutions of a certain form. Since we will deal with linear PDEs, the superposition principle will allow us
More informationQuasi-Particle Dynamics of Linearly Coupled Systems of Nonlinear Schrödinger Equations
Quasi-Particle Dynamics of Linearly Coupled Systems of Nonlinear Schrödinger Equations Michail D. Todorov Faculty of Applied Mathematics and Computer Science Technical University of Sofia, Bulgaria SS25
More information11 Quantum theory: introduction and principles
Part 2: Structure Quantum theory: introduction and principles Solutions to exercises E.b E.2b E.3b E.4b E.5b E.6b Discussion questions A successful theory of black-body radiation must be able to explain
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 informationFourier mode dynamics for NLS Synchronization in fiber lasers arrays
Fourier mode dynamics for NLS Synchronization in fiber lasers arrays Nonlinear Schrodinger workshop Heraklio, May 21 2013 Jean-guy Caputo Laboratoire de Mathématiques INSA de Rouen, France A. Aceves, N.
More informationAsymptotic behaviour of the heat equation in twisted waveguides
Asymptotic behaviour of the heat equation in twisted waveguides Gabriela Malenová Faculty of Nuclear Sciences and Physical Engineering, CTU, Prague Nuclear Physics Institute, AS ČR, Řež Graphs and Spectra,
More informationProper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems. Stefan Volkwein
Proper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems Institute for Mathematics and Scientific Computing, Austria DISC Summerschool 5 Outline of the talk POD and singular value decomposition
More informationPY 351 Modern Physics - Lecture notes, 3
PY 351 Modern Physics - Lecture notes, 3 Copyright by Claudio Rebbi, Boston University, October 2016. These notes cannot be duplicated and distributed without explicit permission of the author. Time dependence
More informationPhysics 486 Discussion 5 Piecewise Potentials
Physics 486 Discussion 5 Piecewise Potentials Problem 1 : Infinite Potential Well Checkpoints 1 Consider the infinite well potential V(x) = 0 for 0 < x < 1 elsewhere. (a) First, think classically. Potential
More informationDifferential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm
Differential Equations 228 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 25 at 2:5pm Instructions: This in-class exam is 5 minutes. No calculators, notes, tables or books. No answer check is
More information4. Supplementary Notes on Time and Space Evolution of a Neutrino Beam
Lecture Notes for Quantum Physics II & III 8.05 & 8.059 Academic Year 1996/1997 4. Supplementary Notes on Time and Space Evolution of a Neutrino Beam c D. Stelitano 1996 As an example of a two-state system
More informationSpecial Instructions:
Be sure that this examination has 20 pages including this cover The University of British Columbia Sessional Examinations - December 2016 Mathematics 257/316 Partial Differential Equations Closed book
More informationBessel functions. The drum problem: Consider the wave equation (with c = 1) in a disc with homogeneous Dirichlet boundary conditions:
Bessel functions The drum problem: Consider the wave equation (with c = 1) in a disc with homogeneous Dirichlet boundary conditions: u = u t, u(r 0,θ,t) = 0, u u(r,θ,0) = f(r,θ), (r,θ,0) = g(r,θ). t (Note
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationPhysics 443, Solutions to PS 1 1
Physics 443, Solutions to PS. Griffiths.9 For Φ(x, t A exp[ a( mx + it], we need that + h Φ(x, t dx. Using the known result of a Gaussian intergral + exp[ ax ]dx /a, we find that: am A h. ( The Schrödinger
More informationMicrolocal analysis and inverse problems Lecture 3 : Carleman estimates
Microlocal analysis and inverse problems ecture 3 : Carleman estimates David Dos Santos Ferreira AGA Université de Paris 13 Monday May 16 Instituto de Ciencias Matemáticas, Madrid David Dos Santos Ferreira
More informationMathematical Tripos Part IB Michaelmas Term Example Sheet 1. Values of some physical constants are given on the supplementary sheet
Mathematical Tripos Part IB Michaelmas Term 2015 Quantum Mechanics Dr. J.M. Evans Example Sheet 1 Values of some physical constants are given on the supplementary sheet 1. Whenasampleofpotassiumisilluminatedwithlightofwavelength3
More informationMath 251 December 14, 2005 Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt
Math 251 December 14, 2005 Final Exam Name Section There are 10 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning of each question
More information16.1. PROBLEM SET I 197
6.. PROBLEM SET I 97 Answers: Problem set I. a In one dimension, the current operator is specified by ĵ = m ψ ˆpψ + ψˆpψ. Applied to the left hand side of the system outside the region of the potential,
More informationBounds 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 informationQualification Exam: Mathematical Methods
Qualification Exam: Mathematical Methods Name:, QEID#41534189: August, 218 Qualification Exam QEID#41534189 2 1 Mathematical Methods I Problem 1. ID:MM-1-2 Solve the differential equation dy + y = sin
More informationFOURIER METHODS AND DISTRIBUTIONS: SOLUTIONS
Centre for Mathematical Sciences Mathematics, Faculty of Science FOURIER METHODS AND DISTRIBUTIONS: SOLUTIONS. We make the Ansatz u(x, y) = ϕ(x)ψ(y) and look for a solution which satisfies the boundary
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationarxiv: v1 [physics.comp-ph] 28 Jan 2008
A new method for studying the vibration of non-homogeneous membranes arxiv:0801.4369v1 [physics.comp-ph] 28 Jan 2008 Abstract Paolo Amore Facultad de Ciencias, Universidad de Colima, Bernal Díaz del Castillo
More information1. Fourier Transform (Continuous time) A finite energy signal is a signal f(t) for which. f(t) 2 dt < Scalar product: f(t)g(t)dt
1. Fourier Transform (Continuous time) 1.1. Signals with finite energy A finite energy signal is a signal f(t) for which Scalar product: f(t) 2 dt < f(t), g(t) = 1 2π f(t)g(t)dt The Hilbert space of all
More informationPHZ 6607 Fall 2004 Homework #4, Due Friday, October 22, 2004
Read Chapters 9, 10 and 20. PHZ 6607 Fall 2004 Homework #4, Due Friday, October 22, 2004 1. The usual metric of four-dimensional flat Minkowski-space in spherical-polar coordinates is ds 2 = dt 2 + dr
More informationExit times of diffusions with incompressible drifts
Exit times of diffusions with incompressible drifts Andrej Zlatoš University of Chicago Joint work with: Gautam Iyer (Carnegie Mellon University) Alexei Novikov (Pennylvania State University) Lenya Ryzhik
More informationPHY 396 K. Problem set #5. Due October 9, 2008.
PHY 396 K. Problem set #5. Due October 9, 2008.. First, an exercise in bosonic commutation relations [â α, â β = 0, [â α, â β = 0, [â α, â β = δ αβ. ( (a Calculate the commutators [â αâ β, â γ, [â αâ β,
More informationApplied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.
Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the
More informationPhysics 137A Quantum Mechanics Fall 2012 Midterm II - Solutions
Physics 37A Quantum Mechanics Fall 0 Midterm II - Solutions These are the solutions to the exam given to Lecture Problem [5 points] Consider a particle with mass m charge q in a simple harmonic oscillator
More informationA Friendly Review of Absorbing Boundary Conditions and Perfectly Matched Layers for Classical and Relativistic Quantum Waves Equations
A Friendly Review of Absorbing Boundary Conditions and Perfectly Matched Layers for Classical and Relativistic Quantum Waves Equations X. Antoine a, E. Lorin b,c and Q. Tang a a Institut Elie Cartan de
More informationTransient Phenomena in Quantum Bound States Subjected to a Sudden Perturbation
Symmetry, Integrability and Geometry: Methods and Applications Vol. (5), Paper 3, 9 pages Transient Phenomena in Quantum Bound States Subjected to a Sudden Perturbation Marcos MOSHINSKY and Emerson SADURNÍ
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 informationProblem 1: A 3-D Spherical Well(10 Points)
Problem : A 3-D Spherical Well( Points) For this problem, consider a particle of mass m in a three-dimensional spherical potential well, V (r), given as, V = r a/2 V = W r > a/2. with W >. All of the following
More informationQuantum Mechanics Solutions
Quantum Mechanics Solutions (a (i f A and B are Hermitian, since (AB = B A = BA, operator AB is Hermitian if and only if A and B commute So, we know that [A,B] = 0, which means that the Hilbert space H
More informationHomework for Math , Fall 2016
Homework for Math 5440 1, Fall 2016 A. Treibergs, Instructor November 22, 2016 Our text is by Walter A. Strauss, Introduction to Partial Differential Equations 2nd ed., Wiley, 2007. Please read the relevant
More informationNumerical Methods 2: Hill s Method
Numerical Methods 2: Hill s Method Bernard Deconinck Department of Applied Mathematics University of Washington bernard@amath.washington.edu http://www.amath.washington.edu/ bernard Stability and Instability
More informationBoundary conditions and estimates for the linearized Navier-Stokes equations on staggered grids
Boundary conditions and estimates for the linearized Navier-Stokes equations on staggered grids Wendy Kress and Jonas Nilsson Department of Scientific Computing Uppsala University Box S-75 4 Uppsala Sweden
More informationCOMPLEX SPHERICAL WAVES AND INVERSE PROBLEMS IN UNBOUNDED DOMAINS
COMPLEX SPHERICAL WAVES AND INVERSE PROBLEMS IN UNBOUNDED DOMAINS MIKKO SALO AND JENN-NAN WANG Abstract. This work is motivated by the inverse conductivity problem of identifying an embedded object in
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationWEYL S LEMMA, ONE OF MANY. Daniel W. Stroock
WEYL S LEMMA, ONE OF MANY Daniel W Stroock Abstract This note is a brief, and somewhat biased, account of the evolution of what people working in PDE s call Weyl s Lemma about the regularity of solutions
More informationPartial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:
Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant
More informationOn Chern-Simons-Schrödinger equations including a vortex point
On Chern-Simons-Schrödinger equations including a vortex point Alessio Pomponio Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari Workshop in Nonlinear PDEs Brussels, September 7
More informationQuantization of Scalar Field
Quantization of Scalar Field Wei Wang 2017.10.12 Wei Wang(SJTU) Lectures on QFT 2017.10.12 1 / 41 Contents 1 From classical theory to quantum theory 2 Quantization of real scalar field 3 Quantization of
More informationThis ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0
Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More informationTransient Schrödinger-Poisson Simulations of a High-Frequency Resonant Tunneling Diode Oscillator
Transient Schrödinger-Poisson Simulations of a High-Frequency Resonant Tunneling Diode Oscillator Jan-Frederik Mennemann Institute for Analysis and Scientific Computing, Vienna University of Technology,
More information1 Distributions (due January 22, 2009)
Distributions (due January 22, 29). The distribution derivative of the locally integrable function ln( x ) is the principal value distribution /x. We know that, φ = lim φ(x) dx. x ɛ x Show that x, φ =
More informationPartial Differential Equations
Partial Differential Equations Xu Chen Assistant Professor United Technologies Engineering Build, Rm. 382 Department of Mechanical Engineering University of Connecticut xchen@engr.uconn.edu Contents 1
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem in a 2D Box Region Let us consider a model problem of parabolic
More informationC. Show your answer in part B agrees with your answer in part A in the limit that the constant c 0.
Problem #1 A. A projectile of mass m is shot vertically in the gravitational field. Its initial velocity is v o. Assuming there is no air resistance, how high does m go? B. Now assume the projectile is
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH College of Informatics and Electronics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MS425 SEMESTER: Autumn 25/6 MODULE TITLE: Applied Analysis DURATION OF EXAMINATION:
More informationTotal Angular Momentum for Hydrogen
Physics 4 Lecture 7 Total Angular Momentum for Hydrogen Lecture 7 Physics 4 Quantum Mechanics I Friday, April th, 008 We have the Hydrogen Hamiltonian for central potential φ(r), we can write: H r = p
More information