Separation of Variables

Size: px
Start display at page:

Download "Separation of Variables"

Transcription

1 Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical implications of Schröinger s equation, which will take up the rest of the semester. So we start with a review of: Linear Algebra, Separation of Variables (SOV), an Probability. There will be no particular completeness to our iscussions at this stage, I want to emphasize those aspects of each of these that will prove useful. Pitfalls an caveats will be aresse as we encounter them in actual problems (where they won t show up, so won t be aresse). Try as one might, it is ifficult to relate these three areas, so we will take each one in turn, starting from the most familiar, an working to the least (that s a matter of taste, so apologies in avance for the wrong orering in anticipation of our section on probability: How many ifferent orerings of these three subjects are there? On average, then, how many people isagree with the current one?). We will begin with separation of variables, a simple technique that you have use recently we will review the basic electrostatic interest in these solutions, an o a few examples of SOV applie to PDE s other than Laplace s equation ( 2 V = 0). 1.1 Separation for the Laplace Problem Separation of variables refers to a class of techniques for probing solutions to partial ifferential equations (PDEs) by turning them into orinary ifferential equations (ODEs). We generally rely on some notion of uniqueness 1 for 1 Uniqueness implying both a solution to the PDE an satisfaction of some bounary conitions. 1 of 11

2 1.1. SEPARATION FOR THE LAPLACE PROBLEM Lecture 1 our PDE then the logic of SOV rests on the notion that fining a solution, even by this somewhat limite technique, suffices. There are two basic flavors of SOV, aitive an multiplicative in either case, the style of argument is the same, so we ll start with the algebraically simpler version Aitive Separation We take the Laplace problem: Fin V (x, y) satisfying 2 V = 0 on some region, with some specifie bounary conitions (for example, in E&M, we woul have an implicit V 0 at spatial infinity conition, an maybe some conucting bounaries where we set V = 0). Working in two imensions (for simplicity), we start with a separation ansatz aitive, in this case: then running this through the Laplacian gives: V (x, y) = V x (x) + V y (y), (1.1) 2 V x x V y = 0. (1.2) y2 Now for the stanar argument: V x (x) epens only on x just as V y (y) epens only on y they cannot be equal unless each is equal to a constant. Suppose we set V x = α, then we must have V y = α. The solution to these ODEs is: V x (x) = 1 2 α x2 + β x + δ V y (y) = 1 2 α y2 + γ y + ρ. (1.3) for arbitrary constants (β, δ, γ, ρ), to be use in imposing the bounary conition(s). The solution, for any α, can be written as: V (x, y) = 1 2 α (x2 y 2 ) + β x + γ y + κ (1.4) As far as electrostatics goes, this solution oes not match our usual notion of potential it oesn t ie at spatial infinity. That s not always a eal-breaker, but it is true that aitive separation is of limite utility in E&M (it plays a much larger role in, for example, Hamilton-Jacobi theory). Nevertheless, the pattern of: ansatz, argument, solution is the same as for the more immeiately useful, an algebraically involve multiplicative separation. 2 of 11

3 1.1. SEPARATION FOR THE LAPLACE PROBLEM Lecture Multiplicative Separation The ifference here is in the starting ansatz rather than aitive single variable functions, we take the multiplicative combination: then Laplace s equation becomes V (x, y) = V x (x) V y (y), (1.5) 2 V x x 2 V 2 V y y + V x = 0, (1.6) y2 an if we ivie by V (x, y) itself, we can write this in a form conucive to making the separation argument: V x + V y = 0. (1.7) V x V y The first term epens only on x, the secon only on y, so for a solution satisfying our functional assumptions, we must have both terms equal to a constant, call it α 2 (a convenient choice, as you will see in a moment) for α complex. Then V x = α 2 = V y (1.8) V x V y will certainly have 2 V = 0. The solutions here are familiar: V x = A e α x + B e α x V y = F e i α y + G e i α y. (1.9) Now we see that for a generic complex number, both V x an V y will be combinations of growing an ecaying exponentials, as well as oscillatory sines an cosines. For example, if α is real, V x is growing an ecaying, while V y is pure sine an cosine. For α = i a with a real, the roles are reverse. In general, we use the bounary conitions to choose a convenient expression for α. Example Consier a pair of infinite, groune conucting sheets separate a istance with a conuctor connecting the two sheets hel at V 0 as shown in Figure 1.1. What is the potential in between the plates? We know that the 3 of 11

