ast 1420 galactic structure and dynamics today disks
play

AST 1420 Galactic Structure and Dynamics Today: disks! NGC 5907 - PowerPoint PPT Presentation

AST 1420 Galactic Structure and Dynamics Today: disks! NGC 5907 M31 Today: disks! Outline Simple gravitational potentials for disks Realistic potentials for disks Orbits in axisymmetric potentials Close-to-circular orbits


  1. AST 1420 Galactic Structure and Dynamics

  2. Today: disks! NGC 5907 M31

  3. Today: disks!

  4. Outline • Simple gravitational potentials for disks • Realistic potentials for disks • Orbits in axisymmetric potentials • Close-to-circular orbits (e.g., stars near the Sun)

  5. Computing the gravitational potential for a disk • Simple: just solve the Poisson equation! • Newton’s theorems don’t apply

  6. Razor-thin disks • ‘Extreme’ assumption: disk is infinitely thin • Integrate Poisson equation • Divergence theorem • For small volume

  7. Example: Kuzmin disk • At any z > 0: point mass at (0,-a) • At any z < 0: point mass at (0,a) • For any z =/= 0 —> density = 0 —> razor-thin • At z=0 ~ Plummer —> non-zero density

  8. Example: Kuzmin disk

  9. Thick disks: Miyamoto- Nagai • Razor-thin is not realistic for orbits in disks • Can thicken Kuzmin by replacing |z| —> √ [z 2 + b 2 ] • Miyamoto-Nagai: • b —> 0: Kuzmin • a —> 0: Plummer (spherical)

  10. Thick disks: Miyamoto- Nagai

  11. Thick disks: Miyamoto- Nagai Potential is much less flattened:

  12. Thick disks: Miyamoto- Nagai • Advantages: • Analytical formula —> fast for, e.g., orbit integration • Can vary b/a to get different ‘thicknesses’ • Disadvantages: • Large R ~ R -3 —> not like a realistic disk, too much density at large R • Similarly, vertical profile not realistic

  13. Miyamoto-Nagai rotation curve

  14. Flattened + flat rotation curve • Spherical model: • Has v c (r) = constant • General flattening strategy: r —> m = √ [R 2 + z 2 /q 2 ] • Has v c (R) = constant at z=0

  15. Flattened logarithmic potential • Potential for q=0.7 is slightly flattened:

  16. Flattened logarithmic potential • But density becomes negative around R=0!

  17. Flattened logarithmic potential: density • When q < 1/ √ 2: density is negative at R=0 • and everywhere between • Thus, you cannot flatten a logarithmic potential too much without creating negative densities • General problem when flattening a potential using r —> m = √ [R 2 + z 2 /q 2 ] [e.g., often used to get flattened DM halos]

  18. Potentials for disk densities

  19. Poisson equation for disks • Write density as • Poisson equation best written in cylindrical coordinates • Take Fourier transform of phi part

  20. Poisson equation for disks • or … • Similar to the Bessel differential equation:

  21. Poisson equation for disks • Solutions of the form

  22. Poisson equation for axisymmetric razor-thin disks • Solution: • With 𝛵

  23. Poisson equation for axisymmetric razor-thin disks • If we can decompose • Then the potential is • Decomposition from the Fourier-Bessel theorem

  24. Circular velocity for axisymmetric razor-thin disks

  25. Example: Mestel disk • Razor-thin disk with surface density • Hankel transform: • Potential:

  26. Example: Mestel disk • Enclosed mass: • Same as for spherical potential! Coincidence!

  27. Example: Exponential disk • Observed disks have exponential(-ish) surface density profiles: • Hankel transform: • Potential:

  28. Example: Exponential disk

  29. Example: Exponential disk

  30. Potentials with finite thickness • In general difficult, but galaxies overall well approximated as • Then each layer gives rise to the ~same potential, but shifted in z • Full potential adds up contribution from each layer

  31. Example: 
 Double- exponential disk • Density for which both radial and vertical profile are exponential • Full potential + forces: one-dimensional numerical integrals

  32. Double- exponential disk: rotation curve • Rotation curve becomes:

  33. Double- exponential disk: rotation curve

  34. Miyamoto-Nagai rotation curve

  35. Orbits in axisymmetric disks

  36. Orbits in axisymmetric disks • Approximate model for disk galaxy: • Flattened axisymmetric disk • Symmetric around z=0 • Often use Miyamoto-Nagai for computational convenience • E.g., galpy ’s Milky-Way model • Miyamoto-Nagai disk with scale length 3 kpc, scale height 280 pc • NFW halo • Spherical bulge with exponential cut-off

  37. Orbits in cylindrical geometry • Lagrangian in cylindrical coordinates • With conjugate momenta • Hamiltonian

  38. Orbits in cylindrical geometry • z component of the angular momentum is conserved • Hamiltonian becomes • with • Effectively a two (four) dimensional system in (R,z) —> meridional plane

  39. Motion in the meridional plane • Equations of motion • Coupled oscillations in R and z • No analytical solutions

  40. Motion in the meridional plane: effective potential

  41. http://astro.utoronto.ca/~bovy/AST1420/orbits/lec5-orbitexample1.html

  42. Motion in the meridional plane

  43. http://astro.utoronto.ca/~bovy/AST1420/orbits/lec5-orbitexample2.html

  44. Motion in the meridional plane

  45. Motion in the meridional plane

  46. Isolating integrals of the motion • Time-independent potential: E conserved —> motion restrained to ϕ < E • Spherical potential: L conserved —> motion restrained to (a) orbital plane (b) vT = | L |/R ( ϕ eff < E) • Motion fills rest of the phase-space

  47. http://astro.utoronto.ca/~bovy/AST1420/orbits/lec5-orbitexample3.html

  48. “Third integral” • Fact that orbits in axisymmetric do not fully explore the area in the meridional plane alllowed by ϕ eff < E means there has to be an additional integral • No known classical integral (like E or L z ) —> non-classical ‘third’ integral • Will come back to this in a few lectures • Focus here on approximate understanding

  49. Separability of disk orbits • Orbits in disks are to a good approximation independent oscillations in R and z (coupling is small) • Will be highly useful to understand equilibrium models of disks • Write the potential as

  50. Motion in the meridional plane

  51. Separability of disk orbits

  52. Separability of disk orbits

  53. Separability of disk orbits

  54. Separability of disk orbits

  55. Motion in the meridional plane

  56. Separability of disk orbits

  57. more typical orbit…

  58. Separability of disk orbits

  59. Separability of disk orbits

  60. Separability of disk orbits

  61. Separability of disk orbits • Orbits in disks are to a good approximation independent oscillations in R and z (coupling is small) • Will be highly useful to understand equilibrium models of disks • Write the potential as • Hamiltonian then splits into two pieces, with separately conserved energies

  62. Epicycle approximation for close-to-circular orbits

  63. Potential for close to circular orbits • Circular orbit is the minimum of the effective potential, can Taylor expand around the minimum • Motion is then explicitly two decoupled oscillators

  64. Frequencies • Frequencies of the oscillations: • +azimuthal Ω = v c (R) / R • Range:

  65. Vertical motion • Solution is sinusoidal • But bad approximation because disks are thin • Vertical decoupling is useful, vertical epicycles not so much

  66. Radial motion • Radial oscillation around guiding-center radius : radius of circular orbit with angular momentum • Azimuthal motion from conservation of angular momentum • Subtracting out motion of guiding center, motion is ellipse: epicycle • Axis ratio

Recommend


More recommend