221A Lecture Notes WKB Method

Save this PDF as:
 WORD  PNG  TXT  JPG

Size: px
Start display at page:

Download "221A Lecture Notes WKB Method"

Transcription

1 A Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using ψ x, t = e is x,t/ h : S [ t ψ = m S i h ] m S + V ψ. Assuming ψ 0, this leds to n eqution S t = m S i h m S + V. 3 Now tking the forml limit h 0, it becomes the sme s the clssicl Hmilgon Jcobi eqution S t = m S + V. It is often sid tht the h 0 limit is indeed the clssicl mechnics becuse it reduces to the Hmilton Jcobi eqution. Clssicl Limit The fct tht the exponent S stisfies the Hmilton Jcobi eqution my not men the clssicl behvior to you right wy. The wy to see it is by forming wve pcket, nd sk the question how does the wve pcket moves s whole. To simplify the discussion, let us consentrte on the cse where the Hmiltonin does not depend explicitly on time, nd we cn seprte the vrible s S x, t = S x, E Et. 5

2 The inverse Legendre trnsformtion sys t = S x, E E. 6 This condition cn be viewed s n implicit eqution to determine the position of the prticle x s function of time t. To form wve pcket, we need to put together wves of slighly different energies. We cn write ψ x, t = degee i S x,e Et/ h, 7 The question here is where the wve pcket is. At most plces, the verging over E mkes the phse oscillte very fst becuse of the smll h, nd there is strong cncelltion. Only t specil points where the phse fctor is sttionry with respect to E, the wve function is sizble. Therefore, the position of the wve pcket is determined by the sttionry condition S x, E S x, E Et = t = 0. 8 E E This is nothing but Eq. 6. In other words, the wve pcket is sizble only where the position stisfies the clssicl eqution of motion. Therefore, the wve pcket follows the clssicl eqution of motion. 3 h Expnsion Now tht we hve rederived clssicl mechnics from quntum mechnics, our interest now turns to quntum mechnics itself. We will use the fct tht h is indeed smll to find useful pproximtion method in quntum mechnics. The Eq. 3 is exct s long s ψ 0. One cn use this eqution, nd consider n expnsion S x, t = S 0 + hs + h S +. 9 This is n expnsion in h, nd hence clled h-expnsion or semi-clssicl expnsion. Plugging in the expnsion into Eq. 3, we find S 0 t S t = m S 0 + V, 0 = [ i S 0 + S m 0 S ],

3 nd similrly for higher terms in h. The leding eqution hs only S 0, nd it is exctly the sme s Hmilton Jcobi eqution. Once you solve these equtions strting from S 0, S, etc, you hve solved the wve function ψ in systemtic expnsion in h. There is one importnt difference between the clssicl limit which reduces to the Hmilton Jcobi eqution nd the ide of expnsion in h. In the clssicl limit, S is of course rel becuse it does not mke sense otherwise. However, in the sense of n expnsion in h, we don t know if S 0 hs to be rel. In fct, in mny interesting cses, S 0 turns out to be complex. We will see these exmples below. WKB Approximtion WKB Approximtion, due to Wentzel, Krmers, nd Brillouin, keeps terms up to O h in the h expnsion. It is used mostly for the time-independent cse, or in other words, for n eigenstte of energy E. Then the wve function hs the ordinry time dependence e iet/ h. We lso restrict ourselves to the one-dimensionl problem. In terms of S, it corresponds to S x, t = S x Et. Therefore only S 0 hs the time dependence S 0 x, t = S 0 x Et, while higher order terms S i = S i x for i 0 do not depend on time. The lowest order term S 0 stisfies the Hmilton Jcobi eqution see Eq. 0 E = m S 0 + V x. 3 The differentil eqution cn be solve immeditely s S 0 x = ± me V x dx = px dx up to n integrtion constnt which cn be determined only fter imposing boundry condition on the wve function. We used the nottion px = ± me V x becuse it is nothing but the momentum of the prticle in the clssicl sense. Once we know S 0, we cn lso solve for S. Strting from Eq., nd using S / t = 0, we find S 0S = is 0, 5 3

