Bertrand s Theorem. October 8, µr 2 + V (r) 0 = dv eff dr. 3 + dv. f (r 0 )

Similar documents
Complex Variables. For ECON 397 Macroeconometrics Steve Cunningham

m = Average Rate of Change (Secant Slope) Example:

Integration of Basic Functions. Session 7 : 9/23 1

SUBJECT:ENGINEERING MATHEMATICS-I SUBJECT CODE :SMT1101 UNIT III FUNCTIONS OF SEVERAL VARIABLES. Jacobians

Section 7.4: Integration of Rational Functions by Partial Fractions

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University

Chapter 6. Inverse Circular Functions and Trigonometric Equations. Section 6.1 Inverse Circular Functions y = 0

Lecture 7 Waveguides. TC 412 Microwave Communications

Figure 1 Probability density function of Wedge copula for c = (best fit to Nominal skew of DRAM case study).

Objectives: We will learn about filters that are carried out in the frequency domain.

FEA Solution Procedure

The Linear Quadratic Regulator

Chapter 9 Flow over Immersed Bodies

Applying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB

3.4-Miscellaneous Equations

Ground Rules. PC1221 Fundamentals of Physics I. Position and Displacement. Average Velocity. Lectures 7 and 8 Motion in Two Dimensions

3. Several Random Variables

FEA Solution Procedure

Unknown Input High Gain Observer for Parametric Fault Detection and Isolation of Dynamical Systems

Pendulum Equations and Low Gain Regime

Higher Maths A1.3 Recurrence Relations - Revision

FEA Solution Procedure

Lecture: Corporate Income Tax

LIGHTWEIGHT STRUCTURES in CIVIL ENGINEERING - CONTEMPORARY PROBLEMS

Second-Order Wave Equation

AMS 212B Perturbation Methods Lecture 05 Copyright by Hongyun Wang, UCSC

Image and Multidimensional Signal Processing

PhysicsAndMathsTutor.com

Chapter 1: Differential Form of Basic Equations

Prandl established a universal velocity profile for flow parallel to the bed given by

OPTI-502 Optical Design and Instrumentation I John E. Greivenkamp Final Exam In Class Page 1/14 Fall, 2017

Formal Methods for Deriving Element Equations

Physics 5153 Classical Mechanics. Solution by Quadrature-1

Math 116 First Midterm October 14, 2009

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.

Lecture: Corporate Income Tax - Unlevered firms

Essentials of optimal control theory in ECON 4140

Optimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications

10.4 Solving Equations in Quadratic Form, Equations Reducible to Quadratics

STEP Support Programme. STEP III Hyperbolic Functions: Solutions

MEG 741 Energy and Variational Methods in Mechanics I

Chem 4501 Introduction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics. Fall Semester Homework Problem Set Number 10 Solutions

FRTN10 Exercise 12. Synthesis by Convex Optimization

APPENDIX B MATRIX NOTATION. The Definition of Matrix Notation is the Definition of Matrix Multiplication B.1 INTRODUCTION

Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS

Linear Strain Triangle and other types of 2D elements. By S. Ziaei Rad

Conditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane

Math 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation:

A Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane

Setting The K Value And Polarization Mode Of The Delta Undulator

Solving a System of Equations

A Note on Johnson, Minkoff and Phillips Algorithm for the Prize-Collecting Steiner Tree Problem

Two identical, flat, square plates are immersed in the flow with velocity U. Compare the drag forces experienced by the SHADED areas.

10.2 Solving Quadratic Equations by Completing the Square

Gravitational Instability of a Nonrotating Galaxy *

2. Find the coordinates of the point where the line tangent to the parabola 2

Einstein Classes, Unit No. 102, 103, Vardhman Ring Road Plaza, Vikas Puri Extn., Outer Ring Road New Delhi , Ph. : ,

Theoretical Fluid Mechanics Turbulent Flow Velocity Profile By James C.Y. Guo, Professor and P.E. Civil Engineering, U. of Colorado at Denver

Physics of the Interstellar and Intergalactic Medium

J. Basic. Appl. Sci. Res., 3(2s) , , TextRoad Publication

Elements of Coordinate System Transformations

PLANETARY ORBITS. According to MATTER (Re-examined) Nainan K. Varghese,

Shooting Method for Ordinary Differential Equations Autar Kaw

III. Demonstration of a seismometer response with amplitude and phase responses at:

Optimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance

SIO 211B, Rudnick. We start with a definition of the Fourier transform! ĝ f of a time series! ( )

MECHANICS OF SOLIDS COMPRESSION MEMBERS TUTORIAL 2 INTERMEDIATE AND SHORT COMPRESSION MEMBERS

6.4 VECTORS AND DOT PRODUCTS

A fundamental inverse problem in geosciences

Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introduction The transmission line equations are given by

