Week 1
Week 1
Week 1
1.1
Introduction/Series
Why this course exists
PHZ3113 was created in the early 1990s in the department of physics as part of a general revamping of the undergraduate physics curriculum. The perception of many faculty at that time was that more and more students were arriving at so-called majors courses like statistical mechanics and quantum mechanics mathematically unprepared to understand lectures and perform exercises. One generally recognized problem was what a physicists would refer to as an impedance mismatch (i.e., diculty of transporting information, energy, current etc. through an interface between two media because of very dierent properties of these media) between the physics and mathematics curricula at UF. As at many universities, mathematics is generally taught from the point of view of a mathematician, i.e. with a view to the formal structure of the subject and emphasizing aspects of interest to majors who might be considering doing research in the eld. Physicists teaching 15 years ago found that students had an appreciation for the beauty of mathematics, and could, under some circumstances, prove some elementary theorems, but were hard pressed to integrate the function sin x, much less solve the dierential equation for a damped harmonic oscillator. Thus PHZ3113 was created to teach students, not just from physics but other scientic or engineering disciplines as well, the nuts and bolts of performing mathematical operations of use to our eld, so that when they encounter a Hamiltonian matrix in PHY4604 (Quantum mechanics), they will know in their gut what a matrix is, how to nd its eigenvalues and eigenvectors, how to invert it, and take its Hermitian conjugate, not just be able to prove that the dual of a monomorphism is an epimorphism (aside: your instructor remembers this statement from his undergrad linear algebra course, but not what it means). At the same time, teaching this material gives us the opportunity to explore aspects of the relationship between physics and mathematics which go well beyond nuts and bolts. Mathematics is sometimes called the language of physics, but can also be thought of as an infrastructure, a sort of highway system which connects dierent branches of physics, enabling insights from one to be exported from one to another. Think of an example you already know: the remarkable set of conclusions which follow from Gausss law in electrostatics, which are rigorously analogous to those derived by Newton for the gravitational force. At the level of laboratory physics, gravity and electromagnetism are completely distinct physical phenomena, yet they are bound together by a common mathematical framework. The elds
1
associated with a mass and a charge both fall o with the inverse square of the distance: Gm 1 q Egrav = 2 ECoul = , (1) r 4 0 r2 where q is the charge, m is the mass, and 1/(4 0 and G are the universal constants of electrostatics and gravity, respectively. But for any vector eld E which falls o like from a point source as E = a/r2 (a is a constant) we have the mathematical result E dA = 4a
S
(2)
for the ux of E through any surface S (recall thats any closed surface S, not just a sphere). Thus we have shown that the gravitational eld inside a uniform spherical shell of mass is exactly zero, by noting that a Gaussian sphere concentric with the center and with radius such that the surface is at the eld point contains zero mass, and that therefore by symmetry (another mathematical idea), the eld must be zero everywhere. Make sure you understand the by symmetry argument; we will use it often. Basically, it appeals to logic. On the Gaussian sphere (dashed line in gure), the eld must point radially because there is no reason for it to point in any other direction. Therefore the integrand of the surface integral is a constant, so E must be zero if the total integral is to be zero. This is true everywhere inside
cha r
ge
,m
E=?
the sphere since the choice of the Gaussian sphere was arbitrary. This little result seems of little practical interest, but is in fact directly related to the more powerful principle of superposition in the theory of gravitation. Note that nothing in our argument really depended on it being a shell of mass, however. It could have been a shell of charge, and the eld E the electric eld.
2
ass
Therefore, without doing any further work, we see that because the mathematics are the same for the electric eld, the result that ECoul should be zero inside a uniform shell of charge follows immediately. Imagine what an enormous help it was to the natural philosophers constructing electromagnetic theory to build on the mathematics already invented by Newton (and Leibniz)! In fact, Newton derived the result above and many others required for his theory of gravitation without the benet of Gausss law, (logical, you say, since Gauss was born 50 years after Newton died). Most of electrostatics was also developed without it. The history is interesting but irrelevant to the point I wish to make, that there are important physical parallels between the theories of gravitation and electrostatics together due to the 1/r2 forces, and that it is the structure of the calculus of vector elds which allows us to see this. (Caveat and question for further thought (HW 1): in E& M you proved using a variant of this approach that the electric eld inside a hollow conductor is zero. Whats the gravitational analog of this statement?) 1.2
Series
Taylor Series
You all know that in most cases functions can be represented by their Taylor expansions, 1 1 f (x) = f (a) + f (a)(x a) + f (a)(x a)2 + f (3) (a)(x a)3 + . . . , 2 6 (3)
innite power series with coecients determined by the derivatives of the functions. We now investigate when this is possible in some more (but not terribly rigorous or complete) detail. Q: is it possible to determine the function described by the polynomial P (x) = (a function of x R) completely for all x if you know all its derivatives at a single point x? Without loss of generality, lets call the point we want to focus on x = 0. The coecients of the Taylor series are
N n n=0 an x
a0 = P (0) ;
a1 = P (1) (0) ;
a2 = P (2) (0)/2 ;
...
(4)
If N is nite, the polynomial P (x) is indeed dened completely at all x merely by knowing all of its derivatives at 0; nothing can go wrong.
3
1.2.2
Convergence
On the other hand, the properties of innite series can be quite dierent from those of nite ones. Consider Taylor series of function Q(x) around x = 0,
Q(x) =
n=0
xn (n) Q (0) n!
(5)
For a given x, the series of coecients might not converge to a nite value, even if the individual terms are decreasing. Standard example is the harmonic series 1 k k , which can be shown to diverge. The gure plots the result of adding up all the terms in the series up to a particular k, as a function of that upper value of k. You can pretty much see that the partial sums function looks like a log function, which will indeed diverge as its argument gets large. On the other hand, if you change the sign of every other term, the result converges, as seen in the last panel of the gure. (There is a very interesting discussion in Boas about the fact that this convergence is delicate, or conditionalif you rearrange the order of the terms you may get a dierent answer!)
On the other hand, we have plenty of examples of simple innite series which converge. Take for examine the series which results when your friends have bought a couple (exactly two) pizzas for a party, but invited too many people, to wit an innite number. Each, being semi-polite, leaves half the remaining pizza for the rest. Then in units of pizzas, the sum is 1 1 S = 1 + + + . . . ... = 2 4
2n .
n=0
(6)
Note that if the fundamental law of conservation of pizza holds, the sum S had better be 2. Since the terms are all positive, if we stop at any nite N , SN = N n < 2. n=0 2 Lets generalize this a bit. Suppose we have the geometric series 1 rN +1 1 rN +1 r = SN = = . 1r 1r 1r n=0
n N
(7)
The result is given as problem 2 (with hint) in Boas chapter 1. Now consider two cases: 1. r < 1: The second (remainder) term vanishes as N , so SN = 1/(1 r). Special case: r = 1/2 and we have our pizza result S = 2. Note for this case, the series converges absolutely, i.e. we can sum the terms in any order. 2. r > 1: You can tell from examining the terms that they are all positive and monotonically increasing. Therefore it is clear that the series diverges in this case as N . For some functions, convergence or divergence of the Taylor series may depend on the value of x, since each x produces a dierent series in general. If the series is known to converge at x, it usually does so over a range of x nearby. For example, take the Taylor expansion of x2 x3 ln(1 + x) = x + + ... 2 3
(8)
You can apply a number of tests (see Boas) to show that this series converges if |x| < 1, called an interval of convergence. The size of this region is called the radius of convergence.
1.2.3 Tests for convergence
1. Comparison test. If you know a series converges, then another series whose terms are bounded by those of the convergent one also converges. 2. Integral test. If the corresponding integral converges, so does the series. If it diverges, so does the series. Ex: 1/n2 . Examine the integral n=1
1
1 n dn = n
2
= 1,
1
(9)
so yes the sum will also converge. Note the lower limit does not matter for the test, since divergence problems arise from the upper limit, except if the integral or sum diverges at the lower limit (had we taken the lower limit to be zero the sum and integral would have both been divergent). 3. Ratio test. Dene ratio of subsequent terms n an+1 an lim = n n (10) (11)
If < 1 series converges. Application to Taylor series. Consider a function with expansion
f (x) =
n=0
bn x n .
(12)
The ratio test says that f (x) converges, i.e. the sum is a valid representation of the function, if
n+1 lim bn+1 x n <1 bn x n bn lim . |x| < n bn+1
(13) (14)
This gives interval of convergence of Taylor series. Ex: f (x) = for |x| < 2.
1.2.4 Remarks
(x)n n=0 2n
1. If two functions have the same Taylor expansion, they are the same function. If f (x) = an xn and g(x) = bn xn , then an =bn and f = g. This is n=0 n=0 not always obvious for simple series, i.e. certainly 1+5 = 2+4 does not imply term by term equality. Intuitive point is that the relation has to be valid for any x, so this constrains coecients. Well prove this later using complex variables. 2. Power series can be dierentiated or integrated term by term, & resulting series converges to derivative/integral of the function, with same radius of convergence.
6
3. Approximations for functions: stopping at nite order. Sometimes we are lazy and want to keep only a nite number of terms in an expansion, & ask what kind of an error we are making. Ex: x3 + O(x5 ) sin x = x 3! 1 3 sin = 0.499, 6 6 6 6 (17) (18)
although the exact solution is 0.5. The remainder of the series makes only an extremely small change, at least near /6. Note this particular series converges lim fast enough, i.e. n a /a , such that the radius of convergence is
n n+1
innite! 4. Nested expansions. Sometimes its easier to expand internal parts of a composite function in power series, rather than take derivatives formally as in Eq. 3. Ex: x exp 1x x 1 = 1+ + 1 x 2! x 1x
2
1 + 3!
x 1x
...
(19)
= 1 + x(1 + x + x2 . . . ) + 3 13 = 1 + x + x2 + x3 . . . . 2 6
1 2 x (1 + x + x2 . . . )2 + . . . (20) 2! (21)