4 which hs solution S x = i S 0 x S 0x dx = i log px + constnt. 6 Therefore the generl solution to the Schrödinger eqution up to this order is ψx, t = e is 0x+ hs x/ h e iet/ h = c ± px exp ī me V x / h dx e iet/ h, 7 nd the overll constnt c is of course undetermined from this nlysis. This solution mkes it immeditely cler tht this pproximtion breks down when px goes to zero. Or in other words, the pproximtion is bd where the clssicl prticle stops nd turns becuse of the potentil. Such points re clled clssicl turning points.. Vlidity of the WKB Approximtion The pproximtion to stop with S in the h expnsion is vlid only when S is much smller thn S 0. Or in other words, if the term with h in Eq. is much smller thn the other terms. In prticulr, we require S h S. 8 In the one-diimensionl time-independent cse discussed bove, this is px h p x. 9 Using the definition of px = ± me V x, we find hdv x/dx. 0 E V xpx Agin we find the sme conclusion: the WKB pproximtion breks down close to the clssicl turning point V x = E e.g., px = 0. For exmple, tke hrmonic oscilltor V x = mω x. The vlidity condition Eq. 0 cn be rewritten s 8 E mω x 3 hω mω x.

5 This inequlity is lwys stisfied exctly t the origin x = 0, but once wy from the origin, it is impossible to stisfy unless E hω. In this sense, we re indeed in the clssicl regime. However, even for lrge E hω, the pproximtion is not vlid close to the clssicl turning points E = mω x. Here is the surprise. The vlidity condition Eq. 0 my be stisfied even in the region where the prticle cnnot enter clssicly E < V x. For exmple with the hrmonic oscilltor gin, the vlidity condition is lwys stisfied for lrge x E/mω for ny vlue of E. In other words, the WKB pproximtion is good wy from the clssicl turning points both where clssicl prticle exists nd where clssicl prticle cnnot exist. This is why the WKB pproximtion is not relly clssicl limit. It pplies lso where physics is truly quntum mechnicl. In the clssiclly forbidden region, the solution Eq. 7 needs to be modified to ψx, t = e is 0x+ hs x/ h e iet/ h = c ± mv x E exp h / mv x Edx e iet/ h, by following the sme steps s in the clssiclly llowed region.. Mtching WKB pproximtion cn be good both in the region E > V x nd the region E < V x but cnnot be good in between the regions close to the clssicl turning point E = V x c. In order to utilize the WKB pproximtion to work out wve functions, we need to somehow overcome this limittion. The stndrd method is to expnd round x c nd solve for the wve function exctly. Then you cn mtch on to the WKB solutions wy from x c to determine the entire wve function. The common method is to pproximte the potentil round the clssicl turning point x c by liner one: V x = V x c + V x c x x c + Ox x c, 3 nd ignore the second order term. By definition V x c = E. Schrödinger eqution round this point is therefore h d m dx + V x E ψ = h d m dx + V x c x x c ψ = 0. 5

6 Using the new vrible u = the differentil eqution simplifes drsticlly to /3 m dv h dx x c x x c, 5 d du u ψ = 0. 6 The solution to this eqution is known s the Airy function Aiu = π 0 dt cos 3 t3 + ut. 7 This cn be checked s follows. By cting the differentil opertor in Eq. 6 on the definition of the Airy function, we find d du u Aiu = π 0 dt d dt sin 3 t3 + ut. 8 The boundry term t t = 0 obviously vnishes. The behvior t t = is trickier. The point is tht the rgument of the sin grows s t 3 nd oscilltes more nd more rpidly s t. Therefore for ny infinitesiml intervl of lrge t, the oscilltion bsiclly cncels the integrnd except for left-over tht goes down s /t. Therefore the boundry term for t cn lso be dropped. It cn be shown using the steepest descent method tht the symptotic behvior of the Airy function smoothly mtches to the WKB solutions. The symptotic behvior is Aiu π u / exp 3 u3/ u 0. 9 u 0 π u / cos 3 u u + π Note first tht, for u 0, the symptotic behvior is Aiu = hmv / h /3 π me V / cos h x c me V dx + π. 30 6

