4.2 Формула Эддингтона

Для заданного сферически-симметричного потенциала (или плотности ν(r)) можно попытаться найти соответствующую эргодическую ФР f(E):

ν(r) =f(E)d3v = 4πf(E)v2dv,v = |v|. (4.29)

Если такая f(E) существует, то она единственна. Но не для каждого распределения плотности можно найти f(E).

Энергию частицы (которая в случае ограниченной системы всегда меньше нуля) и потенциал удобно заменить на относительныю энергию и потенциал. Относительный потенциал определяется следующим образом. Пусть Φ – обычный потенциал, как правило определенный так, что Φ(r) 0 при r , и a – радиус системы (так что за пределы этого радиуса частицы не вылетают). Тогда относительный потенциал определяется как

Ψ(r) = Φ(a) Φ(r). (4.30)

Он оказывается положительно определенным, как и относительная энергия

= Ψ(r) v2 2 . (4.31)

Такое переобозначение естественно и удобно, поскольку функции распределения обращаются в нуль при < 0 (В разделе 3.1.1 мы ввели переменную J, аналогичную 2.) Итак,

ν(r) = 4π02Ψ(r)f()v2dv, (4.32)

или поскольку 2v2dv = vdv2 = 22(Ψ )d, получим

ν(r) 8π = 20Ψf()Ψ d (4.33)

Относительный потенциал Ψ = Ψ(r) – монотонно убывающая функция радиуса. Поэтому можно обратить эту зависимость, т.е. получить r = r(Ψ) и затем рассматривать Ψ как новую независимую переменную (например ν = ν(Ψ)). Дифференцируя по Ψ, найдем

1 8π dν dΨ =0Ψf() 1 Ψ d (4.34)

– уравнение Абеля, решение которого известно:

f() = 1 8π2 d d0 dΨ Ψ dν dΨ. (4.35)

По определению f – вероятность нахождения частицы в точке (x,v), т.е. f() должна быть неотрицательной. Отсюда

0 dΨ Ψ dν dΨ (4.36)

должна быть неубывающей функцией. Очевидно, не для всякой ν это условие выполнено, поэтому не любая модель масс имеет эргодическую ФР. Для моделей Пламмера, изохронной, Hernquist, Jaffe аналитические выражения f() известны.