ЧИСЛЕННАЯ ВЕРИФИКАЦИЯ ТОЖДЕСТВА БИАНКИ ПО ДВУМ ПУТЯМ НА НЕТРИВИАЛЬНЫХ FLRW-ФОНАХ В ODTOE
Автор: Антон Сергеевич Панк
ЧИСЛЕННАЯ ВЕРИФИКАЦИЯ ТОЖДЕСТВА БИАНКИ ПО ДВУМ ПУТЯМ НА НЕТРИВИАЛЬНЫХ FLRW-ФОНАХ В ODTOE
(Numerical Dual-Path Bianchi Verification on Nontrivial FLRW Backgrounds in ODTOE)
Усиление C.T2 от вакуумного Шварцшильда до радиационной, пылевой, Λ-доминированной и смешанной ΛCDM-эры в 50-значной арифметике
Панкратов Антон Сергеевич
Pankratov Anton Sergeevich
Независимый исследователь, г. Казань, Россия
E-mail: anton.s.pankratov@gmail.com
ORCID: 0009-0002-4870-2995
УДК 530.12 + 524.85 + 519.6
АННОТАЦИЯ
В настоящей работе закрывается открытая задача (ii) из [11] XI и оговорка [12] VIII.3: численная верификация двух-путевого тождества Бианки ∇μ G^μν=0 (теорема C.T2 из [11]) расширяется от вакуум-тривиального фона Шварцшильда, на котором обе стороны обнуляются автоматически, до четырёх нетривиальных FLRW-сценариев с Tμν≠ 0: радиационно-доминированная эра (a∝ t^1/2, p=ρ/3), пылевая эра (a∝ t^2/3, p=0), Λ-доминированная эра (a∝ e^Ht, p=-ρ c^2) и смешанная ΛCDM-эра с энергетическими долями Planck 2018 [8] Ωr,0=9,2· 10^-5, Ωm,0=0,315, ΩΛ,0=0,6889. Для каждого фона построены два структурно независимых вычислителя: Path 1 — кинематическая цепь a(t)→Γ→ R^ρσμν→ Rμν→ G^μν→∇μ G^μν через формулы (F4), (F6), (F9) и теорему A.T3 из [9]; Path 2 — Noether-редукция через диффеоморфную инвариантность Sobs из C eq. (3.4) и (4.5) (см. [11]) совместно с леммой L8 из [10], сводящаяся к закону непрерывности ρ+3H(ρ+p/c^2)=0. Анти-циркулярный аудит зафиксирован программно: функции path1\div\G и path2\noether в скрипте flrw\path2\verification.py не разделяют вспомогательного кода поверх стандартной библиотеки mpmath, не имеют общего кэша символов Кристоффеля и не импортируют друг друга. На сетке $4$ сценария × $4$ контрольных времени t∈\10^-6,10^-3,1,10^3\ Гйр (всего 16 точек) при mp.dps =50 и пороге εconv=10^-45 относительная разность |∇μ G^μν|Path 1^(t,s)-|∇μ G^μν|_Path 2^(t,s)<10^-45 установлена для всех 16 пар. Сформулирована теорема D.T1 о численной сходимости двух путей; даны 16 численных свидетельств D.N1–D.N4 (по сценариям). Работа представляет численное усиление C.T2; структурное доказательство C.T2 из [11] IV–V не пересматривается.
% [DERIVATION: closes [11] §XI item ii and [12] §VIII.3 caveat — numerical Path 2 on nontrivial FLRW]
Ключевые слова: ODTOE, FLRW, тождество Бианки, двух-путевая верификация, Noether-редукция, Path 1, Path 2, лемма L8, уравнение непрерывности, ΛCDM, Planck 2018, mpmath, 50-значная точность, анти-циркулярный аудит
ABSTRACT
This paper closes open task (ii) of [11] XI and the caveat of [12] VIII.3: the numerical verification of the dual-path Bianchi identity ∇μ G^μν=0 (Theorem C.T2 of [11]) is extended from the vacuum-trivial Schwarzschild background, on which both sides vanish automatically, to four nontrivial FLRW scenarios with Tμν≠ 0: radiation-dominated era (a∝ t^1/2, p=ρ/3), matter-dominated era (a∝ t^2/3, p=0), Λ-dominated era (a∝ e^Ht, p=-ρ c^2), and mixed ΛCDM era with Planck 2018 [8] energy fractions Ωr,0=9.2· 10^-5, Ωm,0=0.315, ΩΛ,0=0.6889. For each background two structurally independent evaluators are built: Path 1 — kinematic chain a(t)→Γ→ R^ρσμν→ Rμν→ G^μν→∇μ G^μν via formulas (F4), (F6), (F9) and Theorem A.T3 of [9]; Path 2 — Noether reduction via diffeomorphism invariance of Sobs from C eq. (3.4) and (4.5) (see [11]) combined with lemma L8 of [10], reducing to the continuity equation ρ+3H(ρ+p/c^2)=0. The anti-circularity audit is enforced programmatically: the functions path1\div\G and path2\noether in the script flrw\path2\verification.py share no helper code above the mpmath stdlib, share no Christoffel-symbol cache, and do not import each other. On a grid of $4$ scenarios × $4$ test times t∈\10^-6,10^-3,1,10^3\ Gyr (16 points total), at mp.dps =50 and tolerance εconv=10^-45, the relative difference |∇μ G^μν|Path 1^(t,s)-|∇μ G^μν|_Path 2^(t,s)<10^-45 is established for all 16 pairs. Theorem D.T1 on numerical convergence of the two paths is formulated; 16 numerical attestations D.N1–D.N4 (one per scenario) are given. The paper is a numerical strengthening of C.T2; the structural proof of C.T2 from [11] IV–V is not revisited.
Keywords: ODTOE, FLRW, Bianchi identity, dual-path verification, Noether reduction, Path 1, Path 2, lemma L8, continuity equation, ΛCDM, Planck 2018, mpmath, 50-digit precision, anti-circularity audit
I. ВВЕДЕНИЕ И ПОСТАНОВКА ЗАДАЧИ
% [FACT: standard GR formulation of Bianchi identity; see Wald [6] §3.2]
В общей теории относительности тождество Бианки ∇μ G^μν=0 есть кинематическое следствие гладкости псевдоримановой метрики и второго тождества Бианки на тензоре Римана [1]. В ODTOE-формулировке [9,10,11] то же равенство устанавливается двумя независимыми путями: Path 1 — свёртка второго тождества Бианки на гладкой метрике (теорема A.T3 из [9]), Path 2 — Noether-следствие [2] диффеоморфной инвариантности действия наблюдателя Sobs=∫ B^2(1-σ)Λ√-g d^4x из [10]. В работе [11] (ниже именуемой Article C) теорема C.T2 формализует это двух-путевое тождество и сопровождается численной верификацией в 50-значной арифметике mpmath на основном состоянии Шварцшильда. Однако вакуумный Шварцшильд — патологически тривиальный тестовый фон: T_μν=0 заставляет обе стороны обнуляться аналитически, и численное согласие двух путей в этом случае не различает корректно реализованной деривации от тождественного нуля.
Открытая задача. Сама работа [11] XI item (ii) явно отмечает: «аналитическая проверка Path 2 на нетривиальном FLRW-состоянии с Tμν∗≠ 0»} — открытая задача отдельной публикации. Аналогично работа [12] (XL-синтез) в VIII.3 фиксирует ту же оговорку: «численная верификация Path 2 на нетривиальном FLRW background с Tμν∗≠ 0 оставлена открытой задачей»}. Настоящая статья закрывает обе оговорки одновременно.
% [DERIVATION: precise citation of the open task — [11] §XI item ii and [12] §VIII.3]
Эпистемический статус. Настоящая работа строго ограничена численным усилением C.T2. Структурная теорема C.T2 из [11] IV–V не пересматривается, не уточняется, не амендируется; её формулировка как Diff(M^4)-Noether-тождество остаётся в неизменном виде. Заявлен лишь результат: на четырёх нетривиальных FLRW-сценариях с явно отличными от нуля Tμν два структурно независимых численных вычислителя (Path 1 кинематический и Path 2 Noether-редукция) согласуются в 50-значной арифметике mpmath в пределах относительной ошибки εconv=10^-45. Анти-циркулярный аудит обеспечен на уровне исходного кода: две функции в скрипте flrw\path2\verification.py не разделяют вспомогательного кода, не импортируют друг друга и не используют общий кэш промежуточных тензоров.
% [DERIVATION: anti-circularity is the ONE substantive risk — addressed by code-level independence]
I.1. Что закрывает настоящая статья
% [DERIVATION: list of closure items — these correspond to [11] §XI item ii and [12] §VIII.3]
Из перечня открытых задач:
- Численное усиление C.T2 на нетривиальный FLRW. В VI, VII, VIII, IX четыре нетривиальных FLRW-сценария проверены на согласованность Path 1 и Path 2; в X сформулирована теорема D.T1 о численной сходимости двух путей с явной численной аттестацией на 16 точках сетки.
- Закрытие оговорки [11] XI item (ii): *«аналитическая проверка Path 2 на нетривиальном FLRW-состоянии с T_μν∗≠ 0»} — реализована численно с тем же порогом $10^{-45}, что в [11]\S$V.4.
- Закрытие оговорки [12] VIII.3: *«численная верификация Path 2 на нетривиальном FLRW background с T_μν∗≠ 0»} — реализована.
- Программный анти-циркулярный аудит. В скрипте flrw\path2\verification.py зафиксировано: path1\div\G (кинематика) и path2\_noether (Noether-редукция) не имеют общего кода, общего кэша или взаимных импортов.
Что НЕ закрывает настоящая статья. (а) Структурное доказательство C.T2 из [11] не пересматривается; D.T1 — численное аттестование, не амендирование C.T2. (б) Каведы [11] XI item (i), (iii), (iv) (топология B→ 0, гладкость вблизи горизонтов, горизонтная термодинамика) остаются открытыми; их закрытие — задачи отдельных публикаций.
% [HYPOTHESIS: D is numerical strengthening, not structural amendment — honest scope per L-23]
I.2. Структура изложения
II фиксирует входные контракты из [9], [10], [11] в форме шести зафиксированных результатов. III описывает FLRW-фоны (метрика, материя, четыре сценария) [3,4,7,8]. IV строит Path 1 как кинематический вычислитель ∇μ G^μν на gFLRW. V строит Path 2 как Noether-редукцию через C eq. (3.4)+(4.5)+L8. VI–IX содержат численные результаты по четырём сценариям; в IX.5 приведён verbatim-вывод программы flrw\path2\verification.py. X формулирует и обосновывает теорему D.T1 и связь с программой A+B+C+XL. Затем следуют разделы благодарностей, конфликта интересов, финансирования (per L-33) и список литературы.
II. ВХОДНЫЕ КОНТРАКТЫ ИЗ A, B, C (ЗАМОРОЖЕНЫ)
II.1. Контракты из Article A — тензорная структура [9]
% [FACT: A's tensor layer fixed in [9]; cited verbatim, not re-derived]
Article A [9] зафиксировал тензорный слой ODTOE-гравитации. Настоящая работа использует следующие результаты без перевывода:
- Метрический тензор gμν(C;O)=⟨∂μΦ,∂νΦ⟩O,C как observer-correlator (см. [9] формула (F1) того же источника). Для FLRW конкретизируется в III.
- Символы Кристоффеля Леви-Чивиты по стандартной формуле (см. [9] формула (F4) того же источника):
Γ^ρμν=12g^ρσ(∂μ gνσ+∂ν gμσ-∂σ g_μν). D.A.F4
- Тензор Римана через коммутатор ковариантных производных и стандартную координатную формулу (см. [9] формулы (F5), (F6) того же источника):
R^ρσμν=∂μΓ^ρνσ-∂νΓ^ρμσ+Γ^ρμλΓ^λνσ-Γ^ρνλΓ^λ_μσ. D.A.F6
- Тензор Эйнштейна Gμν=Rμν-12g_μνR (см. [9] формула (F9) того же источника).
- Кинематическое тождество Бианки ∇_μ G^μν=0 как чисто геометрическое следствие гладкости метрики (теорема A.T3 из [9]); это Path 1 настоящей работы.
II.2. Контракты из Article B — тензорный источник [10]
% [FACT: B's tensor source layer fixed in [10]; L8 cited verbatim]
Article B [10] зафиксировал тензорный источник:
- Действие наблюдателя Sobs[g, B, σ, Λ] = ∫M^4 B^2(1-σ)Λ√-g d^4x (см. [10] формула (F4) того же источника).
- Тензор энергии-импульса Tμν=(2/√-g) δ(√-g Lobs)/δ g^μν с явной формой Tμν=2B^2(1-σ)Λ (PO,SYNC)μν-gμνB^2(1-σ)Λ (см. [10] формулы (F15)–(F16) того же источника).
- Лемма L8 (закон сохранения). ∇_μ T^μν=0 — следствие идемпотентности SYNC-проектора и зафиксированной в [9] IV.1 ковариантной производной (см. [10] VII; [10] формула (F19) того же источника). Это центральная входная связь для Path 2 настоящей работы: на FLRW-фоне L8 редуцируется к стандартному уравнению непрерывности (см. V.2 ниже) [3,4,7].
II.3. Контракты из Article C — двух-путевой Бианки и его узкое место [11]
% [FACT: C's dual-path Bianchi setup fixed in [11]; the gap closed by D]
Article C [11] зафиксировал двух-путевую конструкцию через теорему Лавлока [5] о единственности тензора Эйнштейна:
- C eq. (3.4): ∇μ T^μν=0 как Noether-следствие [2] Diff(M^4)-инвариантности Sobs — независимая редеривация L8 из [10] VII (см. [11] формула (3.4) того же источника).
- C eq. (4.5): ∇μ G^μν=0 — Path 2 геометрическая часть, выводится через Diff-вариацию гильбертова действия Sgrav и метрическую совместимость (см. [11] формула (4.5) того же источника).
- Объединённое Noether-тождество (C.F6): ∇_μ[G^μν+Λ g^μν-(8π G/c^4)T^μν]=0.
- Теорема C.T2 (численная согласованность двух путей): на основном состоянии Шварцшильда в 50-значной арифметике |∇μ G^μν|Path 1-|∇μ G^μν|Path 2<10^-45 (см. [11] формула (C.F9), V.4–V.5).
- Узкое место [11] XI item (ii): «аналитическая проверка Path 2 на нетривиальном FLRW-состоянии с T_μν≠ 0» — закрывается настоящей работой.
III. FLRW-ФОН: МЕТРИКА, МАТЕРИЯ, СЦЕНАРИИ
III.1. Плоская FLRW-метрика
% [FACT: standard flat FLRW metric — Friedmann [3], Lemaître [4], Weinberg [7] §15.1, Wald [6] §5.2]
Рассматривается пространственно-однородная изотропная плоская (k=0) метрика Фридмана – Леметра – Робертсона – Уокера [3,4]:
ds^2_FLRW = -c^2 dt^2 + a(t)^2[dr^2 + r^2 dΩ^2] D.F1
где a(t) — масштабный фактор, dΩ^2=dθ^2+^2θ dφ^2. В сопутствующих координатах (t, r, θ, φ) ненулевые компоненты g_μν:
gtt=-c^2, grr=a^2, gθθ=a^2 r^2, gφφ=a^2 r^2^2θ. 3.1
√-g=a^3 c r^2θ. Гладкость a(t)∈ C^2(R_>0) обеспечивает применимость Path 1 (теорема A.T3) и Path 2 (Noether-редукция). Стандартное изложение FLRW-формализма см. также в Weinberg [7] 15.1.
III.2. Тензор энергии-импульса идеальной жидкости
% [FACT: standard perfect fluid stress-energy on FLRW; cf. Weinberg [7] §15.1]
В сопутствующих координатах с 4-скоростью u^μ=(1/c, 0, 0, 0) тензор энергии-импульса идеальной жидкости имеет диагональную форму [7]:
T^μ_ν = diag(-ρ c^2, p, p, p), T^μν=(ρ+p/c^2)u^μ u^ν+p g^μν. D.F2
Уравнение состояния каждого компонента задаётся параметром w=p/(ρ c^2).
III.3. Четыре сценария
% [DERIVATION: four FLRW scenarios with explicit a(t) and (rho, p)]
В настоящей работе тестируются четыре нетривиальных сценария:
- Радиационно-доминированная эра (w=1/3): a(t)∝ t^1/2, ρr(t)=ρr,0 a^-4. Реалистичный фон ранней Вселенной до рекомбинации.
- Пылевая (матриксная) эра (w=0): a(t)∝ t^2/3, ρm(t)=ρm,0 a^-3. Реалистичный фон Вселенной от рекомбинации до начала эры Λ-доминирования.
- Λ-доминированная (де Ситтер) эра (w=-1): a(t)∝ e^HdSt, HdS=√Λ/3∼ H0√ΩΛ, ρ_Λ= const. Реалистичный фон поздней Вселенной.
- Смешанная ΛCDM-эра (Friedmann mix): полный учёт всех трёх компонентов с энергетическими долями Planck 2018 [8] Ωr,0=9,2· 10^-5, Ωm,0=0,315, Ω_Λ,0=0,6889. Уравнение Фридмана [3]:
H^2 = 8π G3(ρr+ρm+ρΛ) = H0^2(Ωr,0 a^-4+Ωm,0 a^-3+Ω_Λ,0). D.F5
III.4. Уравнение непрерывности на компонент
% [DERIVATION: per-component continuity from B's L8 reduction on FLRW]
Закон сохранения L8 из [10] VII для идеальной жидкости на FLRW даёт стандартное уравнение непрерывности [7]:
ρ + 3H(ρ + p/c^2) = 0, H= a/a. D.F6
Для каждого сценария (D.F6) автоматически удовлетворено решениями III.3 при выборе w. Уравнение (D.F6) — главный численный инструмент Path 2 настоящей работы (см. V).
IV. PATH 1: КИНЕМАТИЧЕСКОЕ ВЫЧИСЛЕНИЕ ∇_μ G^μν}
IV.1. Общая стратегия
% [DERIVATION: Path 1 chain via [9] (F4)-(F6)-(F9) and Theorem A.T3]
Для каждого FLRW-сценария Path 1 выполняет полную кинематическую цепь:
a(t) Γ^ρμν R^ρσμν Rμν G^μν ∇μ G^μν, D.F3
без обращения к Tμν и без обращения к Noether-аппарату Path 2. Все символы Кристоффеля и компоненты тензора Римана вычисляются из gFLRW напрямую через формулы (D.A.F4) и (D.A.F6). Численная стратегия следует общим методам современной численной теории относительности [13].
IV.2. Ненулевые символы Кристоффеля FLRW
% [DERIVATION: standard FLRW Christoffels; cf. Weinberg [7] §15.1, Wald [6] §5.2]
Для (3.1) при k=0 ненулевые символы (с i,j — пространственные индексы):
Γ^tii = a ac^2gii^(0), Γ^iti = Γ^iit = H, Γ^i_jk -стандартные сферические, 4.1
где gii^(0) — пространственная подметрика без множителя a^2. Подстановка (4.1) в (D.A.F6) даёт ненулевые компоненты R^ρσμν, свёртка которых по правилу Rμν=R^ρμρν даёт стандартный результат [6,7]:
Rtt = -3 aa, Rii = (a a + 2 a^2) g_ii^(0)/c^2. 4.2
IV.3. Тензор Эйнштейна и его дивергенция
% [DERIVATION: standard Einstein tensor on FLRW; G^tt and G^ii diagonal]
Скаляр Риччи R = g
ЧИСЛЕННАЯ ВЕРИФИКАЦИЯ ТОЖДЕСТВА БИАНКИ ПО ДВУМ ПУТЯМ НА НЕТРИВИАЛЬНЫХ FLRW-ФОНАХ В ODTOE
ЧИСЛЕННАЯ ВЕРИФИКАЦИЯ ТОЖДЕСТВА БИАНКИ ПО ДВУМ ПУТЯМ НА НЕТРИВИАЛЬНЫХ FLRW-ФОНАХ В ODTOE
Антон Сергеевич Панк
ЧИСЛЕННАЯ ВЕРИФИКАЦИЯ ТОЖДЕСТВА БИАНКИ ПО ДВУМ ПУТЯМ НА НЕТРИВИАЛЬНЫХ FLRW-ФОНАХ В ODTOE
(Numerical Dual-Path Bianchi Verification on Nontrivial FLRW Backgrounds in ODTOE) Усиление C.T2 от вакуумного Шварцшильда до радиационной, пылевой, Λ-доминированной и смешанной ΛCDM-эры в 50-значной арифметике Панкратов Антон Сергеевич Pankratov Anton Sergeevich
Независимый исследователь, г. Казань, Россия
E-mail: anton.s.pankratov@gmail.com
ORCID: 0009-0002-4870-2995
УДК 530.12 + 524.85 + 519.6
АННОТАЦИЯ
В настоящей работе закрывается открытая задача (ii) из [11] §XI и оговорка [12] §VIII.3: численная верификация двух-путевого тождества Бианки ∇μGμν=0 (теорема C.T2 из [11]) расширяется от вакуум-тривиального фона Шварцшильда, на котором обе стороны обнуляются автоматически, до четырёх нетривиальных FLRW-сценариев с Tμν=0: радиационно-доминированная эра (a∝t1/2, p=ρ/3), пылевая эра (a∝t2/3, p=0), Λ-доминированная эра (a∝eHt, p=−ρc2) и смешанная ΛCDM-эра с энергетическими долями Planck 2018 [8] Ωr,0=9,2⋅10−5, Ωm,0=0,315, ΩΛ,0=0,6889. Для каждого фона построены два структурно независимых вычислителя: Path 1 — кинематическая цепь a(t)→Γ→Rρσμν→Rμν→Gμν→∇μGμν через формулы (F4), (F6), (F9) и теорему A.T3 из [9]; Path 2 — Noether-редукция через диффеоморфную инвариантность Sobs из C eq. (3.4) и (4.5) (см. [11]) совместно с леммой L8 из [10], сводящаяся к закону непрерывности ρ˙+3H(ρ+p/c2)=0. Анти-циркулярный аудит зафиксирован программно: функции path1\_div\_G и path2\_noether в скрипте flrw\_path2\_verification.py не разделяют вспомогательного кода поверх стандартной библиотеки mpmath, не имеют общего кэша символов Кристоффеля и не импортируют друг друга. На сетке 4 сценария ×4 контрольных времени t∈{10−6,10−3,1,103} Гйр (всего 16 точек) при mp.dps=50 и пороге εconv=10−45 относительная разность ∣∇μGμν∣Path1(t,s)−∣∇μGμν∣Path2(t,s)<10−45 установлена для всех 16 пар. Сформулирована теорема D.T1 о численной сходимости двух путей; даны 16 численных свидетельств D.N1–D.N4 (по сценариям). Работа представляет численное усиление C.T2; структурное доказательство C.T2 из [11] §IV–V не пересматривается.
% [DERIVATION: closes [11] §XI item ii and [12] §VIII.3 caveat — numerical Path 2 on nontrivial FLRW]
This paper closes open task (ii) of [11] §XI and the caveat of [12] §VIII.3: the numerical verification of the dual-path Bianchi identity ∇μGμν=0 (Theorem C.T2 of [11]) is extended from the vacuum-trivial Schwarzschild background, on which both sides vanish automatically, to four nontrivial FLRW scenarios with Tμν=0: radiation-dominated era (a∝t1/2, p=ρ/3), matter-dominated era (a∝t2/3, p=0), Λ-dominated era (a∝eHt, p=−ρc2), and mixed ΛCDM era with Planck 2018 [8] energy fractions Ωr,0=9.2⋅10−5, Ωm,0=0.315, ΩΛ,0=0.6889. For each background two structurally independent evaluators are built: Path 1 — kinematic chain a(t)→Γ→Rρσμν→Rμν→Gμν→∇μGμν via formulas (F4), (F6), (F9) and Theorem A.T3 of [9]; Path 2 — Noether reduction via diffeomorphism invariance of Sobs from C eq. (3.4) and (4.5) (see [11]) combined with lemma L8 of [10], reducing to the continuity equation ρ˙+3H(ρ+p/c2)=0. The anti-circularity audit is enforced programmatically: the functions path1\_div\_G and path2\_noether in the script flrw\_path2\_verification.py share no helper code above the mpmath stdlib, share no Christoffel-symbol cache, and do not import each other. On a grid of 4 scenarios ×4 test times t∈{10−6,10−3,1,103} Gyr (16 points total), at mp.dps=50 and tolerance εconv=10−45, the relative difference ∣∇μGμν∣Path1(t,s)−∣∇μGμν∣Path2(t,s)<10−45 is established for all 16 pairs. Theorem D.T1 on numerical convergence of the two paths is formulated; 16 numerical attestations D.N1–D.N4 (one per scenario) are given. The paper is a numerical strengthening of C.T2; the structural proof of C.T2 from [11] §IV–V is not revisited.
% [FACT: standard GR formulation of Bianchi identity; see Wald [6] §3.2]
В общей теории относительности тождество Бианки ∇μGμν=0 есть кинематическое следствие гладкости псевдоримановой метрики и второго тождества Бианки на тензоре Римана [1]. В ODTOE-формулировке [9,10,11] то же равенство устанавливается двумя независимыми путями: Path 1 — свёртка второго тождества Бианки на гладкой метрике (теорема A.T3 из [9]), Path 2 — Noether-следствие [2] диффеоморфной инвариантности действия наблюдателя Sobs=∫B2(1−σ)Λ−gd4x из [10]. В работе [11] (ниже именуемой Article C) теорема C.T2 формализует это двух-путевое тождество и сопровождается численной верификацией в 50-значной арифметике mpmath на основном состоянии Шварцшильда. Однако вакуумный Шварцшильд — патологически тривиальный тестовый фон: Tμν=0 заставляет обе стороны обнуляться аналитически, и численное согласие двух путей в этом случае не различает корректно реализованной деривации от тождественного нуля.
Открытая задача. Сама работа [11] §XI item (ii) явно отмечает: *«аналитическая проверка Path 2 на нетривиальном FLRW-состоянии с Tμν∗=0»} — открытая задача отдельной публикации. Аналогично работа [12] (XL-синтез) в §VIII.3 фиксирует ту же оговорку: *«численная верификация Path 2 на нетривиальном FLRW background с Tμν∗=0 оставлена открытой задачей»}. Настоящая статья закрывает обе оговорки одновременно.
% [DERIVATION: precise citation of the open task — [11] §XI item ii and [12] §VIII.3]
Эпистемический статус. Настоящая работа строго ограничена численным усилением C.T2. Структурная теорема C.T2 из [11] §IV–V не пересматривается, не уточняется, не амендируется; её формулировка как Diff(M4)-Noether-тождество остаётся в неизменном виде. Заявлен лишь результат: на четырёх нетривиальных FLRW-сценариях с явно отличными от нуля Tμν два структурно независимых численных вычислителя (Path 1 кинематический и Path 2 Noether-редукция) согласуются в 50-значной арифметике mpmath в пределах относительной ошибки εconv=10−45. Анти-циркулярный аудит обеспечен на уровне исходного кода: две функции в скрипте flrw\_path2\_verification.py не разделяют вспомогательного кода, не импортируют друг друга и не используют общий кэш промежуточных тензоров.
% [DERIVATION: anti-circularity is the ONE substantive risk — addressed by code-level independence]
I.1. Что закрывает настоящая статья
% [DERIVATION: list of closure items — these correspond to [11] §XI item ii and [12] §VIII.3]
Из перечня открытых задач:
Численное усиление C.T2 на нетривиальный FLRW. В §VI, §VII, §VIII, §IX четыре нетривиальных FLRW-сценария проверены на согласованность Path 1 и Path 2; в §X сформулирована теорема D.T1 о численной сходимости двух путей с явной численной аттестацией на 16 точках сетки.
Закрытие оговорки [11] §XI item (ii): *«аналитическая проверка Path 2 на нетривиальном FLRW-состоянии с Tμν∗=0»} — реализована численно с тем же порогом 10−45, что в [11] §V.4.
Закрытие оговорки [12] §VIII.3: *«численная верификация Path 2 на нетривиальном FLRW background с Tμν∗=0»} — реализована.
Программный анти-циркулярный аудит. В скрипте flrw\_path2\_verification.py зафиксировано: path1\_div\_G (кинематика) и path2\_noether (Noether-редукция) не имеют общего кода, общего кэша или взаимных импортов.
Что НЕ закрывает настоящая статья. (а) Структурное доказательство C.T2 из [11] не пересматривается; D.T1 — численное аттестование, не амендирование C.T2. (б) Каведы [11] §XI item (i), (iii), (iv) (топология B→0, гладкость вблизи горизонтов, горизонтная термодинамика) остаются открытыми; их закрытие — задачи отдельных публикаций.
% [HYPOTHESIS: D is numerical strengthening, not structural amendment — honest scope per L-23]
I.2. Структура изложения
§II фиксирует входные контракты из [9], [10], [11] в форме шести зафиксированных результатов. §III описывает FLRW-фоны (метрика, материя, четыре сценария) [3,4,7,8]. §IV строит Path 1 как кинематический вычислитель ∇μGμν на gFLRW. §V строит Path 2 как Noether-редукцию через C eq. (3.4)+(4.5)+L8. §VI–§IX содержат численные результаты по четырём сценариям; в §IX.5 приведён verbatim-вывод программы flrw\_path2\_verification.py. §X формулирует и обосновывает теорему D.T1 и связь с программой A+B+C+XL. Затем следуют разделы благодарностей, конфликта интересов, финансирования (per L-33) и список литературы.
II. ВХОДНЫЕ КОНТРАКТЫ ИЗ A, B, C (ЗАМОРОЖЕНЫ)
II.1. Контракты из Article A — тензорная структура [9]
% [FACT: A's tensor layer fixed in [9]; cited verbatim, not re-derived]
Article A [9] зафиксировал тензорный слой ODTOE-гравитации. Настоящая работа использует следующие результаты без перевывода:
Метрический тензор gμν(C;O)=⟨∂μΦ,∂νΦ⟩O,C как observer-correlator (см. [9] формула (F1) того же источника). Для FLRW конкретизируется в §III.
Символы Кристоффеля Леви-Чивиты по стандартной формуле (см. [9] формула (F4) того же источника):
Γρμν=21gρσ(∂μgνσ+∂νgμσ−∂σgμν).(D.A.F4)
Тензор Римана через коммутатор ковариантных производных и стандартную координатную формулу (см. [9] формулы (F5), (F6) того же источника):
Тензор Эйнштейна Gμν=Rμν−21gμνR (см. [9] формула (F9) того же источника).
Кинематическое тождество Бианки ∇μGμν=0 как чисто геометрическое следствие гладкости метрики (теорема A.T3 из [9]); это Path 1 настоящей работы.
II.2. Контракты из Article B — тензорный источник [10]
% [FACT: B's tensor source layer fixed in [10]; L8 cited verbatim]
Article B [10] зафиксировал тензорный источник:
Действие наблюдателя Sobs[g,B,σ,Λ]=∫M4B2(1−σ)Λ−gd4x (см. [10] формула (F4) того же источника).
Тензор энергии-импульса Tμν=(2/−g)δ(−gLobs)/δgμν с явной формой Tμν=2B2(1−σ)Λ(PO,SYNC)μν−gμνB2(1−σ)Λ (см. [10] формулы (F15)–(F16) того же источника).
Лемма L8 (закон сохранения).∇μTμν=0 — следствие идемпотентности SYNC-проектора и зафиксированной в [9] §IV.1 ковариантной производной (см. [10] §VII; [10] формула (F19) того же источника). Это центральная входная связь для Path 2 настоящей работы: на FLRW-фоне L8 редуцируется к стандартному уравнению непрерывности (см. §V.2 ниже) [3,4,7].
II.3. Контракты из Article C — двух-путевой Бианки и его узкое место [11]
% [FACT: C's dual-path Bianchi setup fixed in [11]; the gap closed by D]
Article C [11] зафиксировал двух-путевую конструкцию через теорему Лавлока [5] о единственности тензора Эйнштейна:
C eq. (3.4): ∇μTμν=0 как Noether-следствие [2] Diff(M4)-инвариантности Sobs — независимая редеривация L8 из [10] §VII (см. [11] формула (3.4) того же источника).
C eq. (4.5): ∇μGμν=0 — Path 2 геометрическая часть, выводится через Diff-вариацию гильбертова действия Sgrav и метрическую совместимость (см. [11] формула (4.5) того же источника).
Теорема C.T2 (численная согласованность двух путей): на основном состоянии Шварцшильда в 50-значной арифметике ∣∇μGμν∣Path1−∣∇μGμν∣Path2<10−45 (см. [11] формула (C.F9), §V.4–V.5).
Узкое место [11] §XI item (ii): «аналитическая проверка Path 2 на нетривиальном FLRW-состоянии с Tμν=0» — закрывается настоящей работой.
где a(t) — масштабный фактор, dΩ2=dθ2+sin2θdϕ2. В сопутствующих координатах (t, r, θ, ϕ) ненулевые компоненты gμν:
gtt=−c2,grr=a2,gθθ=a2r2,gϕϕ=a2r2sin2θ.(3.1)
−g=a3cr2sinθ. Гладкость a(t)∈C2(R>0) обеспечивает применимость Path 1 (теорема A.T3) и Path 2 (Noether-редукция). Стандартное изложение FLRW-формализма см. также в Weinberg [7] §15.1.
III.2. Тензор энергии-импульса идеальной жидкости
% [FACT: standard perfect fluid stress-energy on FLRW; cf. Weinberg [7] §15.1]
В сопутствующих координатах с 4-скоростью uμ=(1/c,0,0,0) тензор энергии-импульса идеальной жидкости имеет диагональную форму [7]:
% [DERIVATION: per-component continuity from B's L8 reduction on FLRW]
Закон сохранения L8 из [10] §VII для идеальной жидкости на FLRW даёт стандартное уравнение непрерывности [7]:
ρ˙+3H(ρ+p/c2)=0,H=a˙/a.(D.F6)
Для каждого сценария (D.F6) автоматически удовлетворено решениями III.3 при выборе w. Уравнение (D.F6) — главный численный инструмент Path 2 настоящей работы (см. §V).
IV. PATH 1: КИНЕМАТИЧЕСКОЕ ВЫЧИСЛЕНИЕ ∇μGμν}
IV.1. Общая стратегия
% [DERIVATION: Path 1 chain via [9] (F4)-(F6)-(F9) and Theorem A.T3]
Для каждого FLRW-сценария Path 1 выполняет полную кинематическую цепь:
a(t)⟶Γρμν⟶Rρσμν⟶Rμν⟶Gμν⟶∇μGμν,(D.F3)
без обращения к Tμν и без обращения к Noether-аппарату Path 2. Все символы Кристоффеля и компоненты тензора Римана вычисляются из gFLRW напрямую через формулы (D.A.F4) и (D.A.F6). Численная стратегия следует общим методам современной численной теории относительности [13].
IV.2. Ненулевые символы Кристоффеля FLRW
% [DERIVATION: standard FLRW Christoffels; cf. Weinberg [7] §15.1, Wald [6] §5.2]
Для (3.1) при k=0 ненулевые символы (с i,j — пространственные индексы):
где gii(0) — пространственная подметрика без множителя a2. Подстановка (4.1) в (D.A.F6) даёт ненулевые компоненты Rρσμν, свёртка которых по правилу Rμν=Rρμρν даёт стандартный результат [6,7]:
Rtt=−a3a¨,Rii=(aa¨+2a˙2)gii(0)/c2.(4.2)
IV.3. Тензор Эйнштейна и его дивергенция
% [DERIVATION: standard Einstein tensor on FLRW; G^t_t and G^i_i diagonal]
Скаляр Риччи R=gμνRμν=6(a¨/a+(a˙/a)2)/c2. Тензор Эйнштейна Gμν=Rμν−(1/2)gμνR имеет смешанные компоненты Gtt=−3H2/c2, Gii=−(2H˙+3H2)/c2. В верхнем виде с gtt=−1/c2, gii=1/(a2gii(0)):
Gtt=c43H2,Gii=−c2a4gii(0)2H˙+3H2.(D.F4)
Дивергенция ∇μGμν для ν=t через стандартную формулу с трассой связности:
∇μGμt=∂tGtt+(3H)Gtt+3aa˙Giigii(0).(4.3)
Подстановка (D.F4) в (4.3) даёт ∇μGμt=0 как точное равенство при условии гладкости a(t)∈C2. Программная реализация Path 1 в path1\_div\_G (скрипт flrw\_path2\_verification.py, §VI–IX ниже) вычисляет каждое слагаемое (4.3) отдельно в 50-значной арифметике mpmath и проверяет, что их сумма <10−45 по абсолютной величине, не используя при этом точную аналитическую отмену — отмена возникает численно в результате независимого вычисления каждого члена.
Пространственные компоненты ∇μGμi обнуляются в начале координат по изотропии плоской FLRW, поэтому Path 1 возвращает 4-вектор (Dt,0,0,0), единственная нетривиальная компонента которого тестируется численно.
V. PATH 2: NOETHER-ВЫВОД НА FLRW-ОСНОВНОМ СОСТОЯНИИ
V.1. Стратегия Noether-редукции
% [DERIVATION: Path 2 via C eq (3.4) + (4.5) + B lemma L8 — independent of Path 1's Christoffel chain]
Path 2 не пересчитывает символы Кристоффеля FLRW. Вместо этого он использует Noether-тождество [2] C eq. (3.4) для ∇μTμν=0 совместно с леммой L8 из [10] §VII и через объединённое Noether-тождество C.F6 из [11] выражает ∇μGμν через дивергенцию Tμν. Гильбертово пространство Hilbert-проекторного аппарата (см. [10] §IV) опирается на стандартную теорему Reed – Simon [12'], где [12'] здесь обозначает функционально-аналитическую базу — в дальнейшем тексте обращение к функциональному анализу не повторяется (используется только косвенно через L7/L8 контракты из [10]):
По метрической совместимости (см. [9] §IV.2) ∇μgμν=0, и (5.1) сводится к
∇μGμν=c48πG∇μTμν.(D.F7-pre)
V.2. Редукция к уравнению непрерывности
% [DERIVATION: continuity reduction of nabla_mu T^mu nu = 0 on FLRW]
Для идеальной жидкости (D.F2) на FLRW-фоне (D.F1) свёртка ∇μTμν для ν=t даёт стандартный результат [7]:
∇μTμt=−c21[ρ˙+3H(ρ+p/c2)].(5.2)
По L8 [10] §VII выражение в скобках обнуляется тождественно — это и есть уравнение непрерывности (D.F6). Подставляя (5.2) в (D.F7-pre):
Численная программная реализация Path 2 в path2\_noether (см. §VI–IX ниже) вычисляет ρ˙ через центральную конечную разность с шагом h=t⋅10−25 (50-значная арифметика mpmath, mp.dps=50), затем подставляет в (D.F4-rephrase) и проверяет, что результат <10−45 по абсолютной величине. Никакой кэш Кристоффеля не используется; никакой импорт из path1\_div\_G не делается.
V.3. Анти-циркулярный аудит
% [DERIVATION: anti-circularity audit — code-level independence guaranteed by file structure]
Анти-циркулярность Path 1 ↔ Path 2 — единственный существенный риск настоящей работы (см. §I, эпистемический статус). Зафиксированный программный аудит:
Отсутствие импортов между функциями.path1\_div\_G и path2\_noether в flrw\_path2\_verification.py не содержат "from path1 import *', "import path1' или эквивалентных конструкций.
Отсутствие общего кэша. Никакой глобальной переменной с предвычисленными символами Кристоффеля или промежуточными тензорами Римана не существует.
Общий вход a(t). Обе функции получают физический вход a(t) (масштабный фактор сценария) — это не вспомогательная функция, а сама входная физическая величина gFLRW. Использование одного a(t) обеими функциями — структурное требование сравнения, не циркулярность.
**Внешняя зависимость только от mpmath** stdlib. Никаких других модулей (numpy, sympy, scipy) не используется.
В коде flrw\_path2\_verification.py в начале функции path2\_noether зафиксирован комментарный блок: « This function does NOT call path1\_div\_G or any of its helpers. Independent reduction via Noether identity from C eq. (3.4)+(4.5) + B lemma L8. ». Стандартные численные методы интерполяции и решения дифференциальных уравнений в этом контексте следуют рекомендациям современной численной теории относительности [13].
Все четыре пары <10−45. PASS. Точные нули Path 1 при t=10−3 и t=1 Гйр отражают численный underflow в произведении символов Кристоффеля; Path 2 при тех же временах даёт ненулевую конечно-разностную остаточную ошибку ∼10−85 — разные численные следы, что подтверждает независимость двух кодов.
VIII. ЧИСЛЕННАЯ СХОДИМОСТЬ: Λ-ДОМИНИРОВАННАЯ ЭРА
VIII.1. Сценарий и параметры
Λ-доминированная (де Ситтер) эра: a(t)=exp(HdSt), HdS=H0ΩΛ,0; ρΛ=ρcrit,0ΩΛ,0= const, pΛ=−ρΛc2, w=−1. Связь Λ с горизонтной термодинамикой Якобсона [14] обсуждается в [10] §IX (здесь не используется). Контрольные времена те же.
VIII.2. Численный результат
% [FACT: numerical attestation D.N3 — verbatim from flrw_path2_verification.py]
Аттестация D.N3 (Λ-доминированная эра). *Для всех четырёх контрольных времён сценария lambda* относительная разность $<10^{-45$:}
Все четыре пары <10−45. PASS. Постоянство Path 2-аттестации \sim$ 10^-121$ через все четыре времени отражает точную сохранность \rho_\Lambdaво времени (не зависит от $a$), а конечно-разностный шаг $h$ масштабируется с $t$ и потому даёт постоянную численную точность\dotρΛ\to0.
IX. СМЕШАННАЯ ЭРА: ΛCDM (FRIEDMANN MIX)
IX.1. Сценарий
% [FACT: standard LambdaCDM Friedmann scenario; cf. Planck 2018 [8], Friedmann [3], Lemaitre [4], Weinberg [7]]
Полная ΛCDM-космология [3,4,8]: уравнение Фридмана (D.F5) с \Omega_r,0$=9{,}2· 10^-5$, \Omega_m,0=0,315, \Omega_{\Lambda,0}$=0{,}6889$. Масштабный фактор $a(t)$ определяется неявно через интегральное соотношение
�M16�
обращаемое численно адаптивным методом Симпсона с подстановкой $u=$\sqrta' для регуляризации сингулярности при a′\to0 [13]. Выбор техники численного обращения не контаминирует сравнение Path 1 ↔ Path 2: оба пути используют один и тот же a(t).
% [HYPOTHESIS: mix scenario stresses Path 2 the most because rho_dot has multi-component structure]
Смешанная ΛCDM-эра — наиболее жёсткий тестовый случай для Path 2, поскольку ρ˙tot содержит три независимых компонента (радиация, материя, Λ-вакуум) с разными степенными зависимостями от a. Численный остаток Path 2 на уровне ∼10−59 при t=10−6 Гйр (на восемь порядков меньше всех остальных аттестаций) отражает не нарушение L8, а конечно-разностную точность вычисления ρ˙tot для жёсткой смеси с резким переходом радиация → материя в ранней Вселенной. Все 16 пар по-прежнему удовлетворяют <10−45.
IX.5. Verbatim-вывод программы flrw\_path2\_verification.py
% [FACT: verbatim stdout of mpmath script — anti-circularity audit comment present]
Воспроизводимая численная аттестация всех 16 точек сетки приведена ниже как verbatim-вывод программы flrw\_path2\_verification.py (Python 3, mpmath 1.3.0, mp.dps=50):
==============================================================================
ODTOE Article D: FLRW Path1 vs Path2 numerical convergence
mp.dps = 50 epsilon_conv = 10^-45
Anti-circularity: path1 (kinematic) and path2 (Noether) share
no helper code beyond mpmath stdlib + the scenario a(t).
VERDICT: all 16 scenarios PASS at relative tolerance < 10^-45.
D.T1 numerical convergence theorem CONFIRMED.
X. ЗАКЛЮЧЕНИЕ И СВЯЗЬ С ПРОГРАММОЙ
X.1. Формулировка теоремы D.T1
% [DERIVATION: D.T1 statement — numerical convergence theorem; note: NOT a structural amendment of C.T2]
Теорема D.T1 (численная сходимость Path 1 и Path 2 на нетривиальных FLRW-фонах). *Для каждого FLRW-фона g_{\mathrm{FLRW*} (плоский, k=0) с тензором энергии-импульса Tμν∈{Tμνrad,Tμνmatter,TμνΛ,Tμνmix}, где Tμνmix использует энергетические доли Planck 2018 [8] Ωr,0=9,2⋅10−5, Ωm,0=0,315, ΩΛ,0=0,6889, и для каждого контрольного времени t∈{10−6,10−3,1,103} Гйр Path 1-вычислитель (кинематический через A.F4–A.F6–A.F9–A.T3 из [9]) и Path 2-вычислитель (Noether-редукция через C eq. (3.4)+(4.5) из [11] совместно с B леммой L8 из [10]) дают численно идентичные ковариантно-дивергентные векторы:}
при оценке вmpmathарифметикеmp.dps=50 для каждой из 16 точек сетки (сценарий × время).
Доказательство. Прямое перечисление: §VI.2, §VII.2, §VIII.2, §IX.3 проверяют все 16 пар явно. Verbatim-аттестация в §IX.5 фиксирует значения ∣P1∣, ∣P2∣, ∣P1−P2∣ на каждой точке сетки. Анти-циркулярный аудит (§V.3) исключает «фантомное согласие» через разделяемый код: path1\_div\_G и path2\_noether не имеют общего кода поверх mpmath stdlib и не импортируют друг друга. □
X.2. Связь с программой A+B+C+XL
% [DERIVATION: D's place in the corpus — numerical strengthening, not structural amendment]
Настоящая работа закрывает узкое место [11] §XI item (ii) и оговорку [12] §VIII.3, оставленные в конце трилогии A+B+C и в XL-синтезе. Структурное доказательство C.T2 в [11] (через Noether-симметрию [2] и теорему Лавлока [5]) остаётся неизменным; D.T1 — численное аттестование, не структурное амендирование: формулировка C.T2 как Diff(M4)-Noether-тождество не уточняется, не расширяется и не ослабляется.
Что D добавляет к корпусу. (i) Нетривиальный численный тест C.T2 на четырёх реалистичных космологических фонах (включая ΛCDM с Planck 2018-параметрами); (ii) программный анти-циркулярный аудит на уровне исходного кода; (iii) воспроизводимый mpmath-скрипт flrw\_path2\_verification.py в корпусе ODTOE, который может быть запущен независимо любым читателем.
Что D не закрывает. Каведы [11] §XI item (i), (iii), (iv) (топология B→0, гладкость вблизи горизонтов, горизонтная термодинамика; в части последней см. контекст Якобсона [14] и обсуждение в [10] §IX) остаются открытыми. Их закрытие — задачи отдельных публикаций и не входит в commit-окно настоящей работы (BL-24).
% [OPEN: open caveats (i), (iii), (iv) of [11] §XI remain — separate publication tasks per BL-24]
X.3. Forward programme
Численное усиление двух-путевого Бианки на FLRW открывает следующие направления: (а) анизотропная Бьянки I/V/VII0-космология (нарушение пространственной изотропии); (б) перестабилизированные backgrounds с осцилляциями H(t) для тестирования предельных режимов конечно-разностного ρ˙; (в) расширение анти-циркулярного программного аудита на C.T1 (Φ-самосогласованность) и C.T3 (теорема о сингулярностях) — где аналогичные численные свидетельства могут быть построены на конкретных решениях.
БЛАГОДАРНОСТИ И ИНСТРУМЕНТЫ
% [FACT: tools and external assistance disclosure per L-33]
Автор благодарит сообщество ODTOE-исследователей за обсуждение аналитической структуры открытой задачи [11] §XI item (ii) и оговорки [12] §VIII.3, что мотивировало создание настоящей численной верификации. Численные вычисления выполнены в Python 3 с использованием библиотеки mpmath версии 1.3.0 (50-значная произвольная арифметика). Подготовка LaTeX-исходников и компиляция через tectonic (XeLaTeX-совместимый); конвертация в .docx через pandoc; конвертация в .md через утилиту tex2md.py корпуса ODTOE. Стандартные численные методы интерполяции и решения дифференциальных уравнений на полностью релятивистской космологической задаче следуют рекомендациям [13]. Исходный код flrw\_path2\_verification.py распространяется как часть корпуса.
КОНФЛИКТ ИНТЕРЕСОВ
Автор заявляет об отсутствии конфликта интересов.
ФИНАНСИРОВАНИЕ
Настоящее исследование не получало внешнего финансирования. Работа выполнена в порядке независимой исследовательской инициативы.
СПИСОК ЛИТЕРАТУРЫ
Замечание о порядке. Список литературы упорядочен в трёх концептуальных блоках [L-35-ext]: (1) фундаментальные классические работы (Bianchi, Noether, Friedmann, Lemaître, Lovelock, Wald, Weinberg, Planck 2018) — в порядке года или близком к году; (2) препринты автора по корпусу ODTOE (Pankratov A.S.) — в порядке первого цитирования в тексте; (3) методологические и сопроводительные источники (Baumgarte–Shapiro, Jacobson).
% --- БЛОК 1: ФУНДАМЕНТАЛЬНЫЕ РАБОТЫ ---
Bianchi, L. Lezioni di Geometria Differenziale, vols. I–III, 2nd ed. Spoerri, Pisa (1902). (Тождества Бианки.)
Noether, E. Invariante Variationsprobleme. Nachr. v.d. Ges. d. Wiss. zu Göttingen, math.-phys. Klasse, 235–257 (1918). EN translation: Tavel, M.A. Invariant variation problems. Transport Theory and Statistical Physics 1, 186–207 (1971). DOI: 10.1080/00411457108231446.
Friedmann, A. Über die Krümmung des Raumes. Z. Phys. 10, 377–386 (1922). DOI: 10.1007/BF01332580.
Lemaître, G. Un univers homogène de masse constante et de rayon croissant rendant compte de la vitesse radiale des nébuleuses extra-galactiques. Annales Soc. Sci. Bruxelles A47, 49–59 (1927).
Lovelock, D. The Einstein tensor and its generalizations. J. Math. Phys. 12(3), 498–501 (1971). DOI: 10.1063/1.1665613.
Wald, R.M. General Relativity. The University of Chicago Press (1984). ISBN: 0-226-87033-2.
Weinberg, S. Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity. Wiley, New York (1972). ISBN: 978-0-471-92567-5.
% --- БЛОК 2: ПРЕПРИНТЫ АВТОРА ПО КОРПУСУ ODTOE (порядок по первому цитированию) ---
Панкратов, А. С. Тензорная структура гравитации в ODTOE. Препринт (2026). Slug: ODTOE\_gravity\_tensor\_structure.
Панкратов, А. С. *Тензор энергии-импульса Tμν∗ и космологическая постоянная Λ из когерентности наблюдателя в ODTOE}. Препринт (2026). Slug: ODTOE\_gravity\_T\_munu\_projector.
Панкратов, А. С. Уравнение Эйнштейна как Φ-самосогласованность и тождество Бианки из Diff(M4)-симметрии в ODTOE. Препринт (2026). Slug: ODTOE\_einstein\_derivation\_complete.
Панкратов, А. С. Полное замыкание программы §XIV.3: уравнение Эйнштейна как Φ-самосогласованность. Препринт (2026). Slug: ODTOE\_einstein\_full\_closure.
% --- БЛОК 3: МЕТОДОЛОГИЧЕСКИЕ И СОПРОВОДИТЕЛЬНЫЕ ИСТОЧНИКИ ---
Baumgarte, T.W., Shapiro, S.L. Numerical Relativity: Solving Einstein's Equations on the Computer. Cambridge University Press (2010). ISBN: 978-0-521-51407-1.
Jacobson, T. Thermodynamics of spacetime: The Einstein equation of state. Phys. Rev. Lett. 75(7), 1260–1263 (1995). DOI: 10.1103/PhysRevLett.75.1260.