7 Here we used the liner expnsion V x = E +V x c x x c to relte powers of u to me V. This expression is consistent with the WKB solution for prticulr choice of the overll constnt. Similrly, for u 0, Aiu = hmv / h /3 π mv E.3 Bound Sttes / exp h x c mv x Edx. 3 Consider the following sitution. In the region I x <, E < V nd it is clssiclly forbidden. In the region II < x < b, E > V nd the the prticle is clssiclly llowed. But gin in the region III x > b it is clssiclly forbidden. In this cse, we expect bound sttes with discrete energy levels. From the mtching using the Airy function t x = b, we wnt the wve function in the clssicly llowed region II to be s in Eq. 30, ψx = c hmv b/ h /3 / b cos me V x π me V h dx + π. b 3 On the other hnd from the mtching t x =, the wve function in the region II must be ψx = c hmv / h /3 / cos me V x π me V h dx π. 33 c b nd c re normliztion constnts. The minus sign in front of π/ is there becuse the clssiclly llowed region is to the right of the turning point nd hence we need to use Ai u insted of Aiu for mtching. The two behviors obtined from mtching t both sides must be the sme in order to hve consistent wve function for the entire region II. Becuse we do not know the coefficients c,b, which cn in prticulr differ in signs, we hve to require h b me V x dx + π = h me V x dx π nπ. 3 7

8 Further simplifying the eqution we obtin b me V xdx = n + π h. 35 This is reminiscent of Bohr Sommerfeld quntiztion condition tht ppered t the erly stge of quntum mechnics, pdq = nh 36 with h = π h except tht there is n extr contribution / here. Remember tht the Bohr Sommerfeld quntiztion condition ws n d hoc requirement which hppened to give correct energy spectr for the hydrogen-like toms. Nonetheless, this is n importnt eqution becuse it puts constrints on the llowed vlues of E, nd hence we find discrete energy levels. Becuse the WKB pproximtion is better for lrger E, the constrint Eq. 35 is expected to give correct energy levels t lest for highly excited sttes. Let us pply this constrint Eq. 35 for hrmonic oscilltor. The condition then is mx x min me m ω x dx = The l.h.s. is n elementry integrl me mω E = π n + n + π h. 37 nd we find hω. 38 This is precisely the energy levels of hrmonic oscilltor. Of course, obtining the exct energy levels is n ccident for hrmonic oscilltor. In generl, we expect to get pproximtely correct energy levels for high n. Indeed, if we pply the sme constrint for the hydrogen-like toms for l = 0 nd pplying the WKB method to the rdil wve function, we would obtin the correct results if we drop / this is wht Bohr Sommerfeld did, but don t if we keep / s obtined by the WKB method. Wht we obtin from the WKB method is not just the energy levels, but more detils such s the pproximte wve function. To simplify the nottion, let us choose unit system such tht m = ω = h =. Then the wve function in the region II for the hrmonic oscilltor is ψ E x / cos E E x dx π 8. 39

9 Doing the stndrd integrl, we find ψ cos E sin E x / x E + π + x x E E π. 0 Becuse E = n + / in this unit, it further simplifies to [ ] cos E sin x E + x E x E ψ [ E x / ] sin E sin x E + x E x E n even n odd. Tunneling To study tunneling process, we gin hve three regions, clssiclly llowed region I x <, where the prticle initilly exist, clssiclly forbidden region II < x < b, nd clssiclly llowed region III x > b where the prticle tunnels to. We follow the sme mtching procedure s in the bound stte exmple. A difference is tht the mtching t x = b uses different version of the Airy function which rises exponentilly wy from x = b. This version of the Airy function Bi is explined in the note on steepest descent method. The bottom line is tht the mtching t x = is done s ψ hmv / h /3 π me V ψ hmv / h /3 π mv E / cos h / x exp h me V dx + π, x < The mtching t x = b is, on the other hnd, ψ C i hmv b/ h /3 π mv E ψ C i hmv b/ h /3 π me V / / exp + h sin h b b x mv Edx. x > mv Edx, x < b me V dx + π. x > b The overll normliztion fctor C is determined by the requirement tht the behvior of the wve function is consistent for < x < b between two 9

10 mtching procedures. We therefore find C = i V /6 exp h V b b mv x Edx. Compring two clssiclly llowed regions, nd tking dvntge of further normliztion chnge, the mtching reduces to ψ me V cos / h x ψ exp h b mv x Edx me V dx + π, x < me V / sin h In other words, the mplitude in the region x > b due to tunneling from the region x < is suppressed by exp h b mv x Edx. 3 b me V dx + π. Usully people refer to Gmov s trnsmission coefficient exp h b mv x Edx s suppression fctor for the tunneling rte squre of the mplitude. 0

221B Lecture Notes WKB Method

221B Lecture Notes WKB Method Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using

More information

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner

More information

Solutions to Problems in Merzbacher, Quantum Mechanics, Third Edition. Chapter 7