Downloaded 01/04/14 to Redistribution subject to SIAM license or copyright; see

Right Trapezoid Cover for Triangles of Perimeter Two

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1

sin u 5 opp } cos u 5 adj } hyp opposite csc u 5 hyp } sec u 5 hyp } opp Using Inverse Trigonometric Functions

Efficiency Increase and Input Power Decrease of Converted Prototype Pump Performance

OPTI-502 Optical Design and Instrumentation I John E. Greivenkamp Final Exam In Class Page 1/16 Fall, 2015

Change of Variables. (f T) JT. f = U

Chapter 2 Introduction to the Stiffness (Displacement) Method. The Stiffness (Displacement) Method

Vibrational modes of a rotating string

PhysicsAndMathsTutor.com

Mean Value Formulae for Laplace and Heat Equation

Nonparametric Identification and Robust H Controller Synthesis for a Rotational/Translational Actuator

The Theory of Virtual Particles as an Alternative to Special Relativity

Chapter 5 Computational Methods for Mixed Models

Concept of Stress at a Point

Department of Industrial Engineering Statistical Quality Control presented by Dr. Eng. Abed Schokry

Inertial Instability of Arbitrarily Meandering Currents Governed by the Eccentrically Cyclogeostrophic Equation

Lecture 8: September 26

Two-media boundary layer on a flat plate

Electromagnetically induced transparency with a standing-wave drive in the frequency up-conversion regime

Homogeneous Liner Systems with Constant Coefficients

DYNAMICAL LOWER BOUNDS FOR 1D DIRAC OPERATORS. 1. Introduction We consider discrete, resp. continuous, Dirac operators

University of California, Berkeley Physics H7C Fall 1999 (Strovink) SOLUTION TO FINAL EXAMINATION

Constant acceleration 9E

Robust Shortest Path Planning and Semicontractive Dynamic Programming

The Hashemite University Department of Civil Engineering ( ) Dr. Hazim Dwairi 1

Success Center Math Tips

Math 1. 2-hours test May 13, 2017

Pulses on a Struck String

Approximate Solution of Real Definite Integrals Over a Circle in Adaptive Environment

Transcription:

Bertrand s Theorem October 8, Circlar orbits The eective potential, V e = has a minimm or maximm at r i and only i so we mst have = dv e L µr + V r = L µ 3 + dv = L µ 3 r r = L µ 3 At this radis, there is no net radial orce, so that circlar orbits are possible. Sch orbits may be stable or nstable, depending on the sign o the second derivative o the eective potential. Stable orbits occr only i < d V e = 3L µr 4 = 3L µr 4 + d V d r and since r = L, we have µ 3 d r < 3 r r General case Eliminating = L µr, the remaining, radial eqation o motion is µ r L µr 3 r =

Write this as an orbit eqation, sing ṙ = dt = L µr d r = dt = L µr L µr L µr 3 + L d r µr Sbstitting, = L µ r 5 + L µ r 4 d r Now let r = The eqation o motion now becomes = L µr 5 or = L 5 µ = L µ = L µ = µ r L µr 3 r = L µr 5 + L d r µr 4 L µr 3 r = d d r = d + L µr 4 d r d d + L µ = 3 d d d L µr 3 r + L 4 d µ 3 d L µ d L 3 µ r Finally, write the orce in terms o the potential, d + = µ L = dv = d = dv d d L µr 3 r d L 3 µ r dv d

Now we have simply d + = µ dv L d This is the eqation o a iven harmonic oscillator. This striking set o transormations is what happens when people spend 3 years working on a problem. Now consider circlar orbits. This means that does not change with at all, so we have d = and thereore, i we set h = µ dv L d then = h Next, expand and the orce or small pertrbations abot, = + η µ dv L d = h Sbstitting these into the eqation o motion, = h + h η + h η + d η + + η = h + h η + h η + = + h η + h η + d η + [ h ] η = h η + I h <, then the eqation has exponential instead o oscillatory soltions, and the circlar orbits are not stable. For stable orbits, we thereore set λ = h > I we neglect the qaatic and higher terms on the right side o the eqation, with initial condition η = when =, we have soltions η = A sin λ In order to reprodce the initial conditions ater some integer nmber, q, o complete orbits, we reqire so that η = η πq = A sin πqλ qλ = p with p another integer. We see that λ mst be rational, λ = p q 3

Assming the orce is a continos nction o position, h is continos and λ is also continos. Thereore, λ = p q or all nearly circlar orbits. Retrning to the deinition o λ, and sing the constancy o λ, Since the circlar orbit satisies λ = h = d d λ = µ L 3 µ L + µ L = h = µ L = L µ 3 so that λ = µ L 3 L λ = + µ L λ 3 = µ 3 + µ L L µ 3 This mst hold regardless o the vale o, so we may op the sbscript and integrate d d = λ 3 d = λ 3 r r d = λ 3 r r = Ar λ 3 and in order to have stable, pertrbatively closed orbits, the orce law mst be a rational power law. This give h as Now retrn to the ll eqation o motion h = µ A L λ 3 = µa L λ d η + + η = + h η + h η + 6 h η 3 + 4

