Memahami Formulasi Simetris-Hiperbolik dalam Persamaan Medan Einstein

Persamaan medan Einstein berbentuk tensor. Ruas kiri adalah tensor Einstein nan mengandung lengkungan ruang-waktu. Ruas kanan adalah tensor energi-momentum nan mengandung pengedaran massa dan energi. Bentuknya Gμν = (8πG/c⁴) Tμν. Ini adalah sepuluh persamaan diferensial parsial nan tidak linear. Tidak linear artinya medan gravitasi mempengaruhi medan gravitasi itu sendiri. Di alam semesta nyata, tidak ada situasi nan simetri sempurna selain lubang hitam Schwarzschild alias alam semesta FLRW. Untuk tabrakan dua bintang neutron, corak medan gravitasi terus berubah setiap milidetik. Komputer kudu menyelesaikan persamaan ini secara numerik. Tapi masalah besar muncul, sistem persamaan ini tidak dirancang untuk perkembangan waktu. Dalam corak aslinya, hanya ada hambatan (constraints) nan kudu dipenuhi pada setiap lembaran ruang. Tidak ada persamaan perkembangan definitif untuk metrik. Karena itu, fisikawan kudu mengubah corak persamaan Einstein agar bisa dijalankan seperti mesin waktu digital. Inilah awal dari dekomposisi 3+1. Dekomposisi 3+1 Membagi Ruang-Waktu Menjadi Lapisan Dekomposisi 3+1 dimulai dengan menganggap bahwa ruang-waktu empat dimensi bisa dipotong menjadi lembaran-lembaran ruang tiga dimensi nan setiap lembaran punya waktu koordinat konstan. Setiap lembaran punya metrik tiga dimensi γij. Metrik ini mengukur jarak antar titik di dalam lembaran nan sama. Untuk menghubungkan satu lembaran ke lembaran berikutnya, diperlukan dua besaran. Fungsi lapse α menentukan berapa lama waktu proper nan berlalu antara dua lembaran nan berdekatan. Vektor shift βi menentukan gimana koordinat spasial bergeser dari satu lembaran ke lembaran berikutnya. Metrik ruang-waktu empat dimensi ditulis ulang dalam besaran-besaran ini. Bentuk intervalnya: ds² = –α² dt² + γij (dxi + βi dt)(dxj + βj dt). Perhatikan bahwa komponen waktu-waktu dari metrik invers juga berubah. Determinan metrik tiga dimensi γij diberi notasi γ. Besaran ini krusial untuk menghitung kerapatan energi. Para perintis ADM (Arnowitt, Deser, Misner) pada tahun 1959 menurunkan persamaan perkembangan untuk γij dan curvature ekstrinsik Kij. Curvature ekstrinsik ini mengukur gimana lembaran ruang melengkung relatif terhadap ruang-waktu empat dimensi di sekitarnya. Definisi Kij = – (1/(2α)) (∂t γij – ∇i βj – ∇j βi). Tanda negatif di sini adalah konvensi. Sekarang sistem ADM mempunyai dua belas variabel: enam komponen γij dan enam komponen Kij. Ada delapan persamaan hambatan nan kudu dipenuhi di setiap titik ruang. Kendala Hamiltonian: R + K² – Kij Kij = 0. Kendala momentum: ∇j (Kij – γij K) = 0. Sisanya adalah dua belas persamaan evolusi. Sistem ini bentuknya hiperbolik tetapi tidak simetris. Akibatnya, simulasi numerik dengan ADM original luntur dalam waktu singkat. Dua belas persamaan perkembangan dari formalisme ADM original nan diturunkan dari proyeksi tensor Riemann ke lembaran spasial, diantaranya, Turunan waktu dari metrik tiga dimensi γij: ∂t γij = –2α Kij + ∇i βj + ∇j βi. Turunan waktu dari curvature ekstrinsik Kij: ∂t Kij = –∇i ∇j α + α ( Rij – 2 Kik Kkj + K Kij ) + βk ∇k Kij + Kik ∇j βk + Kkj ∇i βk. Di sini indeks i dan j melangkah dari 1 sampai 3, sehingga ada 6 komponen untuk γij dan 6 komponen untuk Kij. Suku Rij adalah tensor Ricci tiga dimensi nan dihitung dari γij dan turunan pertama dan kedua spasialnya. K adalah trace dari Kij ialah K = γij Kij. Simbol ∇i menyatakan turunan kovarian tiga dimensi nan kompatibel dengan γij. Dua belas persamaan ini membentuk sistem perkembangan tertutup lantaran ruas kanan hanya berisi γij, Kij, α, βi, dan turunan spasialnya. Fungsi lapse α dan shift βi tidak mempunyai persamaan perkembangan sendiri dalam ADM murni dimana keduanya bebas dipilih sebagai gauge. Keenam persamaan untuk γij berkarakter kinematik, langsung mendefinisikan Kij dalam corak turunan waktu metrik. Keenam persamaan untuk Kij berkarakter dinamik, berasal dari persamaan Einstein untuk komponen ruang-ruang dari tensor Riemann empat dimensi nan diproyeksikan. Gejala Ketidakstabilan ADM Asli dan Akar Penyebabnya Ketika intelektual mencoba menjalankan ADM original di komputer Cray-1 dan Cray X-MP pada tahun 1980-an dan awal 1990-an, hasilnya buruk. Simulasi lubang hitam Schwarzschild nan semestinya stabil malah meledak setelah beberapa kali putaran waktu. Galat numerik nan mini di awal simulasi membesar secara eksponensial. Analisis matematis menunjukkan bahwa sistem ADM original tidak memenuhi kondisi hiperbolisitas kuat. Ada dua jenis hiperbolisitas nan dikenal dalam kajian persamaan diferensial parsial. Sistem lemah-hiperbolik hanya mempunyai nilai eigen riil dari matriks simbol tetapi tidak mempunyai golongan komplit vektor eigen. Sistem ini tidak stabil terhadap gangguan kecil. Sistem kuat-hiperbolik mempunyai nilai eigen riil dan golongan komplit vektor eigen. Sistem ini stabil dalam ruang Hilbert. Sistem simetris-hiperbolik adalah kasus unik di mana matriks simbol dapat dibuat simetris dengan metrik daya nan konstan positif. Sistem simetris-hiperbolik adalah nan paling stabil secara numerik. ADM original tidak termasuk dalam kategori mana pun lantaran persamaan evolusinya mengandung suku-suku derivatif orde dua nan tidak dapat dipisahkan dengan rapi. Secara teknis, persamaan ADM untuk Kij mengandung turunan spasial dari Kij dan turunan spasial dari γij nan saling mengenai melalui kendala. Ketika hambatan tidak terpenuhi secara eksak akibat galat numerik, sistem menghasilkan mode-mode tak-fisik nan tumbuh tanpa batas. Inilah nan disebut masalah stabilitas constraint-violation. Para intelektual butuh formulasi ulang nan secara otomatis menekan pelanggaran kendala. Formulasi BSSN Sebagai Perbaikan Pertama nan Sukses Pada tahun 1995, Shibata dan Nakamura di Jepang mengusulkan variabel baru. Kemudian Baumgarte dan Shapiro di Amerika Serikat menyadari potensi formulasi tersebut pada tahun 1998. Formulasi ini dikenal sebagai BSSN (Baumgarte–Shapiro–Shibata–Nakamura). Ide dasarnya adalah memisahkan metrik tiga dimensi menjadi aspek konformal dan bagian trace. Variabel baru nan diperkenalkan: ϕ = (1/12) ln(det γij). Faktor ini mengukur skala volume dari metrik. Metrik konformal γ̃ij = e⁻⁴φ γij mempunyai determinan satu. Selanjutnya curvature ekstrinsik dipisahkan menjadi bagian trace K dan bagian traceless Ãij = e⁻⁴φ (Kij – (1/3) γij K). Jumlah variabel menjadi lebih banyak tetapi persamaan perkembangan menjadi lebih sederhana. nan paling penting, BSSN memperkenalkan variabel bantu Γ̃i nan merupakan kontraksi dari turunan metrik konformal. Variabel ini bertindak sebagai peredam numerik. Persamaan perkembangan BSSN tidak lagi mengandung turunan kedua dari metrik secara eksplisit. Banyak suku nan sebelumnya menyebabkan mode tak-stabil sekarang berubah menjadi suku sumber alias suku adveksi. Simulasi lubang hitam tunggal dengan BSSN bisa melangkah hingga waktu koordinat melampaui seribu massa lubang hitam tanpa ledakan. Formulasi BSSN juga bisa menangani lubang hitam nan bergerak. Teknik moving puncture nan dikembangkan tahun 2005 oleh Campanelli, Lousto, Zlochower, dan lain-lain memanfaatkan BSSN dengan pilihan gauge 1+log untuk lapse dan gamma-driver untuk shift. BSSN sendiri berkarakter kuat-hiperbolik tetapi belum sepenuhnya simetris-hiperbolik dalam makna paling ketat. Namun untuk keperluan numerik, sifat kuat-hiperbolik sudah cukup. Formulasi Z4 nan Lebih Radikal dan Simetris Penuh Pada tahun 2003, Bona, Ledvinka, Palenzuela, dan Zacek memperkenalkan pendekatan berbeda. Mereka menambahkan empat variabel baru Zμ nan secara buatan dimasukkan ke dalam persamaan Einstein. Hasilnya adalah sistem Z4. Ide di kembali Z4 sederhana namun jenius. Kendala Hamiltonian dan momentum dalam ADM dapat ditulis sebagai vektor empat dimensi. Dalam Z4, vektor hambatan ini dipromosikan menjadi variabel bergerak nan berevolusi. Persamaan Einstein nan dimodifikasi berbentuk Gμν + ∇μ Zν + ∇ν Zμ – gμν ∇α Zα = 8π Tμν. Tanda dan aspek numerik disesuaikan. Hasilnya, sistem menjadi simetris-hiperbolik untuk pilihan parameter tertentu. Ada jenis Z4 nan disebut Z4c alias Z4 dengan redaman kendala. Parameter κ1 dan κ2 ditambahkan ke persamaan perkembangan Zμ untuk secara eksponensial meluruhkan pelanggaran kendala. Formulasi Z4 lebih unggul dari BSSN dalam simulasi nan sangat lama alias nan melibatkan vakum ekstrim. Misalnya, simulasi lubang hitam dengan rasio massa sangat besar (1000 banding 1) lebih stabil di Z4 daripada BSSN. Analisis spektral dari matriks simbol Z4 menunjukkan bahwa semua nilai eigen adalah riil dan terdapat metrik daya positif nan membikin sistem simetris. Sifat ini menjamin bahwa masalah Cauchy untuk Z4 adalah well-posed dalam ruang Sobolev. Artinya, solusi ada, unik, dan berjuntai kontinu pada info awal. Ini adalah agunan matematis tertinggi nan bisa diberikan untuk persamaan diferensial parsial non-linear. Hingga saat ini, belum ada formulasi lain nan mencapai tingkat simetris-hiperbolik semurni Z4 dengan hambatan damping. Saat ini, kode seperti GRChombo, Einstein Toolkit dengan thorn Z4c, dan IllinoisGRMHD menggunakan jenis Z4 untuk produksi ilmiah. Tiga formulasi sudah kita bahas: ADM original (tidak stabil), BSSN (kuat-hiperbolik, stabil untuk merger), dan Z4 (simetris-hiperbolik penuh, paling stabil secara matematis). Pilihan formulasi tergantung pada masalah bentuk nan disimulasikan. Untuk bintang neutron dengan medan magnet, Z4 lebih aman. Untuk merger lubang hitam cepat, BSSN sudah cukup. Semoga Bermanfaat dan Terima Kasih.