Stellar Dynamics and Structure of Galaxies
|
|
- Peregrine Simon
- 5 years ago
- Views:
Transcription
1 Stellar Dynamics and Structure of Galaxies in a given potential Vasily Belokurov vasily@ast.cam.ac.uk Institute of Astronomy Lent Term / 59
2 1 Collisions Model requirements 2 in spherical 3 4 Orbital periods Example Outline I 2 / 59
3 Galaxies Part II Collisions Do we have to worry about collisions? Collisions Model requirements Globular clusters look densest, so obtain a rough estimate of collision timescale for them 3 / 59
4 Collisions Model requirements Collisions in globular clusters The case of NGC 2808 ρ M pc 3 M 0.8 M. n pc 3 is the star number density. We have σ r 13 km s 1 as the typical 1D speed of a star, so the 3D speed is 3 σ r (= σx 2 + σy 2 + σz 2 ) 20 km s 1. Since M R (see Fluids, or Stars, course notes), have R 0.8R. 4 / 59
5 Collisions Model requirements Collisions in globular clusters The case of NGC 2808 For a collision, need the volume π(2r ) 2 σt coll to contain one star, i.e. n 0 = 1/ ( π(2r ) 2 σt coll ) or (1.1) t coll = 1/ ( 4πR σn 2 ) 0 (1.2) R = 0.8R, n 0 = 10 5 pc 3, σ = 20km/s Putting in the numbers gives t coll yr. So direct collisions between stars are rare, but if you have 10 6 stars then there is a collision every 10 8 years, so they do happen. Note that NGC 2808 is 10 times denser than typical So, for now, ignore collisions, and we are left with stars orbiting in the potential from all the other stars in the system. 5 / 59
6 Collisions Model requirements Model requirements Model (e.g., a globular cluster) just as a self-gravitating collection of objects. Have a gravitational potential well Φ(r), approximately smooth if the number of particles >> 1. Conventionally take Φ( ) = 0. Stars orbit in the potential well, with time per orbit (for a globular cluster) 2R h /σ 10 6 years << age. Remember how to measure age for globular clusters? Stars give rise to Φ(r) by their mass, so for this potential in a steady state could average each star over its orbit to get ρ(r). The key problem is therefore self-consistently building a model which fills in the terms: Φ(r) stellar orbits ρ(r) Φ(r) (1.3) Note that in most observed cases we only have v line of sight (R), so it is even harder to model real systems. Self-consistent = orbits & stellar mass give ρ, which leads to Φ, which supports the orbits used to construct ρ 6 / 59
7 in spherical The law of attraction s of motion and Newtonian gravity GR not needed, since 10 < v < 10 3 km is << c = 3 10 s 5 km s GM rc 2??? The gravitational force per unit mass acting on a body due to a mass M at the origin is f = GM ˆr = GM r 2 r 3 r (1.4) We can write this in terms of a potential Φ, using ( ) 1 = 1 r r 2ˆr (1.5) 7 / 59
8 in spherical The corresponding potential So f = Φ (1.6) where Φ is a scalar, Φ = Φ(r) = GM (1.7) r Hence the potential due to a point mass M at r = r 1 is Φ(r) = GM r r 1 (1.8) 8 / 59
9 Density vs Potential in spherical From Hayashi et al, The shape of the gravitational potential in cold dark matter haloes 9 / 59
10 Galaxies Part II The law of motion in spherical 10 / 59
11 in spherical Particle of constant mass m at position r subject to a force F. : i.e. If F is due to a gravitational potential Φ(r), then d (mṙ) = F (1.9) dt m r = F (1.10) F = mf = m Φ (1.11) The angular momentum about the origin is H = r (mṙ).then dh dt where G is the torque about the origin. = r (m r) + mṙ ṙ = r F G (1.12) 11 / 59
12 in spherical Particle of constant mass m at position r subject to a force F. : i.e. If F is due to a gravitational potential Φ(r), then d (mṙ) = F (1.9) dt m r = F (1.10) F = mf = m Φ (1.11) The angular momentum about the origin is H = r (mṙ).then dh dt where G is the torque about the origin. = r (m r) + mṙ ṙ = r F G (1.12) 11 / 59
13 in spherical The kinetic energy Energy T = 1 2 mṙ. ṙ (1.13) dt dt = mṙ. r = F. ṙ (1.14) If F = m Φ, then dt dt = mṙ. Φ(r) (1.15) But if Φ is independent of t, the rate of change of Φ along an orbit is d dt Φ(r) = Φ. ṙ (1.16) from the chain rule 12 / 59
14 in spherical The kinetic energy Energy T = 1 2 mṙ. ṙ (1.13) dt dt = mṙ. r = F. ṙ (1.14) If F = m Φ, then dt dt = mṙ. Φ(r) (1.15) But if Φ is independent of t, the rate of change of Φ along an orbit is d dt Φ(r) = Φ. ṙ (1.16) from the chain rule 12 / 59
15 in spherical Hence Energy dt dt = m d Φ(r) (1.17) dt m d ( ) 1. ṙ + Φ(r) = 0 (1.18) dt 2ṙ E = 1 2ṙ. ṙ + Φ(r) (1.19) The total energy is constant for a given orbit 13 / 59
16 in spherical Hence Energy dt dt = m d Φ(r) (1.17) dt m d ( ) 1. ṙ + Φ(r) = 0 (1.18) dt 2ṙ E = 1 2ṙ. ṙ + Φ(r) (1.19) The total energy is constant for a given orbit 13 / 59
17 in spherical Φ(r) = Φ( r ) = Φ(r), so f = Φ = ˆr dφ dr. The orbital angular momentum H = mr ṙ, and dh dt in spherical dφ = r mf = m r ˆr = 0. (1.20) dr So the angular momentum per unit mass h = H/m = r ṙ is a constant vector, and is perpendicular to r and ṙ the particle stays in a plane through the origin which is perpendicular to h Check: r h, r + δr = r + ṙδt h since both r and ṙ h, so particle remains in the plane Thus the problem becomes a two-dimensional one to calculate the orbit use 2-D cylindrical coordinates (R, φ, z) at z = 0, or spherical polars (r, θ, φ) with θ = π 2. So, in 2D, use (R, φ) and (r, φ) interchangeably.. 14 / 59
18 in spherical The equation of motion in can be written in radial angular terms, using r = rˆr = rê r + 0ê φ, so r =(r, 0). We know that d dt êr = φê φ (1.21) and d dt êφ = φê r (1.22) ê r = cos(φ)ê x + sin(φ)ê y ê φ = sin(φ)ê x + cos(φ)ê y d dt êr = sin(φ) φê x + cos(φ) φê y d dt êφ = cos(φ) φê x sin(φ) φê y 15 / 59
19 in spherical The equation of motion in can be written in radial angular terms, using r = rˆr = rê r + 0ê φ, so r =(r, 0). We know that d dt êr = φê φ (1.21) and d dt êφ = φê r (1.22) ê r = cos(φ)ê x + sin(φ)ê y ê φ = sin(φ)ê x + cos(φ)ê y d dt êr = sin(φ) φê x + cos(φ) φê y d dt êφ = cos(φ) φê x sin(φ) φê y 15 / 59
20 in spherical Hence and so ṙ = ṙê r + r φê φ (1.23) or ṙ = v =(ṙ, r φ) r = rê r + ṙ φê φ + ṙ φê φ + r φê φ r φ 2 ê r = ( r r φ 2 )ê r + 1 d ( ) r 2 φ ê φ r dt = a = [ r r φ 2, 1 d ( ) r 2 φ ] r dt (1.24) In general f =(f r, f φ ), and then f r = r r φ 2, where the second term ( ) is the centrifugal force, since we are in a rotating frame, and the torque rf φ = d dt r 2 φ (= r f). In a spherical potential f φ = 0, so r 2 φ is constant. 16 / 59
21 in spherical To determine the shape of the orbit we need to remove t from the equations and find r(φ). It is simplest to set u = 1/r, and then from r 2 φ = h obtain Then and φ = hu 2 (1.25) ṙ = 1 u 2 u = 1 u 2 du dφ φ = h du dφ (1.26) r = h d 2 u dφ 2 φ = h 2 u 2 d 2 u dφ 2. (1.27) 17 / 59
22 in spherical So the radial equation of motion becomes r r φ 2 = f r h 2 u 2 d 2 u dφ 2 1 u h2 u 4 = f r (1.28) d 2 u dφ 2 + u = f r h 2 u 2 (1.29) The orbit equation in spherical potential 18 / 59
23 in spherical Since f r is just a function of r (or u) this is an equation for u(φ), i.e. r(φ) - the path of the orbit. Note that it does not give r(t), or φ(t) - you need one of the other equations for those. If we take f r = GM r = GMu 2, then 2 d 2 u dφ 2 + u = GM/h2 (1.30) (which is something you will have seen in the Relativity course). 19 / 59
24 in spherical The solution to this equation is Kepler orbits Solution to the equation of motion l r = lu = 1 + e cos(φ φ 0) (1.31) which you can verify simply by putting it in the differential equation. Then e cos(φ φ 0) l e cos(φ φ 0) l = GM h 2 so l = h 2 /GM and e and φ 0 are constants of integration. 20 / 59
25 in spherical 1 r = 1 + e cos(φ φ 0) l Kepler orbits Bound orbits l Note that if e < 1 then 1/r is never zero, so r is bounded in the range 1+e < r < l 1 e. Also, in all cases the orbit is symmetric about φ = φ 0, so we take φ 0 = 0 as defining the reference line for the angle φ. l is the distance from the origin for φ = ± π 2 (with φ measured relative to φ 0 ). 21 / 59
26 in spherical Kepler orbits Bound orbits We can use different parameters. Knowing that the point of closest approach (perihelion for a planet in orbit around the Sun, periastron for something about a star) is at l/(1 + e) when φ = 0 and the aphelion (or whatever) is at l/(1 e) when φ = π, we can set the distance between these two points (= major axis of the orbit)=2a. Then l 1 + e + l 1 e = 2a l(1 e) + l(1 + e) = 2a(1 e2 ) (1.32) l = a(1 e 2 ) (1.33) r P = a(1 e) is the perihelion distance from the gravitating mass at the origin, and r a = a(1 + e) is the aphelion distance. The distance of the Sun from the midpoint is ae, and the angular momentum h 2 = GMl = GMa(1 e 2 ). 22 / 59
27 in spherical The energy per unit mass E = 1. ṙ + Φ(r) = ṙ 2ṙ 2 r 2 φ 2 GM (1.34) r r p = a(1 e) This is constant along the orbit, so we can evaluate it anywhere convenient - e.g. at perihelion where ṙ = 0. Then φ = h and so rp 2 E = 1 GMa(1 e 2 ) 2 a 2 (1 e) 2 GM a(1 e) = GM [ ( ) e 1 ] a 2 1 e 1 e = GM 2a (1.35) This is < 0 for a bound orbit, and depends only on the semi-major axis a (and not e). 23 / 59
28 in spherical... deduced from observations, and explained by Newtonian theory of gravity. 24 / 59
29 in spherical 1 are ellipses with the Sun at a focus. 2 Planets sweep out equal areas in equal time δa = 1 2 r 2 δφ [= 1 r(rδφ)] (1.36) 2 da dt = 1 2 r 2 h φ = = constant (1.37) 2 Kepler s second law is a consequence of a central force, since this is why h is a constant. 25 / 59
30 in spherical 1 are ellipses with the Sun at a focus. 2 Planets sweep out equal areas in equal time δa = 1 2 r 2 δφ [= 1 r(rδφ)] (1.36) 2 da dt = 1 2 r 2 h φ = = constant (1.37) 2 Kepler s second law is a consequence of a central force, since this is why h is a constant. 25 / 59
31 in spherical 3 (Period) 2 (size of orbit) 3 ( In one period T, the area swept out is A = 1 2 ht = T 0 But A = area of ellipse = πab = πa 2 1 e 2 [ Have π 0 = l2 2 A = = 2π 0 2π 0 2π 0 r dφ r 2 dφ rdr dφ (1 + e cos φ) 2 ) da dt dt 3rd Law l r = lu = 1 + e cos(φ φ 0) dx (a + b cos x) 2 = π a a 2 b 2 a2 b 2 26 / 59
32 in spherical so Since l = a(1 e 2 ) this implies and since b = a 1 e 2, ] A = 2 l2 2 π 1 e e 2 A = πa 2 1 e 2 A = πab 3rd Law def: e = 1 b2 a 2 27 / 59
33 in spherical Therefore T = 2A h = 2πa2 1 e 2 h 2πa 2 1 e = 2 GMa(1 e2 ) 3rd Law since h 2 = GMa(1 e 2 ) a 3 T = 2π GM T 2 a 3 (1.38) where in this case M is the mass of the Sun. Note: Since E = GM 2πGM 2a, the period T =. ( 2E) / 59
34 in spherical What happens to l r = 1 + e cos φ when e 1? If e > 1 then 1 + e cos φ = 0 has solutions φ where r = cos φ = 1/e Then φ φ φ, and, since cos φ is negative, π 2 < φ < π. The orbit is a hyperbola. If e = 1 then the particle just gets to infinity at φ = ±π - it is a parabola. 29 / 59
35 in spherical What happens to l r = 1 + e cos φ when e 1? If e > 1 then 1 + e cos φ = 0 has solutions φ where r = cos φ = 1/e Then φ φ φ, and, since cos φ is negative, π 2 < φ < π. The orbit is a hyperbola. If e = 1 then the particle just gets to infinity at φ = ±π - it is a parabola. 29 / 59
36 in spherical What happens to l r = 1 + e cos φ when e 1? If e > 1 then 1 + e cos φ = 0 has solutions φ where r = cos φ = 1/e Then φ φ φ, and, since cos φ is negative, π 2 < φ < π. The orbit is a hyperbola. If e = 1 then the particle just gets to infinity at φ = ±π - it is a parabola. 29 / 59
37 Kepler orbits in spherical 30 / 59
38 in spherical Energies for these unbound orbits: So, as r E 1 2ṙ 2 E = h 2 2ṙ 2 r 2 GM r r 2 φ = h 31 / 59
39 in spherical Recall d dt of this and since h = r 2 φ l r = 1 + e cos φ l ṙ = e sin φ φ r 2 ṙ = eh l sin φ As r cos φ 1/e E 2ṙ2 1 = 1 e 2 h 2 (1 2 l 2 1e ) 2 = GM 2l (e2 1) (recalling that h 2 = GMl) Thus E > 0 if e > 1 and for parabolic orbits (e = 1) E = / 59
40 in spherical We have seen that in a fixed potential Φ(r) a particle has constant energy E = 1 2ṙ2 + Φ(r) along an orbit. If we adopt the usual convention and take Φ(r) 0 as r, then if at some point r 0 the particle has velocity v 0 such that 1 2 v2 0 + Φ(r 0 ) > 0 then it is able to reach infinity. So at each point r 0 we can define an escape velocity v esc such that v esc = 2Φ(r 0 ) 33 / 59
41 in spherical The escape velocity from the Sun ( ) 1 2GM 2 ( r0 ) 1 2 v esc = = 42.2 a.u. Note: The circular velocity v circ is such that r φ 2 = GM r 2 r φ GM ( r0 ) 1 2 = v circ = = 29.8 r 0 a.u. r 0 From the Solar neighborhood km s 1 km s 1 (= 2π a.u./yr). v esc = 2v circ for a point mass source of the gravitational potential. 34 / 59
42 Galaxies Part II From the Galaxy in spherical 35 / 59
43 Kepler orbits in spherical 36 / 59
44 What we have done so far is assume a potential due to a fixed point mass which we take as being at the origin of our polar coordinates. We now wish to consider a situation in which we have two point masses, M 1 and M 2 both moving under the gravitational attraction of the other. This is a cluster of N stars where N = 2 and we can solve it exactly! Hooray! The potential is no longer fixed at origin Φ(r) = GM 1 r r 1 GM 2 r r 2 37 / 59
45 What we have done so far is assume a potential due to a fixed point mass which we take as being at the origin of our polar coordinates. We now wish to consider a situation in which we have two point masses, M 1 and M 2 both moving under the gravitational attraction of the other. This is a cluster of N stars where N = 2 and we can solve it exactly! Hooray! The potential is no longer fixed at origin Φ(r) = GM 1 r r 1 GM 2 r r 2 37 / 59
46 What we have done so far is assume a potential due to a fixed point mass which we take as being at the origin of our polar coordinates. We now wish to consider a situation in which we have two point masses, M 1 and M 2 both moving under the gravitational attraction of the other. This is a cluster of N stars where N = 2 and we can solve it exactly! Hooray! The potential is no longer fixed at origin Φ(r) = GM 1 r r 1 GM 2 r r 2 37 / 59
47 Or the force acting on star 1, due to star 2 is F 1 = GM 1M 2 r 1 r 2 2 in the direction of r 2 r 1 And by symmetry (or Newton s 3rd law) F 1 = GM 1M 2 r 1 r 2 3 (r 2 r 1 ) F 2 = GM 1M 2 r 1 r 2 3 (r 1 r 2 ) 38 / 59
48 Then we know and where M 1 r 1 = GM 1M 2 d 2 ˆd (1.39) M 2 r 2 = GM ( ) 1M 2 d 2 ˆd is the vector from M 2 to M 1. Using these two we can write for d = r 1 r 2 (1.40) d = r 1 r 2 (1.41) d = G(M 1 + M 2 ) d 2 ˆd (1.42) 39 / 59
49 d = G(M 1 + M 2 ) ˆd which is identical to the equation of motion of a particle subject to a fixed mass M 1 + M 2 at the origin. So we know that the period a T = 2π 3 (1.43) G(M 1 + M 2 ) d 2 where the size (maximum separation) of the relative orbit is 2a. 40 / 59
50 d = G(M 1 + M 2 ) ˆd which is identical to the equation of motion of a particle subject to a fixed mass M 1 + M 2 at the origin. So we know that the period a T = 2π 3 (1.43) G(M 1 + M 2 ) d 2 where the size (maximum separation) of the relative orbit is 2a. 40 / 59
51 If we take the coordinates for the centre of mass r CM = From equations (1.39) and (1.40) we know that M 1 M 2 r 1 + r 2 (1.44) M 1 + M 2 M 1 + M 2 and so or i.e.ṙ CM =constant. M 1 r 1 + M 2 r 2 = 0 (1.45) d dt (M 1ṙ 1 + M 2 ṙ 2 ) = 0 (1.46) (M 1 ṙ 1 + M 2 ṙ 2 ) = constant (1.47) We can choose an inertial frame in which the centre of mass has zero velocity, so might as well do so 41 / 59
52 If we take the coordinates for the centre of mass r CM = From equations (1.39) and (1.40) we know that M 1 M 2 r 1 + r 2 (1.44) M 1 + M 2 M 1 + M 2 and so or i.e.ṙ CM =constant. M 1 r 1 + M 2 r 2 = 0 (1.45) d dt (M 1ṙ 1 + M 2 ṙ 2 ) = 0 (1.46) (M 1 ṙ 1 + M 2 ṙ 2 ) = constant (1.47) We can choose an inertial frame in which the centre of mass has zero velocity, so might as well do so 41 / 59
53 Note that choosing r CM = 0 M 1 r 1 = M 2 r 2, and so r 1 = d + r 2 = d M1 This r 1 = M2 M 1+M 2 d and similarly r 2 = M1 M 1+M 2 d. The angular momentum J (or H if you want) is M 2 r 1 J = M 1 r 1 ṙ 1 + M 2 r 2 ṙ 2 = = M 1 M2 2 (M 1 + M 2 ) 2 d ḋ + M 2M1 2 (M 1 + M 2 ) 2 d ḋ M 1 M 2 d M 1 + M ḋ 2 (1.48) So J = µh (1.49) where µ is the reduced mass, and h is the specific angular momentum. 42 / 59
54 Note that choosing r CM = 0 M 1 r 1 = M 2 r 2, and so r 1 = d + r 2 = d M1 This r 1 = M2 M 1+M 2 d and similarly r 2 = M1 M 1+M 2 d. The angular momentum J (or H if you want) is M 2 r 1 J = M 1 r 1 ṙ 1 + M 2 r 2 ṙ 2 = = M 1 M2 2 (M 1 + M 2 ) 2 d ḋ + M 2M1 2 (M 1 + M 2 ) 2 d ḋ M 1 M 2 d M 1 + M ḋ 2 (1.48) So J = µh (1.49) where µ is the reduced mass, and h is the specific angular momentum. 42 / 59
55 Momentum loss due to mass loss 43 / 59
56 Momentum loss due to Gravitational Radiation 44 / 59
57 Momentum loss due to Gravitational Radiation Question: predict the evolution of the pulsar s orbit. 45 / 59
58 Momentum loss due to Gravitational Radiation Weisberg and Taylor / 59
59 Binary Super-massive Black holes 47 / 59
60 Orbital periods Example Remember the orbit equation? ( d 2 u 1 ) dφ 2 + u = f u h 2 u 2 (1.50) where u 1 r and f r = f for a spherical potential. For f from a gravitational potential, we have ( ) 1 f = dφ = u 2 dφ u dr du (1.51) since gravity is conservative. There are two types of orbit: Unbound: r, u 0 as φ φ Bound: r (and u) oscillate between finite limits. 48 / 59
61 Energy Orbital periods Example If we take (1.50) du dφ : du d 2 u dφ dφ 2 + u du dφ + u2 dφ du h 2 u 2 du dφ = 0 (1.52) [ d ( ) ] 2 1 du + 1 dφ 2 dφ 2 u2 + Φ h 2 = 0 (1.53) and integrating over φ we have 1 2 ( ) 2 du + 1 dφ 2 u2 + Φ h 2 = constant = E h 2 (1.54) 49 / 59
62 Orbital periods Example and using h = r 2 φ E = r 4 φ 2 ( ) 2 du dφ 2 r 2 φ2 + Φ(r) = r 4 ( ) 2 du dt 2 r 2 φ 2 + Φ(r) = r 4 ( ) 2 du 2 dr ṙ r 2 φ2 + Φ(r) E h 2 = 1 2 Energy ( ) 2 du dφ u2 + Φ h 2 = 1 2ṙ r 2 φ2 + Φ(r) (1.55) i.e. we can show that the constant E we introduced is the energy per unit mass. 50 / 59
63 Orbital periods Example Peri and Apo ( ) 2 E h 2 = 1 du 2 dφ u2 + Φ h 2 For bound orbits, the limiting values of u (or r) occur where du dφ = 0, i.e. where u 2 = from (1.54). This has two roots, u 1 = 1 r 1 and u 2 = 1 r 2 For r 1 < r 2, where r 1 is the pericentre, r 2 the apocentre 2E 2Φ(u) h 2 (1.56) this is not obvious, since Φ is not defined, but it can be proved - it is an Example! 51 / 59
64 Orbital periods Example Orbital periods Radial motion The radial period T r is defined as the time to go from r 2 r 1 r 2. Now take (1.55) and re-write: E = 1 2 ṙ r 2 φ 2 + Φ(r) ( ) 2 dr = 2(E Φ(r)) h2 dt r 2 (1.57) where we used h = r 2 φ to eliminate φ So dr dt = ± 2(E Φ(r)) h2 r 2 (1.58) (two signs - ṙ can be either > 0 or < 0, and ṙ = 0 at r 1 & r 2. Then T r = r2 dt r2 odt = 2 r 1 dr dr = 2 dr (1.59) r 1 2(E Φ(r)) h2 r 2 52 / 59
65 Orbital periods Example Orbital periods Radial motion The radial period T r is defined as the time to go from r 2 r 1 r 2. Now take (1.55) and re-write: E = 1 2 ṙ r 2 φ 2 + Φ(r) ( ) 2 dr = 2(E Φ(r)) h2 dt r 2 (1.57) where we used h = r 2 φ to eliminate φ So dr dt = ± 2(E Φ(r)) h2 r 2 (1.58) (two signs - ṙ can be either > 0 or < 0, and ṙ = 0 at r 1 & r 2. Then T r = r2 dt r2 odt = 2 r 1 dr dr = 2 dr (1.59) r 1 2(E Φ(r)) h2 r 2 52 / 59
66 Orbital periods Azimuthal motion Orbital periods Example If travelling from r 2 r 1 r 2 φ is increased by an amount so φ = odφ = 2 r2 φ = 2h r 1 r2 r 1 dφ dr dr = 2 dr r2 r 1 dφ dt dt dr (1.60) dr r 2 2(E Φ(r)) h2 r 2 (1.61) 53 / 59
67 Precession of the orbit Orbital periods Example For a given orbit, the time taken to go around once (i.e. 0 2π) depends in general on where you start, so the azimuthal period is not well defined. Instead use the mean angular velocity ω = φ/t r to obtain a mean azimuthal period T φ, so T φ = 2π/ ω T φ = 2π φ T r is the mean time to go around once. Note that unless φ/2π is a rational number the orbit is not closed. 54 / 59
68 Orbital periods Example Precession of the orbit For Keplerian orbit φ = 2π T r = T φ. In one period T r the apocentre (or pericentre) advances by an angle φ 2π. i.e.the orbit shifts in azimuth at an average rate given by the mean precession rate Ω p = φ 2π T r rad s 1 (1.62) Thus the precession period is T p = 2π = T r Ω φ p 2π 1 (1.63) This precession is in the sense opposite to the rotation of the star In the special case of a Keplerian orbit φ = 2π T φ = T r and Ω p = 0, i.e. orbits are closed and do not precess. Otherwise general orbit is a rosette between r 1 & r 2. This allows us to visualize how we can build a galaxy out of stars on different orbits. 55 / 59
69 Precession of the orbit Orbital periods Example 56 / 59
70 Orbital periods Example T r for the Keplerian case Φ(r) = GM r We have equation (1.59) r2 dr T r = 2 r 1 2(E Φ(r)) h2 r 2 Now r 1 & r 2 are determined from ṙ = 0, i.e. Example 2(E Φ(r)) h2 r 2 = 0 (1.64) 2E + 2GM r r 1 r 2 = h2 2E ; (remember E < 0 for a bound orbit). h2 r 2 = 0 (1.65) r 2 + GM E r h2 2E = 0 (1.66) (r r 1 )(r r 2 ) = 0 (1.67) r 1 + r 2 = GM E (1.68) 57 / 59
71 Orbital periods Example T r for the Keplerian case Φ(r) = GM r We have equation (1.59) r2 dr T r = 2 r 1 2(E Φ(r)) h2 r 2 Now r 1 & r 2 are determined from ṙ = 0, i.e. Example 2(E Φ(r)) h2 r 2 = 0 (1.64) 2E + 2GM r r 1 r 2 = h2 2E ; (remember E < 0 for a bound orbit). h2 r 2 = 0 (1.65) r 2 + GM E r h2 2E = 0 (1.66) (r r 1 )(r r 2 ) = 0 (1.67) r 1 + r 2 = GM E (1.68) 57 / 59
72 Orbital periods Example Rewrite (1.59) as r2 T r = 2 r 1 rdr 2E(r r1 )(r r 2 ) = 2 r2 2 E r 1 Example rdr (r2 r)(r r 1 ) (1.69) if r 1 < r < r 2. This is another of those integrals. If R = a + bx + cx 2 = r 2 + (r 1 + r 2 )r r 1 r 2 and = 4ac b 2 which becomes, using the variables here, = (r 1 r 2 ) 2 then xdx R = b ( ) 1 2cx + b sin 1 R c 2c c for c < 0 and < 0 (See G&R and 2.264). 58 / 59
73 Orbital periods Example The first term is 0 at r 1 and r 2 (R = 0 there), so [ ] [ ( 2 r1 + r T r = 2 sin 1 2r2 + r 1 + r 2 2 E 2 r 1 r 2 = = [ 2 r1 + r 2 2 E 2 2 GM [ π ] 2 E 2( E) 2 ( π 2 ) = 2πGM ( 2E) 3 2 ] [sin 1 (1) sin 1 ( 1)] Example ) ( )] sin 1 2r1 + r 1 + r 2 r 1 r 2 (1.70) 59 / 59
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 informationChapter 8. Orbits. 8.1 Conics
Chapter 8 Orbits 8.1 Conics Conic sections first studied in the abstract by the Greeks are the curves formed by the intersection of a plane with a cone. Ignoring degenerate cases (such as a point, or pairs
More informationUse conserved quantities to reduce number of variables and the equation of motion (EOM)
Physics 106a, Caltech 5 October, 018 Lecture 8: Central Forces Bound States Today we discuss the Kepler problem of the orbital motion of planets and other objects in the gravitational field of the sun.
More informationChapter 13. Gravitation
Chapter 13 Gravitation e = c/a A note about eccentricity For a circle c = 0 à e = 0 a Orbit Examples Mercury has the highest eccentricity of any planet (a) e Mercury = 0.21 Halley s comet has an orbit
More informationChapter 13. Universal Gravitation
Chapter 13 Universal Gravitation Planetary Motion A large amount of data had been collected by 1687. There was no clear understanding of the forces related to these motions. Isaac Newton provided the answer.
More informationCentral Force Problem
Central Force Problem Consider two bodies of masses, say earth and moon, m E and m M moving under the influence of mutual gravitational force of potential V(r). Now Langangian of the system is where, µ
More informationLecture 21 Gravitational and Central Forces
Lecture 21 Gravitational and Central Forces 21.1 Newton s Law of Universal Gravitation According to Newton s Law of Universal Graviation, the force on a particle i of mass m i exerted by a particle j of
More informationLecture XIX: Particle motion exterior to a spherical star
Lecture XIX: Particle motion exterior to a spherical star Christopher M. Hirata Caltech M/C 350-7, Pasadena CA 95, USA Dated: January 8, 0 I. OVERVIEW Our next objective is to consider the motion of test
More informationPotential/density pairs and Gauss s law
Potential/density pairs and Gauss s law We showed last time that the motion of a particle in a cluster will evolve gradually, on the relaxation time scale. This time, however, is much longer than the typical
More informationUniversal Gravitation
Universal Gravitation Newton s Law of Universal Gravitation Every particle in the Universe attracts every other particle with a force that is directly proportional to the product of their masses and inversely
More informationPhysics 351, Spring 2015, Homework #7. Due at start of class, Friday, March 3, 2017
Physics 351, Spring 2015, Homework #7. Due at start of class, Friday, March 3, 2017 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page at positron.hep.upenn.edu/q351
More informationPhysics 351, Spring 2015, Homework #7. Due at start of class, Friday, March 6, 2015
Physics 351, Spring 2015, Homework #7. Due at start of class, Friday, March 6, 2015 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page at positron.hep.upenn.edu/q351
More information( ) 2 1 r S. ( dr) 2 r 2 dφ
General relativity, 4 Orbital motion of small test masses The starting point for analyzing free fall trajectories in the (-space, 1-time) Schwarzschild spacetime is Equation (3) from GR 3: ( dτ ) = 1 r
More informationAST1100 Lecture Notes
AST1100 Lecture Notes 5 The virial theorem 1 The virial theorem We have seen that we can solve the equation of motion for the two-body problem analytically and thus obtain expressions describing the future
More informationF = ma. G mm r 2. S center
In the early 17 th century, Kepler discovered the following three laws of planetary motion: 1. The planets orbit around the sun in an ellipse with the sun at one focus. 2. As the planets orbit around the
More informationToday. Laws of Motion. Conservation Laws. Gravity. tides
Today Laws of Motion Conservation Laws Gravity tides Newton s Laws of Motion Our goals for learning: Newton s three laws of motion Universal Gravity How did Newton change our view of the universe? He realized
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 informationWelcome back to Physics 215
Welcome back to Physics 215 Today s agenda: More rolling without slipping Newtonian gravity Planetary orbits Gravitational Potential Energy Physics 215 Spring 2018 Lecture 13-1 1 Rolling without slipping
More informationBinary star formation
Binary star formation So far we have ignored binary stars. But, most stars are part of binary systems: Solar mass stars: about 2 / 3 are part of binaries Separations from: < 0.1 au > 10 3 au Wide range
More informationStellar Dynamics and Structure of Galaxies
Stellar Dynamics and Structure of Galaxies Gerry Gilmore H47 email: gil@ast.cam.ac.uk Lectures: Monday 12:10-13:00 Wednesday 11:15-12:05 Friday 12:10-13:00 Books: Binney & Tremaine Galactic Dynamics Princeton
More informationQuestion 1: Spherical Pendulum
Question 1: Spherical Pendulum Consider a two-dimensional pendulum of length l with mass M at its end. It is easiest to use spherical coordinates centered at the pivot since the magnitude of the position
More informationL03: Kepler problem & Hamiltonian dynamics
L03: Kepler problem & Hamiltonian dynamics 18.354 Ptolemy circa.85 (Egypt) -165 (Alexandria) Greek geocentric view of the universe Tycho Brahe 1546 (Denmark) - 1601 (Prague) "geo-heliocentric" system last
More information1 Summary of Chapter 2
General Astronomy (9:61) Fall 01 Lecture 7 Notes, September 10, 01 1 Summary of Chapter There are a number of items from Chapter that you should be sure to understand. 1.1 Terminology A number of technical
More informationEinstein s Equations. July 1, 2008
July 1, 2008 Newtonian Gravity I Poisson equation 2 U( x) = 4πGρ( x) U( x) = G d 3 x ρ( x) x x For a spherically symmetric mass distribution of radius R U(r) = 1 r U(r) = 1 r R 0 r 0 r 2 ρ(r )dr for r
More informationFROM NEWTON TO KEPLER. One simple derivation of Kepler s laws from Newton s ones.
italian journal of pure and applied mathematics n. 3 04 (393 400) 393 FROM NEWTON TO KEPLER. One simple derivation of Kepler s laws from Newton s ones. František Mošna Department of Mathematics Technical
More informationWelcome back to Physics 211. Physics 211 Spring 2014 Lecture Gravity
Welcome back to Physics 211 Today s agenda: Newtonian gravity Planetary orbits Gravitational Potential Energy Physics 211 Spring 2014 Lecture 14-1 1 Gravity Before 1687, large amount of data collected
More informationStellar Dynamics and Structure of Galaxies
Stellar Dynamics and Structure of Galaxies. Spherically symmetric objects Vasily Belokurov vasily@ast.cam.ac.uk Institute of Astronomy Lent Term 2016 1 / 21 Outline I 1 2 Globular of galaxies 2 / 21 Why
More informationSpecial Relativity: The laws of physics must be the same in all inertial reference frames.
Special Relativity: The laws of physics must be the same in all inertial reference frames. Inertial Reference Frame: One in which an object is observed to have zero acceleration when no forces act on it
More informationCentral 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
Central force motion/kepler problem This short note summarizes our discussion in the lectures of various aspects of the motion under central force, in particular, the Kepler problem of inverse square-law
More informationNewton s Laws of Motion and Gravity ASTR 2110 Sarazin. Space Shuttle
Newton s Laws of Motion and Gravity ASTR 2110 Sarazin Space Shuttle Discussion Session This Week Friday, September 8, 3-4 pm Shorter Discussion Session (end 3:40), followed by: Intro to Astronomy Department
More informationAstro Lecture 12. Energy and Gravity (Cont d) 13/02/09 Habbal Astro Lecture 12 1
Astro 110-01 Lecture 12 Energy and Gravity (Cont d) 13/02/09 Habbal Astro110-01 Lecture 12 1 Energy due to movement of Kinetic Energy: object E k = ½ m v 2 13/02/09 Habbal Astro110-01 Lecture 12 2 Gravitational
More informationIntroduction to Astronomy
Introduction to Astronomy AST0111-3 (Astronomía) Semester 2014B Prof. Thomas H. Puzia Newton s Laws Big Ball Fail Universal Law of Gravitation Every mass attracts every other mass through a force called
More informationPhysics 115/242 The Kepler Problem
Physics 115/242 The Kepler Problem Peter Young (Dated: April 21, 23) I. INTRODUCTION We consider motion of a planet around the sun, the Kepler problem, see e.g. Garcia, Numerical Methods for Physics, Sec.
More informationdistribution of mass! The rotation curve of the Galaxy ! Stellar relaxation time! Virial theorem! Differential rotation of the stars in the disk
Today in Astronomy 142:! The local standard of rest the Milky Way, continued! Rotation curves and the! Stellar relaxation time! Virial theorem! Differential rotation of the stars in the disk distribution
More informationLecture 15 - Orbit Problems
Lecture 15 - Orbit Problems A Puzzle... The ellipse shown below has one focus at the origin and its major axis lies along the x-axis. The ellipse has a semimajor axis of length a and a semi-minor axis
More information4.3 Conservation Laws in Astronomy
4.3 Conservation Laws in Astronomy Our goals for learning: Why do objects move at constant velocity if no force acts on them? What keeps a planet rotating and orbiting the Sun? Where do objects get their
More informationPotential-Density pair: Worked Problem In a spherical galaxy, the density of matter varies with radius as. a r 2 (r+a) 2, ρ(r) = M 4π
15 Potential-Density pair: Worked Problem In a spherical galaxy, the density of matter varies with radius as where M and a are constants. ρ(r) = M 4π a r 2 (r+a) 2, (a) Verify that the total mass of the
More informationSolution Set Two. 1 Problem #1: Projectile Motion Cartesian Coordinates Polar Coordinates... 3
: Solution Set Two Northwestern University, Classical Mechanics Classical Mechanics, Third Ed.- Goldstein October 7, 2015 Contents 1 Problem #1: Projectile Motion. 2 1.1 Cartesian Coordinates....................................
More informationLecture 9 Chapter 13 Gravitation. Gravitation
Lecture 9 Chapter 13 Gravitation Gravitation UNIVERSAL GRAVITATION For any two masses in the universe: F = Gm 1m 2 r 2 G = a constant evaluated by Henry Cavendish +F -F m 1 m 2 r Two people pass in a hall.
More informationProblem Set 4 Solutions
Problem Set 4 Solutions AY 7b Spring 2012 Problem 1 For the very simple model of the MW where Θ() = 200 km s 1, we know that where Ω = Θ/. Since Θ const, we can rewrite this as v r (l) = (Ω Ω 0 ) sin l
More informationOrbital Motion in Schwarzschild Geometry
Physics 4 Lecture 29 Orbital Motion in Schwarzschild Geometry Lecture 29 Physics 4 Classical Mechanics II November 9th, 2007 We have seen, through the study of the weak field solutions of Einstein s equation
More informationASTRO 310: Galac/c & Extragalac/c Astronomy Prof. Jeff Kenney
ASTRO 310: Galac/c & Extragalac/c Astronomy Prof. Jeff Kenney Class 9 September 26, 2018 Introduc/on to Stellar Dynamics: Poten/al Theory & Mass Distribu/ons & Mo/ons Gravity & Stellar Systems ~Only force
More informationGravitational Efects and the Motion of Stars
Gravitational Efects and the Motion of Stars On the largest scales (galaxy clusters and larger), strong evidence that the dark matter has to be non-baryonic: Abundances of light elements (hydrogen, helium
More informationTwo- and Three-Dimensional Motion (Symon Chapter Three)
Two- and Three-Dimensional Motion (Symon Chapter Three) Physics A300 Fall 2003 Contents 1 Vector Analysis 2 1.1 Definitions..................................... 2 1.1.1 Notation..................................
More informationAim: Understand equilibrium of galaxies
8. Galactic Dynamics Aim: Understand equilibrium of galaxies 1. What are the dominant forces? 2. Can we define some kind of equilibrium? 3. What are the relevant timescales? 4. Do galaxies evolve along
More informationCelestial Mechanics Lecture 10
Celestial Mechanics Lecture 10 ˆ This is the first of two topics which I have added to the curriculum for this term. ˆ We have a surprizing amount of firepower at our disposal to analyze some basic problems
More informationProblem Set 5 Solution
Problem Set 5 Solution Friday, 14 October 011 Physics 111 Proble The Setup For a pair of point masses, and m, interacting via a conservative force directed along the line joining the masses, show that
More informationDynamics of the solar system
Dynamics of the solar system Planets: Wanderer Through the Sky Planets: Wanderer Through the Sky Planets: Wanderer Through the Sky Planets: Wanderer Through the Sky Ecliptic The zodiac Geometry of the
More information!! r r θ! 2 r!φ 2 sinθ = GM / r 2 r!! θ + 2!r θ! r!φ 2 sinθ cosθ = 0 r!! φ sinθ + 2r θ!!φ cosθ + 2!r!φ sinθ = 0
Phys 325 Lecture 9 Tuesday 12 February, 2019 NB Midterm next Thursday 2/21. Sample exam will be posted soon Formula sheet on sample exam will be same as formula sheet on actual exam You will be permitted
More informationChapter 13: universal gravitation
Chapter 13: universal gravitation Newton s Law of Gravitation Weight Gravitational Potential Energy The Motion of Satellites Kepler s Laws and the Motion of Planets Spherical Mass Distributions Apparent
More informationWe have seen how to calculate forces and potentials from the smoothed density ρ.
9. Orbits in stationary Potentials We have seen how to calculate forces and potentials from the smoothed density ρ. Now we can analyse how stars move in this potential. Because two body interactions can
More informationTycho Brahe ( )
Tycho Brahe (1546-1601) Foremost astronomer after the death of Copernicus. King Frederick II of Denmark set him up at Uraniborg, an observatory on the island of Hveen. With new instruments (quadrant),
More informationPhys 7221, Fall 2006: Midterm exam
Phys 7221, Fall 2006: Midterm exam October 20, 2006 Problem 1 (40 pts) Consider a spherical pendulum, a mass m attached to a rod of length l, as a constrained system with r = l, as shown in the figure.
More informationPhysics 12. Unit 5 Circular Motion and Gravitation Part 2
Physics 12 Unit 5 Circular Motion and Gravitation Part 2 1. Newton s law of gravitation We have seen in Physics 11 that the force acting on an object due to gravity is given by a well known formula: F
More information4 Orbits in stationary Potentials (BT 3 to page 107)
-313see http://www.strw.leidenuniv.nl/ franx/college/ mf-sts13-c-1-313see http://www.strw.leidenuniv.nl/ franx/college/ mf-sts13-c Orbits in stationary Potentials (BT 3 to page 17) Now we have seen how
More informationPhys 7221, Fall 2006: Homework # 5
Phys 722, Fal006: Homewor # 5 Gabriela González October, 2006 Prob 3-: Collapse of an orbital system Consier two particles falling into each other ue to gravitational forces, starting from rest at a istance
More informationWelcome back to Physics 215
Welcome back to Physics 215 Today s agenda: Gravity 15-2 1 Current assignments HW#15 due Monday, 12/12 Final Exam, Thursday, Dec. 15 th, 3-5pm in 104N. Two sheets of handwritten notes and a calculator
More informationChapter 4 Making Sense of the Universe: Understanding Motion, Energy, and Gravity. Copyright 2009 Pearson Education, Inc.
Chapter 4 Making Sense of the Universe: Understanding Motion, Energy, and Gravity How do we describe motion? Precise definitions to describe motion: Speed: Rate at which object moves speed = distance time
More informationAy123 Set 1 solutions
Ay13 Set 1 solutions Mia de los Reyes October 18 1. The scale of the Sun a Using the angular radius of the Sun and the radiant flux received at the top of the Earth s atmosphere, calculate the effective
More informationSolutions to PS 2 Physics 201
Solutions to PS Physics 1 1. ke dq E = i (1) r = i = i k eλ = i k eλ = i k eλ k e λ xdx () (x x) (x x )dx (x x ) + x dx () (x x ) x ln + x x + x x (4) x + x ln + x (5) x + x To find the field for x, we
More informationr CM = ir im i i m i m i v i (2) P = i
Physics 121 Test 3 study guide Thisisintendedtobeastudyguideforyourthirdtest, whichcoverschapters 9, 10, 12, and 13. Note that chapter 10 was also covered in test 2 without section 10.7 (elastic collisions),
More informationCopyright 2009, August E. Evrard.
Unless otherwise noted, the content of this course material is licensed under a Creative Commons BY 3.0 License. http://creativecommons.org/licenses/by/3.0/ Copyright 2009, August E. Evrard. You assume
More informationOrbit Characteristics
Orbit Characteristics We have shown that the in the two body problem, the orbit of the satellite about the primary (or vice-versa) is a conic section, with the primary located at the focus of the conic
More informationPhysics 201, Lecture 23
Physics 201, Lecture 23 Today s Topics n Universal Gravitation (Chapter 13) n Review: Newton s Law of Universal Gravitation n Properties of Gravitational Field (13.4) n Gravitational Potential Energy (13.5)
More informationMidterm 3 Thursday April 13th
Welcome back to Physics 215 Today s agenda: rolling friction & review Newtonian gravity Planetary orbits Gravitational Potential Energy Physics 215 Spring 2017 Lecture 13-1 1 Midterm 3 Thursday April 13th
More informationVISUAL PHYSICS ONLINE
VISUAL PHYSICS ONLINE EXCEL SIMULATION MOTION OF SATELLITES DOWNLOAD the MS EXCEL program PA50satellite.xlsx and view the worksheet Display as shown in the figure below. One of the most important questions
More informationAST111 PROBLEM SET 2 SOLUTIONS. RA=02h23m35.65s, DEC=+25d18m42.3s (Epoch J2000).
AST111 PROBLEM SET 2 SOLUTIONS Home work problems 1. Angles on the sky and asteroid motion An asteroid is observed at two different times. The asteroid is located at RA=02h23m35.65s, DEC=+25d18m42.3s (Epoch
More informationKepler, Newton, and laws of motion
Kepler, Newton, and laws of motion First: A Little History Geocentric vs. heliocentric model for solar system (sec. 2.2-2.4)! The only history in this course is this progression: Aristotle (~350 BC) Ptolemy
More informationToday. Events. Energy. Gravity. Homework Due Next time. Practice Exam posted
Today Energy Gravity Events Homework Due Next time Practice Exam posted Autumn is here! Autumnal equinox occurred at 11:09pm last night night and day very nearly equal today days getting shorter Moon is
More informationThe first term involves the cross product of two parallel vectors and so it vanishes. We then get
Physics 3550 Angular Momentum. Relevant Sections in Text: 3.4, 3.5 Angular Momentum You have certainly encountered angular momentum in a previous class. The importance of angular momentum lies principally
More informationAstro 210 Lecture 8 Feb 4, 2011
Astro 210 Lecture 8 Feb 4, 2011 Announcements HW2 due apologies for the erratum HW3 available, due next Friday HW1 Q8 bonus still available register your iclicker; link on course webpage Planetarium: shows
More information5. Universal Laws of Motion
5. Universal Laws of Motion If I have seen farther than others, it is because I have stood on the shoulders of giants. Sir Isaac Newton (164 177) Physicist Image courtesy of NASA/JPL Sir Isaac Newton (164-177)
More information5.1. Accelerated Coordinate Systems:
5.1. Accelerated Coordinate Systems: Recall: Uniformly moving reference frames (e.g. those considered at 'rest' or moving with constant velocity in a straight line) are called inertial reference frames.
More informationGravity and Orbits. Objectives. Clarify a number of basic concepts. Gravity
Gravity and Orbits Objectives Clarify a number of basic concepts Speed vs. velocity Acceleration, and its relation to force Momentum and angular momentum Gravity Understand its basic workings Understand
More informationEinstein s Theory of Gravity. June 10, 2009
June 10, 2009 Newtonian Gravity Poisson equation 2 U( x) = 4πGρ( x) U( x) = G d 3 x ρ( x) x x For a spherically symmetric mass distribution of radius R U(r) = 1 r U(r) = 1 r R 0 r 0 r 2 ρ(r )dr for r >
More informationANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS
ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS A1.1. Kepler s laws Johannes Kepler (1571-1630) discovered the laws of orbital motion, now called Kepler's laws.
More informationLesson 9. Luis Anchordoqui. Physics 168. Tuesday, October 24, 17
Lesson 9 Physics 168 1 Static Equilibrium 2 Conditions for Equilibrium An object with forces acting on it but that is not moving is said to be in equilibrium 3 Conditions for Equilibrium (cont d) First
More informationISIMA lectures on celestial mechanics. 1
ISIMA lectures on celestial mechanics. 1 Scott Tremaine, Institute for Advanced Study July 2014 The roots of solar system dynamics can be traced to two fundamental discoveries by Isaac Newton: first, that
More informationChapter 2 Introduction to Binary Systems
Chapter 2 Introduction to Binary Systems In order to model stars, we must first have a knowledge of their physical properties. In this chapter, we describe how we know the stellar properties that stellar
More informationGravitation. Luis Anchordoqui
Gravitation Kepler's law and Newton's Synthesis The nighttime sky with its myriad stars and shinning planets has always fascinated people on Earth. Towards the end of the XVI century the astronomer Tycho
More informationA Guide to the Next Few Lectures!
Dynamics and how to use the orbits of stars to do interesting things chapter 3 of S+G- parts of Ch 11 of MWB (Mo, van den Bosch, White) READ S&G Ch 3 sec 3.1, 3.2, 3.4 we are skipping over epicycles 1
More informationChapter 1. Basic Concepts. 1.1 Trajectories
Chapter 1 Basic Concepts 1.1 Trajectories We shall be concerned in this course with the motion of particles. Larger bodies will (with a few exceptions) be made up of collections of particles. We will find
More informationA Guide to the Next Few Lectures!
Dynamics and how to use the orbits of stars to do interesting things chapter 3 of S+G- parts of Ch 11 of MWB (Mo, van den Bosch, White) READ S&G Ch 3 sec 3.1, 3.2, 3.4 we are skipping over epicycles 1
More informationChapter 9. Gravitation
Chapter 9 Gravitation 9.1 The Gravitational Force For two particles that have masses m 1 and m 2 and are separated by a distance r, the force has a magnitude given by the same magnitude of force acts on
More informationLecture 22: Gravitational Orbits
Lecture : Gravitational Orbits Astronomers were observing the motion of planets long before Newton s time Some even developed heliocentric models, in which the planets moved around the sun Analysis of
More information9/13/ Describing Motion: Examples from Everyday Life. Chapter 4: Making Sense of the Universe Understanding Motion, Energy, and Gravity
9/13/17 Lecture Outline 4.1 Describing Motion: Examples from Everyday Life Chapter 4: Making Sense of the Universe Understanding Motion, Energy, and Gravity Our goals for learning: How do we describe motion?
More informationAST111, Lecture 1b. Measurements of bodies in the solar system (overview continued) Orbital elements
AST111, Lecture 1b Measurements of bodies in the solar system (overview continued) Orbital elements Planetary properties (continued): Measuring Mass The orbital period of a moon about a planet depends
More informationASTR3007/4007/6007, Class 2: Stellar Masses; the Virial Theorem
ASTR7/47/67, Class : Stellar Masses; the Virial Theorem 4 February In the first class we discussed stars light output and the things than can be derived directly from it luminosity, temperature, and composition.
More informationPost-Keplerian effects in binary systems
Post-Keplerian effects in binary systems Laboratoire Univers et Théories Observatoire de Paris / CNRS The problem of binary pulsar timing (Credit: N. Wex) Some classical tests of General Relativity Gravitational
More informationBlack Holes. Jan Gutowski. King s College London
Black Holes Jan Gutowski King s College London A Very Brief History John Michell and Pierre Simon de Laplace calculated (1784, 1796) that light emitted radially from a sphere of radius R and mass M would
More informationHOMEWORK 10. Applications: special relativity, Newtonian limit, gravitational waves, gravitational lensing, cosmology, 1 black holes
General Relativity 8.96 (Petters, spring 003) HOMEWORK 10. Applications: special relativity, Newtonian limit, gravitational waves, gravitational lensing, cosmology, 1 black holes 1. Special Relativity
More informationF 12. = G m m 1 2 F 21. = G m 1m 2 = F 12. Review: Newton s Law Of Universal Gravitation. Physics 201, Lecture 23. g As Function of Height
Physics 01, Lecture Today s Topics n Universal Gravitation (Chapter 1 n Review: Newton s Law of Universal Gravitation n Properties of Gravitational Field (1.4 n Gravitational Potential Energy (1.5 n Escape
More informationObservational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws
Observational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws Craig Lage New York University - Department of Physics craig.lage@nyu.edu February 24, 2014 1 / 21 Tycho Brahe s Equatorial
More information(b) The period T and the angular frequency ω of uniform rotation are related to the cyclic frequency f as. , ω = 2πf =
PHY 302 K. Solutions for problem set #9. Non-textbook problem #1: (a) Rotation frequency of 1 Hz means one revolution per second, or 60 revolutions per minute (RPM). The pre-lp vinyl disks rotated at 78
More informationThe Restricted 3-Body Problem
The Restricted 3-Body Problem John Bremseth and John Grasel 12/10/2010 Abstract Though the 3-body problem is difficult to solve, it can be modeled if one mass is so small that its effect on the other two
More informationNotes on Planetary Motion
(1) Te motion is planar Notes on Planetary Motion Use 3-dimensional coordinates wit te sun at te origin. Since F = ma and te gravitational pull is in towards te sun, te acceleration A is parallel to te
More informationAy 1 Lecture 2. Starting the Exploration
Ay 1 Lecture 2 Starting the Exploration 2.1 Distances and Scales Some Commonly Used Units Distance: Astronomical unit: the distance from the Earth to the Sun, 1 au = 1.496 10 13 cm ~ 1.5 10 13 cm Light
More informationCelestial Orbits. Adrienne Carter Ottopaskal Rice May 18, 2001
Celestial Orbits Adrienne Carter sillyajc@yahoo.com Ottopaskal Rice ottomanbuski@hotmail.com May 18, 2001 1. Tycho Brache, a Danish astronomer of the late 1500s, had collected large amounts of raw data
More informationNewton s Laws and the Nature of Matter
Newton s Laws and the Nature of Matter The Nature of Matter Democritus (c. 470-380 BCE) posited that matter was composed of atoms Atoms: particles that can not be further subdivided 4 kinds of atoms: earth,
More information