LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS

Size: px
Start display at page:

Download "LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS"

Transcription

1 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October 7, 202 Prof. Aan Guth LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS In these notes I wi describe the separation of variabe technique for soving Lapace s equation, using spherica poar coordinates. The soutions wi invove Legendre poynomias for cases with azimutha symmetry, and more generay they wi invove spherica harmonics. I wi construct these soutions using traceess symmetric tensors, but in Lecture Notes 8 I describe how the soutions in this form reate to the more standard expressions in terms of Legendre poynomias and spherica harmonics. (Logicay Lecture Notes 8 shoud come after these notes, athough they were posted first.) If you are starting from scratch, I think that the traceess symmetric tensor method is the simpest way to understand this mathematica formaism. If you aready know spherica harmonics, I think that you wi find the traceess symmetric tensor approach to be a usefu addition to your arsena of mathematica methods. The symmetric traceess tensor approach is particuary usefu if one needs to extend the formaism beyond what we wi be doing for exampe, there are anaogues of spherica harmonics in higher dimensions, and there are aso vector spherica harmonics that are usefu for expanding vector functions of ange. Vector spherica harmonics are used in the most genera treatments of eectromagnetic radiation, athough we wi not be introducing them in this course. I don t know a reference for the traceess symmetric tensor method, which is the main reason I am writing these notes. For the standard method, imited to the case of azimutha symmetry, our textbook by Griffiths shoud be sufficient. If you woud ike to see an additiona reference on spherica harmonics, which are needed when there is no azimutha symmetry, then I woud recommend J.D. Jackson, Cassica Eectrodynamics, rd Edition (John Wiey & Sons, 999), Sections.,.2,.5, and.6.. LAPLACE S EQUATION IN SPHERICAL COORDINATES: In spherica coordinates, Lapace s equation can be written as 2 ϕ ϕ(r, θ, φ) = r 2 + r2 r r r 2 2 θ ϕ =0, (9.) where the anguar part is given by 2 ϕ 2 ϕ θ ϕ sin θ + sin θ θ θ sin 2 θ φ 2. (9.2)

2 8.07 LECTURE NOTES 9, FALL 202 p. 2 It woud be more ogica to write 2 θ as 2 θ,φ, but it woud be tiresome to write it that way. It is sometimes usefu to rewrite the first term in Eq. (9.) using ) 2 ϕ 2 (r = (rϕ). (9.) r 2 r r r r 2 To use the method of separation of variabes, we can seek a soution of the form Then Lapace s equation can be written as ϕ(r, θ, φ) =R(r)F (θ, ϕ). (9.4) r 2 d 0= 2 dr ϕ = r 2 + RF R dr dr F 2 θ F. (9.5) Since the first term on the right-hand side depends ony on r, and the second term depends ony on θ and φ, the ony way that the equation can be satisfied is if each term is a constant. Thus we can write F 2 θ F = C θ, (9.6) d ( R r 2 d ) = C θ. (9.7) R dr dr 2. THE EXPANSION OF F (θ, φ): We now wish to find the most genera soution to the equation 2 θ F = C θ F, (9.8) which is a rewriting of Eq. (9.6). If such a function F can be found, we say that F is an eigenfunction of the operator 2 θ,witheigenvaue C θ. A function of anges (θ, φ) can equivaenty be thought of as a function of the unit vector nˆ that points in the direction of θ and φ, which can be written expicity as nˆ =sinθ cos φ ê +sinθ sin φ ê 2 +cosθ ê, (9.9) where ê,ê 2,andê can aso be written as ê x,ê y,andê z. I am abeing the unit vectors using the numbers, 2, and when I am thinking about summing over the indices, and otherwise I use x, y, andz. I now caim that the most genera function of (θ, φ) can be written as a power series in nˆ, or more precisey as a power series in the components of nˆ. I wi not prove this,

3 8.07 LECTURE NOTES 9, FALL 202 p. but it is true at east for square-integrabe piece-wise continuous functions F (θ, φ). Such a power series can be written as F (ˆn) =C (0) + C () i nˆi + C ( ij nˆinˆ j C ) i i 2...i nˆi nˆi2...nˆi +..., (9.0) where repeated indices are summed from to (as Cartesian coordinates). Note that C i i 2...i nˆ i nˆi 2...nˆi represents the genera term in the series, where the first three terms correspond to =0, =,and = 2. The indices i, i 2,..., i represent different indices, ike i and j; since they are repeated, they are each summed from to. The coefficients C i i 2...i are caed tensors, and the number of indices is caed the () rank of the tensor. Note that C (0), C i,andc ij are specia cases of tensors, athough they can aso be considered a scaar, a vector, and a matrix. It is possibe to impose restrictions on the coefficients of Eq. (9.0) without actuay restricting what can appear on the right-hand side of the equations. In particuar, we wi insist that ) The tensors C i i 2...i are symmetric under any reordering of the indices: C i i where {j,j 2,...,j } is any permutation of {i,i 2,...,i }. C i i 2...i C 2...i = j j 2...j, (9.) 2) The tensors are traceess, in the sense that if any two indices are set equa to each other and summed, the resut is equa to zero. Since the tensors are aready assumed to be symmetric, it does not matter which indices are summed, so we can choose the ast two: C i i 2...i =0. (9.2) To expain why these restrictions on the C s do not impose any restriction on the righthand side of Eq. (9.0), I wi use the exampe of C ij, but I think you wi be abe to see that that the argument appies to a. The insistence that C ij is symmetric can be seen to make no difference to the right-hand side of Eq. (9.0), because C ij mutipies the symmetric tensor nˆinˆj. Thus,ifC ij had an antisymmetric part, it woud not contribute to the right-hand side of Eq. (9.0). The requirement of traceessness is ess obvious, but suppose that were not traceess. Then we coud write C ij 2jj C ii = λ =0. (9.)

4 8.07 LECTURE NOTES 9, FALL 202 p. 4 We coud then define a new quantity, ij = C ij λδ ij. (9.4) C It foows that C ij is traceess: C ii = C ii λδii = λ λδii =0, (9.5) since δ ii =. The origina term C ij nˆinˆj can then be expressed in terms of C ij : [ ˆ ˆ = C ij n i n j C ij + λδ ij ] nˆinˆ j = C ij nˆinˆj + λ, (9.6) whereweusedthefactthatδ ij nˆinˆj =,sincenˆ is a unit vector. The extra term, λ, can then be absorbed into a redefinition of C (0) : Finay, we can write C (0) = C (0) + λ. (9.7) C (0) + () () (2 C nˆi + ) i C ij nˆ i nˆj = C (0) + C i nˆi + C ij nˆinˆj, (9.8) so we can insist that the tensor that mutipies nˆinˆj be traceess with no restriction on what functions can be expressed in this form. I wi ca the th term of this expansion F (ˆn), so F (ˆn) = i i 2...i i i 2 i C nˆ nˆ...nˆ. (9.9). EVALUATION OF 2 θ F(ˆn): To evauate 2 θ F (ˆn), we are going to take advantage of a convenient trick. Instead of deaing directy with F (ˆn), we wi instead introduce a radia variabe r, usingitto define a coordinate vector r = rn ˆ. (9.20) Foowing the notation of Griffiths (see his Eq. (.9)), I wi denote the coordinates of r by x, y, andz, or in index notation I wi ca them x i.so r = x i ê i = x ê + x 2 ê 2 + x ê. (9.2)

