нам в университете (20 лет назад) на лекциях по ОТО эту формулу выводили как раз таки без псевдятинки.
Напеть сможете?
Пусть, значит, есть система сделанная из обычной материи с тензором энергии импульса

. Она как-то там кувыркается, излучает волны гравитационные, энергия

этой системы из-за этого зависит от времени. Надо найти чему равно

.
Ведём векторнозначную 1-форму:

дуальная ей векторнозначная 3-форма:

Энергия зависит от системы отсчёта. Фиксируем систему отсчёта, вот она:

. Для краткости времениподобный вектор

обозначим

. В выбранной системе отсчёта энергия равна следующему интегралу:

Энергия сохраняется если

вектор Киллинга (

), что демонстрируется следующими выкладками:

Чтобы вычисление

напрямую дало потерю энергии связанную именно с гравитационными волнами надо считать, что в отсутствии этих волн энергия сохраняется

, то есть

вектор Киллинга исходного гравитационного поля

без гравитационной волны.
Пусть теперь на фоне

пробежала слабенькая гравитационная волна


Вектор

был вектором Киллинга для метрики

, а для метрики

он никакой не вектор Киллинга:

Ну, что, пусть фоновая метрика была метрикой Минковского, тогда



Подставляем это дело в формулу (4), в которой заменяем

на

получаем:

значит:

Ну вот почти что и всё, дальше я на словах опишу ибо очень утомительно. Теперь вместо

в формулу (10) надо подставить левую часть линеаризованных уравнений Эйнштена и аккуратно просуммировать по всем поляризациям гравитационной волны. Правая часть формулы (10) распадётся на два члена: первый в виде производной по времени от интеграла по 3D-объёму, а второй в виде интеграла по 2D-поверхности. Первым можно пренебречь (энергия

много больше). За гравитационную волну отвечает второй член.