4 1.1. SEPARATION FOR THE LAPLACE PROBLEM Lecture 1 potential satisfies Laplace s equation in the region between the plates (no charge in there), an the bounary conitions are clear: ( ) 2 π x V (x = 0, y) = V (x =, y) = 0 V (x, y = 0) = V 0 sin. (1.10) ŷ ˆx V 0 sin ( ) 2πx Figure 1.1: Two parallel plates separate by a istance an hel at V = 0, the connecting plate is hel at V 0 sin(2 π x/). We know the solution will be of the form (1.9), but how to choose α an the constants (A, B, F, G)? There is an aitional, implicit bounary conition we like the potential to go to zero in the open spatial irection, y this tells us that we shoul set α = i a for a IR to get growing an ecaying exponentials for the V y (y) solution: V x = A e i a x + B e i a x V y = F e a y + G e a y (1.11) an, furthermore, we shoul set G = 0 to eliminate the growing exponential. By introucing à = (A+B) an B = i (A B), we can write V x (x) in terms of sine an cosine. Our solution so far, then, is: V x = à cos(a x) + B sin(a x) V y = F e a y V = V x V y. (1.12) Now for the rest of the bounary conitions: take V x (x = 0) = 0 that tells us that à = 0. Our solution is now much simpler, the full potential has the 4 of 11

5 1.1. SEPARATION FOR THE LAPLACE PROBLEM Lecture 1 form V = Q sin(a x) e a y (1.13) where we have introuce yet another rewrite of combinations of constants (Q = Ã F ). The bounary conition V (x =, y) = 0 is interesting, this tells us that sin(a ) = 0, not a statement about the obvious constant Q, but rather, a constraint on the allowe wavelengths: sin(a ) = 0 a = n π (1.14) for integer n. Then the only a s that satisfy the bounary conition are a = n π, still an infinite family, but labelle by the integer inex n. So, we have V n (x, y) = Q n sin(n π x/) e n π y/, (1.15) an we put the subscript n here to remin us that any integer n provies a solution. This type of bounary conition, one that limits the allowe wavelength, is what leas to the mathematical representation of quantization when applie to the separable solutions of Schröinger s equation. We have one bounary left because Laplace s equation is a linear PDE, sums of solutions are still solutions, an we can make a general solution out of our n-inexe ones. Each term satisfies Laplace s equation an the bounary conitions at x = 0, an y, so the sum oes as well: V (x, y) = ( n π x ) Q n sin n=1 e n π y. (1.16) As for the final bounary conition: V (x, y = 0) = V 0 sin(2 π x/), we can take Q n = 0 for all n except for n = 2. We know that the solution is unique, so this is certainly vali. But it begs the question, what about a less wellaapte riving potential? That is a question for next time, as it is tie to the ecomposition of functions in sin an cosine bases, a linear algebra issue. The final answer: ( ) 2 π x V (x, y) = V 0 sin e 2 π y. (1.17) This solution is just the potential at the y = 0 plate, ecaying into the region between the plates equal contours are shown in Figure of 11

6 V 0 sin ( ) 2πx Figure 1.2: The potential between the plates for the setup shown in Figure 1.1. There is not much more to say about separation of variables the best approach is to o as many examples as you can get your hans on that makes some of the choices that appear unmotivate at first (why shoul we make the integration constant α 2? Why o we pick the y solution to be growing an ecaying?) more reasonable. 1.2 Separation for Other PDEs We will be applying the SOV technique to Schröinger s equation, an that is not going to be much ifferent than the Laplacian example from E&M. The approach is inepenent of PDE, although the PDE must be linear (not strictly speaking true, but the superposition principle we invoke in writing (1.16) relie on linearity) just to show some of the places where SOV comes up, we can look at the Helmholtz equation an the wave equation. 6 of 11

