Average value of position for the anharmonic oscillator: Classical versus quantum results
|
|
- Tyler Dickerson
- 5 years ago
- Views:
Transcription
1 verage value of position for the anharmonic oscillator: Classical versus quantum results R. W. Robinett Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 682 Receive 8 January 996; accepte June 996 The evaluation of the average value of the position coorinate,, of a particle moving in a harmonic oscillator potential (V()k 2 /2) with a small anharmonic piece (V()k ) is a stanar calculation in classical Newtonian mechanics an statistical mechanics where the problem has relevance to thermal epansion. In each case, the calculation is most easily one using a perturbative epansion. In this note, we perform the same computation of in quantum mechanics using time-inepenent perturbation theory an the laer operator formalism to show how similar results are obtaine. We also inicate how a semiclassical calculation using a classical probability istribution can also be use to obtain the same result. 997 merican ssociation of Physics Teachers. I. INTRODUCTION The simple harmonic oscillator is one of the most often iscusse problems in all of classical, statistical, an quantum mechanics. In aition to its irect connections to important physical systems, one of its most useful features is that at every stage of evelopment, from elementary classical mechanics to quantum fiel theory, it is eactly soluble an so can be easily use as a close-form analytic eample. Many interesting physical phenomena, however, rely irectly on small eviations away from the simple form of the quaratic potential of the harmonic oscillator. s an eample, even introuctory tets 2 4 often stress the fact that the thermal epansion observe in many materials is ue to a lack of symmetry about the minimum in the intermolecular potential. Such a feature of realistic potentials can be moele by the aition of a small anharmonic term of the form V()k to the usual harmonic piece, V()k 2 /2. The average value of position,, in such an oscillator can then be calculate using classical Newtonian mechanics see, e.g., Ref. 5 or statistical mechanics see, e.g., Ref. 6 with the familiar result that increases linearly with amplitue square, energy, or temperature in contrast to the case of the pure oscillator where. Many tets on quantum mechanics consier the effects of small anharmonicities on the energy spectrum of the harmonic oscillator. quaratic term of the form V()a 4 is often iscusse see, e.g., Refs. 7 in the contet of perturbation theory where there is a nontrivial first-orer effect; such calculations are also often use as an illustration of the power of operator methos an raising/lowering operator techniques are use to evaluate the require matri elements. In contrast, the V()k term is presente see, e.g., Refs. as an eample where the first-orer correction vanishes ue to the symmetry of the problem an the perturbation an a secon-orer calculation is require. In a few cases, the first-orer correction to the wave function itself has also been eplicitly presente Refs. an, but no physical information is etracte from it. In this note, we eten these very familiar analyses to evaluate, using perturbation theory an operator techniques, the average value of position for the stanar oscillator problem perturbe by a small cubic anharmonic term an make comparisons with the classical an statistical mechanical results. In aition, we also present a erivation of the same effect using classical probability istributions, P CL (), which are also iscusse in many elementary tets an make a nice connection between the classical Newtonian solutions an the probabilistic quantum mechanical erivation using P QM ()() 2. In Sec. II we very briefly review the classical an statistical mechanical results an their relevance to thermal epansion, while Sec. III presents the main quantum mechanical an semiclassical results. II. CLSSICL RESULTS In classical Newtonian mechanics, a potential energy given by 9 m. J. Phys. 65, March merican ssociation of Physics Teachers 9
2 U 2k 2 k correspons to a force Fkk 2 so that Newton s law reas ẍt 2 t 2 2 t, where k/m as usual. perturbative solution of the form t t t 2 2 t, where (t) cos(t) is the lowest-orer result with amplitue, gives the ifferential equation ẍ t 2 t 2 cos 2 t for the O term. This can be easily solve for (t) using stanar techniques to give t cost cos2t. 6 The average value of position can be efine via t tt tt t using, for eample, as the perio of the unperturbe motion or in the limit that. In either case, the averaging process singles out the c component only, giving 2 2. In orer to make connections with thermal epansion, we note that the total energy of the unperturbe oscillator system is given by E tot k 2 /2. For the one-imensional oscillator with two quaratic egrees of freeom, this energy will also correspon to E2[k B T/2]k B T, where k B is Boltmann s constant. In this way we fin that Eq. 8 implies that k E k k BT, so that the average separation increases with temperature. This last result can be seen more irectly see, e.g., Ref. 6 by using the Boltmann istribution which is proportional to epu(), where /k B T in the thermoynamic average value given by e U e U k 4 e k2 /2 e k2 /2 k k k BT, where we have use the epansion e U e k2 /2 k III. QUNTUM MECHNICL ND SEMICLSSICL RESULTS In quantum mechanics, information on the average value of position is most easily etracte from the position-space stationary state wave functions, in this case labele n (). The average value of position in the unperturbe system is given by n n, 2 which vanishes because of the parity properties of the energy eigenfunctions. The first-orer shift in the energy eigenvalue also vanishes since EE n E n n V n k n n for the same reason. The first-orer wave function is given by where n n jn c jn j, c jn j V n E n E j V jn E n E j. 4 5 Using the anharmonic piece, V()k, as the perturbing potential, we require the values of V jn k j n. 6 To evaluate these, we use the formalism of raising an lowering operators which satisfy  n n n an  n n n in terms of which the position coorinate is given by  Â. The matri elements are given by 9,,2, /2 j n nnn2 j,n n /2 j,n n /2 j,n nn2n j,n One then has for the first-orer wave functions Refs. an n n k /2 nnn2 n 9n /2 n 9n /2 n nn2n n. 2 Using the perturbe wave function, we fin the average value of position to first-orer in is given by n n n   n. 2 We once again make use of raising an lowering operators to fin to O 9 m. J. Phys., Vol. 65, No., March 997 R. W. Robinett 9
3 Fig.. Quantum mechanical soli an classical ot ash probability istributions, P QM () () n () 2 an P CL () vs. The state corresponing to n2 is shown an the vertical ashe lines inicate the locations of the classical turning points at given by E n (n/2)k 2 /2. n  n n  n 6 k /2 n 2 22 so that the average value of position can be written as m 2 n 2 k E n, 2 since km 2. If we then associate the unperturbe energy eigenvalues, E n (), via E n 2k 2, we o inee reprouce the result of Eq. 8, namely, To better visualie how the probability is shifte aroun in response to the anharmonic perturbation to yiel a nonvanishing value of, we first plot in Fig. a stanar image of the quantum mechanical probability ensity P QM () n () 2 soli curve for the case n2. For later comparison, the corresponing classical epression, P CL (), ot ashe curve an the classical turning points vertical ashe lines are also shown. The ifference in probability ensities corresponing to the unperturbe oscillator an the result incluing the effect of the cubic perturbation to first orer given by Eq. 2 is given by P QM n 2 n 2 n nnn2 n 9n /2 n 9n /2 n nn2n n. 26 This is seen to be an o function of for all values of n an is plotte in Fig. 2 for the same n2 case. The classical connection is better realie by suitably averaging locally over the obvious rapi oscillations in the quantum case, an in Fig. we show the same probability ensity, but now Fig. 2. The ifference in probability ensities efine in Eq. 26, P QM () vs between the first-orer wave function () n () of Eq. 2 an the eroth orer result () n (). We use the same value of n2 as in Fig.. binne into finite intervals. The shift in probability ensity from left to right near both the classical turning points is quite similar to what woul be epecte from classical arguments. One final way to erive these results is to use semiclassical probability istributions. These are often iscusse in elementary tets on quantum mechanics see, e.g., Ref. 4, where the unperturbe harmonic oscillator is the stanar eample as in Fig., but which can be applie to many other systems as well. 5 We recall that the average value of a function of the position coorinate can obtaine by generaliing Eq. 7 to yiel f f tt 2 f f PCL, v 27 where we have use the connection v()/t an efine the classical probability istribution via Fig.. The ifference in probability ensities, P QM (), between the eroth- an first-orer results vs integrate over small bins in the position coorinate; we use the same value of n2 as in Figs. an 2. Note that the probability has been shifte slightly from left to right near both classical turning points as epecte from a classical picture. 92 m. J. Phys., Vol. 65, No., March 997 R. W. Robinett 92
4 P CL 2 v. 28 The integral is now taken between the classical turning points, i.e., in the range (, ), which is over only half the perio, /2. The local spee is relate to the potential energy functions via E v 2 U. For the harmonic oscillator, the stanar result is P CL 2 2, which satisfies PCL 29 as it shoul; this is the classical result reprouce in Fig.. For the perturbe oscillator, the turning points are etermine by the roots of the cubic equation 2k 2 E 2 k 2 k or The two roots which represent the actual turning points of the classical motion in the slightly moifie well are given by an, where an has the series epansion The thir etraneous root, which is ue to the eventual turnover of the potential cause by the cubic piece, is given by /22 2 O ( 2 ) an properly is pushe to infinity when. With these changes in turning points an potential, we can write the epectation value of the position coorinate in the form v v I R, 5 where we have efine R()I()/ in the epectation that the integral giving information on the shift in the average value of position will have an epansion in starting with. Our first inclination might be to attempt to epan both require integrals using the stanar calculus result for the ifferentiation of efinite integrals, namely, b b f, f b, f a, a a b a f,. 6 The singular nature of the integrans near the classical turning points, however, makes such an epansion far from obviously convergent. To eplicitly check the behavior of I for, we perform the necessary integrals numerically in this case, using MTHEMTIC for several values of an the results are shown in Table I. For smaller values of, the numerical integration results are not reliable an even for. we are somewhat in the noise. Over this range, however, it certainly seems that the integral is well escribe by the epansion I/2 an this implies, of course, that the epectation value of position satisfies 2 2 2, 7 as epecte from all our earlier results. Using values of in the range.., it also seems likely that the net term in the epansion is 9 2 /4. We can attempt to etract these results in a more analytic fashion using some rather formal properties of the integrals involve. For eample, we can argue that the new classical probability istribution can actually be written in the form P CL with the same normaliation to orer or as before. This can be justifie by calculating? PCL y y 2 2y FF2FO 2 O 2, 9 where an y an we have assume a series epansion for F; one can perform this normaliation integral 9 m. J. Phys., Vol. 65, No., March 997 R. W. Robinett 9
5 Table I. Results of the numerical integration of Eq. 5 for various values of. I//2 R/2/ numerically as mentione above for several small values of an fit to a lower-power polynomial in an one inee fins the result 2 O 4. This implies that the normaliation of P CL () is unchange in this case, at least by a first-orer change in the potential; this is similar to the quantum result where first-orer changes in () only change the normaliation of () 2 to secon orer. Using this notation, we fin that the average value of position in the perturbe oscillator is given by PCL GG 2 G, y y 2 2y y 4 where we also assume that G has a series epansion. We can etract information on G by writing G H H The first integral can be performe yieling H , 42 since efines the classical turning point where the integral vanishes. The secon integral can be evaluate with since this term is alreay eplicitly of orer ; this results in H This implies that G()H 2 ()/4 an from Eq. 4 we fin , 44 as epecte. Once again, one can perform this integral numerically an fit to powers of an this result is confirme. The assumption that we can systematically epan all of these semiclassical quantities in powers of is, perhaps, only justifie to leaing orer in as that correspons, in some sense, to the use of first-orer perturbation theory. Higherorer terms in, which woul correspon to secon-orer an higher perturbation theory woul arise from truly quantum mechanical effects coupling to ifferent states, weighte by energy ifferences an cannot necessarily be epecte to be reprouce by purely classical arguments. CKNOWLEDGMENTS We are grateful for conversations with G. Fleming, an to P. Gol for proviing us with a copy of Ref... B. Pippar, The Physics of Vibration Cambrige U.P., Cambrige, 978, Vol. ; Cambrige U.P., Cambrige, 98, Vol D. Halliay, R. Resnick, an J. Walker, Funamentals of Physics Wiley, New York, 99, 4th e., p. 54. H. D. Young, University Physics ison-wesley, Reaing, M, 992, 8th e., p R.. Serway, Physics for Scientists an Engineers Sauners, Philaelphia, 992, r e., upate version, p C. Kittel, W. D. Knight, an M. Ruerman, Berkeley Physics Course, Mechanics Eucational Services Inc., McGraw-Hill, New York, 965, pp C. Kittel, Introuction to Soli State Physics Wiley, New York, 97, 4th e., p D. S. Saon, Elementary Quantum Mechanics Holen-Day, San Francisco, 968, pp R. L. Liboff, Introuctory Quantum Mechanics ison-wesley, Reaing, 99, 2n e., pp T.-Y. Wu, Quantum Mechanics Worl Scientific, Singapore, 986, pp C. Cohen-Tannouji, B. Diu, an F. Laloë, Quantum Mechanics Wiley, New York, 977, pp L. Pauling an E. B. Wilson, Introuction to Quantum Mechanics McGraw Hill, New York, 95, pp E. E. nerson, Moern Physics an Quantum Mechanics Sauners, Philaelphia, 97, pp C. S. Johnson, Jr. an L. G. Peerson, Problems an Solutions in Quantum Chemistry an Physics Dover, New York, 986, pp See Ref. 8, pp R. W. Robinett, Quantum an classical probability istributions for position an momentum, m. J. Phys. 6, m. J. Phys., Vol. 65, No., March 997 R. W. Robinett 94
II. First variation of functionals
II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent
More informationd 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 informationCalculus in the AP Physics C Course The Derivative
Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.
More informationCalculus 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 informationensembles 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 informationTable 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 informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationLagrangian 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 informationOptimization of Geometries by Energy Minimization
Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.
More information1 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 informationSummary: Differentiation
Techniques of Differentiation. Inverse Trigonometric functions The basic formulas (available in MF5 are: Summary: Differentiation ( sin ( cos The basic formula can be generalize as follows: Note: ( sin
More informationGeneralization 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 informationIntroduction 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 informationIntegration Review. May 11, 2013
Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In
More informationSeparation of Variables
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
More information3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes
Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More information05 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 informationThe 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 informationA. Incorrect! The letter t does not appear in the expression of the given integral
AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationLeChatelier Dynamics
LeChatelier Dynamics Robert Gilmore Physics Department, Drexel University, Philaelphia, Pennsylvania 1914, USA (Date: June 12, 28, Levine Birthay Party: To be submitte.) Dynamics of the relaxation of a
More informationEuler 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 informationUnit #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 informationfv = 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 informationSemiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom
PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,
More informationThe 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 informationExtinction, σ/area. Energy (ev) D = 20 nm. t = 1.5 t 0. t = t 0
Extinction, σ/area 1.5 1.0 t = t 0 t = 0.7 t 0 t = t 0 t = 1.3 t 0 t = 1.5 t 0 0.7 0.9 1.1 Energy (ev) = 20 nm t 1.3 Supplementary Figure 1: Plasmon epenence on isk thickness. We show classical calculations
More information4. Important theorems in quantum mechanics
TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this
More informationWJEC Core 2 Integration. Section 1: Introduction to integration
WJEC Core Integration Section : Introuction to integration Notes an Eamples These notes contain subsections on: Reversing ifferentiation The rule for integrating n Fining the arbitrary constant Reversing
More informationand 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 informationApplications of First Order Equations
Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more
More informationThe Sokhotski-Plemelj Formula
hysics 25 Winter 208 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),
More informationMATH2231-Differentiation (2)
-Differentiation () The Beginnings of Calculus The prime occasion from which arose my iscovery of the metho of the Characteristic Triangle, an other things of the same sort, happene at a time when I ha
More informationLecture 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 informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More information12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes
Differentiation of vectors 12.5 Introuction The area known as vector calculus is use to moel mathematically a vast range of engineering phenomena incluing electrostatics, electromagnetic fiels, air flow
More informationCalculus 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 informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationSchrö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 informationQuantum 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 informationLecture 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 informationMath 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 informationThe 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 informationPHYS 414 Problem Set 2: Turtles all the way down
PHYS 414 Problem Set 2: Turtles all the way own This problem set explores the common structure of ynamical theories in statistical physics as you pass from one length an time scale to another. Brownian
More informationAnalytic Scaling Formulas for Crossed Laser Acceleration in Vacuum
October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945
More informationVectors 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 informationImplicit 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 informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationIntroduction 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 informationMath 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 informationLinear and quadratic approximation
Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function
More informationQuantum 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 informationHomework 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 informationPlanar sheath and presheath
5/11/1 Flui-Poisson System Planar sheath an presheath 1 Planar sheath an presheath A plasma between plane parallel walls evelops a positive potential which equalizes the rate of loss of electrons an ions.
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationLinear 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 informationAPPROXIMATE 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 informationAssignment 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 informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More informationAN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A
AN INTRODUCTION TO AIRCRAFT WIN FLUTTER Revision A By Tom Irvine Email: tomirvine@aol.com January 8, 000 Introuction Certain aircraft wings have experience violent oscillations uring high spee flight.
More informationConservation 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 informationLagrangian 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 informationChapter 9 Method of Weighted Residuals
Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe
More informationA Model of Electron-Positron Pair Formation
Volume PROGRESS IN PHYSICS January, 8 A Moel of Electron-Positron Pair Formation Bo Lehnert Alfvén Laboratory, Royal Institute of Technology, S-44 Stockholm, Sween E-mail: Bo.Lehnert@ee.kth.se The elementary
More informationarxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003
Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep
More informationTHE 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 informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationSensors & 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 informationThermal 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 informationarxiv:nlin/ v1 [nlin.cd] 21 Mar 2002
Entropy prouction of iffusion in spatially perioic eterministic systems arxiv:nlin/0203046v [nlin.cd] 2 Mar 2002 J. R. Dorfman, P. Gaspar, 2 an T. Gilbert 3 Department of Physics an Institute for Physical
More informationTrigonometric Functions
72 Chapter 4 Trigonometric Functions 4 Trigonometric Functions To efine the raian measurement system, we consier the unit circle in the y-plane: (cos,) A y (,0) B So far we have use only algebraic functions
More informationThe Sokhotski-Plemelj Formula
hysics 24 Winter 207 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),
More informationOn the uniform tiling with electrical resistors
On the uniform tiling with electrical resistors M. Q. Owaiat Department of Physics, Al-Hussein Bin Talal University, Ma an,7, Joran E-mail: Owaiat@ahu.eu.jo Abstract We calculate the effective resistance
More informationMath 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 information18 EVEN MORE CALCULUS
8 EVEN MORE CALCULUS Chapter 8 Even More Calculus Objectives After stuing this chapter you shoul be able to ifferentiate an integrate basic trigonometric functions; unerstan how to calculate rates of change;
More informationQubit channels that achieve capacity with two states
Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March
More informationJUST 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 informationTrack Initialization from Incomplete Measurements
Track Initialiation from Incomplete Measurements Christian R. Berger, Martina Daun an Wolfgang Koch Department of Electrical an Computer Engineering, University of Connecticut, Storrs, Connecticut 6269,
More informationChapter 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 informationTMA 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 informationExamining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing
Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June
More informationKinetic Energy Is Important in the Nanoscale World
Kinetic Energy Is Important in the Nanoscale Worl Frank Riou Department of Chemistry College of St. Beneict & St. John's University St. Joseph, MN 56374 Most eplanations of atomic an molecular phenomena
More informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More information12.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 information3.2 Differentiability
Section 3 Differentiability 09 3 Differentiability What you will learn about How f (a) Might Fail to Eist Differentiability Implies Local Linearity Numerical Derivatives on a Calculator Differentiability
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
More informationmodel considered before, but the prey obey logistic growth in the absence of predators. In
5.2. First Orer Systems of Differential Equations. Phase Portraits an Linearity. Section Objective(s): Moifie Preator-Prey Moel. Graphical Representations of Solutions. Phase Portraits. Vector Fiels an
More informationAPPPHYS 217 Thursday 8 April 2010
APPPHYS 7 Thursay 8 April A&M example 6: The ouble integrator Consier the motion of a point particle in D with the applie force as a control input This is simply Newton s equation F ma with F u : t q q
More informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationWitten s Proof of Morse Inequalities
Witten s Proof of Morse Inequalities by Igor Prokhorenkov Let M be a smooth, compact, oriente manifol with imension n. A Morse function is a smooth function f : M R such that all of its critical points
More informationN may be the number of photon interactions in a given volume or the number of radioactive disintegrations in a given time. Expectation value:
DESCRIPTION OF IONIZING RADIATION FIELDS (Chapter 1 p5-19) Stochastic Vs Non-Stochastic Description Raiation interaction is always stochastic to some egree. o entities whether photons or charge particles
More informationON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract
More informationFree 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 informationStrength Analysis of CFRP Composite Material Considering Multiple Fracture Modes
5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation
More informationinvolve: 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 informationConstruction 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 informationState observers and recursive filters in classical feedback control theory
State observers an recursive filters in classical feeback control theory State-feeback control example: secon-orer system Consier the riven secon-orer system q q q u x q x q x x x x Here u coul represent
More informationTopic 2.3: The Geometry of Derivatives of Vector Functions
BSU Math 275 Notes Topic 2.3: The Geometry of Derivatives of Vector Functions Textbook Sections: 13.2 From the Toolbox (what you nee from previous classes): Be able to compute erivatives scalar-value functions
More information