Designer basis functions for potentials in galactic dynamics Prasenjit Saha Canadian Institute for Theoretical Astrophysics McLennan Labs, University
|
|
- Ralf Barrett
- 5 years ago
- Views:
Transcription
1 Designer basis functions for potentials in galactic dynamics Prasenjit Saha Canadian Institute for Theoretical Astrophysics McLennan Labs, University of Toronto, 60 St. George St., Toronto M5S 1A7, Canada Summary: Existing methods for solving Poisson's equation using basis functions can be made considerably more exible, yet simpler, by using non-orthogonal basis functions. The rst part of this paper describes how such basis sets can be constructed. The second part gives explicit formulas (up to quadratures and matrix inversions) for basis functions being spherical harmonics times radial functions; two basis sets in the literature appear as special cases. To appear in Mon. Not. R. astr. Soc. 1
2 1 THE GENERAL PROBLEM Suppose we are doing a problem that requires approximate solutions of Poisson's equation. This may be a particle simulation, a gas simulation, or a normal mode analysis of a stellar dynamical system. Suppose also that we have a crude idea of the behaviour of the potential (such as asymptotic small-r and large-r forms) and further that is very smooth in particular we may want to smooth over the graininess in a simulation. Then it would be sensible to express the potential as ha `prior' form W i hpolynomial in some variable(s)i: (1:1) or as a sum of such terms. The prior form W is thus a zeroth order approximation to. This paper suggests how we can solve eciently for such forms of. 1.1 Known bases For some special W, there are elegant basis sets in the literature which allow one to solve Poisson's equation r 2 = 4 (1:2) in a simple and ecient way. These basis sets n are biorthonormal, i.e., m r 2 n d 3 x = mn : (1:3) Assuming the m are complete, we have = X n c n n ; c n = 1 4 n d 3 x: (1:4) Such basis sets are typically eigenfunctions of hsome weight functioni r 2 ; biorthogonality and completeness follow since r 2 is Hermitian. Since spherical harmonics are eigenfunctions of the angular part of r 2, one obvious form is lmp = Y lm (; ) F lp (r): (1:5) The basis sets given by Clutton-Brock (1973; hereafter CB73), Polyachenko & Shukhman (1981), and Hernquist & Ostriker (1992; hereafter HO92) are of the type (1.5), with F lp being eigenfunctions of the radial part of r 2 (apart from some weight function). For analogous examples in disc systems see Clutton-Brock (1972) and Kalnajs (1976). These bases are useful, but they are inexible in their choice of the weight function (or prior form) W ; they are also highly non-trivial to discover, so a worker wishing to tailor a basis set to a particular application has little hope of nding a new eigenfunction set with a well-matched W. Generality can be gained by dropping the requirement that the n be eigenfunctions of r 2 and simply constructing a biorthonormal set by a Gramm-Schmidt method. In Saha (1991) I adopted a basis of the form (1.5), with F lp not being eigenfunctions of anything in particular, but made biorthogonal by Gramm-Schmidt. Robijn (1992) uses a similar method, but with spheroidal harmonics rather than spherical harmonics for the angular part. 2
3 1.2 Non-biorthogonal bases In fact the biorthogonality requirement can be dropped as well. If we take any linearly independent and complete set n, we can expand = X n c n n (1:6a) where the expansion coecients c n satisfy X n M mn c n = 4 m d 3 x; where M mn = m r2 n d 3 x (1:6b) and M mn must be non-singular. Evaluating and inverting M mn needs to done only once, which will typically be a trivial overhead compared to computing R n d 3 x. The proposal of this paper is to take n of the form (1.1) and then use (1.6) to nd the expansion coecients. The polynomial factor guarantees linear independence of the n. Completeness in the strict sense does not follow, (but as illustrated later in this paper) a weaker form of completeness will hold that seems adequate for applications. 1.3 Discussion The second part of this paper is devoted to the details of the spherical-harmonic-timesradial-functions form (1.5). It turns out that the bases of CB73 and HO92 strictly speaking, linearly transformed versions of them appear as special cases. In two dimensions the problem is more dicult and are not related by a two-dimensional operator and it is not clear whether the scheme of this paper could be adapted to thin discs. For discs of nite thickness, however, bases using cylindrical coordinates seem feasible and appropriate. It is also worth mentioning that the technique suggested here is not specic to solving Poisson's equation. In particular, we could expand the momentum density p(x) from a sample of points in terms of some basis functions f n (x), thus: p = X n a n f n (1:7a) where the a n satisfy X n M mn a n = vf m d 3 x; where now M mn = f mf n d 3 x: (1:7b) Now if we combine such an expansion with a Poisson solver (1.6), the ratio 4 P m a nf n Pn c nr 2 n (1:8) gives a smoothed velocity eld. 3
4 2 HARMONIC EXPANSIONS 2.1 Formulas We assume the potential may be expanded as Then the c lmp are given by where S lm p = X lmp c lmp Y lm (; )F lp (r): (2:1) 4 pq 1 0 X q pq c lmq = S lm p (2:2) Y lm(; )F lp r 2 d 3 x (2:3a) F lp d dr r2 d dr? l(l + 1) F lq dr: (2:3b) We now do two things. First we assume a separate prior form W l for each multipole component, and factor this out, thus: pq =?1 F lp (r) = W l (r) U lp (r): (2:4) Then we change variables from r to some u having a nite domain, (say?1 to 1). This changes (2.3b) to 1 ( r 2 W l U lp r 0 W l d 2 + 2r 2 W 0 2 l + r 2? rr00 d W r 02 l + r 0h r 2 W 00 l + 2rW 0 l? l(l + 1)W li ) U lq (2:5a) where r 0 dr ; r00 d2 r 2 ; W 0 dw l l dr ; W 00 l d2 W l dr : 2 (2:5b) The cartesian components of the acceleration can be evaluated from the derivatives of thus: 10 1 x xz y r rr 2? CB C R 2 r? 0 B y r z r yz rr 2? 1 r x R 2 0 C A B sin C A (2:6) where R means p x 2 + y 2. This of course will require values of the Legendre polynomials and their derivatives, which are best got from their recursion relations. 4
5 2.2 Special cases CB73 takes W l = r l (1 + r 2 ) l+ 1 2?1 and u = r2? 1 r : (2:7) Inserting (2.7) into the form (2.5) and laboriously collecting terms yields: 1 Mpq lm = 2?(2l+1) (1? u 2 ) l+ 2 1 Ulp (1? u 2 ) d2 d? (2l + 3)u 2? 1 4 (2l + 1)(2l + 3) U lq : (2:8) HO92 take which gives: W l = pq = 2?(4l+3) 1?1 rl and u = r? 1 (1 + r) 2l+1 r + 1 : (2:9) (1? u 2 ) 2l+1 U lp (1? u 2 ) d2? 2(l + 1)(2l + 1) U lq : d? 4(l + 1)u 2 (2:10) It turns out that in both cases the eigenfunctions of the operator within the curly brackets in (2.10) and (2.8) are known functions ultraspherical polynomials and the integrals express the orthogonality relation for these eigenfunctions. Thus, if U lp is taken to be the appropriate ultraspherical polynomial then the matrix Mpq lm is diagonal. Now, for our purpose we don't need Mpq lm to be diagonal. Therefore, instead of taking U lp to be an p-th degree ultraspherical polynomial, we may take U lp = u p. This is a linearly independent basis related to the ultraspherical polynomials by a linear transformation, and will therefore give the same results for. What other choices of W l would be apt for studying (say) elliptical galaxies? One possibility would be W l = r l (1 + r)?2l ; =?1 r (1 + r)?? 1 ; (2:11a) with chosen depending on the problem at hand. The density corresponding to is = (4)?1 (1 + )r?2 (1 + r)??2 ; (2:11b) which is adjustable from a point mass ( =?1), through a Jae model ( = 0), to a Hernquist model ( = 1), and on to models with nite central density. 5
6 2.3 Completeness questions If we are trying to approximate some potential multipole moment G l (r) using our basis functions, can we make the error jg l j in our approximation arbitrarily small by including enough basis functions? Completeness in this sense does not hold, but the following weaker form of completeness does apply. Any continuous function f(u) with a nite domain can be approximated by a polynomial Q(u) in the sense that jf(u)? Q(u)j 2 (2:12) can be made arbitrarily small by making the degree of Q high enough. Since our basis consists of W l times polynomials in u, we can make jgl j 2 W 2 l dr (2:13) dr arbitrarily small. This type of argument could, of course, be applied in other coordinate systems as well. 2.4 Implementation I have written some C functions implementing the formulas of this section (and also some extra code necessary to test these). The chosen prior form is specied through user-supplied functions for r(u), u(r), dr=, W l (r), dw l =dr, and the integral in (2.5a). The code will then expand to the desired harmonic and polynomial order from a sample of masses and positions, and supply the potential and accelerations. There is also some code to t r 2 v to basis functions again with some chosen prior form and thus nd a smoothed velocity eld. (Since may become singular in some model systems, it seems better to expand r 2 v rather than v in basis functions.) Programs and documentation are available to anyone interested. ACKNOWLEDGMENTS I am grateful to James Binney, M. Slim Fayache, Lars Hernquist, Kathryn Johnston, and Scott Tremaine for discussions and correspondence. I also thank the Natural Sciences and Engineering Research Council of Canada for nancial support through an operating grant, and the Physics Department, Oxford University for their hospitality. REFERENCES Clutton-Brock, M., Astrophys. & Space Sci., 16, 101. Clutton-Brock, M., Astrophys. & Space Sci., 23, 55. Hernquist, L. & J. P. Ostriker, Astrophys. J., 386, 375. Kalnajs, A. J., Astrophys. J., 205, 745. Polyachenko, V. L. & I. G. Shukhman, Astron. h., 58, 933, English translation Sov. Astr., 25, 533. Robijn, F., Leiden preprint. Saha, P., Mon. Not. R. astr. Soc., 248,
High-accuracy minimum relaxation N-body simulations using orthogonal series force computation
University of Massachusetts Amherst ScholarWorks@UMass Amherst Astronomy Department Faculty Publication Series Astronomy 1996 High-accuracy minimum relaxation N-body simulations using orthogonal series
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
More informationAngular momentum. Quantum mechanics. Orbital angular momentum
Angular momentum 1 Orbital angular momentum Consider a particle described by the Cartesian coordinates (x, y, z r and their conjugate momenta (p x, p y, p z p. The classical definition of the orbital angular
More informationGalactic Astronomy 2016
10 Potential theory The study of dynamics is fundamental to any understanding of galaxies. The following few chapters provide a brief summary of the most important results. More details can be found in
More informationSpherical Coordinates and Legendre Functions
Spherical Coordinates and Legendre Functions Spherical coordinates Let s adopt the notation for spherical coordinates that is standard in physics: φ = longitude or azimuth, θ = colatitude ( π 2 latitude)
More information1 r 2 sin 2 θ. This must be the case as we can see by the following argument + L2
PHYS 4 3. The momentum operator in three dimensions is p = i Therefore the momentum-squared operator is [ p 2 = 2 2 = 2 r 2 ) + r 2 r r r 2 sin θ We notice that this can be written as sin θ ) + θ θ r 2
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 21 Quantum Mechanics in Three Dimensions Lecture 21 Physics 342 Quantum Mechanics I Monday, March 22nd, 21 We are used to the temporal separation that gives, for example, the timeindependent
More informationLinear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space
Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) Contents 1 Vector Spaces 1 1.1 The Formal Denition of a Vector Space.................................. 1 1.2 Subspaces...................................................
More informationSUMMARY PHYSICS 707 Electrostatics. E(x) = 4πρ(x) and E(x) = 0 (1)
SUMMARY PHYSICS 707 Electrostatics The basic differential equations of electrostatics are E(x) = 4πρ(x) and E(x) = 0 (1) where E(x) is the electric field and ρ(x) is the electric charge density. The field
More informationwhich implies that we can take solutions which are simultaneous eigen functions of
Module 1 : Quantum Mechanics Chapter 6 : Quantum mechanics in 3-D Quantum mechanics in 3-D For most physical systems, the dynamics is in 3-D. The solutions to the general 3-d problem are quite complicated,
More informationLegendre Polynomials and Angular Momentum
University of Connecticut DigitalCommons@UConn Chemistry Education Materials Department of Chemistry August 006 Legendre Polynomials and Angular Momentum Carl W. David University of Connecticut, Carl.David@uconn.edu
More information5.1 Circular Velocities and Rotation Curves
Chapter 5 otation Curves 5.1 Circular Velocities and otation Curves The circular velocity v circ is the velocity that a star in a galaxy must have to maintain a circular orbit at a specified distance from
More informationarxiv: v1 [gr-qc] 26 Jun 2008
Chaotic Motion Around Stellar Objects with Octupolar Deformation: Newtonian and Post Newtonian Approaches Javier Ramos-Caro, Framsol López-Suspes, and Guillermo A. González Grupo de Investigación en Relatividad
More informationwhich arises when we compute the orthogonal projection of a vector y in a subspace with an orthogonal basis. Hence assume that P y = A ij = x j, x i
MODULE 6 Topics: Gram-Schmidt orthogonalization process We begin by observing that if the vectors {x j } N are mutually orthogonal in an inner product space V then they are necessarily linearly independent.
More information2 where (u v w) = u 0 denote the deviations from the mean ow in a Cartesian coordinate system (x y z). We adopt a right-handed coordinate system (e.g.
SOLAR MESOGRANULATION CORRELATION OBSERVATIONS PROBING THE HELICITY G. R UDIGER Astrophysikalisches Institut Potsdam, An der Sternwarte 16, D-14482 Potsdam, Germany A. BRANDENBURG Department of Mathematics,
More informationElectrodynamics I Midterm - Part A - Closed Book KSU 2005/10/17 Electro Dynamic
Electrodynamics I Midterm - Part A - Closed Book KSU 5//7 Name Electro Dynamic. () Write Gauss Law in differential form. E( r) =ρ( r)/ɛ, or D = ρ, E= electricfield,ρ=volume charge density, ɛ =permittivity
More informationComputation of the scattering amplitude in the spheroidal coordinates
Computation of the scattering amplitude in the spheroidal coordinates Takuya MINE Kyoto Institute of Technology 12 October 2015 Lab Seminar at Kochi University of Technology Takuya MINE (KIT) Spheroidal
More informationMultipole moments. November 9, 2015
Multipole moments November 9, 5 The far field expansion Suppose we have a localized charge distribution, confined to a region near the origin with r < R. Then for values of r > R, the electric field must
More informationarxiv: v1 [gr-qc] 2 Feb 2015
arxiv:1502.00424v1 [gr-qc] 2 Feb 2015 Valiente Kroon s obstructions to smoothness at infinity James Grant Department of Mathematics, University of Surrey, Paul Tod Mathematical Institute University of
More informationPhysics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I
Physics 342 Lecture 23 Radial Separation Lecture 23 Physics 342 Quantum Mechanics I Friday, March 26th, 2010 We begin our spherical solutions with the simplest possible case zero potential. Aside from
More information8.1 The hydrogen atom solutions
8.1 The hydrogen atom solutions Slides: Video 8.1.1 Separating for the radial equation Text reference: Quantum Mechanics for Scientists and Engineers Section 10.4 (up to Solution of the hydrogen radial
More informationN-BODY SIMULATIONS OF DISC GALAXIES CAN SHED LIGHT ON THE DARK-MATTER PROBLEM ALESSANDRO B. ROMEO Onsala Space Observatory, Chalmers University of Tec
N-BODY SIMULATIONS OF DISC GALAXIES CAN SHED LIGHT ON THE DARK-MATTER PROBLEM ALESSANDRO B. ROMEO Onsala Space Observatory, Chalmers University of Technology, S-43992 Onsala, Sweden E-Mail: romeo@oso.chalmers.se
More informationContents. 2.1 Vectors in R n. Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v. 2.50) 2 Vector Spaces
Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v 250) Contents 2 Vector Spaces 1 21 Vectors in R n 1 22 The Formal Denition of a Vector Space 4 23 Subspaces 6 24 Linear Combinations and
More information4.1 Eigenvalues, Eigenvectors, and The Characteristic Polynomial
Linear Algebra (part 4): Eigenvalues, Diagonalization, and the Jordan Form (by Evan Dummit, 27, v ) Contents 4 Eigenvalues, Diagonalization, and the Jordan Canonical Form 4 Eigenvalues, Eigenvectors, and
More information20 The Hydrogen Atom. Ze2 r R (20.1) H( r, R) = h2 2m 2 r h2 2M 2 R
20 The Hydrogen Atom 1. We want to solve the time independent Schrödinger Equation for the hydrogen atom. 2. There are two particles in the system, an electron and a nucleus, and so we can write the Hamiltonian
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More informationarxiv: v1 [astro-ph] 6 Sep 2007
Mon. Not. R. Astron. Soc. 000, 1 22 (2007) Printed 5 February 2008 (MN LATEX style file v2.2) Appropriate SCF basis sets for orbital studies of galaxies and a quantum-mechanical method to compute them
More informationAngular Momentum. Classically the orbital angular momentum with respect to a fixed origin is. L = r p. = yp z. L x. zp y L y. = zp x. xpz L z.
Angular momentum is an important concept in quantum theory, necessary for analyzing motion in 3D as well as intrinsic properties such as spin Classically the orbital angular momentum with respect to a
More informationSOLUTIONS TO SELECTED PROBLEMS FROM ASSIGNMENTS 3, 4
SOLUTIONS TO SELECTED POBLEMS FOM ASSIGNMENTS 3, 4 Problem 5 from Assignment 3 Statement. Let be an n-dimensional bounded domain with smooth boundary. Show that the eigenvalues of the Laplacian on with
More informationCHAPTER 3 POTENTIALS 10/13/2016. Outlines. 1. Laplace s equation. 2. The Method of Images. 3. Separation of Variables. 4. Multipole Expansion
CHAPTER 3 POTENTIALS Lee Chow Department of Physics University of Central Florida Orlando, FL 32816 Outlines 1. Laplace s equation 2. The Method of Images 3. Separation of Variables 4. Multipole Expansion
More information1 Vectors. Notes for Bindel, Spring 2017 Numerical Analysis (CS 4220)
Notes for 2017-01-30 Most of mathematics is best learned by doing. Linear algebra is no exception. You have had a previous class in which you learned the basics of linear algebra, and you will have plenty
More information1 A complete Fourier series solution
Math 128 Notes 13 In this last set of notes I will try to tie up some loose ends. 1 A complete Fourier series solution First here is an example of the full solution of a pde by Fourier series. Consider
More informationChemistry 532 Practice Final Exam Fall 2012 Solutions
Chemistry 53 Practice Final Exam Fall Solutions x e ax dx π a 3/ ; π sin 3 xdx 4 3 π cos nx dx π; sin θ cos θ + K x n e ax dx n! a n+ ; r r r r ˆL h r ˆL z h i φ ˆL x i hsin φ + cot θ cos φ θ φ ) ˆLy i
More informationCurvilinear coordinates
C Curvilinear coordinates The distance between two points Euclidean space takes the simplest form (2-4) in Cartesian coordinates. The geometry of concrete physical problems may make non-cartesian coordinates
More informationSeparation of Variables in Polar and Spherical Coordinates
Separation of Variables in Polar and Spherical Coordinates Polar Coordinates Suppose we are given the potential on the inside surface of an infinitely long cylindrical cavity, and we want to find the potential
More informationElectromagnetic Theory I
Electromagnetic Theory I Final Examination 18 December 2009, 12:30-2:30 pm Instructions: Answer the following 10 questions, each of which is worth 10 points. Explain your reasoning in each case. Use SI
More informationRelevant sections from AMATH 351 Course Notes (Wainwright): Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls):
Lecture 5 Series solutions to DEs Relevant sections from AMATH 35 Course Notes (Wainwright):.4. Relevant sections from AMATH 35 Course Notes (Poulin and Ingalls): 2.-2.3 As mentioned earlier in this course,
More informationTime part of the equation can be separated by substituting independent equation
Lecture 9 Schrödinger Equation in 3D and Angular Momentum Operator In this section we will construct 3D Schrödinger equation and we give some simple examples. In this course we will consider problems where
More informationPractical in Numerical Astronomy, SS 2012 LECTURE 9
Practical in Numerical Astronomy, SS 01 Elliptic partial differential equations. Poisson solvers. LECTURE 9 1. Gravity force and the equations of hydrodynamics. Poisson equation versus Poisson integral.
More informationLecture #1. Review. Postulates of quantum mechanics (1-3) Postulate 1
L1.P1 Lecture #1 Review Postulates of quantum mechanics (1-3) Postulate 1 The state of a system at any instant of time may be represented by a wave function which is continuous and differentiable. Specifically,
More informationJackson 6.4 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jackson 6.4 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell PROBLEM: A uniformly magnetized and conducting sphere of radius R and total magnetic moment m = 4πMR 3
More informationEigenfunctions on the surface of a sphere. In spherical coordinates, the Laplacian is. u = u rr + 2 r u r + 1 r 2. sin 2 (θ) + 1
Eigenfunctions on the surface of a sphere In spherical coordinates, the Laplacian is u = u rr + 2 r u r + 1 r 2 [ uφφ sin 2 (θ) + 1 sin θ (sin θ u θ) θ ]. Eigenfunctions on the surface of a sphere In spherical
More informationSolutions to Laplace s Equations- II
Solutions to Laplace s Equations- II Lecture 15: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Laplace s Equation in Spherical Coordinates : In spherical coordinates
More informationThe two body problem involves a pair of particles with masses m 1 and m 2 described by a Lagrangian of the form:
Physics 3550, Fall 2011 Two Body, Central-Force Problem Relevant Sections in Text: 8.1 8.7 Two Body, Central-Force Problem Introduction. I have already mentioned the two body central force problem several
More informationQuantum Orbits. Quantum Theory for the Computer Age Unit 9. Diving orbit. Caustic. for KE/PE =R=-3/8. for KE/PE =R=-3/8. p"...
W.G. Harter Coulomb Obits 6-1 Quantum Theory for the Computer Age Unit 9 Caustic for KE/PE =R=-3/8 F p' p g r p"... P F' F P Diving orbit T" T T' Contact Pt. for KE/PE =R=-3/8 Quantum Orbits W.G. Harter
More informationRoots of Unity, Cyclotomic Polynomials and Applications
Swiss Mathematical Olympiad smo osm Roots of Unity, Cyclotomic Polynomials and Applications The task to be done here is to give an introduction to the topics in the title. This paper is neither complete
More informationNow derive both sides of above equation with respect to ψ: (known as Abel integral equation) the above equation is invertible thanks to the identity:
Change of variable in the integral: Now derive both sides of above equation with respect to ψ: (known as Abel integral equation) the above equation is invertible thanks to the identity: so that one gets
More information1 Ordinary points and singular points
Math 70 honors, Fall, 008 Notes 8 More on series solutions, and an introduction to \orthogonal polynomials" Homework at end Revised, /4. Some changes and additions starting on page 7. Ordinary points and
More information13 Spherical Coordinates
Utah State University DigitalCommons@USU Foundations of Wave Phenomena Library Digital Monographs 8-204 3 Spherical Coordinates Charles G. Torre Department of Physics, Utah State University, Charles.Torre@usu.edu
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationQuantum Theory of Angular Momentum and Atomic Structure
Quantum Theory of Angular Momentum and Atomic Structure VBS/MRC Angular Momentum 0 Motivation...the questions Whence the periodic table? Concepts in Materials Science I VBS/MRC Angular Momentum 1 Motivation...the
More informationSpezielle Funktionen 3 Orthogonale Polynome
Spezielle Funktionen 3 Orthogonale Polynome M. Gutting 28. Mai 2015 3.5 Application of the Legendre Polynomials in Electrostatics Let x R 3 and ρ(x) be a charge distribution with total charge R 3 ρ(x)dx.
More information3. Quantum Mechanics in 3D
3. Quantum Mechanics in 3D 3.1 Introduction Last time, we derived the time dependent Schrödinger equation, starting from three basic postulates: 1) The time evolution of a state can be expressed as a unitary
More informationSpiral Structure Formed in a Pair of Interacting Galaxies
J. Astrophys. Astr. (1993) 14, 19 35 Spiral Structure Formed in a Pair of Interacting Galaxies Ch. L. Vozikis & Ν. D.Caranicolas Department of Physics, Section of Astrophysics, Astronomy and Mechanics,
More informationQuantum Mechanics-I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras. Lecture - 21 Square-Integrable Functions
Quantum Mechanics-I Prof. Dr. S. Lakshmi Bala Department of Physics Indian Institute of Technology, Madras Lecture - 21 Square-Integrable Functions (Refer Slide Time: 00:06) (Refer Slide Time: 00:14) We
More informationCollisionless Boltzmann Eq (Vlasov eq)" S+G sec 3.4! Collisionless Boltzmann Eq S&G 3.4!
Collisionless Boltzmann Eq (Vlasov eq)" S+G sec 3.4 When considering the structure of galaxies cannot follow each individual star (10 11 of them), Consider instead stellar density and velocity distributions.
More informationLinear Algebra. and
Instructions Please answer the six problems on your own paper. These are essay questions: you should write in complete sentences. 1. Are the two matrices 1 2 2 1 3 5 2 7 and 1 1 1 4 4 2 5 5 2 row equivalent?
More information(a) Determine the general solution for φ(ρ) near ρ = 0 for arbitary values E. (b) Show that the regular solution at ρ = 0 has the series expansion
Problem 1. Curious Wave Functions The eigenfunctions of a D7 brane in a curved geometry lead to the following eigenvalue equation of the Sturm Liouville type ρ ρ 3 ρ φ n (ρ) = E n w(ρ)φ n (ρ) w(ρ) = where
More informationPHYS 404 Lecture 1: Legendre Functions
PHYS 404 Lecture 1: Legendre Functions Dr. Vasileios Lempesis PHYS 404 - LECTURE 1 DR. V. LEMPESIS 1 Legendre Functions physical justification Legendre functions or Legendre polynomials are the solutions
More informationBMT 2016 Orthogonal Polynomials 12 March Welcome to the power round! This year s topic is the theory of orthogonal polynomials.
Power Round Welcome to the power round! This year s topic is the theory of orthogonal polynomials. I. You should order your papers with the answer sheet on top, and you should number papers addressing
More informationTransformed E&M I homework. Multipole Expansion (Griffiths Chapter 3)
Transformed E&M I homework Multipole Expansion (Griffiths Chapter 3) Multipole Expansion Question 1. Multipole moment of charged wire CALCULATION; EXPANSION (U. Nauenberg, HW3, solutions available) A charge
More information+ y = 1 : the polynomial
Notes on Basic Ideas of Spherical Harmonics In the representation of wavefields (solutions of the wave equation) one of the natural considerations that arise along the lines of Huygens Principle is the
More informationLegendre s Equation. PHYS Southern Illinois University. October 13, 2016
PHYS 500 - Southern Illinois University October 13, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 13, 2016 1 / 10 The Laplacian in Spherical Coordinates The Laplacian is given
More informationbe made, in which the basis functions are localised in real space. Examples include truncated Gaussian orbitals, orbitals based on pseudoatomic wave f
Localised spherical-wave basis set for O(N ) total-energy pseudopotential calculations P.D. Haynes and M.C. Payne Theory of Condensed Matter, Cavendish Laboratory, Madingley Road, Cambridge, CB3 0HE, U.K.
More informationSpecial Functions of Mathematical Physics
Arnold F. Nikiforov Vasilii B. Uvarov Special Functions of Mathematical Physics A Unified Introduction with Applications Translated from the Russian by Ralph P. Boas 1988 Birkhäuser Basel Boston Table
More informationMore On Carbon Monoxide
More On Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results Jerry Gilfoyle The Configurations of CO 1 / 26 More On Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results Jerry Gilfoyle The Configurations
More informationExpansion of 1/r potential in Legendre polynomials
Expansion of 1/r potential in Legendre polynomials In electrostatics and gravitation, we see scalar potentials of the form V = K d Take d = R r = R 2 2Rr cos θ + r 2 = R 1 2 r R cos θ + r R )2 Use h =
More informationMidterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015
Midterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015 The test lasts 1 hour and 15 minutes. No documents are allowed. The use of a calculator, cell phone or other equivalent electronic
More informationCHEM-UA 127: Advanced General Chemistry I
1 CHEM-UA 127: Advanced General Chemistry I Notes for Lecture 11 Nowthatwehaveintroducedthebasicconceptsofquantummechanics, wecanstarttoapplythese conceptsto build up matter, starting from its most elementary
More informationby Ron Buckmire Mathematics Department Occidental College Los Angeles, CA , U.S.A.
A new nite-dierence scheme for singular dierential equations in cylinical or spherical coordinates by Ron Buckmire Mathematics Department Occidental College Los Angeles, CA 90041-3338, U.S.A. December
More informationLecture 9: Molecular integral. Integrals of the Hamiltonian matrix over Gaussian-type orbitals
Lecture 9: Molecular integral evaluation Integrals of the Hamiltonian matrix over Gaussian-type orbitals Gaussian-type orbitals The de-facto standard for electronic-structure calculations is to use Gaussian-type
More informationA Fourier-Bessel Expansion for Solving Radial Schrödinger Equation in Two Dimensions
A Fourier-Bessel Expansion for Solving Radial Schrödinger Equation in Two Dimensions H. TAŞELI and A. ZAFER Department of Mathematics, Middle East Technical University, 06531 Ankara, Turkey ABSTRACT The
More information1 Reduced Mass Coordinates
Coulomb Potential Radial Wavefunctions R. M. Suter April 4, 205 Reduced Mass Coordinates In classical mechanics (and quantum) problems involving several particles, it is convenient to separate the motion
More informationSection 5.2 Series Solution Near Ordinary Point
DE Section 5.2 Series Solution Near Ordinary Point Page 1 of 5 Section 5.2 Series Solution Near Ordinary Point We are interested in second order homogeneous linear differential equations with variable
More informationEigenvectors and Hermitian Operators
7 71 Eigenvalues and Eigenvectors Basic Definitions Let L be a linear operator on some given vector space V A scalar λ and a nonzero vector v are referred to, respectively, as an eigenvalue and corresponding
More informationMODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione
MODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione, Univ. di Roma Tor Vergata, via di Tor Vergata 11,
More informationSpherical Harmonics on S 2
7 August 00 1 Spherical Harmonics on 1 The Laplace-Beltrami Operator In what follows, we describe points on using the parametrization x = cos ϕ sin θ, y = sin ϕ sin θ, z = cos θ, where θ is the colatitude
More informationDynamical friction, galaxy merging, and radial-orbit instability in MOND
Dynamical friction, galaxy merging, and radial-orbit instability in MOND Carlo Nipoti Dipartimento di Astronomia Università di Bologna Strasbourg, 29 June 2010 Outline MOND and Equivalent Newtonian Systems
More information1 Commutators (10 pts)
Final Exam Solutions 37A Fall 0 I. Siddiqi / E. Dodds Commutators 0 pts) ) Consider the operator  = Ĵx Ĵ y + ĴyĴx where J i represents the total angular momentum in the ith direction. a) Express both
More informationQuantum Continuum Mechanics for Many-Body Systems
Workshop on High Performance and Parallel Computing Methods and Algorithms for Materials Defects, 9-13 February 2015 Quantum Continuum Mechanics for Many-Body Systems J. Tao 1,2, X. Gao 1,3, G. Vignale
More information1 Astronomical Observation
Structural Unification of Newtonian and Rational Gravity Jin He Wuhan FutureSpace Scientific Corporation Limited, Wuhan, Hubei 430074, China E-mail: mathnob@yahoo.com Abstract If there were no electromagnetic
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationMATH 31 - ADDITIONAL PRACTICE PROBLEMS FOR FINAL
MATH 3 - ADDITIONAL PRACTICE PROBLEMS FOR FINAL MAIN TOPICS FOR THE FINAL EXAM:. Vectors. Dot product. Cross product. Geometric applications. 2. Row reduction. Null space, column space, row space, left
More informationSchrödinger equation for central potentials
Chapter 2 Schrödinger equation for central potentials In this chapter we will extend the concepts and methods introduced in the previous chapter ifor a one-dimenional problem to a specific and very important
More informationOn prediction and density estimation Peter McCullagh University of Chicago December 2004
On prediction and density estimation Peter McCullagh University of Chicago December 2004 Summary Having observed the initial segment of a random sequence, subsequent values may be predicted by calculating
More informationAn Overview of Fast Multipole Methods
An Overview of Fast Multipole Methods Alexander T. Ihler April 23, 2004 Sing, Muse, of the low-rank approximations Of spherical harmonics and O(N) computation... Abstract We present some background and
More information(ii) Observational Geomagnetism. Lecture 5: Spherical harmonic field models
(ii) Observational Geomagnetism Lecture 5: Spherical harmonic field models Lecture 5: Spherical harmonic field models 5.1 Introduction 5.2 How to represent functions on a spherical surface 5.3 Spherical
More informationClassical Field Theory: Electrostatics-Magnetostatics
Classical Field Theory: Electrostatics-Magnetostatics April 27, 2010 1 1 J.D.Jackson, Classical Electrodynamics, 2nd Edition, Section 1-5 Electrostatics The behavior of an electrostatic field can be described
More informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationGravitational Potential Energy of Interpenetrating Spherical Galaxies in Hernquist s Model
J. Astrophys. Astr. (1997) 18, 23 32 Gravitational Potential Energy of Interpenetrating Spherical Galaxies in Hernquist s Model Κ. S. V. S. Narasimhan, K. S. Sastry* & Saleh Mohammed Alladin Inter University,
More informationFrame Diagonalization of Matrices
Frame Diagonalization of Matrices Fumiko Futamura Mathematics and Computer Science Department Southwestern University 00 E University Ave Georgetown, Texas 78626 U.S.A. Phone: + (52) 863-98 Fax: + (52)
More informationFun With Carbon Monoxide. p. 1/2
Fun With Carbon Monoxide p. 1/2 p. 1/2 Fun With Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results p. 1/2 Fun With Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results C V (J/K-mole) 35 30 25
More informationEnergy spectrum for a short-range 1/r singular potential with a nonorbital barrier using the asymptotic iteration method
Energy spectrum for a short-range 1/r singular potential with a nonorbital barrier using the asymptotic iteration method A. J. Sous 1 and A. D. Alhaidari 1 Al-Quds Open University, Tulkarm, Palestine Saudi
More informationQuantum Physics II (8.05) Fall 2002 Assignment 11
Quantum Physics II (8.05) Fall 00 Assignment 11 Readings Most of the reading needed for this problem set was already given on Problem Set 9. The new readings are: Phase shifts are discussed in Cohen-Tannoudji
More informationFast Multipole Methods for The Laplace Equation. Outline
Fast Multipole Methods for The Laplace Equation Ramani Duraiswami Nail Gumerov Outline 3D Laplace equation and Coulomb potentials Multipole and local expansions Special functions Legendre polynomials Associated
More informationKadath: a spectral solver and its application to black hole spacetimes
Kadath: a spectral solver and its application to black hole spacetimes Philippe Grandclément Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris F-92195 Meudon, France philippe.grandclement@obspm.fr
More informationMethod of Frobenius. General Considerations. L. Nielsen, Ph.D. Dierential Equations, Fall Department of Mathematics, Creighton University
Method of Frobenius General Considerations L. Nielsen, Ph.D. Department of Mathematics, Creighton University Dierential Equations, Fall 2008 Outline 1 The Dierential Equation and Assumptions 2 3 Main Theorem
More informationBondi-Sachs Formulation of General Relativity (GR) and the Vertices of the Null Cones
Bondi-Sachs Formulation of General Relativity (GR) and the Vertices of the Null Cones Thomas Mädler Observatoire de Paris/Meudon, Max Planck Institut for Astrophysics Sept 10, 2012 - IAP seminar Astrophysical
More informationSpherical Wave Approach (SWA) for the electromagnetic modelling of 3D GPR scenarios. Nicola Tedeschi
Spherical Wave Approach (SWA) for the electromagnetic modelling of 3D GPR scenarios Nicola Tedeschi Outline Solutions of the vector wave equation Spherical wave approach Spherical reflected and transmitted
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 9, February 8, 2006
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer Lecture 9, February 8, 2006 The Harmonic Oscillator Consider a diatomic molecule. Such a molecule
More information