Solutions to Problems in Merzbacher, Quantum Mechanics, Third Edition. Chapter 7 Solutions to Problems in Merzbcher, Quntum Mechnics, Third Edition Homer Reid April 5, 200 Chpter 7 Before strting on these problems I found it useful to review how the WKB pproimtion works in the first

More information

Name Solutions to Test 3 November 8, 2017

Name Solutions to Test 3 November 8, 2017 Nme Solutions to Test 3 November 8, 07 This test consists of three prts. Plese note tht in prts II nd III, you cn skip one question of those offered. Some possibly useful formuls cn be found below. Brrier

More information

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1 3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =

More information

Best Approximation. Chapter The General Case

Best Approximation. Chapter The General Case Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given

More information

Improper Integrals, and Differential Equations

Improper Integrals, and Differential Equations Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted

More information

Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon

Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon Phys463.nb 49 7 Energy Bnds Ref: textbook, Chpter 7 Q: Why re there insultors nd conductors? Q: Wht will hppen when n electron moves in crystl? In the previous chpter, we discussed free electron gses,

More information

Quantum Mechanics Qualifying Exam - August 2016 Notes and Instructions

Quantum Mechanics Qualifying Exam - August 2016 Notes and Instructions Quntum Mechnics Qulifying Exm - August 016 Notes nd Instructions There re 6 problems. Attempt them ll s prtil credit will be given. Write on only one side of the pper for your solutions. Write your lis

More information

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1 Tor Kjellsson Stockholm University Chpter 5 5. Strting with the following informtion: R = m r + m r m + m, r = r r we wnt to derive: µ = m m m + m r = R + µ m r, r = R µ m r 3 = µ m R + r, = µ m R r. 4

More information

p(t) dt + i 1 re it ireit dt =

p(t) dt + i 1 re it ireit dt = Note: This mteril is contined in Kreyszig, Chpter 13. Complex integrtion We will define integrls of complex functions long curves in C. (This is bit similr to [relvlued] line integrls P dx + Q dy in R2.)

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Statistical Physics I Spring Term Solutions to Problem Set #1

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Statistical Physics I Spring Term Solutions to Problem Set #1 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Deprtment 8.044 Sttisticl Physics I Spring Term 03 Problem : Doping Semiconductor Solutions to Problem Set # ) Mentlly integrte the function p(x) given in

More information

c n φ n (x), 0 < x < L, (1) n=1

