Derivasi Presesi Perihelion
DERIVASI RUMUS PRESESI PERIHELION
I. PENDAHULUAN
Dokumen ini bertujuan menurunkan secara sistematis rumus presesi perihelion partikel uji yang mengorbit objek pusat simetris bola dalam kerangka Medan Prasion. Derivasi dilakukan dengan menggunakan metrik yang dibangun dari dua parameter fundamental:- Medan Prasion (ᛟ) – medan fundamental yang merupakan manifestasi realitas pra-spasial dalam ruang-waktu emergent. Medan ini bervariasi terhadap posisi dan menjadi sumber seluruh dinamika fisika.
- Signature Modulasi (❀) – parameter tak berdimensi yang mengkarakterisasi kekuatan modulasi medan ᛟ oleh suatu objek. Nilainya terbatas pada rentang: 0 ≤ ❀ ≤ 2 dengan ❀ = 0 berarti tidak ada modulasi (objek tidak ada), dan ❀ = 2 merupakan batas maksimum yang diizinkan oleh prinsip kausalitas.
II. METRIK MEDAN PRASION DENGAN SIGNATURE MODULASI ❀
A. Definisi Medan Prasion untuk Objek Simetris Bola
Untuk objek masif dengan simetri bola sempurna, berjari-jari R, dan Signature Modulasi ❀, medan Prasion pada jarak radial r dari pusat objek (untuk r ≥ R) didefinisikan sebagai: \boxed{ᛟ(r) = \left(1 + \frac{❀ R}{2r}\right)^{-1/2}} \tag{1} Penjelasan dan sifat-sifat:- ᛟ(r) adalah besaran tak berdimensi yang merepresentasikan "intensitas keberadaan" medan Prasion pada titik tersebut.
- Nilai ᛟ = 1 di tak hingga (r → ∞) sebagai referensi dasar (void kosmik).
- Nilai ᛟ < 1 di sekitar objek, menunjukkan bahwa potensialitas telah teraktualisasi menjadi materi.
- Signature Modulasi ❀ mengkarakterisasi kekuatan modulasi: semakin besar ❀, semakin dalam "cekungan" medan ᛟ.
- Parameter R adalah radius objek, yang menjadi skala panjang karakteristik.
- Fungsi ᛟ(r) bersifat monoton naik terhadap r (turunannya positif) dan selalu positif untuk semua r > 0.
B. Metrik Ruang-Waktu Emergent
Ruang-waktu fisika muncul sebagai manifestasi dari distribusi medan ᛟ. Metrik yang konsisten dengan definisi ᛟ sebagai rasio waktu proper terhadap waktu koordinat (dτ = ᛟ ,dt untuk pengamat yang diam terhadap objek) adalah: \boxed{ds^2 = -c^2 ᛟ(r)^2 dt^2 + ᛟ(r)^{-2} dr^2 + r^2 d\Omega^2} \tag{3} dengan:- c = kecepatan cahaya dalam vakum,
- dΩ² = dθ² + sin²θ dφ² = elemen sudut pada bola satuan.
C. Komponen Metrik dan Notasi
Untuk memudahkan penurunan selanjutnya, kita definisikan komponen metrik yang relevan untuk gerak pada bidang ekuator (θ = π/2): g_{tt} = -c^2 A(r), \quad g_{rr} = B(r), \quad g_{\phi\phi} = r^2 \tag{5} dengan: A(r) = ᛟ(r)² = (1 + ❀R/2r)⁻¹ (6)B(r) = ᛟ(r)⁻² = 1 + ❀R/2r (7) Perhatikan hubungan sederhana namun penting antara kedua komponen: \boxed{A(r) \cdot B(r) = 1} \tag{8} Sifat ini akan menyederhanakan banyak perhitungan dalam penurunan persamaan orbit.
D. Ekspansi Deret untuk Analisis Medan Lemah
Untuk keperluan analisis medan lemah (yakni ketika ❀R ≪ a , dengan a adalah sumbu semi-mayor orbit), kita memerlukan ekspansi deret A(r) dalam kebalikan jarak 1/r. Ekspansi hingga orde yang diperlukan (orde 1/r² atau dalam variabel u = 1/r nanti hingga u²) adalah: A(r) = \left(1 + \frac{❀ R}{2r}\right)^{-1} = 1 - \frac{❀ R}{2r} + \frac{❀^2 R^2}{4r^2} - \frac{❀^3 R^3}{8r^3} + \cdots \tag{9} Sedangkan B(r) sudah dalam bentuk linear dalam 1/r sehingga tidak perlu diekspansi lebih lanjut: B(r) = 1 + \frac{❀ R}{2r} \quad \text{(bentuk eksak)} \tag{10} Dalam penurunan nanti, kita akan bekerja dengan variabel u = 1/r. Ekspansi A dalam u menjadi: A(u) = 1 - \frac{❀ R}{2} u + \frac{❀^2 R^2}{4} u^2 + \mathcal{O}(u^3) \tag{11} Turunan A terhadap u yang akan diperlukan adalah: \frac{dA}{du} = -\frac{❀ R}{2} + \frac{❀^2 R^2}{2} u + \mathcal{O}(u^2) \tag{12} Dengan fondasi metrik ini, kita siap untuk menurunkan persamaan geodesik dan menghitung presesi perihelion pada bagian-bagian selanjutnya.III. PERSAMAAN GEODESIK DAN KONSTANTA GERAK
A. Lagrangian Partikel Uji
Untuk partikel uji bermassa yang bergerak dalam ruang‑waktu yang dimetrikkan oleh (4), lintasannya merupakan kurva geodesik. Kita gunakan waktu proper τ sebagai parameter kurva. Lagrangian partikel (per satuan massa) dapat ditulis sebagai: \mathcal{L} = \frac{1}{2} g_{\mu\nu} \frac{dx^\mu}{d\tau} \frac{dx^\nu}{d\tau}. Dengan metrik (4) dan dengan membatasi gerak pada bidang ekuator (θ = π/2, sehingga θ̇=0 dan sin²θ = 1), Lagrangian menjadi: \mathcal{L} = \frac{1}{2} \left[ -c^2 A(r) \dot{t}^2 + B(r) \dot{r}^2 + r^2 \dot{\phi}^2 \right], \tag{13} di mana titik menyatakan turunan terhadap τ: ṫ = dt/dτ, ṙ = dr/dτ, φ̇ = dφ/dτ.B. Konstanta Gerak dari Simetri
Metrik (4) tidak bergantung secara eksplisit pada koordinat t dan φ. Hal ini mengindikasikan adanya dua vektor Killing, ∂t dan ∂φ, yang masing‑masing menghasilkan konstanta gerak sepanjang geodesik. · Energi per satuan massa (dari simetri waktu): \frac{\partial \mathcal{L}}{\partial \dot{t}} = -c^2 A(r) \dot{t} \equiv -E \quad \text{(konstan)}. Tanda negatif dipilih agar E positif untuk partikel yang bergerak maju dalam waktu. Dengan demikian: \boxed{E = c^2 A(r) \dot{t}} \tag{14} · Momentum sudut per satuan massa (dari simetri azimuth): \frac{\partial \mathcal{L}}{\partial \dot{\phi}} = r^2 \dot{\phi} \equiv h \quad \text{(konstan)}. Maka: \boxed{h = r^2 \dot{\phi}} \tag{15} Kedua konstanta ini akan sangat berguna dalam mereduksi persamaan gerak.C. Normalisasi Metrik
Untuk partikel bermassa, vektor kecepatan memenuhi kondisi normalisasi: g_{\mu\nu} \frac{dx^\mu}{d\tau} \frac{dx^\nu}{d\tau} = -c^2. Substitusi komponen metrik dan pembatasan gerak pada bidang ekuator memberikan: -c^2 A(r) \dot{t}^2 + B(r) \dot{r}^2 + r^2 \dot{\phi}^2 = -c^2. \tag{16}D. Persamaan Radial
Kita sekarang menggunakan konstanta gerak (14) dan (15) untuk mengeliminasi ṫ dan φ̇ dari persamaan normalisasi (16). Dari (14): ṫ = E / (c² A(r)). Dari (15): φ̇ = h / r². Substitusikan ke (16): -c^2 A(r) \left( \frac{E}{c^2 A(r)} \right)^2 + B(r) \dot{r}^2 + r^2 \left( \frac{h}{r^2} \right)^2 = -c^2. Sederhanakan suku pertama: -c^2 A(r) \cdot \frac{E^2}{c^4 A(r)^2} = -\frac{E^2}{c^2 A(r)}. Maka persamaan menjadi: -\frac{E^2}{c^2 A(r)} + B(r) \dot{r}^2 + \frac{h^2}{r^2} = -c^2. Pindahkan suku yang tidak mengandung ṙ ke ruas kanan: B(r) \dot{r}^2 = \frac{E^2}{c^2 A(r)} - \frac{h^2}{r^2} - c^2. Kita dapat menulis ulang sebagai: \boxed{B(r) \dot{r}^2 = \frac{E^2}{c^2 A(r)} - \left(c^2 + \frac{h^2}{r^2}\right)} \tag{17} Persamaan (17) adalah persamaan radial yang mengatur gerak partikel dalam koordinat r dan waktu proper τ. Sisi kanan hanya bergantung pada r (melalui A(r) dan suku h²/r²) serta konstanta E dan h. Ini merupakan langkah awal untuk mendapatkan bentuk orbit r(φ).IV. PERSAMAAN RADIAL DALAM VARIABEL u = 1/r
A. Transformasi ke Variabel u
Untuk mempelajari bentuk orbit, lebih mudah menggunakan variabel u = 1/r sebagai fungsi sudut φ. Hal ini karena φ adalah koordinat siklus yang baik dan kita memiliki kekekalan momentum sudut. Pertama, kita nyatakan ṙ dalam turunan terhadap φ: \dot{r} = \frac{dr}{d\tau} = \frac{dr}{d\phi} \frac{d\phi}{d\tau} = \frac{dr}{d\phi} \dot{\phi}. Dari (15), φ̇ = h/r² = h u². Jadi: \dot{r} = \frac{dr}{d\phi} \cdot h u^2. Selanjutnya, kita ubah dr/dφ menjadi turunan u terhadap φ. Karena r = 1/u, maka: \frac{dr}{d\phi} = -\frac{1}{u^2} \frac{du}{d\phi}. Sehingga: \dot{r} = -\frac{1}{u^2} \frac{du}{d\phi} \cdot h u^2 = -h \frac{du}{d\phi}. \tag{18} Hubungan (18) sangat penting: ia menghubungkan kecepatan radial dengan turunan u terhadap φ.B. Substitusi ke Persamaan Radial
Kita substitusikan (18) ke dalam persamaan radial (17). Perhatikan bahwa B(r) = B(u) dan A(r) = A(u) adalah fungsi dari u melalui r = 1/u. Dari (18), ṙ² = h² (du/dφ)². Maka (17) menjadi: B(u) \cdot h^2 \left(\frac{du}{d\phi}\right)^2 = \frac{E^2}{c^2 A(u)} - \left(c^2 + h^2 u^2\right). Bagi kedua ruas dengan h² (diasumsikan h ≠ 0 untuk orbit non-radial): B(u) \left(\frac{du}{d\phi}\right)^2 = \frac{E^2}{c^2 h^2 A(u)} - \frac{c^2}{h^2} - u^2. \tag{19} Ini adalah persamaan diferensial orde pertama untuk u(φ). Namun, untuk mendapatkan persamaan orde kedua yang lebih mudah dianalisis, kita akan mendiferensiasikannya terhadap φ.C. Penurunan Persamaan Orde Kedua
Tulis (19) sebagai: \left(\frac{du}{d\phi}\right)^2 = \frac{1}{B(u)} \left[ \frac{E^2}{c^2 h^2 A(u)} - \frac{c^2}{h^2} - u^2 \right]. \tag{20} Misalkan F(u) menyatakan ruas kanan. Maka: \left(\frac{du}{d\phi}\right)^2 = F(u). Diferensialkan kedua ruas terhadap φ: 2 \frac{du}{d\phi} \frac{d^2u}{d\phi^2} = \frac{dF}{du} \frac{du}{d\phi}. Untuk titik-titik di mana du/dφ ≠ 0, kita dapat membagi dengan 2 \, du/dφ: \frac{d^2u}{d\phi^2} = \frac{1}{2} \frac{dF}{du}. \tag{21} Sekarang kita hitung dF/du dari ekspresi F(u). Ingat bahwa A(u) dan B(u) adalah fungsi yang diketahui, dan hubungan A(u) B(u) = 1 akan menyederhanakan perhitungan. Tulis: F(u) = \frac{1}{B(u)} \left[ \frac{E^2}{c^2 h^2 A(u)} - \frac{c^2}{h^2} - u^2 \right]. Karena B(u) = 1/A(u), maka 1/B(u) = A(u). Jadi: F(u) = A(u) \left[ \frac{E^2}{c^2 h^2 A(u)} - \frac{c^2}{h^2} - u^2 \right] = \frac{E^2}{c^2 h^2} - A(u) \left( \frac{c^2}{h^2} + u^2 \right). \tag{22} Bentuk ini jauh lebih sederhana! Perhatikan bahwa suku pertama konstan terhadap u, sehingga turunannya mudah. Sekarang kita diferensialkan F(u) terhadap u: \frac{dF}{du} = -\frac{dA}{du} \left( \frac{c^2}{h^2} + u^2 \right) - A(u) \cdot 2u. Maka dari (21): \frac{d^2u}{d\phi^2} = \frac{1}{2} \left[ -\frac{dA}{du} \left( \frac{c^2}{h^2} + u^2 \right) - 2u A(u) \right]. Atau: \boxed{\frac{d^2u}{d\phi^2} = -\frac{1}{2} \frac{dA}{du} \left( \frac{c^2}{h^2} + u^2 \right) - u A(u)} \tag{23} Persamaan (23) adalah persamaan diferensial orde kedua yang mengatur bentuk orbit u(φ). Ini adalah persamaan eksak dalam kerangka metrik (4). semua bergantung pada fungsi A(u) yang diberikan oleh (6) dalam variabel u: A(u) = \left(1 + \frac{❀ R}{2} u\right)^{-1}. \tag{24}D. Bentuk dalam Parameter Prasion
Untuk kejelasan, kita tulis ulang (23) dengan A(u) dari (24). Namun, untuk analisis selanjutnya kita akan menggunakan ekspansi deret karena kita tertarik pada orbit yang hampir Newtonian (medan lemah). Ekspansi A(u) dan turunannya telah disiapkan di bagian II. Persamaan (23) akan menjadi dasar bagi seluruh perhitungan presesi. Pada bagian berikut, kita akan melakukan ekspansi untuk medan lemah (❀R ≪ a) dan menyederhanakannya hingga orde yang relevan untuk menangkap efek presesi.V. PERSAMAAN DIFERENSIAL ORDE DUA DAN EKSPANSI MEDAN LEMAH
A. Persamaan Eksak untuk Orbit
Kita telah memperoleh persamaan diferensial orde kedua yang mengatur bentuk orbit u(φ) dari partikel uji: \frac{d^2u}{d\phi^2} = -\frac{1}{2} \frac{dA}{du} \left( \frac{c^2}{h^2} + u^2 \right) - u A(u), \tag{23} dengan A(u) adalah komponen metrik yang dalam kerangka Medan Prasion diberikan oleh: A(u) = \left(1 + \frac{❀ R}{2} u\right)^{-1}. \tag{24} Persamaan (23) bersifat eksak dan berlaku untuk seluruh rentang nilai u yang sesuai dengan r ≥ R. Namun, untuk memperoleh solusi yang dapat diinterpretasikan sebagai orbit planet, kita perlu melakukan ekspansi yang sesuai dengan kondisi medan lemah, yaitu ketika pengaruh modulasi relatif kecil terhadap gerak Newtonian.B. Kondisi Medan Lemah
Medan lemah didefinisikan oleh kondisi: \frac{❀ R}{2} u \ll 1 \quad \text{atau} \quad ❀ R \ll r. Untuk orbit planet, r berorde sumbu semi-mayor a, sehingga kondisi ini setara dengan: ❀ R \ll a. Dalam kondisi ini, kita dapat memperlakukan suku yang mengandung ❀ sebagai koreksi kecil terhadap gerak Newtonian. Ekspansi deret A(u) dan turunannya hingga orde yang diperlukan menjadi sah.C. Ekspansi Deret A(u) dan Turunannya
Dari (24), kita ekspansikan A(u) dalam deret pangkat u (yang setara dengan deret dalam 1/r): A(u) = \left(1 + \frac{❀ R}{2} u\right)^{-1} = 1 - \frac{❀ R}{2} u + \frac{❀^2 R^2}{4} u^2 - \frac{❀^3 R^3}{8} u^3 + \cdots \tag{25} Turunan pertama terhadap u diperlukan dalam persamaan (23): \frac{dA}{du} = -\frac{❀ R}{2} \left(1 + \frac{❀ R}{2} u\right)^{-2} = -\frac{❀ R}{2} + \frac{❀^2 R^2}{2} u - \frac{3❀^3 R^3}{8} u^2 + \cdots \tag{26} Untuk analisis presesi hingga orde pertama dalam ❀ (yang akan memberikan koreksi relativistik), kita perlu mempertahankan:- A(u) hingga orde u² (karena suku u² akan berkontribusi pada persamaan melalui perkalian dengan u dan melalui suku dA/du).
- dA/du hingga orde u (karena akan dikalikan dengan c²/h² + u²).
D. Substitusi ke Persamaan (23)
Kita substitusikan ekspansi (25) dan (26) ke dalam (23). Pertama, hitung suku pertama di ruas kanan: -\frac{1}{2} \frac{dA}{du} \left( \frac{c^2}{h^2} + u^2 \right) = -\frac{1}{2} \left( -\frac{❀ R}{2} + \frac{❀^2 R^2}{2} u + \cdots \right) \left( \frac{c^2}{h^2} + u^2 \right). Kalikan dan pertahankan suku hingga orde ❀ u² (dengan mengingat bahwa u sendiri bisa kecil, tetapi kita ingin semua suku yang memberikan kontribusi pada persamaan akhir): = \frac{❀ R}{4} \left( \frac{c^2}{h^2} + u^2 \right) - \frac{❀^2 R^2}{4} u \left( \frac{c^2}{h^2} + u^2 \right) + \cdots Ekspansi lebih lanjut: = \frac{❀ R}{4} \cdot \frac{c^2}{h^2} + \frac{❀ R}{4} u^2 - \frac{❀^2 R^2}{4} \cdot \frac{c^2}{h^2} u - \frac{❀^2 R^2}{4} u^3 + \cdots \tag{27} Suku kedua di ruas kanan (23) adalah: - u A(u) = -u \left( 1 - \frac{❀ R}{2} u + \frac{❀^2 R^2}{4} u^2 - \cdots \right) = -u + \frac{❀ R}{2} u^2 - \frac{❀^2 R^2}{4} u^3 + \cdots \tag{28}E. Penggabungan dan Penyederhanaan
Sekarang kita jumlahkan (27) dan (28) untuk mendapatkan ruas kanan lengkap persamaan (23): \frac{d^2u}{d\phi^2} = \left[ \frac{❀ R}{4} \cdot \frac{c^2}{h^2} + \frac{❀ R}{4} u^2 - \frac{❀^2 R^2}{4} \cdot \frac{c^2}{h^2} u - \frac{❀^2 R^2}{4} u^3 \right] + \left[ -u + \frac{❀ R}{2} u^2 - \frac{❀^2 R^2}{4} u^3 \right] + \cdots Kelompokkan suku-suku sejenis:- Suku konstan (tidak bergantung u): ❀R/4 ⋅ c²/h²
- Suku linear dalam u: -u - (❀²R²/4) ⋅ (c²/h²)u
- Suku kuadrat dalam u: ❀R/4 u² + ❀R/2 u² = 3❀R/4 u²
- Suku kubik dalam u: -❀²R²/4 u³ - ❀²R²/4 u³ = -❀²R²/2 u³
F. Pemertahanan Suku yang Relevan
Untuk analisis presesi hingga orde pertama dalam ❀ (yang merupakan koreksi relativistik terendah), kita hanya perlu mempertahankan suku-suku yang:- Berorde ❀ (yaitu linear dalam ❀)
- Memberikan kontribusi pada pergeseran fase orbit.
G. Interpretasi Persamaan (30)
Persamaan (30) memiliki struktur sebagai berikut: \frac{d^2u}{d\phi^2} + u = \alpha + \beta u^2, dengan: \alpha = \frac{c^2 ❀ R}{4h^2}, \quad \beta = \frac{3❀ R}{4}. Jika suku β u² tidak ada, persamaan ini akan menghasilkan solusi elips Newtonian dengan parameter p = 1/α (setelah disesuaikan). Suku β u² adalah koreksi relativistik yang akan menyebabkan orbit tidak lagi tertutup sempurna, melainkan mengalami presesi.VI. SOLUSI ORDE NOL DAN PERSAMAAN KOREKSI
A. Solusi Orde Nol (Newtonian)
Kita mulai dengan menyelesaikan persamaan (30) dengan mengabaikan suku non-linear u² untuk mendapatkan aproksimasi orde nol. Persamaan orde nol adalah: \frac{d^2u_0}{d\phi^2} + u_0 = \frac{c^2 ❀ R}{4h^2}. \tag{31} Ini adalah persamaan diferensial linear orde dua dengan koefisien konstan dan suku konstan di ruas kanan. Solusi umumnya terdiri dari solusi homogen ditambah solusi partikular.- Solusi homogen: u_h = C cos(φ - φ₀).
- Solusi partikular: Karena ruas kanan konstan, kita coba u_p = konstan. Substitusi u_p = K ke (31) memberikan K = c² ❀ R / 4h².
B. Parameter Orbit Newtonian
Kita definisikan parameter p (semi-latus rektum) sebagai: p = \frac{4h^2}{c^2 ❀ R}. \tag{33} Maka solusi (32) menjadi: u_0(\phi) = \frac{1}{p} \left(1 + e \cos \phi\right). \tag{34} Untuk orbit elips, hubungan antara p, sumbu semi-mayor a, dan eksentrisitas e adalah: p = a(1 - e^2). \tag{35} Hubungan ini akan digunakan nanti untuk menyatakan hasil akhir dalam elemen orbit yang umum digunakan.C. Persamaan untuk Koreksi Orde Pertama
Sekarang kita mencari solusi yang lebih akurat dengan memperhitungkan suku u² dalam (30). Kita tulis: u(\phi) = u_0(\phi) + u_1(\phi), dengan u_1 adalah koreksi kecil (|u_1| ≪ |u_0|). Substitusikan ke persamaan (30): \frac{d^2(u_0 + u_1)}{d\phi^2} + (u_0 + u_1) = \frac{c^2 ❀ R}{4h^2} + \frac{3❀ R}{4} (u_0 + u_1)^2. Karena u_0 memenuhi (31), maka d²u₀/dφ² + u₀ = c² ❀ R / 4h². Kedua ruas saling menghilangkan, menyisakan: \frac{d^2u_1}{d\phi^2} + u_1 = \frac{3❀ R}{4} (u_0^2 + 2u_0 u_1 + u_1^2). Dalam aproksimasi linear, kita abaikan suku yang mengandung u₁² dan juga suku u₀ u₁ karena produk dua besaran kecil. Namun, kita harus berhati-hati: meskipun u₁ kecil, perkalian u₀ u₁ bisa seorde dengan suku lain? Untuk konsistensi hingga orde pertama dalam parameter kecil (❀), kita hanya perlu mempertahankan suku u₀² di ruas kanan. Suku u₀ u₁ berorde ❀ dikali u₁, sehingga lebih tinggi ordenya. Jadi: \frac{d^2u_1}{d\phi^2} + u_1 = \frac{3❀ R}{4} u_0^2. \tag{36}D. Ekspansi u₀²
Dari (34), kita hitung u₀²: u_0^2(\phi) = \frac{1}{p^2} \left(1 + e \cos \phi\right)^2 = \frac{1}{p^2} \left(1 + 2e \cos \phi + e^2 \cos^2 \phi\right). Gunakan identitas cos² φ = 1/2 (1 + cos 2φ): u_0^2 = \frac{1}{p^2} \left[1 + 2e \cos \phi + \frac{e^2}{2} (1 + \cos 2\phi)\right] = \frac{1}{p^2} \left[ \left(1 + \frac{e^2}{2}\right) + 2e \cos \phi + \frac{e^2}{2} \cos 2\phi \right]. \tag{37} Jadi ruas kanan persamaan (36) menjadi: \frac{3❀ R}{4} u_0^2 = \frac{3❀ R}{4p^2} \left[ \left(1 + \frac{e^2}{2}\right) + 2e \cos \phi + \frac{e^2}{2} \cos 2\phi \right]. \tag{38}E. Identifikasi Suku Resonan
Persamaan (36) adalah osilator harmonik terpaksa dengan gaya penggerak periodik. Solusi homogennya adalah cos φ dan sin φ. Jika gaya penggerak mengandung suku yang sebanding dengan cos φ (atau sin φ), akan terjadi resonansi yang menghasilkan solusi dengan amplitudo yang tumbuh linear terhadap φ. Inilah yang menyebabkan pergeseran fase sekuler (presesi). Dari (38), suku yang mengandung cos φ adalah: \frac{3❀ R}{4p^2} \cdot 2e \cos \phi = \frac{3❀ R e}{2p^2} \cos \phi. \tag{39} Suku ini adalah penyebab utama presesi. Suku konstan dan suku cos 2φ hanya menghasilkan osilasi periodik yang tidak berkontribusi pada pergeseran fase jangka panjang.F. Solusi Partikular untuk Suku Resonan
Kita selesaikan persamaan: \frac{d^2u_1}{d\phi^2} + u_1 = \frac{3❀ R e}{2p^2} \cos \phi. \tag{40} Karena cos φ adalah solusi homogen, kita gunakan metode koefisien tak tentu dengan bentuk u₁ₚ = D φ sin φ (bentuk yang umum untuk kasus resonansi). Substitusikan ke (40): \frac{d}{d\phi}(D \phi \sin \phi) = D (\sin \phi + \phi \cos \phi), \frac{d^2}{d\phi^2}(D \phi \sin \phi) = D (\cos \phi + \cos \phi - \phi \sin \phi) = D (2 \cos \phi - \phi \sin \phi). Maka: \frac{d^2u_{1p}}{d\phi^2} + u_{1p} = D (2 \cos \phi - \phi \sin \phi) + D \phi \sin \phi = 2D \cos \phi. Agar sama dengan ruas kanan (40), kita perlukan: 2D = \frac{3❀ R e}{2p^2} \quad \Rightarrow \quad D = \frac{3❀ R e}{4p^2}. Jadi solusi partikular untuk bagian resonan adalah: u_{1,\text{res}} = \frac{3❀ R e}{4p^2} \, \phi \sin \phi. \tag{41}G. Solusi Lengkap dan Interpretasi Fisis
Solusi lengkap hingga orde pertama adalah u(φ) = u₀(φ) + u₁,res(φ) (kita abaikan bagian non-resonan karena hanya memberikan osilasi kecil). Jadi: u(\phi) = \frac{1}{p} \left(1 + e \cos \phi\right) + \frac{3❀ R e}{4p^2} \, \phi \sin \phi. Faktorkan 1/p : u(\phi) = \frac{1}{p} \left[1 + e \cos \phi + \frac{3❀ R}{4p} e \, \phi \sin \phi \right]. \tag{42} Untuk 3❀ R / 4p ≪ 1, kita dapat menulis ulang suku dalam kurung sebagai bentuk kosinus dengan fase yang bergeser. Gunakan identitas trigonometri untuk sudut kecil: \cos\left((1-\delta)\phi\right) \approx \cos \phi + \delta \phi \sin \phi, dengan δ ≪ 1. Dengan membandingkan, kita lihat bahwa: e \cos\left((1-\delta)\phi\right) \approx e \cos \phi + e \delta \phi \sin \phi. Agar sesuai dengan (42), kita perlukan: e \delta = \frac{3❀ R}{4p} e \quad \Rightarrow \quad \delta = \frac{3❀ R}{4p}. Dengan demikian, solusi (42) dapat ditulis sebagai: u(\phi) \approx \frac{1}{p} \left[1 + e \cos\left((1-\delta)\phi\right)\right]. \tag{43}H. Pergeseran Perihelion
Bentuk (43) menunjukkan bahwa sudut φ yang diperlukan untuk kembali ke nilai u yang sama (yaitu kembali ke jarak radial yang sama) tidak lagi 2π, melainkan φ harus bertambah hingga cos kembali ke nilai awalnya. Perihelion terjadi ketika cos((1-δ)φ) = 1, yaitu ketika: (1-\delta)\phi = 2\pi n \quad \Rightarrow \quad \phi = \frac{2\pi n}{1-\delta} \approx 2\pi n (1 + \delta). Jadi setelah satu putaran (n=1), sudut yang ditempuh adalah 2π(1 + δ). Kelebihan sudut ini adalah presesi per revolusi: \Delta\phi = 2\pi \delta = 2\pi \cdot \frac{3❀ R}{4p}. Sederhanakan: \Delta\phi = \frac{3\pi ❀ R}{2p}. \tag{44}I. Rumus Akhir dalam Elemen Orbit
Dengan menggunakan hubungan p = a(1 - e²) dari (35), kita peroleh rumus final presesi perihelion: \boxed{\Delta\phi = \frac{3\pi ❀ R}{2a(1-e^2)} \quad \text{(radian per revolusi)}} \tag{45} Rumus ini hanya bergantung pada Signature Modulasi ❀, radius objek pusat R, serta elemen orbit partikel uji (a dan e). Tidak terdapat konstanta fisika lain selain yang sudah terdefinisi dalam kerangka Medan Prasion.VII. DISKUSI
A. Makna Fisis Rumus Presesi
Rumus presesi perihelion yang diperoleh: \Delta\phi = \frac{3\pi ❀ R}{2a(1-e^2)} memiliki beberapa interpretasi penting dalam kerangka Medan Prasion:- Ketergantungan pada Signature Modulasi ❀: Presesi berbanding lurus dengan ❀, yang merupakan ukuran seberapa kuat objek pusat memodulasi medan Prasion di sekitarnya. Semakin besar ❀, semakin dalam "cekungan" medan ᛟ yang diciptakan objek, dan semakin besar efek relativistik yang teramati sebagai presesi.
- Peran radius R: Kehadiran R dalam pembilang menunjukkan bahwa presesi bergantung pada ukuran fisik objek pusat, bukan hanya "kekuatan" modulasinya. Dua objek dengan ❀ sama tetapi radius berbeda akan menghasilkan presesi yang berbeda untuk orbit dengan parameter yang sama.
- Ketergantungan pada elemen orbit: Bentuk 1/[a(1 - e²)] = 1/p menunjukkan bahwa presesi lebih besar untuk orbit yang lebih kecil (a kecil) dan lebih eksentrik (e mendekati 1). Ini sesuai dengan intuisi fisika: efek relativistik paling kuat di dekat perihelion, sehingga orbit dengan eksentrisitas tinggi mengalami presesi lebih besar.
B. Batas Fundamental ❀ ≤ 2
Salah satu hasil terpenting dari kerangka Medan Prasion adalah keberadaan batas atas Signature Modulasi: 0 \le ❀ \le 2. Batas ini berasal dari prinsip kausalitas yang melarang kecepatan lepas melebihi kecepatan cahaya. Implikasinya terhadap presesi perihelion adalah adanya batas atas maksimum untuk setiap sistem: \Delta\phi_{\text{max}} = \frac{3\pi (2) R}{2a(1-e^2)} = \frac{3\pi R}{a(1-e^2)}. Ini berarti untuk objek pusat dengan radius R, tidak mungkin mengamati presesi yang melebihi nilai ini, apa pun konfigurasi internal objek tersebut. Batasan ini dapat diuji secara observasional, terutama untuk objek kompak seperti bintang neutron atau kandidat lubang hitam.C. Perbandingan dengan Metrik Linear
Sebelumnya, derivasi presesi juga telah dilakukan menggunakan metrik linear dengan medan Prasion ᛟ(r) = 1 - ❀R/4r. Hasilnya identik dengan (45) pada orde pertama. Perbedaan antara kedua metrik hanya muncul pada orde yang lebih tinggi (suku ❀² dan seterusnya). Ini berarti:- Untuk medan lemah (❀ ≪ 1), kedua metrik memberikan prediksi yang sama dan konsisten dengan observasi tata surya.
- Untuk medan kuat (❀ mendekati 2), metrik baru yang digunakan dalam derivasi ini lebih dapat diandalkan karena dibangun secara eksak dan konsisten dengan batas kecepatan lepas.
D. Implikasi untuk Objek Kompak
Untuk objek dengan ❀ mendekati 2, misalnya bintang neutron atau kandidat lubang hitam, presesi perihelion partikel yang mengorbit sangat dekat dengan permukaan dapat mendekati batas maksimum. Sebagai ilustrasi, ambil objek dengan radius R dan orbit dengan a ≈ R (orbit sangat dekat permukaan) dan e = 0 (melingkar). Maka: \Delta\phi_{\text{max}} \approx \frac{3\pi R}{R} = 3\pi \text{ radian} \approx 540^\circ per revolusi. Ini berarti orbit tidak lagi berupa elips yang mengalami presesi kecil, melainkan lintasan yang sangat berbeda dari deskripsi Newtonian. Pada kondisi ini, analisis orde pertama mungkin tidak lagi memadai dan diperlukan penanganan eksak dari persamaan (23).E. Keterujian Teori
Rumus presesi (45) memberikan beberapa prediksi yang dapat diuji:- Proporsionalitas dengan 1/[a(1-e²)]: Untuk planet-planet dalam satu sistem (misal tata surya), rasio presesi antara dua planet harus sama dengan rasio 1/[a(1-e²)] masing-masing. Ini sudah terkonfirmasi oleh mekanika langit.
- Batas atas presesi: Untuk objek dengan radius R yang diketahui, tidak mungkin mengamati presesi melebihi 3πR / a(1 - e²). Ini dapat diuji dengan mengamati orbit di sekitar bintang neutron atau objek kompak lainnya.
- Konsistensi internal: Nilai ❀ yang dihitung dari presesi untuk suatu objek harus konsisten dengan nilai ❀ yang diperoleh dari fenomena lain seperti pembelokan cahaya atau dilasi waktu. Jika pengukuran independen ini memberikan nilai ❀ yang berbeda, teori perlu dikaji ulang.
F. Keterbatasan Derivasi
Perlu dicatat bahwa derivasi ini memiliki beberapa keterbatasan:- Aproksimasi medan lemah: Kita mengasumsikan ❀R ≪ a sehingga ekspansi deret valid. Untuk objek sangat kompak dengan orbit sangat dekat, aproksimasi ini mungkin tidak cukup.
- Abaikan suku orde tinggi: Suku ❀² dan seterusnya diabaikan. Untuk presisi sangat tinggi atau untuk ❀ tidak terlalu kecil, koreksi ini perlu diperhitungkan.
- Objek pusat statis: Kita mengasumsikan objek pusat simetris bola dan statis. Efek rotasi (frame dragging) dan multipol tidak diperhitungkan.
- Partikel uji: Partikel dianggap tidak mempengaruhi medan (test particle). Untuk sistem biner dengan massa sebanding, diperlukan perlakuan lebih lanjut.
VIII. KESIMPULAN
A. Ringkasan Hasil
Telah diturunkan secara sistematis rumus presesi perihelion partikel uji yang mengorbit objek pusat simetris bola dalam kerangka Medan Prasion dengan Signature Modulasi ❀. Derivasi dimulai dari metrik eksak: ds^2 = -c^2 \left(1 + \frac{❀ R}{2r}\right)^{-1} dt^2 + \left(1 + \frac{❀ R}{2r}\right) dr^2 + r^2 d\Omega^2, yang konsisten dengan definisi medan Prasion ᛟ(r) = (1 + ❀R/(2r))⁻¹/². Langkah-langkah derivasi meliputi:- Perumusan persamaan geodesik dan identifikasi konstanta gerak (E dan h).
- Penurunan persamaan radial dan transformasi ke variabel u = 1/r.
- Diferensiasi untuk memperoleh persamaan orbit orde dua (23).
- Ekspansi deret untuk medan lemah (❀R ≪ a) menghasilkan persamaan (30): d²u/dφ² + u = c²❀R/4h² + 3❀R/4 u²
- Solusi orde nol memberikan orbit elips Newtonian: u₀ = 1/p (1 + e cos φ) dengan p = 4h²/(c² ❀R)
- Analisis perturbatif untuk koreksi orde pertama menghasilkan suku resonan yang menyebabkan pergeseran fase.
- Interpretasi pergeseran fase sebagai presesi perihelion dengan rumus akhir: Δφ = 3π❀R / 2a(1 - e²) (radian per revolusi)
B. Posisi dalam Kerangka Medan Prasion
Derivasi ini melengkapi pemahaman tentang bagaimana fenomena gravitasi muncul sebagai manifestasi dinamis dari medan Prasion. Presesi perihelion, yang dalam relativitas umum dijelaskan sebagai efek kelengkungan ruang-waktu, di sini muncul sebagai konsekuensi alami dari non-linearitas medan ᛟ yang termanifestasi dalam suku u² pada persamaan orbit. Dengan demikian, derivasi presesi perihelion ini tidak hanya merupakan latihan matematis, tetapi juga membuka jalan untuk pemahaman yang lebih tentang hubungan antara medan fundamental Prasion dan fenomena yang teramati.