and expand in a Forier series, We also need powers o η. Keeping p to third order, η = η + η cos λ + η cos λ + η 3 cos 3λ + η + η cos λ + η cos λ + η 3 cos 3λ + = η + η η cos λ + η η cos λ + η η 3 cos 3λ +η η cos λ + η cos λ + η η cos λ cos λ +η η cos λ + η η cos λ cos λ + η η 3 cos 3λ η + η cos λ + η cos λ + η 3 cos 3λ + 3 = η 3 η + η η cos λ + η η cos λ + η η 3 cos 3λ + and sing addition ormlas Thereore, +η η cos λ + η η cos λ + η η η cos λ cos λ +η η η cos λ + η η η cos λ cos λ +η η 3 cos 3λ +η η cos λ + η η cos λ + η η η cos λ cos λ +η η cos λ + η 3 cos 3 λ +η η η cos λ cos λ +η η cos λ + η η η cos λ cos λ +η η η cos λ cos λ +η η 3 cos 3λ cos λ = + cos λ cos λ cos λ = [cos λ + λ + cos λ λ] = [cos 3λ + cos λ] cos 3 λ = cos λ + cos λ = cos λ + 4 cos 3λ + cos λ 4 = 3 4 cos λ + cos 3λ 4 η + η cos λ + η cos λ + η 3 cos 3λ + = η + η + η η + η η cos λ + η η + η cos λ + η η 3 + η η cos 3λ η + η cos λ + η cos λ + η 3 cos 3λ + 3 = η 3 + 3η η cos λ + 3η η cos λ + 3η η 3 cos 3λ +3η η cos λ + 6η η η cos λ cos λ + η 3 cos 3 λ = η 3 + 3ηη cos λ + 3η η cos λ + 3η η 3 cos 3λ +3η η + cos λ + 6η η η [cos 3λ + cos λ] +η 3 3 4 cos λ + 4 η3 cos 3λ = η 3 + 3 η η + 3η η + 3η η η + 3 4 η3 cos λ 5

Sbstitting and expanding, + 3ηη + 3 η η cos λ + 3η η 3 + 4 η3 + 3η η η cos 3λ d η + + η = + h η + h η + 6 h η 3 + η λ cos λ 4λ η cos λ 9λ η 3 cos 3λ + + η = h η + η cos λ + η cos λ + η 3 cos 3λ + Eqating like terms, + h η + η cos λ + η cos λ + η 3 cos 3λ + + 6 h η + η cos λ + η cos λ + η 3 cos 3λ + 3 η = h η + h η + h η + 6 h η 3 + 6 h 3 η η [ η λ + η cos λ = h η + h η η + η η + 6 h 3ηη + 3η η η + 3 ] 4 η3 cos λ [ 4λ η cos λ + η cos λ = h η + h η η + η + 6 h 3ηη + 3 ] η η cos λ [ 9λ η 3 cos 3λ + η 3 cos 3λ = h η 3 + h η η 3 + η η + 6 h 3η η 3 + 4 ] η3 + 3η η η cos 3λ Now we keep only the lowest order terms rom each eqation, i.e., we have η η η and so on. We set h = λ and solve = λ η + h η + 6 h η 3 + 4 h η + 4 h η or, smmarizing, = λ η + 4 h η η = 4λ h η = h η η + η η + 8 η3 h = h λ h η 3 λ h η 3 + 8 η3 h = η 3 4λ 5 h + 3λ h η = λ h η η 3 = [ 8λ h η η + ] 4 η3 h η = 4λ h η = 4λ η3 5 h + 3λ h η = η 3 = 8λ λ h η [ h η η + 4 η3 h 6 ]

where η is constant. Now, sing we may ind the derivatives, h = µa L λ so we may solve the eqality λ = h = λ µa L λ µa L λ = h = λ λ µa L λ = λ λ h = + λ λ λ = η 3 4λ 5 h + 3λ h = λ 5λ 4 λ + 3λ + λ λ λ = 5λ λ + 3λ + λ λ = 8λ λ 4 + λ 6 = λ λ 4 λ and we see that the only candidate power laws, r = Ar λ 3, or stable, closed orbits are: r = Ar 3 r = Ar r = Ar The irst o these yields only perectly circlar orbits, so the only nontrivial cases are the inverse sqare law and Hookes law. This condition only shows that these power laws are necessary. That they are sicient to prodce closed orbits reqires solving or their orbits exactly. 7