c n φ n (x), 0 < x < L, (1) n=1 SECTION : Fourier Series. MATH4. In section 4, we will study method clled Seprtion of Vribles for finding exct solutions to certin clss of prtil differentil equtions (PDEs. To do this, it will be necessry

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties

More information

Quantum Physics III (8.06) Spring 2005 Solution Set 5

Quantum Physics III (8.06) Spring 2005 Solution Set 5 Quntum Physics III (8.06 Spring 005 Solution Set 5 Mrch 8, 004. The frctionl quntum Hll effect (5 points As we increse the flux going through the solenoid, we increse the mgnetic field, nd thus the vector

More information

The Regulated and Riemann Integrals

The Regulated and Riemann Integrals Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue

More information

Physics 137A - Quantum Mechanics - Spring 2018 Midterm 1. Mathematical Formulas

Physics 137A - Quantum Mechanics - Spring 2018 Midterm 1. Mathematical Formulas Copyright c 8 by Austin J. Hedemn Physics 7A - Quntum Mechnics - Spring 8 Midterm Mondy, Februry 6, 6:-8: PM You hve two hours, thirty minutes for the exm. All nswers should be written in blue book. You

More information

ODE: Existence and Uniqueness of a Solution

ODE: Existence and Uniqueness of a Solution Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct

More information

Lecture 14: Quadrature

Lecture 14: Quadrature Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl

More information

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1 Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution

More information

Problem Set 3 Solutions

Problem Set 3 Solutions Chemistry 36 Dr Jen M Stndrd Problem Set 3 Solutions 1 Verify for the prticle in one-dimensionl box by explicit integrtion tht the wvefunction ψ ( x) π x is normlized To verify tht ψ ( x) is normlized,

More information

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

1.2. Linear Variable Coefficient Equations. y + b ! = a y + b  Remark: The case b = 0 and a non-constant can be solved with the same idea as above. 1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt

More information

Abstract inner product spaces

Abstract inner product spaces WEEK 4 Abstrct inner product spces Definition An inner product spce is vector spce V over the rel field R equipped with rule for multiplying vectors, such tht the product of two vectors is sclr, nd the

More information

Continuous Quantum Systems

Continuous Quantum Systems Chpter 8 Continuous Quntum Systems 8.1 The wvefunction So fr, we hve been tlking bout finite dimensionl Hilbert spces: if our system hs k qubits, then our Hilbert spce hs n dimensions, nd is equivlent

More information

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as Improper Integrls Two different types of integrls cn qulify s improper. The first type of improper integrl (which we will refer to s Type I) involves evluting n integrl over n infinite region. In the grph

More information

Math& 152 Section Integration by Parts

Math& 152 Section Integration by Parts Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible

More information

Physics 202H - Introductory Quantum Physics I Homework #08 - Solutions Fall 2004 Due 5:01 PM, Monday 2004/11/15

Physics 202H - Introductory Quantum Physics I Homework #08 - Solutions Fall 2004 Due 5:01 PM, Monday 2004/11/15 Physics H - Introductory Quntum Physics I Homework #8 - Solutions Fll 4 Due 5:1 PM, Mondy 4/11/15 [55 points totl] Journl questions. Briefly shre your thoughts on the following questions: Of the mteril

More information

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl

More information

The Wave Equation I. MA 436 Kurt Bryan

The Wave Equation I. MA 436 Kurt Bryan 1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string

More information

Math 8 Winter 2015 Applications of Integration

Math 8 Winter 2015 Applications of Integration Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl

More information

7.2 The Definite Integral

7.2 The Definite Integral 7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where

More information

Quantum Physics I (8.04) Spring 2016 Assignment 8

Quantum Physics I (8.04) Spring 2016 Assignment 8 Quntum Physics I (8.04) Spring 206 Assignment 8 MIT Physics Deprtment Due Fridy, April 22, 206 April 3, 206 2:00 noon Problem Set 8 Reding: Griffiths, pges 73-76, 8-82 (on scttering sttes). Ohnin, Chpter

More information

Quantum Physics II (8.05) Fall 2013 Assignment 2

Quantum Physics II (8.05) Fall 2013 Assignment 2 Quntum Physics II (8.05) Fll 2013 Assignment 2 Msschusetts Institute of Technology Physics Deprtment Due Fridy September 20, 2013 September 13, 2013 3:00 pm Suggested Reding Continued from lst week: 1.

More information

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019 ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS MATH00030 SEMESTER 208/209 DR. ANTHONY BROWN 7.. Introduction to Integrtion. 7. Integrl Clculus As ws the cse with the chpter on differentil

More information

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS. THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem

More information

Chapter 0. What is the Lebesgue integral about?

Chapter 0. What is the Lebesgue integral about? Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous

More information

1 The Lagrange interpolation formula

1 The Lagrange interpolation formula Notes on Qudrture 1 The Lgrnge interpoltion formul We briefly recll the Lgrnge interpoltion formul. The strting point is collection of N + 1 rel points (x 0, y 0 ), (x 1, y 1 ),..., (x N, y N ), with x

More information

Math 1B, lecture 4: Error bounds for numerical methods

Math 1B, lecture 4: Error bounds for numerical methods Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the

More information

13: Diffusion in 2 Energy Groups

13: Diffusion in 2 Energy Groups 3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups

More information

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007 A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus

More information

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0) 1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this

More information

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one

More information

Review of basic calculus

Review of basic calculus Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below

More information

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1 MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further

More information

Physics 215 Quantum Mechanics 1 Assignment 2

Physics 215 Quantum Mechanics 1 Assignment 2 Physics 15 Quntum Mechnics 1 Assignment Logn A. Morrison Jnury, 16 Problem 1 Clculte p nd p on the Gussin wve pcket α whose wve function is x α = 1 ikx x 1/4 d 1 Solution Recll tht where ψx = x ψ. Additionlly,

More information

Review of Calculus, cont d

Review of Calculus, cont d Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some

More information

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2010 Homework Assignment 4; Due at 5p.m. on 2/01/10

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2010 Homework Assignment 4; Due at 5p.m. on 2/01/10 University of Wshington Deprtment of Chemistry Chemistry 45 Winter Qurter Homework Assignment 4; Due t 5p.m. on // We lerned tht the Hmiltonin for the quntized hrmonic oscilltor is ˆ d κ H. You cn obtin

More information

Classical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011

Classical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011 Clssicl Mechnics From Moleculr to Con/nuum Physics I WS 11/12 Emilino Ippoli/ October, 2011 Wednesdy, October 12, 2011 Review Mthemtics... Physics Bsic thermodynmics Temperture, idel gs, kinetic gs theory,

More information

AP Calculus Multiple Choice: BC Edition Solutions

AP Calculus Multiple Choice: BC Edition Solutions AP Clculus Multiple Choice: BC Edition Solutions J. Slon Mrch 8, 04 ) 0 dx ( x) is A) B) C) D) E) Divergent This function inside the integrl hs verticl symptotes t x =, nd the integrl bounds contin this