7 1.2.1 Helmholtz Equation Suppose we consier the same problem as above, but for the Helmholtz PDE: ( 2 + µ 2 ) V = 0 ( ) 2 π x V (x, y = 0, ) = 0 V (x, y = 0) = V 0 sin V (x, y ) 0, (1.18) we are imagining a slight moification to electrostatics (actually amounts to a theory of electrostatics with a massive photon, but forget about that for now). First we ll fin the ODE solutions to the Helmholtz PDE via V (x, y) = V x (x) V y (y): V x + V y + µ 2 = 0. (1.19) V x V y We see immeiately that we will again get oscillatory an exponential solutions, since we want the solution in the increasing y irection to ecay, we ll call V y V y = α 2, with α real. With the constant µ 2 in place, we on t nee to set the x term equal to α 2, we have aitional freeom. We know that the bounary conitions for x are easily impose for sines, so set V x V x = β 2 for β real. Then the constraint we must impose to solve the Helmholtz PDE is: β 2 + α 2 + µ 2 = 0. Imposing the bounary conitions at x = 0, gives back β = n π as before, so that ( n π x ) V x = A sin. (1.20) Now, for the V y equation, we are given µ, an we have β 2 = n 2 π 2 / 2, so the growing an ecaying exponentials have the form: V y = F e n 2 π 2 / 2 µ 2 y + G e n 2 π 2 / 2 µ 2 y. (1.21) There are interesting elements to this solution, most notably, a cutoff in n for which V y is actually oscillatory (n < µ π ) ignoring that for the moment, if we assume the above terms for V y are inee growing an ecaying for n = 2, we have the solution: ( ) 2 π x V (x, y) = V 0 sin e 4 π 2 / 2 µ 2 y, (1.22) which just has a ifferent funamental ecay length when compare to the Laplace form for electrostatics. 7 of 11

8 1.2.2 The Wave Equation Separation of variables is not limite to purely spatial problems we can mix in time. For the wave equation in one spatial imension (call it x), efine in a meium with funamental spee v, we have 1 v 2 2 u(x, t) t u(x, t) x 2 = 0. (1.23) Take a multiplicative ansatz: u(x, t) = u x (x) u t (t), then by the now familiar proceure of inputting an iviing by u: u t an again, setting both terms equal to a constant: u t v 2 + u x = 0 (1.24) u t u x v 2 = α 2 u x = α 2. (1.25) u t u x What is interesting here is that both the spatial solution u x an the temporal u t can be oscillatory that is not true for the Laplace equation separable solutions we saw above. Without worrying about bounary conitions, we have the easy solution: u t = A cos(α v t) u x = B cos(α x) (1.26) so we can write u(x, t) = u 0 cos(α x) cos(α v t) of course, just as obvious is the solution û(x, t) = u 0 sin(α x) sin(α v t). Since the wave equation is linear, the sum of these two solutions is also a solution, an it is an instructive one: u(x, t) + û(x, t) = u 0 (cos(α x) cos(α v t) + u 0 sin(α x) sin(α v t)) = u 0 cos(α (x v t)), (1.27) precisely a plane wave. Of course, in general, we woul be given spatial bounary conitions an/or an initial waveform at t = Poisson s Equation The utility of separation is not limite to source-free equations consier Poisson s equation for electrostatic potential in the presence of source charge ensity ρ: 2 V = ρ ɛ 0. (1.28) 8 of 11

9 We must be given a ρ that is itself appropriately separable. The simplest possible such istribution woul be a constant, but we can imagine more interesting charge istributions. As an example, suppose we are in spherical coorinates, an we have a spherically symmetric charge ensity: ρ(r, θ, φ) = ρ(r), epening only on our istance from the origin. Poisson s equation becomes, in these variables ( 1 r 2 r ( r 2 V r ) + 1 r 2 sin θ ( sin θ V ) + θ θ 1 r 2 sin 2 θ 2 ) V φ 2 V = ρ(r). ɛ 0 (1.29) The separation ansatz now becomes V (r, θ, φ) = V r (r) V θ (θ) V φ (φ), an we can see what will happen to satisfy the functional epenence assumption, we will have constant V θ an V φ, leaving us with just the raial ODE: 1 r 2 r ( r 2 V r ) ρ(r) =. (1.30) ɛ 0 Suppose we take ρ(r) = ρ 0 a constant, then the solution to the above is: V r (r) = α + β r ρ 0 r 2 6 ɛ 0. (1.31) Notice that the first two terms are just solutions to 2 V = 0 (an overall constant an the potential outsie a spherically symmetric istribution of charge) if we are insie a uniformly charge sphere, with r = 0 inclue in the omain, then we must set β = 0. The constant α, of course, can be set to zero, an we are left with the usual potential for a uniformly charge sphere. 9 of 11

