**Information**

The methods used in stellar dynamics are adopted from theoretical physics. The student is expected to have a fundamental understanding of statistical physics, (vector) calculus, electrodynamics, classical mechanics and should be able to program a computer

Teacher: Simon Portegies Zwart

Assistant: Moein Mosleh

**Time & Place:**

Each friday from 13:45 to 15:30 in HL427

**Exam**

Example Exam

- The example exam can be found here

- The REAL exam can be found here

#### Lecture information

**Lecture 1**

Subjects:

- Gravitational Force

- The 3 laws of Newton

-Energy conservation

Assignment:

- At which mass does gravity become important. (How many protons are required before they become self gravitating.)

Paper:

A family of models for Spherical Stellar Systems

**Lecture 2**

Subjects:

- Potential Theory

- Newtons first and second theorm

Assignment:

- Given p(r) = R / (r^2(r+R^2)) where R = constant

- Whhat is normalisation p(r)
- What is M(r)
- What is phi(r)
- What can say about circular speed at:
- r << R (deep in the potential well)
- r >> R (far away from the center)

Paper:

Standardised Units and Time Scales (Heggie et. al.)

**Lecture 3**

Subjects:

- What to do with Potentials

- Keplers problem (Orbits)

- Equation of motion (in circular coordinates)

- Angular Momentum

Assignment:

- Given F = -(GM / r^{2}) , r_{dot} = phi_{dot}(dr/dphi) , (1/r) = C cos (phi - phi_{0}) + GM/L^{2}

- Demonstrate that: sin (phi - phi
_{0}) = (Lr_{dot}) / eGM - Give an expression for V∞ (t -> ∞ ) = V∞(e, L)
- Give an expression for lcord = a(1-e
^{2})- Is it true?
- What is special about it?

Paper:

Ejection of Hypervelocity Stars By The (Binary) Block Hole In The Galactic Center (section 1)

**Lecture 4**

Subjects:

- Motion of the center of mass

- Deflection Angle

- Encounter time

- Coulomb logarithm

Paper:

Discovery of an Unbound Hypervelocity Star in the Milky Way Halo

Assignment:

Read the following paper:

Mass Segregation, Relaxtion And the Coulomb Logarithm in N-body Systems

- Now measure the Coulomb logarithm in a computer simulation.

Make use of a N-body code, examples can be found, for example, on www.nbabel.org

Input files can be found here: http://nbabel.org/files/input/SD/SD.zip (Note that these files are in the format as descripbed on nbabel.org )

**Lecture 5**

Subjects:

- Dynamical friction

Paper:

Gravitational focusing and shielding of meteoroid streams

**Lecture 6**

Subjects:

- Distribution functions

- Continuity function

- Collisionless Boltzman Equation

Paper:

Self-consistent determinations of the total amount of matter near the sun

**Reading material for lecture 6 and 7**

(Read in the order as written down here)

B&T first edition:

- 8.2.1
- 1.E.a
- 4.1
- 4.2
- 2.2.3
- 4.2.1b
- 4.2.1d

B&T second edition:

- 7.2.2
- 7.2.3
- 7.2.4
- F.1.1
- F.1.2
- 4.1
- 4.1.2
- 2.3.3

**Lecture 7**

Subjects:

- Collisionless Boltzman Equation

- Jeans equation

- Distribution functions

Assignments: Error propegation in the Oort limit

- Rewrite the Jeans equations for a Cylindrical symmetric system

- Use the Poisson equations for a flattend stellar system to eliminate the potentioal term in the Jeans equations

- Calculate the Oort limit by introducing an error ' y' in the observed distance to a tracer populations

- How does the energy propagate in the calculation of the Oort limit?

Paper:

The kinematics and dynamics of the galactic globular cluster system

**Lecture 8**

Subjects:

- Sphericla Stellar Systems

- Plummer

- Distribution functions

Assignments: None

Paper: None

**Lecture 9**

Subjects:

Assignments:

Stability of Lagrange points and core collapse

In the original 1987 edition of BT, chapt. 2 (p. 140), the authors

proved that the Lagrange points L4 and L5 of the rotating logarithmic

potential are always stable. Accoeding to Pfenniger (1990, A&A 230,

55) this proof was incorrect according to was modified in the new

edition of BT.

1. What was wrong with the original proof?

2. Settle any doubts we still might have by numerically solving an

initial-value problem of Newton's equation of motion in the

corotating frame of reference for the appropriately chosen

potential and rotation speed to demonstrate the instability around

L4 and L5. Use one of the N-body codes from nbabel.org to settle

the issue. Explain how you select initial conditions and

demonstrate what orbits are stable, and when they become unstable.

Start with the stability around L1, L2 and L3, and then construct

initial conditions for L4 and L5.

3. Integrate a 1000 particle Plummer sphere using one of the Nbabel

N-body codes until core collapse. Does the code crash?

4. Reproducing fig.8.1 (BT87) using this simulation.

Explain what happens.

Paper:

**Lecture 10**

Subjects:

- Binaries

Papers:

The evolution of a primordial binary population in a globular cluster

Monte Carlo Simulations of Globular Cluster Evolution. III. Primordial Binary Interactions

**Added on 10 -december.**

Page 140 of the old book. Used for assignment of lecture 9. Click here for the page