More information

( ) 2. ( ) is the Fourier transform of! ( x). ( ) ( ) ( ) = Ae i kx"#t ( ) = 1 2" ( )"( x,t) PC 3101 Quantum Mechanics Section 1

( ) 2. ( ) is the Fourier transform of! ( x). ( ) ( ) ( ) = Ae i kx#t ( ) = 1 2 ( )( x,t) PC 3101 Quantum Mechanics Section 1 1. 1D Schrödinger Eqution G chpters 3-4. 1.1 the Free Prticle V 0 "( x,t) i = 2 t 2m x,t = Ae i kxt "( x,t) x 2 where = k 2 2m. Normliztion must hppen: 2 x,t = 1 Here, however: " A 2 dx " " As this integrl

More information

Massachusetts Institute of Technology Quantum Mechanics I (8.04) Spring 2005 Solutions to Problem Set 6

Massachusetts Institute of Technology Quantum Mechanics I (8.04) Spring 2005 Solutions to Problem Set 6 Msschusetts Institute of Technology Quntum Mechnics I (8.) Spring 5 Solutions to Problem Set 6 By Kit Mtn. Prctice with delt functions ( points) The Dirc delt function my be defined s such tht () (b) 3

More information

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a). The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples

More information

Method of stationary phase

Method of stationary phase Physics 4 Spring 16 Method of sttionry phse Lecture notes by M. G. Rozmn Lst modified: April 13, 16 There is n immedite generliztion of the Lplce integrls f t)e φt) dt 1) which we obtin by llowing the

More information

Lecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature

Lecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 6: Singulr Integrls, Open Qudrture rules, nd Guss Qudrture (Compiled 6 August 7) In this lecture we discuss the

More information

Practice Problems Solution

Practice Problems Solution Prctice Problems Solution Problem Consier D Simple Hrmonic Oscilltor escribe by the Hmiltonin Ĥ ˆp m + mwˆx Recll the rte of chnge of the expecttion of quntum mechnicl opertor t A ī A, H] + h A t. Let

More information

Review of Gaussian Quadrature method

Review of Gaussian Quadrature method Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge

More information

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response

More information

Lecture 3. In this lecture, we will discuss algorithms for solving systems of linear equations.

Lecture 3. In this lecture, we will discuss algorithms for solving systems of linear equations. Lecture 3 3 Solving liner equtions In this lecture we will discuss lgorithms for solving systems of liner equtions Multiplictive identity Let us restrict ourselves to considering squre mtrices since one

More information

Definition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim

Definition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim Mth 9 Course Summry/Study Guide Fll, 2005 [1] Limits Definition of Limit: We sy tht L is the limit of f(x) s x pproches if f(x) gets closer nd closer to L s x gets closer nd closer to. We write lim f(x)

More information

Numerical Analysis: Trapezoidal and Simpson s Rule

Numerical Analysis: Trapezoidal and Simpson s Rule nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =

More information

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying

W. We shall do so one by one, starting with I 1, and we shall do it greedily, trying Vitli covers 1 Definition. A Vitli cover of set E R is set V of closed intervls with positive length so tht, for every δ > 0 nd every x E, there is some I V with λ(i ) < δ nd x I. 2 Lemm (Vitli covering)

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

1 The Riemann Integral

