Proof of Bertrand's Theorem*

Similar documents
A = 6561 times greater. B. 81 times greater. C. equally strong. D. 1/81 as great. E. (1/81) 2 = 1/6561 as great Pearson Education, Inc.

Chapter 1. Introduction to Nonlinear Space Plasma Physics

Question 1: Spherical Pendulum

Energy Wave Equations: Correction Factors

The Two -Body Central Force Problem

Central force motion/kepler problem. 1 Reducing 2-body motion to effective 1-body, that too with 2 d.o.f and 1st order differential equations

Physics Mechanics. Lecture 32 Oscillations II

QUARK EFFECT ON H-ATOM SPECTRA(1S)

Problem Set 5 Solution

The first term involves the cross product of two parallel vectors and so it vanishes. We then get

Physics 115/242 The Kepler Problem

Chapter 9 Circular Motion Dynamics

MATHEMATICAL MODELLING, MECHANICS AND MOD- ELLING MTHA4004Y

Q1. (a) Explain what is meant by the gravitational potential at a point in a gravitational field (2)

Problem 1: Spin 1 2. particles (10 points)

Fig. 1. On a sphere, geodesics are simply great circles (minimum distance). From

On logarithmic ratio quantities

Stability of circular orbits

Extending the Patched-Conic Approximation to the Restricted Four-Body Problem

01 Harmonic Oscillations

Physics 351, Spring 2015, Homework #7. Due at start of class, Friday, March 3, 2017

Physics 2020: Sample Problems for Exam 1

TAYLOR AND MACLAURIN SERIES

Physics 351, Spring 2015, Homework #7. Due at start of class, Friday, March 6, 2015

Time part of the equation can be separated by substituting independent equation

The 3 dimensional Schrödinger Equation

Lecture 1 Notes: 06 / 27. The first part of this class will primarily cover oscillating systems (harmonic oscillators and waves).

Tutorial on obtaining Taylor Series Approximations without differentiation

Lecture 6 Scattering theory Partial Wave Analysis. SS2011: Introduction to Nuclear and Particle Physics, Part 2 2

Explanation: The escape velocity and the orbital velocity for a satellite are given by

Harmonic Oscillator - Model Systems

II. Universal Gravitation - Newton 4th Law

Sample paper 1. Question 1. What is the dimensional formula of torque? A. MLT -2 B. MT -2 C. ML 2 T -2 D. MLT -1 E. ML 3 T -2.

Session 6: Analytical Approximations for Low Thrust Maneuvers

Polynomials and Rational Functions. Quadratic Equations and Inequalities. Remainder and Factor Theorems. Rational Root Theorem

The Restricted 3-Body Problem

LECTURE 10: REVIEW OF POWER SERIES. 1. Motivation

Perturbations to the Lunar Orbit

Q1. (a) Explain what is meant by the gravitational potential at a point in a gravitational field (2)

Physics 351, Spring 2017, Homework #3. Due at start of class, Friday, February 3, 2017

9.2 GRAVITATIONAL FIELD, POTENTIAL, AND ENERGY 9.4 ORBITAL MOTION HW/Study Packet

Physics 106a, Caltech 13 November, Lecture 13: Action, Hamilton-Jacobi Theory. Action-Angle Variables

Computational Problem: Keplerian Orbits

Department of Physics, Princeton University. Graduate Preliminary Examination Part I. Thursday, May 7, :00 am - 12:00 noon

Kinetic Theory. Motivation - Relaxation Processes Violent Relaxation Thermodynamics of self-gravitating system

Dual Program Level 1 Physics Course

Theory of mean motion resonances.

Quantum Mechanics for Scientists and Engineers

Brenda Rubenstein (Physics PhD)

Atomic Theory and Atomic structure. Part A. d) Distance of electrons from the nucleus

1 Supplementary Figure

PHYS 385 Lecture 1 - Introduction 1-1

What is the energy - Two laws of motion c.s. fields

2m 2 Ze2. , where δ. ) 2 l,n is the quantum defect (of order one but larger

ORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points)

Topic Outline for Algebra 2 & and Trigonometry One Year Program

GRAVITATIONAL LENSING BY A COLLECTION OF OBJECTS WITH RADIAL DENSITIES. 1. Introduction

MIT Weakly Nonlinear Things: Oscillators.

S9: Momentum Spread Effects and Bending S9A: Formulation