10 Homework Reaing: Griffiths Introuction to Electroynamics, pp Problem 1.1 Solve Laplace s equation 2 V = 0 using separation: V (x, y) = V x (x) V y (y) for the following two-imensional bounary conitions: V (0, y) = V (, y) = 0, V (x, ) = 0 an V (x, 0) = V 0 sin ( ) 2 π x. y (V = 0) (V = 0) (V = 0) x ( ) 2 π x V 0 sin Physically, we are escribing the electrostatic potential insie a square that has groune lines on three sies, an a spatially-varying potential on the fourth sie. Problem 1.2 The heat equation in one (spatial) imension reas 1 k u(x, t) t = 2 u(x, t) x 2, (1.32) where u(x, t) is the temperature of an object with thermal conuctivity k (a material property). For a complete solution, we must specify two bounary conitions (since this equation is secon orer in x) an an initial conition (the value of u(x, t = 0), for example). 10 of 11

11 a. Using separation of variables: u(x, t) = u x (x) u t (t), an separation constant α 2, write own the separate heat equation an solve for u x an u t. b. Suppose we have a material that goes from x = 0 to x = we set the temperature on both ens of the material to zero i.e. u(0, t) = u(, t) = 0. Initially, at time t = 0, we have temperature: u(x, 0) = u 0 sin ( ) π x. Using your separation solution, fin the particular solution for this set of bounary an initial conitions. u(x, 0) u 0 ( π x ) u 0 sin x 11 of 11

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7. Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency

More information

12.11 Laplace s Equation in Cylindrical and

12.11 Laplace s Equation in Cylindrical and SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

Partial Differential Equations

Partial Differential Equations Chapter Partial Differential Equations. Introuction Have solve orinary ifferential equations, i.e. ones where there is one inepenent an one epenent variable. Only orinary ifferentiation is therefore involve.

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

Solving the Schrödinger Equation for the 1 Electron Atom (Hydrogen-Like)

Solving the Schrödinger Equation for the 1 Electron Atom (Hydrogen-Like) Stockton Univeristy Chemistry Program, School of Natural Sciences an Mathematics 101 Vera King Farris Dr, Galloway, NJ CHEM 340: Physical Chemistry II Solving the Schröinger Equation for the 1 Electron

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Lecture 2 - First order linear PDEs and PDEs from physics

Lecture 2 - First order linear PDEs and PDEs from physics 18.15 - Introuction to PEs, Fall 004 Prof. Gigliola Staffilani Lecture - First orer linear PEs an PEs from physics I mentione in the first class some basic PEs of first an secon orer. Toay we illustrate

More information

PDE Notes, Lecture #11

PDE Notes, Lecture #11 PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

More information

Second Major Solution Q1. The three capacitors in the figure have an equivalent capacitance of 2.77 µf. What is C 2?

Second Major Solution Q1. The three capacitors in the figure have an equivalent capacitance of 2.77 µf. What is C 2? Secon Major Solution Q1. The three capacitors in the figure have an equivalent capacitance of.77 µf. What is C? C 4.0 µf.0 µf A) 7 µf B) µf C) 4 µf D) 3 µf E) 6 µf Q. When the potential ifference across

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

5-4 Electrostatic Boundary Value Problems

5-4 Electrostatic Boundary Value Problems 11/8/4 Section 54 Electrostatic Bounary Value Problems blank 1/ 5-4 Electrostatic Bounary Value Problems Reaing Assignment: pp. 149-157 Q: A: We must solve ifferential equations, an apply bounary conitions

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

The Ehrenfest Theorems

The Ehrenfest Theorems The Ehrenfest Theorems Robert Gilmore Classical Preliminaries A classical system with n egrees of freeom is escribe by n secon orer orinary ifferential equations on the configuration space (n inepenent

More information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1 Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of

More information

6 Wave equation in spherical polar coordinates

6 Wave equation in spherical polar coordinates 6 Wave equation in spherical polar coorinates We now look at solving problems involving the Laplacian in spherical polar coorinates. The angular epenence of the solutions will be escribe by spherical harmonics.

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2 Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

Homework 7 Due 18 November at 6:00 pm

Homework 7 Due 18 November at 6:00 pm Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine

More information

Many problems in physics, engineering, and chemistry fall in a general class of equations of the form. d dx. d dx

Many problems in physics, engineering, and chemistry fall in a general class of equations of the form. d dx. d dx Math 53 Notes on turm-liouville equations Many problems in physics, engineering, an chemistry fall in a general class of equations of the form w(x)p(x) u ] + (q(x) λ) u = w(x) on an interval a, b], plus

More information

1 Boas, p. 643, problem (b)

