Зовнішній розв`язок Шварцшильда
Знайдемо розв`язок
релятивістського рівняння гравітаційного поля у разі, коли джерелом поля є
центрально симетрична статична маса. Розв`язок
шукатимемо за межами області простору, зайнятою цією масою. Оскільки компоненти
тензора енергії-імпульсу у цьому разі дорівнюють нулю, то рівняння поля має
вигляд
Rij = 0.
Розв`язок цього рівняння шукатимемо у вигляді
квадрату інтервалу між двома нескінченно близькими подіями – точками чотиривимірного
простору, що мають три просторові координати і час. В якості першої специфічної
межової умови оберемо таку: на дуже великих відстанях від джерела поля, де гравітаційне
поле практично відсутнє, квадрат цього інтервалу збігається з квадратом
інтервалу у релятивістській механіці у відсутності гравітаційного поля
(ds)2 = - (dr)2 – r2
(dθ)2 - r2 sin2(θ) (dφ)2 + (dt)2.
Ця умова означає, що такою ж має залишатись структура
інтервалу і на будь-якій відстані від джерела поля. Умова сферичної симетрії
задачі, пов`язана із сферичною симетричністю джерела поля і його нерухомістю,
дозволяє стверджувати, що гравітаційне поле не випливатиме на кутові доданки
інтервалу. Інші два доданки – перший і четвертий матимуть додаткові змінні
множники, що прямуватимуть до 1 при нескінченному віддалені від джерела поля.
Такі міркування дозволяють записати інтервал у присутності поля у вигляді
(ds)2 = - exp(α) (dr)2 – r2 (dθ)2 - r2 sin2(θ) (dφ)2 + exp(β) (dt)2.
Тут α = α(r), β = β(r), що знову ж таки зумовлене
статичністю поля і його сферичною симетричністю. Якщо ввести стандартні
позначення релятивістської теорії гравітації: r=x1, θ=x2, φ=x3, t=x4, а також використати запис інтервали через
метричний тензор, то квадрат інтервалу можна записати так
(ds)2 = g11 (dx1)2
+ g22
(dx2)2
+ g33
(dx3)2
+ g44
(dx4)2.
Співставляючи два вирази для інтервалу легко бачити, що коваріантні
компоненти визначаються так
g11
= - exp(α);
g22
= - r2; g33 = - r2 sin2(θ);
g44
= exp(β);
gij
= 0, i≠j.
Алгоритм знаходження контраваріантних компонент
метричного тензора – це алгоритм знаходження оберненої матриці. Для його
реалізації нам потрібний визначник метричного тензора, який для діагональної
матриці дорівнює
g = g11 g22 g33 g44
= - r4 sin2(θ) exp(α + β).
Відповідні алгебраїчні
доповнення матимуть вигляд
A11 = g22 g33 g44; A22 = g11 g33
g44; A33 = g11 g22 g44; A11 = g11 g22
g33.
Контраваріантні компоненти метричного тензора тепер
будуть такими
g11 = A11 / g = - exp(-α); g22 = A22 / g = - 1 / r2;
g33 = A33 / g = - 1 / [r2
sin2(θ); g44 = A44
/ g = exp(-β).
Вирази для ненульових символів Кристофеля будуть
такими
Γijj
= - gii
(∂gjj/∂xi) / 2; Γiij = - gii (∂gii/∂xj) / 2; Γiii = - gii (∂gii/∂xi) / 2.
Після підстановки сюди виразів для компонент
метричного тензора, матимемо
Γ111
= α`/ 2; Γ122
= - r
exp(-
α);
Γ133
= - r
sin2(θ) exp(- α); Γ144
= β exp(β
- α) / 2;
Γ233 =
- sin(θ) cos(θ);
Γ221
=
1 / r; Γ331
= 1 / r; Γ441=
β`/ 2; Γ332=
cot(θ).
Тепер діагональні
компоненти тензора Річі визначатимуться так
R11 = β``/ 2 – α` β`/ 4 + (β`)2/
4 – α` / r;
R22 = exp(- α) [1 + r (β` - α`) / 2] – 1;
R33 = exp(- α) sin2(θ) [1 + r (β` - α`) / 2] – sin2(θ);
R44 = exp(β - α) [- β``/ 2 + α` β`/ 4 - (β`)2/ 4 – β`/ r].
Прирівнюючи ці компоненти нулю, отримуємо чотири рівняння
щодо невідомих функцій α і β. Видно, що
друге і третє рівняння збігаються, отже достатньо записати лише три наступні
рівняння
β``/ 2 – α` β`/ 4 + (β`)2/ 4 – α` / r = 0;
exp(- α) [1 + r (β` -
α`)
/ 2] – 1 = 0;
-
β``/
2 + α`
β`/
4 - (β`)2/ 4 – β`/ r = 0.
Лише два з цих трьох рівнянь незалежні. В якості цих двох
рівнянь доцільно взяти суму першого і третього і друге рівняння, яке можна
спростити використанням першого рівняння
α` + β` = 0;
exp(- α) (1 - r α`) – 1 = 0.
Після інтегрування
першого рівняння матимемо наступний загальний інтеграл
α + β = С,
де С – довільна стала. Її можна знайти з межової умови на
нескінченості, тобто з умови того, що на нескінченості квадрат інтервалу з
урахуванням гравітаційного поля переходить у квадрат інтервалу при його
відсутності (α=0, β=0). Отже, ця стала дорівнює нулю і
α + β = 0.
Друге рівняння системи
можна записати у такому спрощеному вигляді
d [r exp(- α)] / dr = 1.
Це рівняння легко
інтегрується і його загальний інтеграл має вигляд
r exp(- α) = r + D,
де D – нова довільна стала. Її можна знайти ще з
однієї межової умови, яка полягає у тому, що на достатньо великій, але
скінченій, відстані від джерела поля гравітаційαне поле має визначатись відомим
з теорії Ньютона класичним виразом. Ця умова має вигляд
g44 / 2
= - m / r + const.
Два отримані нами загальні інтеграли дозволяють записати
наступні вирази для відповідних компонент метричного тензора
g11 = - exp(α) = - 1 / (1 + D / r);
g44 = exp(β) = 1 + D / r.
Межова умова тепер має вигляд
g44
/ 2 = 1 + D / r = - m / r + const.
Звідси
D=-2m. Остаточно вираз для квадрату інтервалу
тепер буде таким
(ds)2 = - (1 – 2 m / r)-1 (dr)2 – r2 (dθ)2 - r2 sin2(θ) (dφ)2 + (1 – 2 m / r) (dt)2.
Це і є зовнішній розв`язок Шварцшильда. Замість одного Ньютонівського потенціалу гравітаційного поля
U(r) = m / r + const
ми маємо чотири ненульові компоненти метричного тензора, які відіграють роль гравітаційних потенціалів,
g11 = - 1 / (1 - 2 m / r); g22 = - r2; g33 = - r2 sin2(θ); g44 = 1 - 2 m / r
Цей розв`язок має особливу точку при виконанні умови
1
– 2 m
/ r
= 0.
Відповідне значення радіус-вектора називається
гравітаційним радіусом тіла, що є джерелом гравітаційного поля
rg = 2 m = 2 γ M / c2.
Для Сонця гравітаційний радіус дорівнює 3 км. Якщо вся маса тіла, що є джерелом гравітаційного поля,
зосереджена у сфері меншого радіусу, ніж гравітаційний, то геометрія
простору-часу всередині цієї сфери принципово відрізняється від геометрії за її
межами. Світло, що випромінюється всередині такої сфери, не зможе подолати
тяжіння такої маси і вирватись за межі такої сфери. Такий об`єкт називається чорною дірою. Такі чорні
діри знаходяться в центрах галактик, забезпечуючи їх видиму структуру.
Коментарі
Дописати коментар