A new class of sum rules for products of Bessel functions. G. Bevilacqua, 1, a) V. Biancalana, 1 Y. Dancheva, 1 L. Moi, 1 2, b) and T.
|
|
- Lorin Lucas
- 6 years ago
- Views:
Transcription
1 A new class of sum rules for products of Bessel functions G. Bevilacqua,, a V. Biancalana, Y. Dancheva, L. Moi,, b and T. Mansour CNISM and Dipartimento di Fisica, Università di Siena, Via Roma 56, 5300 Siena, Italy Department of Mathematics, University of Haifa, 3905 Haifa, Israel (Dated: February 3, 0 We derive a new class of sum rules for products of the Bessel functions of first kind. Using standard algebraic manipulations we extend some of the well known properties of J n. Some physical applications of the results are also discussed. A comparison with the Newberger[J. Math. Phys. 3 (98 78] sum rules is performed on a typical example. PACS numbers: 0.30.Lt, 78.0.Bh a Electronic mail: bevilacqua@unisi.it b Electronic mail: toufik@math.haifa.ac.il
2 I. INTRODUCTION Bessel functions of first kind (J n are among the most ubiquitous special functions of the mathematical physics. Their properties are described in several monographs 4, in the encyclopedic Abramowitz & Stegun 5 as well as in every advanced mathematics textbook. The generating function identity e z (w w J n (z w n ( is the key formula used in discussing the frequency modulation (FM of the electric field of a laser. In fact a pure sinusoidal FM reads as E(t = E 0 ( e i(ω L t+m sin(ω t + c.c., ( where ω L is the laser carrier frequency, Ω is the modulation frequency, M the modulation index and MΩ the modulation depth. Using ( with w = e i Ω t one finds the well-known sidebands structure E(t = E 0 J n (M e i(ω L+nΩt + c.c. (3 In some applications, for instance studying the response of a two-level system to a FM field, one has to deal with expressions like A s = J n (MJ n s (M, (4 γ + i n Ω where γ represents a resonance linewidth. Even if the sum could be done exactly using the Newberger sum rule 6, the result is quite involved (the details are given in Section III A. Alternatively, in the physical interesting limit ɛ = MΩ/γ < (see below for more details, a geometric expansion gives A s = γ J n (MJ n s (M reducing the problem to the evaluation of sums of the form B k,s + [ iω γ n ( Ω γ n + O(ɛ 3 ] (5 n k J n (MJ n s (M (6
3 for integer values of k 0 and s. To the authors best knowledge these sums are not addressed in the literature in the general case although special values are known, for instance B 0,0 = from of [5]. Finding a closed expression for B k,s is the main goal of the present work. The paper is organized as follows. In Section II the main results are derived and discussed, while in Section III some applications are pointed out. Finally Section IV contains the conclusions. II. RESULTS We start the derivation of the main result by rewriting ( in a slight different form f(y, θ e i y sin θ = J n (y e i n θ. (7 Deriving the left-hand side k times one gets f (k (y, θ k f(y, θ θ k = p k (y, θf(y, θ (8 where p k (y, θ is a polynomial in cos θ and sin θ and can be easily recast in the form p k (y, θ = C k,n (y e i n θ, (9 with C k,n (y polynomials in y. An explicit form for these coefficients is worked out in Appendix A. Putting all together one finds ( f (k (y, θ = C k,q (y J n q (y e i n θ (0 q= On the other hand, deriving the right-hand side of (7 one obtains f (k (y, θ = (i n k J n (y e i n θ. ( An explicit comparison of ( with (0 is performed in Appendix B. Now consider the expression F (y, y, θ f (k (y, θf(y, θ. A first form can be obtained using (9 and (7 as follows F (y, y, θ = p k (y, θ e i (y +y sin θ ( ( = C k,n (y e i n θ = q= ( m= m= C k,q m (y J m (y + y 3 J m (y + y e i m θ e i q θ. (
4 In the meantime an equivalent formula can be written using ( ( ( F (y, y, θ = (i n k J n (y e i n θ J m (y e i m θ = ( q= m= (i n k J n (y J q n (y e i q θ, (3 and, given the completeness of the set {e i q θ, q = 0, ±, ±,...}, the coefficients in both expansions must be equal C k,q m (y J m (y + y = m= (i n k J n (y J q n (y. (4 This is the central result of the paper and can be considered as a generalization of the J n addition formula [5, formula 9..75], which, in fact, is obtained for k = 0 noticing that C 0,n = δ n,0 (see Appendix A. Besides, with a proper choice of y and y different sum rules can be obtained. In fact, substituting y = y = y in (4 and remembering that J n (0 = δ n,0 one finds C k,q (y = = = (i n k J n (y J q n ( y (i n k J n (y ( q n J q n (y (i n k J n (y J n q (y. This gives an answer for the initial problem, namely and the B k,s for the first k values can be read from Appendix A. (5 B k,s = i k C k,s (6 The companion case y = y = y is worked out with similar algebraic elaborations leading to the formula ( n n k J n (y J n q (y = ( q i k m= C k,q m (y J m (y. (7 A. Generalizations Many generalizations are possible once one dives in the world of generalized Bessel functions 4. Let us discuss some cases. 4
5 Consider the function e i (x cos θ+y sin θ J n (x, y e i n θ, (8 which arises in the context of general sinusoidal modulation. The J n (x, y are a kind of generalized Bessel functions, which can be given in terms of usual Bessel functions by means of simple algebraic manipulations (see also equation (.0 of Dattoli 4 J n (x, y i q J q (x J n q (y. (9 q= The same steps as before lead to the sum rule n J [ n (x, y Jn q (x, y] (y + i x δq, + (y i x δ q,. (0 Furthermore k-order derivatives lead to sum rules with n k which the interested reader can easily work out. The second generalization is achieved considering e i (y sin θ+y sin θ Ĵ n (y, y e i n θ. ( These Ĵn(y, y generalized Bessel functions can be expressed as Ĵ n (y, y J q (y J n q (y, ( q= and are real function for real values of y i (see chapter of Dattoli 4 for further details. By means of the same device one gets n Ĵn(y, y Ĵn s(y, y y (δ s, + δ s, + y (δ s, + δ s,. (3 Also in this case higher derivatives lead to other sum rules. In a complete general setting, one could consider the problem of a FM signal with no special properties, i.e. e i φ(t = G n e i n Ω t (4 where the only requirement is that φ(t is a real and periodic well-behaved function (a more precise statement is given below φ(t = 5 φ n e i n Ω t. (5
6 The G n coefficients are related in a very complicated way to multiple sums of products of ordinary Bessel functions, but the explicit form is not important here. They satisfy the conservation of energy sum rule G n G n s = δ s,0. (6 Taking the time derivative of (4 and doing the same steps as above we get n G n G n s = i s φ s. (7 Equation (7 holds provided that φ(t can be expanded as a meaningful Fourier series, i.e. φ(t = i n Ω φ n e i n Ω t (8 must be convergent. As before, higher order derivatives will give higher order sum rules. III. APPLICATIONS An example is considered where formula (6 comes handy. Let us model an optical resonance as a damped harmonic oscillator of unit mass, forced by the electric field of the light beam z + γż + ω 0z = fe +iωt, (9 where, as usual, the physical oscillator displacement is x(t = R(z(t and all other symbols have a clear meaning. Such a simple and analytical model avoids all the complications tied with the solution and let us concentrate on the sum rule. In fact, after the transient, the solution settles into the well-known z(t = and the work done by the light field per time unit is dw = f ( dt 4 (eiωt + e iωt f e +iωt ω 0 ω + iγω, (30 iω e iωt ω 0 ω + iγω + iω e iωt. (3 ω0 ω iγω The experimental apparatus usually averages out high frequencies and thus the measured signal is dw ( dt = f I ω ω0 ω + iγω 6 = f ω γ (ω ω 0 + ω γ, (3
7 which, for nearly resonant light ω = ω 0 + δ, δ ω 0, becomes dw dt = f [ γ + + ( 3 ( + + O ω 0 ω 0 ], (33 where = δ/γ. It is easily seen that the leading term shows the well-known Lorentzian form. Replacing the forcing term with a frequency modulated one e i ω t e i (ω t+m sin Ωt = + J n (M e i ωn t (34 in the right-hand side of (9 we get the exact result [ dw ( dt = f I ] ω n J n (M J n s (M e i s Ω t ω0 ωn + iγω n s= (35 where ω n = ω 0 + δ + n Ω. To proceed any further it is necessary to simplify the expression F = ω n /(ω 0 ω n + iγω n. In the optical range it is not restrictive to assume ω 0 γ, δ getting F = [ + δ + nω i γ (δ + nω δ + nω i γ Considering only the leading term in F we are faced in (35 with the sums ( ] + O. (36 ω 0 ω0 iγ J n (MJ n s (M + i + n i Ω/γ (37 which are very similar to (4. To proceed any further with the geometric expansion of the denominator, we notice that the product J n (M J n s (M shows the well-known bridge structure and it is exponentially small when n > N MAX M. So if N MAX Ω/γ < N MAX Ω/γ + i < (38 and the sum can be elaborated as follows iγ J n (MJ n s (M + i + n i Ω/γ iγ = iγ iγ N MAX n= N MAX N MAX J n (MJ n s(m + i + n i Ω/γ J n (MJ n s(m + i n= N MAX J n (MJ n s (M + i [ + nɛ + n ɛ + O(ɛ 3 ] [ + nɛ + n ɛ + O(ɛ 3 ], (39 7
8 where ɛ = i (Ω/γ/(+i and the last row is justified because of the exponentially small nature of terms included. The inner sum in (35 can now worked out by means of the sum rules developed above, and, after some straightforward algebra, one finds dw dt = f { γ + + MΩ cos Ωt γ ( + + ( MΩ ( 3 sin Ωt M γ ( ( } MΩ 3 ( + cos Ωt +..., γ ( + 3 (40 where the terms in parentheses, besides the d.c. component, represent the higher harmonic modulated absorption. Notice how the dependence of the (MΩ/γ n term is related to the real and imaginary part of n-th derivative of /( + i. It seems worth to stress that the validity of (4 is limited to modulation depths smaller than the linewidth as pointed out in (38. Usually in typical experimental conditions one is interested in the first harmonic signal, catching it, for instance, by means of a lock-in amplifier. The above formula shows that it is more convenient to use the biggest value of the modulation depth M Ω compatible with (38. More generally, using an arbitrary but limited modulation signal as in (5 one finds dw dt = f = f { γ + + γ { γ + + σ γ } ( + φ(t +... ω(t ( + σ +... where in the last row the instantaneous frequency ω(t = φ(t is introduced and the modulation depth σ = ω MAX ω min is put in evidence to compare with previous formulae. }. (4 A. Comparison with Newberger sum rule As stated in the Introduction, the quantity A s can be summed by mean of Newberger 6 sum rule. The result is A s = ( s γ π γ/ω sinh(π γ/ω J s i γ/ω(m J i γ/ω (M, s 0 (4 8
9 where one has to deal with complex-order Bessel functions. Even if this is an exact result, the physics is obscured by the formula complexity. For instance it is not easy to extract the behaviour when η = Ω/γ is small. Instead using (5 with s = we quickly found A = M γ i η M γ η + M γ ( + 3 M /4 i η 3 + O(η 4 (43 To proceed further using (4 the product of the Bessel functions must be developed 5 (formula 9..4, and, making use of the properties of Euler Gamma function 5 (formula 6..3, one gets A s = ( s γ ( s ( M ( M s 4 k (s + k! (s + k! k! k + p i γ/ω k=0 p= which, to the authors best knowledge, is an original elaboration. k p= p + (γ/ω (44 This formula is still complicated, but better suited to extract asymptotic behaviour. In fact for s = we have A = M ( γ iγ/ω 3M 4 iγ/ω + (γ/ω +... = M γ i η M γ η + M (45 γ ( + 3 M /4 i η 3 + O(η 4. IV. CONCLUSIONS We have derived some new classes of sum rules and recursion relations obeyed by the Bessel functions of first kind. The results have broad application in the physics of modulated excitation for instance in the case of light-matter interaction. A comparison of our results with an elaborated form of the Newberger sum rule is also given. Appendix A: Determination of C k,n The starting point is the recursion obeyed by the p k polynomials. In fact, it is easily seen that p 0, and assuming (8 then follows ( p k+ (y, θ = θ + i y cos θ p k (y, θ. Projecting this equation as shown in (9 one obtains (A C k+,n = i n C k,n + i y (C k,n+ + C k,n with C 0,n = δ n,0. (A 9
10 To the authors best knowledge this recursion cannot be solved in term of known special functions, even if the homogeneous right-hand side, namely i n C k,n + i y (C k,n+ + C k,n = 0, is satisfied by the J n (y itself. The first few loops give i C,n = y (δ n,+ + δ n, i C,n = y 4 (δ n,0 + δ n,+ + δ n, + y (δ n,+ δ n, i C 3 3,n = y3 8 (δ n,+3 + δ n, y 4 (δ n,+ δ n, + ( 3 y y (δ n,+ + δ n, i C 4 4,n = y4 6 (δ n,+4 + δ n, 4 + 3y3 4 (δ n,+3 δ n, 3 + ( y y (δ n,+ + δ n, + 4 ( 3y 3 (δ n,+ δ n, y ( 3y y 4 δ n,0 The recursion is straightforward, but tedious. A simplification is achieved thinking (A as a matrix vector product and implementing it in any computer algebra system. Alternatively an approach based on the Faà di Bruno formula can be developed. The Faà di Bruno formula generalizes the chain rule to higher derivative and can be stated as follows d k d x k f(g(x = {m j } k! m! m! m k! f (m + +m k (g(x k ( g (j mj (x. j! j= (A3 where the sum is over all k-tuples of non-negative integers m,..., m k satisfying the constraint m + m + 3 m k m k = k. Applying that device to our problem we find i y sin θ p k = e dk y sin θ ei d θk = {m j } k! j m j! j! m j (i ym k (sin (j θ m j, j= (A4 where m = j m j. Next using sin (j θ = ( (j cos θ and sin (j θ = ( (j sin θ we 0
11 get p k = {m j } k! j m j! j! m j (i ym ( φ sin a θ cos b θ (A5 where φ = j 0 (m +4 j + m +4 j+, a = j m j and b = j m j = m a. Finally expressing the trigonometric functions as sin a θ cos b θ = m i a a r=0 b ( ( a b ( r r r r =0 e i (m r r θ (A6 and substituting in (A5 we find C k,n = {m j } k!i b ( φ (y/ m j m j! j! m j a ( ( a ( r b m n r r r=0 (A7 with m n m, where we follow the convention that the binomial coefficients are zero if the lower index is not integer or negative or larger than the upper one. Appendix B: Some new recursion relations By comparing the two different forms (0 and ( of the same quantity, a relation between the coefficients results q k J q (y i k C k,n(y J q n (y. The k = 0 case is the identity J q (y = J q (y. For k = we obtain the well-known 5 (formula 9..7 first row three term recursion relation (B q J q (y = y [ J q+ (y + J q (y ]. (B The interesting and, to the authors best knowledge, new relations are found for k. In fact for k = one gets ( y q J q (y = y [ Jq+ (y J q (y ] y [ Jq+ (y + J q (y ]. 4 For higher values of k other relations are easily obtained using the coefficients of Appendix A. (B3
12 REFERENCES A. Gray and G. B. Mathews, A treatise on Bessel functions and their applications to physics (Macmillan and Co., 895. G. N. Watson, A treatise on the theory of Bessel functions (Cambridge University Press, 9. 3 B. G. Korenev, Bessel functions and their applications (Taylor and Francis, G. Dattoli and A. Torre, Theory and Applications of Generalized Bessel Functions (Aracne Editrice, Rome, Milton Abramowitz and Irene A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, ninth dover printing, tenth gpo printing ed. (Dover, New York, 964 ISBN Barry S. Newberger, New sum rule for products of bessel functions with application to plasma physics, J. Math. Phys. 3, 78 (98.
Notes on the Periodically Forced Harmonic Oscillator
Notes on the Periodically orced Harmonic Oscillator Warren Weckesser Math 38 - Differential Equations 1 The Periodically orced Harmonic Oscillator. By periodically forced harmonic oscillator, we mean the
More informationBESSEL FUNCTIONS APPENDIX D
APPENDIX D BESSEL FUNCTIONS D.1 INTRODUCTION Bessel functions are not classified as one of the elementary functions in mathematics; however, Bessel functions appear in the solution of many physical problems
More informationLab10: FM Spectra and VCO
Lab10: FM Spectra and VCO Prepared by: Keyur Desai Dept. of Electrical Engineering Michigan State University ECE458 Lab 10 What is FM? A type of analog modulation Remember a common strategy in analog modulation?
More informationBessel function - Wikipedia, the free encyclopedia
Bessel function - Wikipedia, the free encyclopedia Bessel function Page 1 of 9 From Wikipedia, the free encyclopedia In mathematics, Bessel functions, first defined by the mathematician Daniel Bernoulli
More informationChapter 2: Complex numbers
Chapter 2: Complex numbers Complex numbers are commonplace in physics and engineering. In particular, complex numbers enable us to simplify equations and/or more easily find solutions to equations. We
More informationAnswers to Problem Set Number MIT (Fall 2005).
Answers to Problem Set Number 6. 18.305 MIT (Fall 2005). D. Margetis and R. Rosales (MIT, Math. Dept., Cambridge, MA 02139). December 12, 2005. Course TA: Nikos Savva, MIT, Dept. of Mathematics, Cambridge,
More informationAtomic cross sections
Chapter 12 Atomic cross sections The probability that an absorber (atom of a given species in a given excitation state and ionziation level) will interact with an incident photon of wavelength λ is quantified
More information221B Lecture Notes on Resonances in Classical Mechanics
1B Lecture Notes on Resonances in Classical Mechanics 1 Harmonic Oscillators Harmonic oscillators appear in many different contexts in classical mechanics. Examples include: spring, pendulum (with a small
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More informationComplex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,
Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )
More information1 Pushing your Friend on a Swing
Massachusetts Institute of Technology MITES 017 Physics III Lecture 05: Driven Oscillations In these notes, we derive the properties of both an undamped and damped harmonic oscillator under the influence
More informationMultiple scale methods
Multiple scale methods G. Pedersen MEK 3100/4100, Spring 2006 March 13, 2006 1 Background Many physical problems involve more than one temporal or spatial scale. One important example is the boundary layer
More informationTime Response of Systems
Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =
More informationSolution Set 2 Phys 4510 Optics Fall 2014
Solution Set Phys 4510 Optics Fall 014 Due date: Tu, September 16, in class Scoring rubric 4 points/sub-problem, total: 40 points 3: Small mistake in calculation or formula : Correct formula but calculation
More information1 (2n)! (-1)n (θ) 2n
Complex Numbers and Algebra The real numbers are complete for the operations addition, subtraction, multiplication, and division, or more suggestively, for the operations of addition and multiplication
More informationDamped Harmonic Oscillator
Damped Harmonic Oscillator Wednesday, 23 October 213 A simple harmonic oscillator subject to linear damping may oscillate with exponential decay, or it may decay biexponentially without oscillating, or
More informationMathematical Physics
Mathematical Physics MP205 Vibrations and Waves Lecturer: Office: Lecture 9-10 Dr. Jiří Vala Room 1.9, Mathema
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,
More informationNew asymptotic expansion for the Γ (z) function.
New asymptotic expansion for the Γ z function. Gergő Nemes Institute of Mathematics, Eötvös Loránd University 7 Budapest, Hungary September 4, 007 Published in Stan s Library, Volume II, 3 Dec 007. Link:
More informationFourier Series and Fourier Transforms
Fourier Series and Fourier Transforms EECS2 (6.082), MIT Fall 2006 Lectures 2 and 3 Fourier Series From your differential equations course, 18.03, you know Fourier s expression representing a T -periodic
More informationTheory of bifurcation amplifiers utilizing the nonlinear dynamical response of an optically damped mechanical oscillator
Theory of bifurcation amplifiers utilizing the nonlinear dynamical response of an optically damped mechanical oscillator Research on optomechanical systems is of relevance to gravitational wave detection
More informationVibrations and waves: revision. Martin Dove Queen Mary University of London
Vibrations and waves: revision Martin Dove Queen Mary University of London Form of the examination Part A = 50%, 10 short questions, no options Part B = 50%, Answer questions from a choice of 4 Total exam
More information4. The interaction of light with matter
4. The interaction of light with matter The propagation of light through chemical materials is described by a wave equation similar to the one that describes light travel in a vacuum (free space). Again,
More informationMacroscopic dielectric theory
Macroscopic dielectric theory Maxwellʼs equations E = 1 c E =4πρ B t B = 4π c J + 1 c B = E t In a medium it is convenient to explicitly introduce induced charges and currents E = 1 B c t D =4πρ H = 4π
More informationProblem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension
105A Practice Final Solutions March 13, 01 William Kelly Problem 1: Lagrangians and Conserved Quantities Consider the following action for a particle of mass m moving in one dimension S = dtl = mc dt 1
More informationMATH 251 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam
MATH 51 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam A collection of previous exams could be found at the coordinator s web: http://www.math.psu.edu/tseng/class/m51samples.html
More informationAnswers to Problem Set Number MIT (Fall 2005).
Answers to Problem Set Number 5. 18.305 MIT Fall 2005). D. Margetis and R. Rosales MIT, Math. Dept., Cambridge, MA 02139). November 23, 2005. Course TA: Nikos Savva, MIT, Dept. of Mathematics, Cambridge,
More informationNORMAL MODES, WAVE MOTION AND THE WAVE EQUATION. Professor G.G.Ross. Oxford University Hilary Term 2009
NORMAL MODES, WAVE MOTION AND THE WAVE EQUATION Professor G.G.Ross Oxford University Hilary Term 009 This course of twelve lectures covers material for the paper CP4: Differential Equations, Waves and
More information01 Harmonic Oscillations
Utah State University DigitalCommons@USU Foundations of Wave Phenomena Library Digital Monographs 8-2014 01 Harmonic Oscillations Charles G. Torre Department of Physics, Utah State University, Charles.Torre@usu.edu
More information10. Operators and the Exponential Response Formula
52 10. Operators and the Exponential Response Formula 10.1. Operators. Operators are to functions as functions are to numbers. An operator takes a function, does something to it, and returns this modified
More informationOutline of parts 1 and 2
to Harmonic Loading http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March, 6 Outline of parts and of an Oscillator
More informationA.1 Introduction & Rant
Scott Hughes 1 March 004 A.1 Introduction & Rant Massachusetts Institute of Technology Department of Physics 8.0 Spring 004 Supplementary notes: Complex numbers Complex numbers are numbers that consist
More informationReview of Linear Time-Invariant Network Analysis
D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x
More informationMulti Degrees of Freedom Systems
Multi Degrees of Freedom Systems MDOF s http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March 9, 07 Outline, a System
More informationWave Phenomena Physics 15c. Lecture 2 Damped Oscillators Driven Oscillators
Wave Phenomena Physics 15c Lecture Damped Oscillators Driven Oscillators What We Did Last Time Analyzed a simple harmonic oscillator The equation of motion: The general solution: Studied the solution m
More informationDiscrete Simulation of Power Law Noise
Discrete Simulation of Power Law Noise Neil Ashby 1,2 1 University of Colorado, Boulder, CO 80309-0390 USA 2 National Institute of Standards and Technology, Boulder, CO 80305 USA ashby@boulder.nist.gov
More informationThe formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d
Solving ODEs using Fourier Transforms The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d 2 dx + q f (x) R(x)
More informationDerivation of the General Propagation Equation
Derivation of the General Propagation Equation Phys 477/577: Ultrafast and Nonlinear Optics, F. Ö. Ilday, Bilkent University February 25, 26 1 1 Derivation of the Wave Equation from Maxwell s Equations
More informationA plane autonomous system is a pair of simultaneous first-order differential equations,
Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium
More informationLab 1: Damped, Driven Harmonic Oscillator
1 Introduction Lab 1: Damped, Driven Harmonic Oscillator The purpose of this experiment is to study the resonant properties of a driven, damped harmonic oscillator. This type of motion is characteristic
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 9 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,
More informationLab 1: damped, driven harmonic oscillator
Lab 1: damped, driven harmonic oscillator 1 Introduction The purpose of this experiment is to study the resonant properties of a driven, damped harmonic oscillator. This type of motion is characteristic
More information( ) f (k) = FT (R(x)) = R(k)
Solving ODEs using Fourier Transforms The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d 2 dx + q f (x) = R(x)
More information23.6. The Complex Form. Introduction. Prerequisites. Learning Outcomes
he Complex Form 3.6 Introduction In this Section we show how a Fourier series can be expressed more concisely if we introduce the complex number i where i =. By utilising the Euler relation: e iθ cos θ
More informationMath Assignment 5
Math 2280 - Assignment 5 Dylan Zwick Fall 2013 Section 3.4-1, 5, 18, 21 Section 3.5-1, 11, 23, 28, 35, 47, 56 Section 3.6-1, 2, 9, 17, 24 1 Section 3.4 - Mechanical Vibrations 3.4.1 - Determine the period
More informationBessel Functions Michael Taylor. Lecture Notes for Math 524
Bessel Functions Michael Taylor Lecture Notes for Math 54 Contents 1. Introduction. Conversion to first order systems 3. The Bessel functions J ν 4. The Bessel functions Y ν 5. Relations between J ν and
More informationConvergence of sequences and series
Convergence of sequences and series A sequence f is a map from N the positive integers to a set. We often write the map outputs as f n rather than f(n). Often we just list the outputs in order and leave
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationLecture 3: Central Limit Theorem
Lecture 3: Central Limit Theorem Scribe: Jacy Bird (Division of Engineering and Applied Sciences, Harvard) February 8, 003 The goal of today s lecture is to investigate the asymptotic behavior of P N (
More information1 Simple Harmonic Oscillator
Massachusetts Institute of Technology MITES 2017 Physics III Lectures 02 and 03: Simple Harmonic Oscillator, Classical Pendulum, and General Oscillations In these notes, we introduce simple harmonic oscillator
More informationMechanics IV: Oscillations
Mechanics IV: Oscillations Chapter 4 of Morin covers oscillations, including damped and driven oscillators in detail. Also see chapter 10 of Kleppner and Kolenkow. For more on normal modes, see any book
More informationThe Interaction of Light and Matter: α and n
The Interaction of Light and Matter: α and n The interaction of light and matter is what makes life interesting. Everything we see is the result of this interaction. Why is light absorbed or transmitted
More informationHandout 11: AC circuit. AC generator
Handout : AC circuit AC generator Figure compares the voltage across the directcurrent (DC) generator and that across the alternatingcurrent (AC) generator For DC generator, the voltage is constant For
More informationINFINITE SEQUENCES AND SERIES
11 INFINITE SEQUENCES AND SERIES INFINITE SEQUENCES AND SERIES In section 11.9, we were able to find power series representations for a certain restricted class of functions. INFINITE SEQUENCES AND SERIES
More informationDamped Oscillation Solution
Lecture 19 (Chapter 7): Energy Damping, s 1 OverDamped Oscillation Solution Damped Oscillation Solution The last case has β 2 ω 2 0 > 0. In this case we define another real frequency ω 2 = β 2 ω 2 0. In
More informationCOMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE
BIT 6-85//41-11 $16., Vol. 4, No. 1, pp. 11 118 c Swets & Zeitlinger COMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE WALTER GAUTSCHI Department of Computer Sciences,
More informationLecture 6: Differential Equations Describing Vibrations
Lecture 6: Differential Equations Describing Vibrations In Chapter 3 of the Benson textbook, we will look at how various types of musical instruments produce sound, focusing on issues like how the construction
More informationLaplace Transform of Spherical Bessel Functions
Laplace Transform of Spherical Bessel Functions arxiv:math-ph/01000v 18 Jan 00 A. Ludu Department of Chemistry and Physics, Northwestern State University, Natchitoches, LA 71497 R. F. O Connell Department
More informationElectronic resonances in broadband stimulated Raman spectroscopy: Supplementary Information
Electronic resonances in broadband stimulated Raman spectroscopy: Supplementary Information G. Batignani,2, E. Pontecorvo, G. Giovannetti, C. Ferrante, G. Fumero, T. Scopigno,3 Dipartimento di Fisica,
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationLINEAR RESPONSE THEORY
MIT Department of Chemistry 5.74, Spring 5: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff p. 8 LINEAR RESPONSE THEORY We have statistically described the time-dependent behavior
More informationMATH1190 CALCULUS 1 - NOTES AND AFTERNOTES
MATH90 CALCULUS - NOTES AND AFTERNOTES DR. JOSIP DERADO. Historical background Newton approach - from physics to calculus. Instantaneous velocity. Leibniz approach - from geometry to calculus Calculus
More informationCHEE 319 Tutorial 3 Solutions. 1. Using partial fraction expansions, find the causal function f whose Laplace transform. F (s) F (s) = C 1 s + C 2
CHEE 39 Tutorial 3 Solutions. Using partial fraction expansions, find the causal function f whose Laplace transform is given by: F (s) 0 f(t)e st dt (.) F (s) = s(s+) ; Solution: Note that the polynomial
More informationDIFFRACTION OF PLANE SH WAVES BY A CIRCULAR CAVITY IN QUARTER-INFINITE MEDIUM
11 th International Conference on Vibration Problems Z. Dimitrovová et al. (eds.) Lisbon, Portugal, 9-12 September 2013 DIFFRACTION OF PLANE SH WAVES BY A CIRCULAR CAVITY IN QUARTER-INFINITE MEDIUM Hasan
More informationSOME INEQUALITIES FOR THE q-digamma FUNCTION
Volume 10 (009), Issue 1, Article 1, 8 pp SOME INEQUALITIES FOR THE -DIGAMMA FUNCTION TOUFIK MANSOUR AND ARMEND SH SHABANI DEPARTMENT OF MATHEMATICS UNIVERSITY OF HAIFA 31905 HAIFA, ISRAEL toufik@mathhaifaacil
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More information2. If the values for f(x) can be made as close as we like to L by choosing arbitrarily large. lim
Limits at Infinity and Horizontal Asymptotes As we prepare to practice graphing functions, we should consider one last piece of information about a function that will be helpful in drawing its graph the
More informationThe Phasor Analysis Method For Harmonically Forced Linear Systems
The Phasor Analysis Method For Harmonically Forced Linear Systems Daniel S. Stutts, Ph.D. April 4, 1999 Revised: 10-15-010, 9-1-011 1 Introduction One of the most common tasks in vibration analysis is
More informationProblem Set 1 October 30, 2017
1. e π can be calculated from e x = x n. But that s not a good approximation method. The n! reason is that π is not small compared to 1. If you want O(0.1) accuracy, then the last term you need to include
More informationB2.III Revision notes: quantum physics
B.III Revision notes: quantum physics Dr D.M.Lucas, TT 0 These notes give a summary of most of the Quantum part of this course, to complement Prof. Ewart s notes on Atomic Structure, and Prof. Hooker s
More informationL = 1 2 a(q) q2 V (q).
Physics 3550 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6, 11.1 11.3 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More information3 + 4i 2 + 3i. 3 4i Fig 1b
The introduction of complex numbers in the 16th century was a natural step in a sequence of extensions of the positive integers, starting with the introduction of negative numbers (to solve equations of
More informationSecond Order Systems
Second Order Systems independent energy storage elements => Resonance: inertance & capacitance trade energy, kinetic to potential Example: Automobile Suspension x z vertical motions suspension spring shock
More informationSupplemental Notes on Complex Numbers, Complex Impedance, RLC Circuits, and Resonance
Supplemental Notes on Complex Numbers, Complex Impedance, RLC Circuits, and Resonance Complex numbers Complex numbers are expressions of the form z = a + ib, where both a and b are real numbers, and i
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall 2007
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering.4 Dynamics and Control II Fall 7 Problem Set #9 Solution Posted: Sunday, Dec., 7. The.4 Tower system. The system parameters are
More informationSolutions for homework 5
1 Section 4.3 Solutions for homework 5 17. The following equation has repeated, real, characteristic roots. Find the general solution. y 4y + 4y = 0. The characteristic equation is λ 4λ + 4 = 0 which has
More informationDifferential Equations
Electricity and Magnetism I (P331) M. R. Shepherd October 14, 2008 Differential Equations The purpose of this note is to provide some supplementary background on differential equations. The problems discussed
More informationSelected Topics in Physics a lecture course for 1st year students by W.B. von Schlippe Spring Semester 2007
Selected Topics in Physics a lecture course for st year students by W.B. von Schlippe Spring Semester 7 Lecture : Oscillations simple harmonic oscillations; coupled oscillations; beats; damped oscillations;
More informationC. Complex Numbers. 1. Complex arithmetic.
C. Complex Numbers. Complex arithmetic. Most people think that complex numbers arose from attempts to solve quadratic equations, but actually it was in connection with cubic equations they first appeared.
More informationHigher Order Averaging : periodic solutions, linear systems and an application
Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,
More informationMath 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
More information37. f(t) sin 2t cos 2t 38. f(t) cos 2 t. 39. f(t) sin(4t 5) 40.
28 CHAPTER 7 THE LAPLACE TRANSFORM EXERCISES 7 In Problems 8 use Definition 7 to find {f(t)} 2 3 4 5 6 7 8 9 f (t),, f (t) 4,, f (t) t,, f (t) 2t,, f (t) sin t,, f (t), cos t, t t t 2 t 2 t t t t t t t
More informationLecture 38: Equations of Rigid-Body Motion
Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can
More informationMATH 113: ELEMENTARY CALCULUS
MATH 3: ELEMENTARY CALCULUS Please check www.ualberta.ca/ zhiyongz for notes updation! 6. Rates of Change and Limits A fundamental philosophical truth is that everything changes. In physics, the change
More informationChapter 3: Capacitors, Inductors, and Complex Impedance
hapter 3: apacitors, Inductors, and omplex Impedance In this chapter we introduce the concept of complex resistance, or impedance, by studying two reactive circuit elements, the capacitor and the inductor.
More informationNon-Linear Plasma Wave Decay to Longer Wavelength
Non-Linear Plasma Wave Decay to Longer Wavelength F. Anderegg 1, a), M. Affolter 1, A. Ashourvan 1, D.H.E. Dubin 1, F. Valentini 2 and C.F. Driscoll 1 1 University of California San Diego Physics Department
More informationOutline. Introduction
Outline Periodic and Non-periodic Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Milano March 5, 4 A periodic loading is characterized by the identity p(t) = p(t + T ) where T is the period
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,
More informationThe Harmonic Oscillator
The Harmonic Oscillator Math 4: Ordinary Differential Equations Chris Meyer May 3, 008 Introduction The harmonic oscillator is a common model used in physics because of the wide range of problems it can
More informationTHE FOURIER TRANSFORM (Fourier series for a function whose period is very, very long) Reading: Main 11.3
THE FOURIER TRANSFORM (Fourier series for a function whose period is very, very long) Reading: Main 11.3 Any periodic function f(t) can be written as a Fourier Series a 0 2 + a n cos( nωt) + b n sin n
More informationLecture 34. Fourier Transforms
Lecture 34 Fourier Transforms In this section, we introduce the Fourier transform, a method of analyzing the frequency content of functions that are no longer τ-periodic, but which are defined over the
More informationNotes on Special Functions
Spring 25 1 Notes on Special Functions Francis J. Narcowich Department of Mathematics Texas A&M University College Station, TX 77843-3368 Introduction These notes are for our classes on special functions.
More information1. (10 points) Find the general solution to the following second-order differential equation:
Math 307A, Winter 014 Midterm Solutions Page 1 of 8 1. (10 points) Find the general solution to the following second-order differential equation: 4y 1y + 9y = 9t. To find the general solution to this nonhomogeneous
More informationTopics Covered in Calculus BC
Topics Covered in Calculus BC Calculus BC Correlation 5 A Functions, Graphs, and Limits 1. Analysis of graphs 2. Limits or functions (including one sides limits) a. An intuitive understanding of the limiting
More informationModule 9 : Infinite Series, Tests of Convergence, Absolute and Conditional Convergence, Taylor and Maclaurin Series
Module 9 : Infinite Series, Tests of Convergence, Absolute and Conditional Convergence, Taylor and Maclaurin Series Lecture 27 : Series of functions [Section 271] Objectives In this section you will learn
More informationter. on Can we get a still better result? Yes, by making the rectangles still smaller. As we make the rectangles smaller and smaller, the
Area and Tangent Problem Calculus is motivated by two main problems. The first is the area problem. It is a well known result that the area of a rectangle with length l and width w is given by A = wl.
More information1 Simple Harmonic Oscillator
Physics 1a Waves Lecture 3 Caltech, 10/09/18 1 Simple Harmonic Oscillator 1.4 General properties of Simple Harmonic Oscillator 1.4.4 Superposition of two independent SHO Suppose we have two SHOs described
More informationarxiv:quant-ph/ v1 20 Apr 1995
Combinatorial Computation of Clebsch-Gordan Coefficients Klaus Schertler and Markus H. Thoma Institut für Theoretische Physik, Universität Giessen, 3539 Giessen, Germany (February, 008 The addition of
More information2.3 Damping, phases and all that
2.3. DAMPING, PHASES AND ALL THAT 107 2.3 Damping, phases and all that If we imagine taking our idealized mass on a spring and dunking it in water or, more dramatically, in molasses), then there will be
More information