1 The Riemann Integral The Riemnn Integrl. An exmple leding to the notion of integrl (res) We know how to find (i.e. define) the re of rectngle (bse height), tringle ( (sum of res of tringles). But how do we find/define n re

More information

1.9 C 2 inner variations

1.9 C 2 inner variations 46 CHAPTER 1. INDIRECT METHODS 1.9 C 2 inner vritions So fr, we hve restricted ttention to liner vritions. These re vritions of the form vx; ǫ = ux + ǫφx where φ is in some liner perturbtion clss P, for

More information

Physics Graduate Prelim exam

Physics Graduate Prelim exam Physics Grdute Prelim exm Fll 2008 Instructions: This exm hs 3 sections: Mechnics, EM nd Quntum. There re 3 problems in ech section You re required to solve 2 from ech section. Show ll work. This exm is

More information

Aike ikx Bike ikx. = 2k. solving for. A = k iκ

Aike ikx Bike ikx. = 2k. solving for. A = k iκ LULEÅ UNIVERSITY OF TECHNOLOGY Division of Physics Solution to written exm in Quntum Physics F0047T Exmintion dte: 06-03-5 The solutions re just suggestions. They my contin severl lterntive routes.. Sme/similr

More information

Lecture 1. Functional series. Pointwise and uniform convergence.

Lecture 1. Functional series. Pointwise and uniform convergence. 1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is

More information

Green s functions. f(t) =

Green s functions. f(t) = Consider the 2nd order liner inhomogeneous ODE Green s functions d 2 u 2 + k(t)du + p(t)u(t) = f(t). Of course, in prctice we ll only del with the two prticulr types of 2nd order ODEs we discussed lst

More information

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2009

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2009 University of Wshington Deprtment of Chemistry Chemistry Winter Qurter 9 Homework Assignment ; Due t pm on //9 6., 6., 6., 8., 8. 6. The wve function in question is: ψ u cu ( ψs ψsb * cu ( ψs ψsb cu (

More information

Do the one-dimensional kinetic energy and momentum operators commute? If not, what operator does their commutator represent?

Do the one-dimensional kinetic energy and momentum operators commute? If not, what operator does their commutator represent? 1 Problem 1 Do the one-dimensionl kinetic energy nd momentum opertors commute? If not, wht opertor does their commuttor represent? KE ˆ h m d ˆP i h d 1.1 Solution This question requires clculting the

More information

Polynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230

Polynomial Approximations for the Natural Logarithm and Arctangent Functions. Math 230 Polynomil Approimtions for the Nturl Logrithm nd Arctngent Functions Mth 23 You recll from first semester clculus how one cn use the derivtive to find n eqution for the tngent line to function t given

More information

SYDE 112, LECTURES 3 & 4: The Fundamental Theorem of Calculus

SYDE 112, LECTURES 3 & 4: The Fundamental Theorem of Calculus SYDE 112, LECTURES & 4: The Fundmentl Theorem of Clculus So fr we hve introduced two new concepts in this course: ntidifferentition nd Riemnn sums. It turns out tht these quntities re relted, but it is

More information

Bernoulli Numbers Jeff Morton

Bernoulli Numbers Jeff Morton Bernoulli Numbers Jeff Morton. We re interested in the opertor e t k d k t k, which is to sy k tk. Applying this to some function f E to get e t f d k k tk d k f f + d k k tk dk f, we note tht since f

More information

x = b a n x 2 e x dx. cdx = c(b a), where c is any constant. a b

x = b a n x 2 e x dx. cdx = c(b a), where c is any constant. a b CHAPTER 5. INTEGRALS 61 where nd x = b n x i = 1 (x i 1 + x i ) = midpoint of [x i 1, x i ]. Problem 168 (Exercise 1, pge 377). Use the Midpoint Rule with the n = 4 to pproximte 5 1 x e x dx. Some quick

More information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn

More information

3.4 Numerical integration

3.4 Numerical integration 3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,

More information

UNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3

UNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3 UNIFORM CONVERGENCE Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3 Suppose f n : Ω R or f n : Ω C is sequence of rel or complex functions, nd f n f s n in some sense. Furthermore,

More information

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004 Advnced Clculus: MATH 410 Notes on Integrls nd Integrbility Professor Dvid Levermore 17 October 2004 1. Definite Integrls In this section we revisit the definite integrl tht you were introduced to when

More information

4 7x =250; 5 3x =500; Read section 3.3, 3.4 Announcements: Bell Ringer: Use your calculator to solve

4 7x =250; 5 3x =500; Read section 3.3, 3.4 Announcements: Bell Ringer: Use your calculator to solve Dte: 3/14/13 Objective: SWBAT pply properties of exponentil functions nd will pply properties of rithms. Bell Ringer: Use your clcultor to solve 4 7x =250; 5 3x =500; HW Requests: Properties of Log Equtions

More information

Math 124A October 04, 2011

Math 124A October 04, 2011 Mth 4A October 04, 0 Viktor Grigoryn 4 Vibrtions nd het flow In this lecture we will derive the wve nd het equtions from physicl principles. These re second order constnt coefficient liner PEs, which model

More information

Lecture 19: Continuous Least Squares Approximation

Lecture 19: Continuous Least Squares Approximation Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for

More information

MAA 4212 Improper Integrals

MAA 4212 Improper Integrals Notes by Dvid Groisser, Copyright c 1995; revised 2002, 2009, 2014 MAA 4212 Improper Integrls The Riemnn integrl, while perfectly well-defined, is too restrictive for mny purposes; there re functions which

More information

Theoretical foundations of Gaussian quadrature

Theoretical foundations of Gaussian quadrature Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of

More information

Continuous Random Variables

Continuous Random Variables STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 217 Néhémy Lim Continuous Rndom Vribles Nottion. The indictor function of set S is rel-vlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht

More information

1 Probability Density Functions

1 Probability Density Functions Lis Yn CS 9 Continuous Distributions Lecture Notes #9 July 6, 28 Bsed on chpter by Chris Piech So fr, ll rndom vribles we hve seen hve been discrete. In ll the cses we hve seen in CS 9, this ment tht our

More information

8 Laplace s Method and Local Limit Theorems

8 Laplace s Method and Local Limit Theorems 8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved

More information

INTRODUCTION TO INTEGRATION

INTRODUCTION TO INTEGRATION INTRODUCTION TO INTEGRATION 5.1 Ares nd Distnces Assume f(x) 0 on the intervl [, b]. Let A be the re under the grph of f(x). b We will obtin n pproximtion of A in the following three steps. STEP 1: Divide

More information

p-adic Egyptian Fractions

p-adic Egyptian Fractions p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction

More information

APPROXIMATE INTEGRATION

APPROXIMATE INTEGRATION APPROXIMATE INTEGRATION. Introduction We hve seen tht there re functions whose nti-derivtives cnnot be expressed in closed form. For these resons ny definite integrl involving these integrnds cnnot be

More information

Numerical integration

Numerical integration 2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter

More information

Calculus of Variations

Calculus of Variations Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function

More information

Week 10: Line Integrals

Week 10: Line Integrals Week 10: Line Integrls Introduction In this finl week we return to prmetrised curves nd consider integrtion long such curves. We lredy sw this in Week 2 when we integrted long curve to find its length.

More information

Chapter 3 The Schrödinger Equation and a Particle in a Box

Chapter 3 The Schrödinger Equation and a Particle in a Box Chpter 3 The Schrödinger Eqution nd Prticle in Bo Bckground: We re finlly ble to introduce the Schrödinger eqution nd the first quntum mechnicl model prticle in bo. This eqution is the bsis of quntum mechnics

More information

Section 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40

Section 5.1 #7, 10, 16, 21, 25; Section 5.2 #8, 9, 15, 20, 27, 30; Section 5.3 #4, 6, 9, 13, 16, 28, 31; Section 5.4 #7, 18, 21, 23, 25, 29, 40 Mth B Prof. Audrey Terrs HW # Solutions by Alex Eustis Due Tuesdy, Oct. 9 Section 5. #7,, 6,, 5; Section 5. #8, 9, 5,, 7, 3; Section 5.3 #4, 6, 9, 3, 6, 8, 3; Section 5.4 #7, 8,, 3, 5, 9, 4 5..7 Since

More information

Sufficient condition on noise correlations for scalable quantum computing

Sufficient condition on noise correlations for scalable quantum computing Sufficient condition on noise correltions for sclble quntum computing John Presill, 2 Februry 202 Is quntum computing sclble? The ccurcy threshold theorem for quntum computtion estblishes tht sclbility

More information

Chapter 28. Fourier Series An Eigenvalue Problem.

Chapter 28. Fourier Series An Eigenvalue Problem. Chpter 28 Fourier Series Every time I close my eyes The noise inside me mplifies I cn t escpe I relive every moment of the dy Every misstep I hve mde Finds wy it cn invde My every thought And this is why

More information

1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ),

1. Gauss-Jacobi quadrature and Legendre polynomials. p(t)w(t)dt, p {p(x 0 ),...p(x n )} p(t)w(t)dt = w k p(x k ), 1. Guss-Jcobi qudrture nd Legendre polynomils Simpson s rule for evluting n integrl f(t)dt gives the correct nswer with error of bout O(n 4 ) (with constnt tht depends on f, in prticulr, it depends on

More information