Symmetries 2 - Rotations in Space

Astromechanics. 6. Changing Orbits

Select/Special Topics in Atomic Physics Prof. P.C. Deshmukh Department Of Physics Indian Institute of Technology, Madras

Two-Body Problem. Central Potential. 1D Motion

d 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.

Physics 139B Solutions to Homework Set 4 Fall 2009

1 Chapter 2 Perform arithmetic operations with polynomial expressions containing rational coefficients 2-2, 2-3, 2-4

On the Disintegration of Very Small Satellites

1 Reduced Mass Coordinates

Forces: Explained and Derived by Energy Wave Equations

THE COMPUTATION OF FLUID-INDUCED FORCES ON CENTRIFUGAL IMPELLERS ROTATING AND WHIRLING IN A VOLUTE CASING. N s Specific speed (Q 1 2 =(gh) 3 4 )

Simplified Analytical Model of a Six-Degree-of-Freedom Large-Gap Magnetic Suspension System

A half submerged metal sphere (UIC comprehensive

Step index planar waveguide

Newton s Gravitational Law

A space probe to Jupiter

Physics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I

Mechanics Departmental Exam Last updated November 2013

Module 24: Angular Momentum of a Point Particle

Kinematics of the Solar Neighborhood

PHYS 1444 Section 02. Lecture #3

0n Stokes' Problem for a Non-Newtonian Fluid. K. R. Rajagopal, Pittsburgh, Pennsylvania, and T. Y. Na, Dearborn, Michigan. (Received August 30, 1982)

ν ηˆαˆβ The inverse transformation matrices are computed similarly:

Pade approximants and noise: rational functions

11.10a Taylor and Maclaurin Series

Lecture 1. Hydrodynamic Stability F. H. Busse NotesbyA.Alexakis&E.Evstatiev

Module 12: Work and the Scalar Product

cos 6 λ m sin 2 λ m Mirror Point latitude Equatorial Pitch Angle Figure 5.1: Mirror point latitude as function of equatorial pitch angle.

In deriving this we ve used the fact that the specific angular momentum

Basic Physics. What We Covered Last Class. Remaining Topics. Center of Gravity and Mass. Sun Earth System. PHYS 1411 Introduction to Astronomy

THE EDUCARE (SIROHI CLASSES) TEST SERIES 2018

Problem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension

February 18, In the parallel RLC circuit shown, R = Ω, L = mh and C = µf. The source has V 0. = 20.0 V and f = Hz.

GENERAL PRINCIPLES. Statics, Units, Calculations & problem Solving Students will be able to:

9.3 Worked Examples Circular Motion

The Relation of Relativistic Energy to Particle Wavelength

ME 012 Engineering Dynamics

Class XII_Delhi_Physics_Set-1

221B Lecture Notes Scattering Theory II

Angle recap. Angular position: Angular displacement: s. Angular velocity: Angular Acceleration:

HW5 Physics 311 Mechanics

Transcription:

APPENDIX A Proof of Bertr's Theorem* The orbit equation under a conservative central force, Eq. (3-34), may be written where The condition for a circular orbit of radius ro = uil, Eq. (3-41), now takes the form u0 = J(uo). (A-3) In addition, of course, the energy must satisfy the condition of Eq. (3-42). If the energy is slightly above that needed for circularity, the potential is such that the motion is stable, then u will remain bounded vary only slightly from u, J(u) can be expressed in terms of the first term in a Taylor series expansion about J(uo): As is customary, the derivative appearing in Eq. (A-4) is a shorth symbol for the derivative of J with respect to u evaluated at u = u,. If the difference u - uo is represented by x, the orbit equation for motion in the vicinity of the circularity conditions is then * See Section 3-6.

602 APPENDIX A where In order for x to describe a bounded stable oscillation,.p2 must be positive definite. From the definition, Eq. (A-2), we have In view of the circularity conditions, Eq. (3-41) or Eq. (A-3), it then follows that where, in addition to the convention empioyed for the derivatives, fo sts for f (llu,). Hence P2 is given by which is the same as Eq. (3-46), the stability condition, P2 > 0, thus reduces to Eq. (3-43). By a suitable choice of origin of 6, the solution to Eq. (A-5) for P2 positive definite can be written (cf. Eq. (3-45)). In order for the orbit to remain closed when the energy angular momentum are thus slightly disturbed from circularity, the quantity P must be a rational number. We are concerned with finding force laws such that for a wide range of initial conditions, i.e., for a wide range of uo, the orbits deviating slightly from circularity remain closed. Under these circumstances, as argued in the main text, fl must have the same value over the entire domain of u,, Eq. (A-7) can be looked on as a differential equation for f (llu) or f (r). The desired force law must therefore conform to a dependence on r given in Eq. (3-48): where k is some constant fi is a rational number. The force law of Eq. (A-8) still permits a wide variety of behavior for the force. However we seek more stringent conditions on the force law, by requiring that even when the deviations from circularity are considerable the orbit remain closed. We must therefore at least deal with deviations of u from uo so large that we must keep more terms than the first in the Taylor series expansion of J(u). Equation (A-4) may therefore be replaced by

PROOF OF BERTRAND'S THEOREM 603 where it is understood the derivatives are evaluated at u = u,. In terms of this expansion of J(u) the orbit equation becomes We seek to find the nature of the source law such that even when the deviation from the circular orbit, x, is large enough that the terms on the right cannot be neglected, the solution to Eq. (A-10) still represents a closed orbit. For small perturbations from circularity we know x has the behavior described by Eq. (A-8), which represents the fundamental term in a Fourier expansion in terms of PO. We seek therefore a closed-orbit solution by including a few more terms in the Fourier expansion : x = a, + a, cos PO + a, cos 2PO + a, cos 3P8. (A-1 1) The amplitudes a, a, must be of smaller magnitude than a, because they vanish faster than a, as circularity is approached. As will be seen a, must be of even lower order of magnitude than a, or a,, which is why terms in cos 4PO so on can be neglected. Consequently in the x2 term on the right of Eq. (A-10) the factors in cos 3p6' are dropped, in the x3 term only factors in cos PO are kept. In evaluating the right-h side powers products of the cosine functions are reduced by means of such identities as cos 3p0 = i(3 cos PO + cos 380). Consistently keeping terms through the order of a: in this manner, Eq. (A-10) with the solution (A-11) can be reduced to P2ao - 3p2a2 cos 2/38-8P3a3 cos 3pe =-J"+ a: cos (A-12) 4 I + -J"cos~PO a: + cos 3pO. 4 24 In order for the solution to be valid, the coefficient of each cosine term must separately vanish, leading to four conditions on the amplitude the derivatives of J: (A- 13a) (A- 13b)

604 APPENDIX A (A- 13d) It should be remembered that we already have shown, on the basis of slight deviations from circularity, that for closed orbits the force law must have the form of Eq. (A-8) or that J(u) is given by (A- 14) Keeping in mind the circularity condition, Eq. (A-3), the various derivatives at uo can be evaluated as J"' = -P2(1 - P2)(1 + P2) u; (A- 15a) Equations (A-13a, b) thus say that ao/a, a2/a, are of the order of al/uo, which is by supposition a small number. Further, Eq. (A- 13d) shows that a3/al is of the order of (a,/~,)~, which justifies the earlier statement that a, is oflower order of magnitude than a, or a,. Equation (A-13c) is a condition on P only, the condition that in fact is the principal conclusion of Bertr's theorem. Substituting Eqs. (A-13a, b) (A-15) into Eq. (A-13c) yields the condition P2(1 - P2)(4 - P2) = 0. (A- 16) For deviations from a circular orbit, that is, p # 0, the only solutions are (A- 1 7a) P2 = 4, f (r) = - kr. (A-17b) Thus the only two possible force laws consistent with the solution are either the gravitational inverse-square law or Hooke's law! We started out with orbits that were circular. These are possible for all attractive force laws over a wide range of 1 E, whose values in turn fix the orbital radius. The requirement that the circular orbit be stable for all radii already restricts the form of the force law through the inequality condition of P2 > 0 (Eq. (3-48)). If we further seek force laws such that orbits that deviate only slightly from a circular orbit are still closed, no matter what the radius of the

PROOF OF BERTRAND'S THEOREM 605 reference orbit, then the force laws are restricted to the discrete set given by Eq. (A-8) with rational values of P. In order for the orbits to remain closed for larger deviations from circularity, no matter what the initial conditions of the reference orbit, only two of these rational values are permitted: = 1 IPI = 2. Since we know that these attractive force laws do in fact give closed orbits for all E 1 leading to bounded motion, they must be the only force laws leading to closed orbits for all bounded motion.