5.1 Эпициклическое приближение

Рассмотрим в первом приближении характер движения звезд по почти круговым орбитам.

Радиус круговой орбиты в экваториальной плоскости Rg определяется минимумом эффективного потенциала Φeff = Φ(R,z) + Lz22R2, т.е.

Φ(R,z) R (Rg,0) = Lz2 Rg3 = Ω2(R g)Rg. (5.2)

Ω(R) – угловая частота вращения частицы на круговой орбите.

Пусть x = R Rg – отклонение радиуса звезды от радиуса круговой орбиты Rg, отвечающей заданному значению углового момента Lz (Rg = Rg(Lz)). Разложим эффективный потенциал Φeff в ряд Тейлора:

Φeff(Rg,0) + 1 22Φeff(R,z) R2 (Rg,0)x2 + 1 22Φeff(R,z) z2 (Rg,0)z2 + 𝒪(xz2) (5.3)

(член, пропорциональный xz равен нулю ввиду предполагаемой симметрии потенциала относительно плоскости z = 0). Определим частоту радиальных (или эпициклических) колебаний:

κ2(R g) 2Φeff(R,z) R2 (Rg,0)(Rg,0), (5.4)

где

κ2 = 4Ω2 + RdΩ2 dR R=Rg. (5.5)

Частота вертикальных колебаний, соответственно, есть:

ν2(R g) = 2Φeff(R,z) z2 (Rg,0). (5.6)

Звезды совершают гармонические колебания:

= κ2x,z̈ = ν2z (5.7)

В центре галактики угловая частота вращения Ω примерно постоянна (твердотельное вращение), поэтому κ 2Ω. Далее Ω убывает, но обычно не быстрее, чем в кеплеровском потенциале ( Ω R32), который дает κ = Ω. Поэтому

Ω κ 2Ω (5.8)

В окрестности Солнца, согласно наблюдениям κΩ = 1.35 ± 0.05

Уравнение (5.7) для движения перпендикулярно плоскости справедливо для достаточно малых z, при которых ρ = const, т.е. для z 300 пк. Однако, поскольку большинство звезд поднимается над экваториальной плоскостью выше 300 пк, эпициклическое приближение неприменимо к движению поперек диска галактики.

Рассмотрим движение звезд параллельно экваториальной плоскости. Примем решение уравнения (5.7) в виде x = Xcos(κt + α). Тогда для движения по азимуту имеем:

φ̇ = Lz R2 = Lz Rg21 + x Rg2 Ω g1 2x Rg (5.9)

Интегрирование по времени дает

φ = φ0 + Ω(Rg)t γ X Rg sin(κt + α),γ 2Ω κ (5.10)

Траектория движение звезды представляет собой эллипс с полуосями (X,Y ), центр которого вращается по окружности радиусом Rg с угловой частотой Ω(Rg). Соотношение полуосей эллипса YX = γ 1. Вращение звезды по эллипсу происходит в сторону, противоположную движению звезды вокруг центра.