5 8.07 LECTURE NOTES 9, FALL 202 p. 5 Then, given any F (ˆn) of the form given in Eq. (9.9), we can define a function F ( r )by F ( r )=C i i 2...i x i x i 2...x i = r F (ˆn). (9.22) Note that C i i...i isthesamerank traceess symmetric tensor used to define F 2 (ˆn), but we are defining F ( r ) by mutipying C i i...i by x i x i2...x i and then summing over 2 indices, instead of mutipying by nˆ i nˆi2...nˆi and then summing. Now we can make use of Eq. (9.), which reates the fu Lapacian 2 to the anguar Lapacian, 2 θ. We wi find that in this case the fu Lapacian and the radia derivative piece of Eq. (9.) wi both be simpe, so we wi be abe to determine the anguar Lapacian by evauating the other terms in Eq. (9.). To evauate 2 F ( r ), we start with = 0. Ceary 2 F 0( r) = 2 C (0) =0, (9.2) since the derivative of a constant vanishes. Simiary for =, 2 F 2 () ( r) = C i x i =0, (9.24) since the first derivative produces a constant, so the second derivative vanishes. The first nontrivia case is =2: 2 F ( r )= 2 2 C ij x i x j = C ij (x i x j ) x m x m (9.25) = C ij [δ im x j + δ jm x i ] x m =2 C ij [ δ im δ jm ]=2C ii =0, where in the ast step we used the a-important fact that C ij is traceess. We wi ook at one more case, and then I hope it wi be cear that it generaizes. For =, 2 F r 2 () = C ijk x i x j x k = () C ijk (x i x j x k ) x m x m ) = C () ijk (δ im x j x k + (terms that symmetrize in ijk) x m ) = C () ijk (δ im δ jm x k + (terms that symmetrize in ijk) ) = C () ijk (δ ij x k + (terms that symmetrize in ijk) = () C iik x k + (terms that symmetrize in ijk) =0, (9.26)

6 8.07 LECTURE NOTES 9, FALL 202 p. 6 where again it is the traceessness of C i i 2...i that caused the term to vanish. Thinking about the genera term, one can see that after the derivatives are cacuated, there are 2 factors of x i that remain, but there are sti indices on C i i 2...i. Since a indices are summed, there are aways two indices on C i i 2...i which are contracted (i.e., set equa to each other) and summed, which causes the resut to vanish by the traceessness condition. The bottom ine, then, is that 2 F ( r) =0fora. (9.27) To see what this says about 2 θ F (ˆn), reca that F ( r) =r F (ˆn). Using Eq. (9.), we can write 2 2 F 0= F ( r) (r) = r + 2 θ r F ( r) 2 r r r 2 ( d ) = r 2 dr F (ˆn)+ r r 2 dr dr r 2 2 θ F (ˆn) = r 2[ ] ( +)F (ˆn)+ 2 θ F (ˆn), and therefore 2 θ F (ˆn) = ( +)F (ˆn). (9.28) ( ) Thus, we have found the eigenfunctions F (ˆ) = n C i i 2...i nˆi nˆi2...nˆi and eigenvaues ( +) of the differentia operator 2 θ. This is a very usefu resut! 4. GENERAL SOLUTION TO LAPLACE S EQUATION IN SPHERICAL COORDINATES: Now that we know the eigenfunctions of 2 θ, we can return to the soution to Lapace s equation by the separation of variabes in spherica coordinates. We now know that in Eq. (9.6), the ony aowed vaues of C θ are ( +), where is an integer. Thus, Eq. (9.7) becomes ( d r 2 dr ) = ( +)R, (9.29) dr dr and we can ook for soutions by trying R(r) =r p. We find consistency provided that p(p +)=( +), (9.0)

7 8.07 LECTURE NOTES 9, FALL 202 p. 7 which is a quadratic equation with two roots, p = and p = ( + ). Since we found two soutions to a second order inear differentia equation, we know that any soution can be written as a inear sum of these two. Thus we can write R (r) =A r B +, (9.) r + aowing for different vaues of A and B for each. The most genera soution to Lapace s equation, in spherica coordinates, can then be written as Φ( r )= B A r + =0 r + C i i 2... i nˆi nˆi 2...nˆi, (9.2) where the A s and B s are arbitrary constants, and each C i i 2...i is an arbitrary traceess symmetric tensor. In Lecture Notes 8 it is shown expicity that the th term here, when compared with the standard expansion in the spherica harmonic functions Y m (θ, φ), corresponds to the sum of a terms with the same, but for a m. TheY m s are defined for integer vaues of m from to, sothereare2 + terms for each vaue of. 5. COUNTING THE NUMBER LINEARLY INDEPENDENT TRACELESS SYMMETRIC TENSORS: Using the fact that C i i 2...i is traceess and symmetric, we can determine how many ineary independent such tenso rs exist, for any given. In other words, how many rea constants are needed to parameterize the most genera traceess symmetric tensor of rank? We begin by cacuating N sym, the number of ineary independent symmetric tensors of rank. Note that we are postponing the consideration of traceessness. For grounding, we can start by saying that it takes one number to specify S (0),arank0 symmetric tensor, because it is just a number. (I am using S for symmetric tensors, () whie reserving C for traceess symmetric tensors.) It takes numbers to specify S i, () () () since the vaues S, S 2,andS can each be specified independenty. For S ij, however, we see the constraints of symmetry: S ij has to equa S ji, so there are fewer (2 (2 than 9 independent vaues; there are 6, which can be taken to be S, 2, ) ) S S, S 22, S 2, S,with S 2, S,and S 2 determined by symmetry. Since the order of the indices does not matter, we can aways ist the indices in ascending order (as I did for

8 8.07 LECTURE NOTES 9, FALL 202 p. 8 S ij ) and then each independent entry wi occur once. When the indices are so written in ascending order, then the index vaues are competey determined if we just specify how many indices are equa to, how many are equa to 2, and how many are equa to. If it heps, we can imagine three hats, abeed, 2, and, and indistinguishabe bas, representing the indices of a rank tensor. There is then a : correspondence between independent tensor eements and the different ways that the bas can be put into the hats. For exampe, if there are 9 bas, with in the first hat, 2 in the second, and 4 in the third, then this arrangement of bas corresponds to S (9) 22.Sonowwejusthave to figure out how many different ways we can put indistinguishabe bas into hats. One way of counting the bas-in-hats probem is to imagine first abeing each ba with a number, so they are no onger indistinguishabe. We aso introduce 2 dividers, where 2 is one ess than the number of hats. Initiay we wi aso assign numbers to the 2 dividers. Thinking of the bas and dividers together, we have +2 distinguishabe objects. We can imagine isting them in a possibe orderings, and with + 2 distinguishabe objects there are ( + 2)! orderings. For each ordering there is an equivaent bas-inhats assignment. The bas to the eft of the eft-most divider are assigned to hat, the bas between the two dividers are assigned to hat 2, and the bas to the right of the right-most divider are assigned to hat. We have of course overcounted, since many different orderings of our + 2 objects wi ead to the same number of bas in each hat. However, we can see exacty by how much we have overcounted. We can re-order the bas without changing the number of bas in each hat, and we can interchange the two dividers. So we have overcounted by a factor of 2!. Thus, ( +2)! N sym = = ( +)( +2). (9.) 2! 2 It is easiy checked that this gives N sym (0) =, N sym () =, and N sym = 6, consistent with the exampes we started with. To impose traceessness, we require that our traceess tensors aso satisfy S i...i 2 jj =0. (9.4) How many conditions is this? One can see that the tensor on the eft is a symmetric tensor of rank 2, with free indices i...i 2. The number of conditions is then Nsym( 2). The number of ineary indendent traceess symmetric tensors is then given by N traceess sym =N sym N sym ( 2) = ( +)( +2) ( ) 2 2 = 2 +. (9.5) So the correspondence with the standard Y m s is consistent, as it woud have to be. The spherica harmonic expansion is just a rewriting of Eq. (9.2), with a particuar choice of basis for the 2 + independent traceess symmetric tensors of rank.

9 8.07 LECTURE NOTES 9, FALL 202 p SPECIAL CASE: AZIMUTHAL SYMMETRY: Azimutha symmetry means symmetry under rotation about an axis, which we wi take to be the z-axis. Equivaenty, we can say that a probem is azimuthay symmetric if nothing depends on the coordinate φ. To speciaize the genera expansion (9.0) for a function of (θ, φ) to the azimuthay symmetric case, we need to construct traceess symmetric tensors which are invariant under rotations about the z-axis. One straightforward way to this is to buid the traceess symmetric tensor from the vector ẑ, the unit vector in the z direction. Note that the x-, y-, and z- componentsofẑ are 0, 0, and, respectivey, so ẑ i = δ i. [You may have noticed an inconsistency in my notation, as earier (e.g., Eq. (9.9)) Iusedê or ê z for the unit vector in the z direction. In this case my inconsistency was intentiona, with two motivations. First, previousy we never needed a notation for the components of a unit basis vector, but here we wi. It woud be a rea pain to write (ê z ) i. Second, in Lecture Notes 8 I describe a convenient way to construct a basis for the traceess symmetric tensors, which invoves the use of a basis for vectors consisting of ẑ and two compex vectors û + and û. So, the use of ẑ rather than ê z wi remind us that we are thinking about the (û +, û, ẑ) basis, rather than the (ê x, ê y, ê z ) basis.] Arank tensor can be constructed from ẑ simpy by taking the product ẑ i ẑ i2...ẑ i. This is symmetric, and can be made traceess by extracting the traceess part. Extracting the traceess part means subtracting terms proportiona to one or more Kronecker δ- functions in such a way that the resut is traceess. It gets rather compicated to describe how this can be done for a genera symmetric tensor of arbitrary rank, so I wi just iustrate it by exampe. Tensors of rank 0and (i.e., scaars and vectors) are by definition traceess. I wi use cury brackets {...} to denote the traceess symmetric part of... Thus, { } =, { ẑ i } =ẑ i. (9.6) But for rank 2, the trace of ẑ i ẑ j is equa to ẑ i ẑ i =ẑ ẑ =. But we can subtract a constant times δ ij so that the resut is traceess: { ẑ i ẑ j } =ẑ i ẑ j δ ij. (9.7) The coefficient is /, because the trace of δ ij is δ ii =. (9.8)

10 8.07 LECTURE NOTES 9, FALL 202 p. 0 For rank, ẑ i ẑ j ẑ k has trace ẑ i ẑ i ẑ k =ẑ k, but we can make it traceess with a subtraction { ẑ i ẑ j ẑ k } =ẑ i ẑ j ẑ k 5( ẑi δ jk +ẑ j δ ik +ẑ k δ ij ). (9.9) The subtraction must of course be symmetrized, as shown, since we are trying to construct a traceess symmetric tensor. To verify that /5 is the right coefficient to make the expression traceess, we can take its trace. Since it is symmetric we can sum over any pair of indices. I wi choose to sum over i and j: δ ij { ẑ i ẑ j ẑ k } =ẑ i ẑ i ẑ k ( 5 z ˆ i δ ) ik +ẑ i δ ik +ẑ k δ ii (9.40) =ẑ k (ẑ k +ẑ k +ẑ k )=0. 5 For rank 4 there is the option of subtracting terms with either one or two Kronecker δ-functions, and both are needed to give a traceess resut. We can start with arbitrary coefficients, and see what they have to be: ( { ẑ i ẑ j ẑ k ẑ m } =ẑ i ẑ j ẑ k ẑ m + c ẑi ẑ j δ km +ẑ i ẑ k δ mj +ẑ i ẑ m δ jk +ẑ j ẑ k δ im ) (9.4) +ẑ j ẑ m δ ik +ẑ k ẑ m δ ij + c2 δij δ km + δ ik δ jm + δ im δ jk. Within each set of parentheses, the terms are chosen to make the expression symmetric in i, j, k, andm. If we cacuate the trace over i and j, we find ( δ ij { ẑ i ẑ j ẑ k ẑ m } =ẑ k ẑ m + c δkm +ẑ k ẑ m +ẑ k ẑ m +ẑ k ẑ m ) +ẑ k ẑ m +ẑ k ẑ m + c2 δkm + δ km + δ km (9.42) = +7c ẑ k ẑ m + c +5c 2 δ km. For the expression to vanish for a k and m, the c = /7 andc 2 =/5. Thus, two terms must vanish separatey, so { ẑiẑ j ẑ k ẑ m } =ẑ i ẑ j ẑ k ẑ m ( 7 ẑ i ẑ j δ km +ẑ i ẑ k δ mj +ẑ i ẑ m δ jk +ẑ j ẑ k δ im ) ( ) +ẑ j ẑ m δ ik +ẑ k ẑ m δ ij + δ ij δ km + δ ik δ jm + δ im δ jk. 5 (9.4) It can be shown that any traceess symmetric tensor of rank that is invariant under rotations about the z-axis is proportiona to { ẑ i...ẑ i }. (To see this, note that in Lecture Notes 8 we construct a compete (2 + )-dimensiona basis for the traceess symmetry tensors of rank. They depend on the azimutha ange φ as z m (φ) e imφ, with m taking integer vaues from to. Since these functions of φ are orthogona in

11 8.07 LECTURE NOTES 9, FALL 202 p. the sense that 2π zm (φ)z 0 m (φ)dφ =2πδ m m, any traceess symmetric tensor of rank that is independent of φ must proportiona to the m = 0basis tensor.) Since Eq. (9.0) tes us how to expand an arbitrary function of θ and φ in terms of traceess symmetric tensors, we can now say that functions of θ aone (i.e., azimuthay symmetric functions of nˆ) can be expanded as F (θ) =c 0 + c { ẑ i }nˆi + c 2 { ẑ i ẑ j } nˆinˆj c { ẑ i...ẑ i } nˆi...nˆi +..., (9.44) where the c s are constants. This corresponds to what is standardy caed an expansion in Legendre poynomias. In Lecture Notes 8 I show exacty how to reate these terms to the standard conventions for normaizing the Legendre poynomias, but we can see here exacty what these functions are. Using ẑ nˆ =cosθ, wehave { } = { ẑ i } nˆi =cosθ 2 { ẑiẑ j } nˆ inˆj =cos θ (9.45) { ẑ i ẑ j ẑ k } nˆinˆjnˆk =cos θ cos θ { ẑ i ẑ j ẑ k ẑ m } nˆ inˆ jnˆknˆm =cos θ cos θ Up to a normaization convention described in Lecture Notes 8, these are the Legendre poynomias P (cos θ). 7. THE MULTIPOLE EXPANSION: The most genera soution to Lapace s equation, in spherica coordinates, was given as Eq. (9.2). We now wish to appy that resut to a common situation: suppose we have a charged object, and we wish to describe the potentia outside of the object. Let s say for definiteness that the charge of the object is entirey contained within a sphere of radius R, centered at the origin. In that case Lapace s equation wi hod for a r>r,sothere shoud be a soution of the form of Eq. (9.2) that is vaid throughout this region. At infinity the potentia of a ocaized charge distribution wi aways approach a constant, which we can take to be zero, we can see that the A coefficients that appear in Eq. (9.2) must a vanish. The B factors can be absorbed into the definition of C i...i,sowecan write the expansion as Φ( r )= r + =0 C i...i n ˆi...nˆi. (9.46)

12 8.07 LECTURE NOTES 9, FALL 202 p. 2 Since each successive term comes with an extra factor of /r, at arge distances the sum is dominated by the first term or maybe the first few terms. A the information about the charge distribution of the object is contained in the C i...i, so knowedge of the first few C i the ob...i is enough to describe the fied at arge distances, no matter how compicated ject. The first few terms of this series have specia names: the = 0term is the monopoe term, the = term is the dipoe, the = 2 term is the quadrupoe, and the =term is the octupoe. If we want to cacuate the C i...i in terms of the charge distribution, we can start with the genera equation for the potentia of an arbitrary charge distribution: ρ( r ) V ( r) = d x. (9.47) 4πɛ0 r r The mutipoe expansion can then be derived by expanding / r r in a power series in r. I begin by doing it as Griffiths does, which gives the simpest but not the most usefu form of the mutipoe expansion. Griffiths rewrote the denominator as = =, (9.48) r r r 2 + r 2 2 r r r 2 + r 2 2rr cos θ where r and r are the engths of the vectors r and r, respectivey, and θ is the ange between these vectors. Next he used the fact that the Legendre poynomias can be defined by the generating function g(x, λ) =, (9.49) +λ 2 2λx which means that the Legendre poynomias P (x) can be obtained by expanding g(x, λ) in a power series in λ: g(x, λ) = = λ P(x). (9.50) + λ 2 2λx =0 Eq. (9.50) is sometimes taken as the definition of the Legendre poynomias, and sometimes it is derived from another definition. In any case, if we accept Eq. (9.50) as vaid, then r = = P (cos θ). (9.5) r r r + r 2 r r 2 r cos θ r =0 r

13 8.07 LECTURE NOTES 9, FALL 202 p. Inserting this reation into Eq. (9.47), we find V ( r) = 4 πɛ 0 r + =0 r ρ( r )P(cos θ )d x. (9.52) This is the easiest way that I know to show that there is an expansion of V ( r) inpowersof /r, but the compication is that cos θ appears inside the integra. If we coud impement Eq. (9.46), we woud be abe to cacuate (or maybe measure) a sma number of the quantities C i...i, and then we woud be abe to evauate V ( r ) at arge distances in any direction. To use Eq. (9.52) directy, however, one woud have to repeat the integration for every direction of r. Griffiths works around this probem by massaging the formua to extract the monopoe and dipoe terms, and in Probem of Probem Set 5 you had the opportunity to carry this out for the quadrupoe and octupoe terms. The standard method of improving Eq. (9.52) is to use spherica harmonics, but here I wi derive the equivaent reations using the traceess symmetric tensor approach. Instead of expanding / r r in powers of r, we wi think of it as a function of three variabes the components r i of r, and we wi expand it as a Tayor series in variabes. To make the formaism cear, I wi define the function f( r ). (9.5) r r The function can then be expanded in a power series using the standard muti-variabe Tayor expansion: f 2 f f( r )=f( 0)+ x i + x i x +..., (9.54) x j i r x x = 0 2! i j where the repeated indices are summed. To separate the anguar behavior, we write r = 0 x i = r nˆ i, (9.55) so Eq. (9.54) becomes f( r )=f( 0)+r f x i r 2 2 f + nˆ n ˆ ! x x i j r =0 i j r = 0 (9.56) The notation can now be simpified by noting that since f is a function of r r,the derivatives with respect to x i can be repaced by derivatives with respect to x i with a change of sign: f = = =. (9.57) r r x i r r r x = 0 i r x i x i r = 0

14 8.07 LECTURE NOTES 9, FALL 202 p. 4 This aows us to write the derivatives in the expansion (9.56) much more simpy. The th derivative is found by repeating the above operation times: f =( ). (9.58) x i... x i r = x i... x i r 0 Combining Eqs. (9.58) with (9.56), we can write ( ) r ( ) = nˆ r r! x... x r i i =0 i...nˆ i. (9.59) Note that the quantity in parentheses in the equation above is traceess, because ( ) = 2 x i x i x i+2... x i r x i+2... x i r (9.60) = 2 =0, xi+2... x i r because 2 (/ r ) = 0except at r = 0. So we can see the traceess symmetric tensor formaism emerging. To evauate this quantity, we wi work out the first severa terms unti we recognize the pattern. We write r rn ˆ, (9.6) and adopt the abbreviation i. (9.62) x i It is usefu to start by evauating the derivatives of the basic quantities r and nˆi: / r = i(x j x j ) 2 xi i = ( xkx k ) /2 i (x j x j )= 2xjδ ij = =ˆn i, 2 2r r ( xj ) δ ij nˆ = = x r = δ nˆ nˆ. r r r 2 r i j i j i ij i j (9.6) It is then straightforward to show that i = nˆi, r r 2 i j = r { r nˆinˆj }, (9.64) 5 i j k = { nˆiˆjnˆk n }, r r4

15 8.07 LECTURE NOTES 9, FALL 202 p. 5 where {}denotes the traceess symmetric part, and the reevant cases are shown expicity in Eqs. (9.6) (9.9). It becomes cear that the genera formua, which can be proven by induction, is ( ) (2 )!! = { nˆii...nˆi x i... x r + }, (9.65) r where i (2)! (2 )!! (2 )(2 )(2 5)...=, with ( )!!. (9.66) 2! Inserting this resut into Eq. (9.59), we find = (2 )!! r { nˆi...nˆ r r i! r } + nˆi...nˆ i. (9.67) =0 One can write this more symmetricay by writing )!! r = (2 nˆ...nˆ nˆ...nˆ r r { i i! r }{ i i }, (9.67) + =0 since { nˆ i...nˆ i } differs from nˆ i...nˆ i by terms proportiona to Kronecker δ-functions, which vanish when summed with the traceess tensor { nˆi...nˆi }. Starting with Eq. (9.67), one can if one wishes drop the cury brackets around either factor (but not both!). Inserting this expression for / r r into Eq. (9.47), we have the fina resut V ( r )= 4πɛ0 C r + i...i n =0 ˆi...nˆi, (9.68) where (2 C i...i = )!! ρ( r ) { r i!... r i } d x. (9.69) Note that the coefficient in the above expression can aso be written as (2 )!! (2)! =. (9.70)! 2 (!) 2

16 8.07 LECTURE NOTES 9, FALL 202 p. 6 For purposes of iustration, I wi write out the first two terms the monopoe and dipoe terms in a bit more detai. The monopoe term can be written as Q V mono ( r) =, (9.7) 4πɛ 0 r where Q = C (0) = ρ( r )d x. (9.72) The dipoe term is where p nˆ V dip ( r) =, (9.7) 4πɛ 0 r 2 p i = C () i = ρ( r ) r i d x. (9.74)

17 MIT OpenCourseWare Eectromagnetism II Fa 202 For information about citing these materias or our Terms of Use, visit:

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October, 202 Prof. Aan Guth LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

More information

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

$, (2.1) n="# #. (2.2)

$, (2.1) n=# #. (2.2) Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n. Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

More information

Legendre Polynomials - Lecture 8

Legendre Polynomials - Lecture 8 Legendre Poynomias - Lecture 8 Introduction In spherica coordinates the separation of variabes for the function of the poar ange resuts in Legendre s equation when the soution is independent of the azimutha

More information

David Eigen. MA112 Final Paper. May 10, 2002

David Eigen. MA112 Final Paper. May 10, 2002 David Eigen MA112 Fina Paper May 1, 22 The Schrodinger equation describes the position of an eectron as a wave. The wave function Ψ(t, x is interpreted as a probabiity density for the position of the eectron.

More information

Physics 116C Helmholtz s and Laplace s Equations in Spherical Polar Coordinates: Spherical Harmonics and Spherical Bessel Functions

Physics 116C Helmholtz s and Laplace s Equations in Spherical Polar Coordinates: Spherical Harmonics and Spherical Bessel Functions Physics 116C Hemhotz s an Lapace s Equations in Spherica Poar Coorinates: Spherica Harmonics an Spherica Besse Functions Peter Young Date: October 28, 2013) I. HELMHOLTZ S EQUATION As iscusse in cass,

More information

Jackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Jackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell Jackson 4.10 Homework Probem Soution Dr. Christopher S. Baird University of Massachusetts Lowe PROBLEM: Two concentric conducting spheres of inner and outer radii a and b, respectivey, carry charges ±.

More information

Physics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27

Physics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27 Physics 55 Fa 7 Homework Assignment #5 Soutions Textook proems: Ch. 3: 3.3, 3.7, 3.6, 3.7 3.3 Sove for the potentia in Proem 3., using the appropriate Green function otained in the text, and verify that

More information

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

More information

Radiation Fields. Lecture 12

Radiation Fields. Lecture 12 Radiation Fieds Lecture 12 1 Mutipoe expansion Separate Maxwe s equations into two sets of equations, each set separatey invoving either the eectric or the magnetic fied. After remova of the time dependence

More information

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization

More information

More Scattering: the Partial Wave Expansion

More Scattering: the Partial Wave Expansion More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction

More information

Lecture 8 February 18, 2010

Lecture 8 February 18, 2010 Sources of Eectromagnetic Fieds Lecture 8 February 18, 2010 We now start to discuss radiation in free space. We wi reorder the materia of Chapter 9, bringing sections 6 7 up front. We wi aso cover some

More information

Physics 505 Fall Homework Assignment #4 Solutions

Physics 505 Fall Homework Assignment #4 Solutions Physics 505 Fa 2005 Homework Assignment #4 Soutions Textbook probems: Ch. 3: 3.4, 3.6, 3.9, 3.0 3.4 The surface of a hoow conducting sphere of inner radius a is divided into an even number of equa segments

More information

Formulas for Angular-Momentum Barrier Factors Version II

Formulas for Angular-Momentum Barrier Factors Version II BNL PREPRINT BNL-QGS-06-101 brfactor1.tex Formuas for Anguar-Momentum Barrier Factors Version II S. U. Chung Physics Department, Brookhaven Nationa Laboratory, Upton, NY 11973 March 19, 2015 abstract A

More information

Higher dimensional PDEs and multidimensional eigenvalue problems

Higher dimensional PDEs and multidimensional eigenvalue problems Higher dimensiona PEs and mutidimensiona eigenvaue probems 1 Probems with three independent variabes Consider the prototypica equations u t = u (iffusion) u tt = u (W ave) u zz = u (Lapace) where u = u

More information

FFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection

FFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection FFTs in Graphics and Vision Spherica Convoution and Axia Symmetry Detection Outine Math Review Symmetry Genera Convoution Spherica Convoution Axia Symmetry Detection Math Review Symmetry: Given a unitary

More information

Math 124B January 31, 2012

Math 124B January 31, 2012 Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat

More information

Physics 506 Winter 2006 Homework Assignment #6 Solutions

Physics 506 Winter 2006 Homework Assignment #6 Solutions Physics 506 Winter 006 Homework Assignment #6 Soutions Textbook probems: Ch. 10: 10., 10.3, 10.7, 10.10 10. Eectromagnetic radiation with eiptic poarization, described (in the notation of Section 7. by

More information

Physics 566: Quantum Optics Quantization of the Electromagnetic Field

Physics 566: Quantum Optics Quantization of the Electromagnetic Field Physics 566: Quantum Optics Quantization of the Eectromagnetic Fied Maxwe's Equations and Gauge invariance In ecture we earned how to quantize a one dimensiona scaar fied corresponding to vibrations on

More information

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed

More information

14 Separation of Variables Method

14 Separation of Variables Method 14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt

More information

221B Lecture Notes Notes on Spherical Bessel Functions

221B Lecture Notes Notes on Spherical Bessel Functions Definitions B Lecture Notes Notes on Spherica Besse Functions We woud ike to sove the free Schrödinger equation [ h d r R(r) = h k R(r). () m r dr r m R(r) is the radia wave function ψ( x) = R(r)Y m (θ,

More information

C. Fourier Sine Series Overview

C. Fourier Sine Series Overview 12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a

More information

4 Separation of Variables

4 Separation of Variables 4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE

More information

V.B The Cluster Expansion

V.B The Cluster Expansion V.B The Custer Expansion For short range interactions, speciay with a hard core, it is much better to repace the expansion parameter V( q ) by f(q ) = exp ( βv( q )) 1, which is obtained by summing over

More information

4 1-D Boundary Value Problems Heat Equation

4 1-D Boundary Value Problems Heat Equation 4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x

More information

Applied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation

Applied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation 22.101 Appied Nucear Physics (Fa 2006) Lecture 7 (10/2/06) Overview of Cross Section Cacuation References P. Roman, Advanced Quantum Theory (Addison-Wesey, Reading, 1965), Chap 3. A. Foderaro, The Eements

More information

Math 124B January 17, 2012

Math 124B January 17, 2012 Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia

More information

Module 22: Simple Harmonic Oscillation and Torque

Module 22: Simple Harmonic Oscillation and Torque Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque

More information

Section 6: Magnetostatics

Section 6: Magnetostatics agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The

More information

2M2. Fourier Series Prof Bill Lionheart

2M2. Fourier Series Prof Bill Lionheart M. Fourier Series Prof Bi Lionheart 1. The Fourier series of the periodic function f(x) with period has the form f(x) = a 0 + ( a n cos πnx + b n sin πnx ). Here the rea numbers a n, b n are caed the Fourier

More information

FOURIER SERIES ON ANY INTERVAL

FOURIER SERIES ON ANY INTERVAL FOURIER SERIES ON ANY INTERVAL Overview We have spent considerabe time earning how to compute Fourier series for functions that have a period of 2p on the interva (-p,p). We have aso seen how Fourier series

More information

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t. Strauss PDEs e: Section 5.6 - Exercise Page 1 of 1 Exercise For probem (1, compete the cacuation of the series in case j(t = and h(t = e t. Soution With j(t = and h(t = e t, probem (1 on page 147 becomes

More information

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

Physics 235 Chapter 8. Chapter 8 Central-Force Motion Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion

More information

Lecture 17 - The Secrets we have Swept Under the Rug

Lecture 17 - The Secrets we have Swept Under the Rug Lecture 17 - The Secrets we have Swept Under the Rug Today s ectures examines some of the uirky features of eectrostatics that we have negected up unti this point A Puzze... Let s go back to the basics

More information

HYDROGEN ATOM SELECTION RULES TRANSITION RATES

HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOING PHYSICS WITH MATLAB QUANTUM PHYSICS Ian Cooper Schoo of Physics, University of Sydney ian.cooper@sydney.edu.au HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS

More information

MA 201: Partial Differential Equations Lecture - 10

MA 201: Partial Differential Equations Lecture - 10 MA 201: Partia Differentia Equations Lecture - 10 Separation of Variabes, One dimensiona Wave Equation Initia Boundary Vaue Probem (IBVP) Reca: A physica probem governed by a PDE may contain both boundary

More information

VI.G Exact free energy of the Square Lattice Ising model

VI.G Exact free energy of the Square Lattice Ising model VI.G Exact free energy of the Square Lattice Ising mode As indicated in eq.(vi.35), the Ising partition function is reated to a sum S, over coections of paths on the attice. The aowed graphs for a square

More information

V.B The Cluster Expansion

V.B The Cluster Expansion V.B The Custer Expansion For short range interactions, speciay with a hard core, it is much better to repace the expansion parameter V( q ) by f( q ) = exp ( βv( q )), which is obtained by summing over

More information

Assignment 7 Due Tuessday, March 29, 2016

Assignment 7 Due Tuessday, March 29, 2016 Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges

More information

Week 6 Lectures, Math 6451, Tanveer

Week 6 Lectures, Math 6451, Tanveer Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n

More information

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18 Quantum Mechanica Modes of Vibration and Rotation of Moecues Chapter 18 Moecuar Energy Transationa Vibrationa Rotationa Eectronic Moecuar Motions Vibrations of Moecues: Mode approximates moecues to atoms

More information

Midterm 2 Review. Drew Rollins

Midterm 2 Review. Drew Rollins Midterm 2 Review Drew Roins 1 Centra Potentias and Spherica Coordinates 1.1 separation of variabes Soving centra force probems in physics (physica systems described by two objects with a force between

More information

Problem set 6 The Perron Frobenius theorem.

Problem set 6 The Perron Frobenius theorem. Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator

More information

An approximate method for solving the inverse scattering problem with fixed-energy data

An approximate method for solving the inverse scattering problem with fixed-energy data J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999

More information

XSAT of linear CNF formulas

XSAT of linear CNF formulas XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

(Refer Slide Time: 2:34) L C V

(Refer Slide Time: 2:34) L C V Microwave Integrated Circuits Professor Jayanta Mukherjee Department of Eectrica Engineering Indian Intitute of Technoogy Bombay Modue 1 Lecture No 2 Refection Coefficient, SWR, Smith Chart. Heo wecome

More information

A Brief Introduction to Markov Chains and Hidden Markov Models

A Brief Introduction to Markov Chains and Hidden Markov Models A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,

More information

6 Wave Equation on an Interval: Separation of Variables

6 Wave Equation on an Interval: Separation of Variables 6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.

More information

Some Measures for Asymmetry of Distributions

Some Measures for Asymmetry of Distributions Some Measures for Asymmetry of Distributions Georgi N. Boshnakov First version: 31 January 2006 Research Report No. 5, 2006, Probabiity and Statistics Group Schoo of Mathematics, The University of Manchester

More information

Lecture 6: Moderately Large Deflection Theory of Beams

Lecture 6: Moderately Large Deflection Theory of Beams Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey

More information

Math 220B - Summer 2003 Homework 1 Solutions

Math 220B - Summer 2003 Homework 1 Solutions Math 0B - Summer 003 Homework Soutions Consider the eigenvaue probem { X = λx 0 < x < X satisfies symmetric BCs x = 0, Suppose f(x)f (x) x=b x=a 0 for a rea-vaued functions f(x) which satisfy the boundary

More information

Joel Broida UCSD Fall 2009 Phys 130B QM II. Homework Set 2 DO ALL WORK BY HAND IN OTHER WORDS, DON T USE MATHEMAT- ICA OR ANY CALCULATORS.

Joel Broida UCSD Fall 2009 Phys 130B QM II. Homework Set 2 DO ALL WORK BY HAND IN OTHER WORDS, DON T USE MATHEMAT- ICA OR ANY CALCULATORS. Joe Broida UCSD Fa 009 Phys 30B QM II Homework Set DO ALL WORK BY HAND IN OTHER WORDS, DON T USE MATHEMAT- ICA OR ANY CALCULATORS. You may need to use one or more of these: Y 0 0 = 4π Y 0 = 3 4π cos Y

More information

Introduction to LMTO method

Introduction to LMTO method 1 Introduction to MTO method 24 February 2011; V172 P.Ravindran, FME-course on Ab initio Modeing of soar ce Materias 24 February 2011 Introduction to MTO method Ab initio Eectronic Structure Cacuations

More information

Lecture Note 3: Stationary Iterative Methods

Lecture Note 3: Stationary Iterative Methods MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or

More information

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p

More information

1D Heat Propagation Problems

1D Heat Propagation Problems Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2

More information

TRANSFORMATION OF REAL SPHERICAL HARMONICS UNDER ROTATIONS

TRANSFORMATION OF REAL SPHERICAL HARMONICS UNDER ROTATIONS Vo. 39 (008) ACTA PHYSICA POLONICA B No 8 TRANSFORMATION OF REAL SPHERICAL HARMONICS UNDER ROTATIONS Zbigniew Romanowski Interdiscipinary Centre for Materias Modeing Pawińskiego 5a, 0-106 Warsaw, Poand

More information

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION J. Korean Math. Soc. 46 2009, No. 2, pp. 281 294 ORHOGONAL MLI-WAVELES FROM MARIX FACORIZAION Hongying Xiao Abstract. Accuracy of the scaing function is very crucia in waveet theory, or correspondingy,

More information

ANISOTROPIES OF THE MICROWAVE BACKGROUND

ANISOTROPIES OF THE MICROWAVE BACKGROUND ANISOTROPIES OF THE MICROWAVE BACKGROUND The Universe just before recombination is a very tighty couped fuid, due to the arge eectromagnetic Thomson cross section. Photons scatter off charged partices

More information

Cryptanalysis of PKP: A New Approach

Cryptanalysis of PKP: A New Approach Cryptanaysis of PKP: A New Approach Éiane Jaumes and Antoine Joux DCSSI 18, rue du Dr. Zamenhoff F-92131 Issy-es-Mx Cedex France eiane.jaumes@wanadoo.fr Antoine.Joux@ens.fr Abstract. Quite recenty, in

More information

Lecture 6 Povh Krane Enge Williams Properties of 2-nucleon potential

Lecture 6 Povh Krane Enge Williams Properties of 2-nucleon potential Lecture 6 Povh Krane Enge Wiiams Properties of -nuceon potentia 16.1 4.4 3.6 9.9 Meson Theory of Nucear potentia 4.5 3.11 9.10 I recommend Eisberg and Resnik notes as distributed Probems, Lecture 6 1 Consider

More information

Solution Set Seven. 1 Goldstein Components of Torque Along Principal Axes Components of Torque Along Cartesian Axes...

Solution Set Seven. 1 Goldstein Components of Torque Along Principal Axes Components of Torque Along Cartesian Axes... : Soution Set Seven Northwestern University, Cassica Mechanics Cassica Mechanics, Third Ed.- Godstein November 8, 25 Contents Godstein 5.8. 2. Components of Torque Aong Principa Axes.......................

More information

Bohr s atomic model. 1 Ze 2 = mv2. n 2 Z

Bohr s atomic model. 1 Ze 2 = mv2. n 2 Z Bohr s atomic mode Another interesting success of the so-caed od quantum theory is expaining atomic spectra of hydrogen and hydrogen-ike atoms. The eectromagnetic radiation emitted by free atoms is concentrated

More information

arxiv:nlin/ v2 [nlin.cd] 30 Jan 2006

arxiv:nlin/ v2 [nlin.cd] 30 Jan 2006 expansions in semicassica theories for systems with smooth potentias and discrete symmetries Hoger Cartarius, Jörg Main, and Günter Wunner arxiv:nin/0510051v [nin.cd] 30 Jan 006 1. Institut für Theoretische

More information

Srednicki Chapter 51

Srednicki Chapter 51 Srednici Chapter 51 QFT Probems & Soutions A. George September 7, 13 Srednici 51.1. Derive the fermion-oop correction to the scaar proagator by woring through equation 5., and show that it has an extra

More information

CONGRUENCES. 1. History

CONGRUENCES. 1. History CONGRUENCES HAO BILLY LEE Abstract. These are notes I created for a seminar tak, foowing the papers of On the -adic Representations and Congruences for Coefficients of Moduar Forms by Swinnerton-Dyer and

More information

Vector Spherical Harmonics

Vector Spherical Harmonics Vector Spherica Haronics Lecture Introduction Previousy we have seen that the Lapacian operator is different when operating on a vector or a scaar function. We avoided this probe by etting the Lapacian

More information

Homework #04 Answers and Hints (MATH4052 Partial Differential Equations)

Homework #04 Answers and Hints (MATH4052 Partial Differential Equations) Homework #4 Answers and Hints (MATH452 Partia Differentia Equations) Probem 1 (Page 89, Q2) Consider a meta rod ( < x < ), insuated aong its sides but not at its ends, which is initiay at temperature =

More information

This shoud be true for any, so the equations are equivaent to the Schrodingerequation (H E) : () Another important property of the (2) functiona is, t

This shoud be true for any, so the equations are equivaent to the Schrodingerequation (H E) : () Another important property of the (2) functiona is, t I. Variationa method Introduction Usefu for the determination of energies and wavefunctions of different atomic states. H the time-independent Hamitonian of the system E n eigenenergies Ψ n - eigenfunctions

More information

Nuclear Size and Density

Nuclear Size and Density Nucear Size and Density How does the imited range of the nucear force affect the size and density of the nucei? Assume a Vecro ba mode, each having radius r, voume V = 4/3π r 3. Then the voume of the entire

More information

Week 5 Lectures, Math 6451, Tanveer

Week 5 Lectures, Math 6451, Tanveer Week 5 Lectures, Math 651, Tanveer 1 Separation of variabe method The method of separation of variabe is a suitabe technique for determining soutions to inear PDEs, usuay with constant coefficients, when

More information

<C 2 2. λ 2 l. λ 1 l 1 < C 1

<C 2 2. λ 2 l. λ 1 l 1 < C 1 Teecommunication Network Contro and Management (EE E694) Prof. A. A. Lazar Notes for the ecture of 7/Feb/95 by Huayan Wang (this document was ast LaT E X-ed on May 9,995) Queueing Primer for Muticass Optima

More information

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After

More information

In Coulomb gauge, the vector potential is then given by

In Coulomb gauge, the vector potential is then given by Physics 505 Fa 007 Homework Assignment #8 Soutions Textbook probems: Ch. 5: 5.13, 5.14, 5.15, 5.16 5.13 A sphere of raius a carries a uniform surface-charge istribution σ. The sphere is rotate about a

More information

Identites and properties for associated Legendre functions

Identites and properties for associated Legendre functions Identites and properties for associated Legendre functions DBW This note is a persona note with a persona history; it arose out off y incapacity to find references on the internet that prove reations that

More information

Electromagnetism Spring 2018, NYU

Electromagnetism Spring 2018, NYU Eectromagnetism Spring 08, NYU March 6, 08 Time-dependent fieds We now consider the two phenomena missing from the static fied case: Faraday s Law of induction and Maxwe s dispacement current. Faraday

More information

MAT 167: Advanced Linear Algebra

MAT 167: Advanced Linear Algebra < Proem 1 (15 pts) MAT 167: Advanced Linear Agera Fina Exam Soutions (a) (5 pts) State the definition of a unitary matrix and expain the difference etween an orthogona matrix and an unitary matrix. Soution:

More information

Partial permutation decoding for MacDonald codes

Partial permutation decoding for MacDonald codes Partia permutation decoding for MacDonad codes J.D. Key Department of Mathematics and Appied Mathematics University of the Western Cape 7535 Bevie, South Africa P. Seneviratne Department of Mathematics

More information

FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 2017

FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 2017 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II November 5, 207 Prof. Alan Guth FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 207 A few items below are marked

More information

8 Digifl'.11 Cth:uits and devices

8 Digifl'.11 Cth:uits and devices 8 Digif'. Cth:uits and devices 8. Introduction In anaog eectronics, votage is a continuous variabe. This is usefu because most physica quantities we encounter are continuous: sound eves, ight intensity,

More information

11.1 One-dimensional Helmholtz Equation

11.1 One-dimensional Helmholtz Equation Chapter Green s Functions. One-dimensiona Hemhotz Equation Suppose we have a string driven by an externa force, periodic with frequency ω. The differentia equation here f is some prescribed function) 2

More information

Theory and implementation behind: Universal surface creation - smallest unitcell

Theory and implementation behind: Universal surface creation - smallest unitcell Teory and impementation beind: Universa surface creation - smaest unitce Bjare Brin Buus, Jaob Howat & Tomas Bigaard September 15, 218 1 Construction of surface sabs Te aim for tis part of te project is

More information

Supporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers

Supporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers Supporting Information for Suppressing Kein tunneing in graphene using a one-dimensiona array of ocaized scatterers Jamie D Was, and Danie Hadad Department of Chemistry, University of Miami, Cora Gabes,

More information

The Group Structure on a Smooth Tropical Cubic

The Group Structure on a Smooth Tropical Cubic The Group Structure on a Smooth Tropica Cubic Ethan Lake Apri 20, 2015 Abstract Just as in in cassica agebraic geometry, it is possibe to define a group aw on a smooth tropica cubic curve. In this note,

More information

(f) is called a nearly holomorphic modular form of weight k + 2r as in [5].

(f) is called a nearly holomorphic modular form of weight k + 2r as in [5]. PRODUCTS OF NEARLY HOLOMORPHIC EIGENFORMS JEFFREY BEYERL, KEVIN JAMES, CATHERINE TRENTACOSTE, AND HUI XUE Abstract. We prove that the product of two neary hoomorphic Hece eigenforms is again a Hece eigenform

More information

A SIMPLIFIED DESIGN OF MULTIDIMENSIONAL TRANSFER FUNCTION MODELS

A SIMPLIFIED DESIGN OF MULTIDIMENSIONAL TRANSFER FUNCTION MODELS A SIPLIFIED DESIGN OF ULTIDIENSIONAL TRANSFER FUNCTION ODELS Stefan Petrausch, Rudof Rabenstein utimedia Communications and Signa Procesg, University of Erangen-Nuremberg, Cauerstr. 7, 958 Erangen, GERANY

More information

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z Chapter W3 Mechanica Systems II Introduction This companion website chapter anayzes the foowing topics in connection to the printed-book Chapter 3: Lumped-parameter inertia fractions of basic compiant

More information

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA ON THE SYMMETRY OF THE POWER INE CHANNE T.C. Banwe, S. Gai {bct, sgai}@research.tecordia.com Tecordia Technoogies, Inc., 445 South Street, Morristown, NJ 07960, USA Abstract The indoor power ine network

More information

Cluster modelling. Collisions. Stellar Dynamics & Structure of Galaxies handout #2. Just as a self-gravitating collection of objects.

Cluster modelling. Collisions. Stellar Dynamics & Structure of Galaxies handout #2. Just as a self-gravitating collection of objects. Stear Dynamics & Structure of Gaaxies handout # Custer modeing Just as a sef-gravitating coection of objects. Coisions Do we have to worry about coisions? Gobuar custers ook densest, so obtain a rough

More information

Estimating the Power Spectrum of the Cosmic Microwave Background

Estimating the Power Spectrum of the Cosmic Microwave Background Estimating the Power Spectrum of the Cosmic Microwave Background J. R. Bond 1,A.H.Jaffe 2,andL.Knox 1 1 Canadian Institute for Theoretica Astrophysics, Toronto, O M5S 3H8, CAADA 2 Center for Partice Astrophysics,

More information

c 2007 Society for Industrial and Applied Mathematics

c 2007 Society for Industrial and Applied Mathematics SIAM REVIEW Vo. 49,No. 1,pp. 111 1 c 7 Society for Industria and Appied Mathematics Domino Waves C. J. Efthimiou M. D. Johnson Abstract. Motivated by a proposa of Daykin [Probem 71-19*, SIAM Rev., 13 (1971),

More information

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 2: The Moment Equations

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 2: The Moment Equations .615, MHD Theory of Fusion ystems Prof. Freidberg Lecture : The Moment Equations Botzmann-Maxwe Equations 1. Reca that the genera couped Botzmann-Maxwe equations can be written as f q + v + E + v B f =

More information

1.2 Partial Wave Analysis

1.2 Partial Wave Analysis February, 205 Lecture X.2 Partia Wave Anaysis We have described scattering in terms of an incoming pane wave, a momentum eigenet, and and outgoing spherica wave, aso with definite momentum. We now consider

More information

17 Lecture 17: Recombination and Dark Matter Production

17 Lecture 17: Recombination and Dark Matter Production PYS 652: Astrophysics 88 17 Lecture 17: Recombination and Dark Matter Production New ideas pass through three periods: It can t be done. It probaby can be done, but it s not worth doing. I knew it was

More information