1 Boas, p. 643, problem (b) Physics 6C Solutions to Homework Set #6 Fall Boas, p. 643, problem 3.5-3b Fin the steay-state temperature istribution in a soli cyliner of height H an raius a if the top an curve surfaces are hel at an

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

Unit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule

Unit #6 - Families of Functions, Taylor Polynomials, l Hopital s Rule Unit # - Families of Functions, Taylor Polynomials, l Hopital s Rule Some problems an solutions selecte or aapte from Hughes-Hallett Calculus. Critical Points. Consier the function f) = 54 +. b) a) Fin

More information

1.4.3 Elementary solutions to Laplace s equation in the spherical coordinates (Axially symmetric cases) (Griffiths 3.3.2)

1.4.3 Elementary solutions to Laplace s equation in the spherical coordinates (Axially symmetric cases) (Griffiths 3.3.2) 1.4.3 Elementary solutions to Laplace s equation in the spherical coorinates (Axially symmetric cases) (Griffiths 3.3.) In the spherical coorinates (r, θ, φ), the Laplace s equation takes the following

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

there is no special reason why the value of y should be fixed at y = 0.3. Any y such that

there is no special reason why the value of y should be fixed at y = 0.3. Any y such that 25. More on bivariate functions: partial erivatives integrals Although we sai that the graph of photosynthesis versus temperature in Lecture 16 is like a hill, in the real worl hills are three-imensional

More information

ELECTRON DIFFRACTION

ELECTRON DIFFRACTION ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

ABCD42BEF F2 F8 5 4D658 CC89

ABCD42BEF F2 F8 5 4D658 CC89 ABCD BEF F F D CC Vetri Velan GSI, Physics 7B Miterm 2: Problem Solution. Outsie sphere, E looks like a point charge. E = The total charge on the sphere is Q sphere = ρ 4 3 πr3 Thus, outsie the sphere,

More information

Math 1B, lecture 8: Integration by parts

Math 1B, lecture 8: Integration by parts Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores

More information

Convergence of Random Walks

Convergence of Random Walks Chapter 16 Convergence of Ranom Walks This lecture examines the convergence of ranom walks to the Wiener process. This is very important both physically an statistically, an illustrates the utility of

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Bair, faculty.uml.eu/cbair University of Massachusetts Lowell 1. Pre-Einstein Relativity - Einstein i not invent the concept of relativity,

More information

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0.

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0. Engineering Mathematics 2 26 February 2014 Limits of functions Consier the function 1 f() = 1. The omain of this function is R + \ {1}. The function is not efine at 1. What happens when is close to 1?

More information

Electromagnetism Answers to Problem Set 3 Spring Jackson Prob. 2.1: Charge above a grounded plane. (z d)

Electromagnetism Answers to Problem Set 3 Spring Jackson Prob. 2.1: Charge above a grounded plane. (z d) Electromagnetism 76 Answers to Problem Set 3 Spring 6. Jackson Prob..: Charge above a groune plane (a) Surface charge ensity E z (ρ, z) = φ z = q [ (z ) [ρ + (z ) (z + ) 3/ [ρ + (z + ) 3/ Evaluate E z

More information

1 Heisenberg Representation

1 Heisenberg Representation 1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.

More information

Introduction to variational calculus: Lecture notes 1

Introduction to variational calculus: Lecture notes 1 October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational

More information

Math 1271 Solutions for Fall 2005 Final Exam

Math 1271 Solutions for Fall 2005 Final Exam Math 7 Solutions for Fall 5 Final Eam ) Since the equation + y = e y cannot be rearrange algebraically in orer to write y as an eplicit function of, we must instea ifferentiate this relation implicitly

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

Some Examples. Uniform motion. Poisson processes on the real line

Some Examples. Uniform motion. Poisson processes on the real line Some Examples Our immeiate goal is to see some examples of Lévy processes, an/or infinitely-ivisible laws on. Uniform motion Choose an fix a nonranom an efine X := for all (1) Then, {X } is a [nonranom]

More information

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

Department of Physics University of Maryland College Park, Maryland. Fall 2005 Final Exam Dec. 16, u 2 dt )2, L = m u 2 d θ, ( d θ

Department of Physics University of Maryland College Park, Maryland. Fall 2005 Final Exam Dec. 16, u 2 dt )2, L = m u 2 d θ, ( d θ Department of Physics University of arylan College Park, arylan PHYSICS 4 Prof. S. J. Gates Fall 5 Final Exam Dec. 6, 5 This is a OPEN book examination. Rea the entire examination before you begin to work.

More information

Lagrangian and Hamiltonian Dynamics

Lagrangian and Hamiltonian Dynamics Lagrangian an Hamiltonian Dynamics Volker Perlick (Lancaster University) Lecture 1 The Passage from Newtonian to Lagrangian Dynamics (Cockcroft Institute, 22 February 2010) Subjects covere Lecture 2: Discussion

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

SYDE 112, LECTURE 1: Review & Antidifferentiation

SYDE 112, LECTURE 1: Review & Antidifferentiation SYDE 112, LECTURE 1: Review & Antiifferentiation 1 Course Information For a etaile breakown of the course content an available resources, see the Course Outline. Other relevant information for this section

More information

PARALLEL-PLATE CAPACITATOR

PARALLEL-PLATE CAPACITATOR Physics Department Electric an Magnetism Laboratory PARALLEL-PLATE CAPACITATOR 1. Goal. The goal of this practice is the stuy of the electric fiel an electric potential insie a parallelplate capacitor.

More information

Physics 115C Homework 4

Physics 115C Homework 4 Physics 115C Homework 4 Problem 1 a In the Heisenberg picture, the ynamical equation is the Heisenberg equation of motion: for any operator Q H, we have Q H = 1 t i [Q H,H]+ Q H t where the partial erivative

More information

fv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n

fv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n Chapter 11 Rossby waves Supplemental reaing: Pelosky 1 (1979), sections 3.1 3 11.1 Shallow water equations When consiering the general problem of linearize oscillations in a static, arbitrarily stratifie

More information

Armenian Transformation Equations For Relativity

Armenian Transformation Equations For Relativity Armenian Transformation Equations For Relativity Robert Nazaryan an Haik Nazaryan 00th Anniversary of the Special Relativity Theory This Research one in Armenia 968-988, Translate from the Armenian Manuscript

More information

Vectors in two dimensions

Vectors in two dimensions Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication

More information

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1 Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems

More information

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012 CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration

More information

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems

Construction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu

More information

Introduction to Markov Processes

Introduction to Markov Processes Introuction to Markov Processes Connexions moule m44014 Zzis law Gustav) Meglicki, Jr Office of the VP for Information Technology Iniana University RCS: Section-2.tex,v 1.24 2012/12/21 18:03:08 gustav

More information

The Three-dimensional Schödinger Equation

The Three-dimensional Schödinger Equation The Three-imensional Schöinger Equation R. L. Herman November 7, 016 Schröinger Equation in Spherical Coorinates We seek to solve the Schröinger equation with spherical symmetry using the metho of separation

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity AP Physics Multiple Choice Practice Electrostatics 1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity. A soli conucting sphere is given a positive charge Q.

More information

Exam #2, Electrostatics

Exam #2, Electrostatics Exam #2, Electrostatics Prof. Maurik Holtrop Department of Physics PHYS 408 University of New Hampshire March 27 th, 2003 Name: Stuent # NOTE: There are 5 questions. You have until 9 pm to finish. You

More information

Chapter 2 Derivatives

Chapter 2 Derivatives Chapter Derivatives Section. An Intuitive Introuction to Derivatives Consier a function: Slope function: Derivative, f ' For each, the slope of f is the height of f ' Where f has a horizontal tangent line,

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Calculus Class Notes for the Combined Calculus and Physics Course Semester I

Calculus Class Notes for the Combined Calculus and Physics Course Semester I Calculus Class Notes for the Combine Calculus an Physics Course Semester I Kelly Black December 14, 2001 Support provie by the National Science Founation - NSF-DUE-9752485 1 Section 0 2 Contents 1 Average

More information

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const.

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const. G4003 Avance Mechanics 1 The Noether theorem We alreay saw that if q is a cyclic variable, the associate conjugate momentum is conserve, q = 0 p q = const. (1) This is the simplest incarnation of Noether

More information

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS

THE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS

More information

Solution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010

Solution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010 NTNU Page of 6 Institutt for fysikk Fakultet for fysikk, informatikk og matematikk This solution consists of 6 pages. Solution to the exam in TFY423 STATISTICAL PHYSICS Wenesay ecember, 2 Problem. Particles

More information

constant ρ constant θ

constant ρ constant θ Answers to Problem Set Number 3 for 8.04. MIT (Fall 999) Contents Roolfo R. Rosales Λ Boris Schlittgen y Zhaohui Zhang ẓ October 9, 999 Problems from the book by Saff an Snier. 2. Problem 04 in section

More information

PH 132 Exam 1 Spring Student Name. Student Number. Lab/Recitation Section Number (11,,36)

PH 132 Exam 1 Spring Student Name. Student Number. Lab/Recitation Section Number (11,,36) PH 13 Exam 1 Spring 010 Stuent Name Stuent Number ab/ecitation Section Number (11,,36) Instructions: 1. Fill out all of the information requeste above. Write your name on each page.. Clearly inicate your

More information

SOLUTIONS for Homework #3

SOLUTIONS for Homework #3 SOLUTIONS for Hoework #3 1. In the potential of given for there is no unboun states. Boun states have positive energies E n labele by an integer n. For each energy level E, two syetrically locate classical

More information

Generalization of the persistent random walk to dimensions greater than 1

Generalization of the persistent random walk to dimensions greater than 1 PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the

More information

Review of Differentiation and Integration for Ordinary Differential Equations

Review of Differentiation and Integration for Ordinary Differential Equations Schreyer Fall 208 Review of Differentiation an Integration for Orinary Differential Equations In this course you will be expecte to be able to ifferentiate an integrate quickly an accurately. Many stuents

More information

February 21 Math 1190 sec. 63 Spring 2017

February 21 Math 1190 sec. 63 Spring 2017 February 21 Math 1190 sec. 63 Spring 2017 Chapter 2: Derivatives Let s recall the efinitions an erivative rules we have so far: Let s assume that y = f (x) is a function with c in it s omain. The erivative

More information

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com

More information

involve: 1. Treatment of a decaying particle. 2. Superposition of states with different masses.

involve: 1. Treatment of a decaying particle. 2. Superposition of states with different masses. Physics 195a Course Notes The K 0 : An Interesting Example of a Two-State System 021029 F. Porter 1 Introuction An example of a two-state system is consiere. involve: 1. Treatment of a ecaying particle.

More information

A Second Time Dimension, Hidden in Plain Sight

A Second Time Dimension, Hidden in Plain Sight A Secon Time Dimension, Hien in Plain Sight Brett A Collins. In this paper I postulate the existence of a secon time imension, making five imensions, three space imensions an two time imensions. I will

More information

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0.

x f(x) x f(x) approaching 1 approaching 0.5 approaching 1 approaching 0. Engineering Mathematics 2 26 February 2014 Limits of functions Consier the function 1 f() = 1. The omain of this function is R + \ {1}. The function is not efine at 1. What happens when is close to 1?

More information

MA 2232 Lecture 08 - Review of Log and Exponential Functions and Exponential Growth

MA 2232 Lecture 08 - Review of Log and Exponential Functions and Exponential Growth MA 2232 Lecture 08 - Review of Log an Exponential Functions an Exponential Growth Friay, February 2, 2018. Objectives: Review log an exponential functions, their erivative an integration formulas. Exponential

More information

Laplace s Equation in Cylindrical Coordinates and Bessel s Equation (II)

Laplace s Equation in Cylindrical Coordinates and Bessel s Equation (II) Laplace s Equation in Cylinrical Coorinates an Bessel s Equation (II Qualitative properties of Bessel functions of first an secon kin In the last lecture we foun the expression for the general solution

More information

Quantum mechanical approaches to the virial

Quantum mechanical approaches to the virial Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from

More information

The electrodynamics of rotating electrons.

The electrodynamics of rotating electrons. Die Elektroynamik es rotierenen Elektrons, Zeit. f. Phys. 37 (196), 43-6. The electroynamics of rotating electrons. By J. Frenkel 1 ) in Leningra. (Receive on May 196) Translate by D. H. Delphenich The

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

and from it produce the action integral whose variation we set to zero:

and from it produce the action integral whose variation we set to zero: Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

Sensors & Transducers 2015 by IFSA Publishing, S. L.

Sensors & Transducers 2015 by IFSA Publishing, S. L. Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity

More information

Physics Courseware Electromagnetism

Physics Courseware Electromagnetism Phsics Courseware Electromagnetism Electric potential Problem.- a) Fin the electric potential at points P, P an P prouce b the three charges Q, Q an Q. b) Are there an points where the electric potential

More information

Conservation Laws. Chapter Conservation of Energy

Conservation Laws. Chapter Conservation of Energy 20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action

More information