Mengkaji Batas Perubahan Aliran Energi Cahaya pada Cakram Benda Hitam Supermasif

Di sekitar lubang hitam supermasif, ada wilayah nan suhunya mencapai miliaran derajat. Materi berputar membentuk piringan, medan magnet terpilin, dan foton terperangkap alias melesat bebas. Nama wilayah itu adalah piringan akresi (accretion disk). Piringan akresi terbentuk ketika materi dari bintang pendamping alias awan gas jatuh menuju lubang hitam. Karena materi mempunyai momentum sudut, dia tidak bisa jatuh langsung secara radial, melainkan membentuk struktur pipih nan berputar mengelilingi lubang hitam. Piringan ini terbagi menjadi beberapa zona. Bagian terluar disebut piringan luar (outer disk), temperaturnya relatif rendah sekitar ribuan hingga jutaan kelvin, didominasi atom netral dan molekul. Semakin ke dalam, gravitasi lubang hitam memampatkan dan memanaskan materi. Di bagian tengah, suhu mencapai jutaan kelvin, materi terionisasi penuh menjadi plasma, dan medan magnet terpilin sempurna lantaran pengaruh dinamo. Di bagian paling dalam, tepat sebelum horizon peristiwa, suhu dapat melonjak hingga miliaran kelvin. Di sinilah foton dengan daya sinar-X dan sinar gamma diproduksi secara masif. Wilayah terdalam ini juga tempat terjadinya transisi optik di bagian dalam sekali, plasma sangat panas dan rapat sehingga optik tebal, foton terperangkap dan berdifusi keluar seperti dalam oven. Sedikit lebih ke luar, kerapatan turun drastis, foton tiba-tiba bisa melesat bebas tanpa tumbukan, menciptakan apa nan disebut korona panas (hot corona). Batas antara kedua area ini nan menjadi konsentrasi utama dalam Radiation-GRMHD, lantaran di sinilah fluks radiasi berubah secara ekstrem dalam jarak nan sangat pendek.

GRMHD, ialah campuran persamaan kontinuitas, persamaan Euler relativistik, persamaan induksi magnet, dan persamaan Einstein untuk metrik latar. Namun plasma ini juga memancarkan dan menyerap foton. Maka perlu ditambahkan persamaan transpor radiasi. Di sinilah masalah dimulai. Persamaan transpor radiasi corak paling umum adalah persamaan transfer radiative relativistik dalam ruang-waktu melengkung. Bentuknya adalah turunan kovarian dari intensitas spesifik terhadap frekuensi. Jika diintegralkan terhadap momen, muncul sistem tertutup nan disebut momen Eddington. Sistem ini terdiri dari daya radiasi, fluks radiasi, dan tekanan radiasi. Kopling dengan GRMHD terjadi melalui suku sumber pada sisi kanan persamaan energi-momentum total. Secara matematis, tensor energi-momentum total adalah jumlah tensor untuk materi plus tensor untuk radiasi. Divergensi dari tensor total nol, tetapi divergensi masing-masing tidak nol, saling berlawanan tanda. Inilah inti dari masalah fluks-konservatif. Tidak boleh ada daya nan lenyap dari campuran sistem, meskipun terjadi transfer antara materi dan radiasi. Memahami Batas Optik Tebal dan Tipis dalam Bilangan Tak Dimensi Parameter esensial dalam masalah ini adalah kedalaman optik, dinotasikan dengan huruf Yunani tau. Definisi matematisnya adalah integral dari koefisien opasitas sepanjang garis pandang. Opasitas sendiri adalah jumlah dari koefisien absorpsi dan hamburan. Dalam sistem nan kita bahas, opasitas tergantung pada densitas plasma, suhu, dan gelombang foton. Untuk menyederhanakan, dalam pendekatan momen digunakan opasitas rata-rata frekuensi, sering disebut opasitas Rosseland untuk wilayah tebal dan opasitas Planck untuk wilayah tipis. Batas optik tebal terjadi ketika tau jauh lebih besar dari satu. Pada kondisi ini, foton mengalami banyak hubungan sebelum keluar. Jarak tempuh bebas rata-rata foton sangat mini dibandingkan skala panjang sistem. Akibatnya, pengedaran perspektif intensitas radiasi nyaris isotropik. Fluks radiasi sebanding dengan gradien daya radiasi, dengan konstanta proporsionalitas satu per tiga opasitas. Ini adalah norma Fick relativistik. Batas optik tipis terjadi ketika tau jauh lebih mini dari satu. Foton nyaris tidak pernah bertumbukan. Distribusi perspektif sangat anisotropik, terpusat pada arah tertentu. Fluks radiasi nyaris sama dengan daya radiasi dikalikan kecepatan cahaya. Antara kedua pemisah ini tidak ada transisi tajam secara fisis, tetapi secara numerik terjadi singularitas lantaran turunan fluks terhadap posisi berubah nilai dari sangat mini menjadi sangat besar dalam rentang spasial nan sangat sempit. Lebar lapisan transisi ini sebanding dengan jalur bebas rata-rata foton lokal. Untuk piringan akresi lubang hitam, jalur bebas ini bisa bervariasi dari beberapa meter hingga ribuan kilometer dalam jarak nan sama. Inilah sebabnya simulasi konvensional gagal. Formulasi Asimtotik Batas Singular Fluks-Konservatif Untuk mengatasi masalah di atas, digunakan teknik ekspansi asimtotik. Ide dasarnya adalah memisahkan domain menjadi tiga zona. Zona dalam dengan tau besar, area luar dengan tau kecil, dan lapisan pemisah di antaranya. Di area dalam, daya radiasi diekspansi dalam deret pangkat parameter epsilon, ialah kebalikan dari tau. Suku utama memenuhi persamaan difusi. Suku berikutnya memberikan koreksi nan melibatkan turunan orde lebih tinggi. Di area luar, ekspansi menggunakan epsilon itu sendiri sebagai parameter kecil, tetapi persamaan nan muncul adalah persamaan adveksi murni untuk fluks radiasi. Di lapisan batas, variabel posisi direskalakan dengan epsilon, sehingga persamaan nan muncul adalah persamaan difusi-reaksi dalam skala mikro. Solusi di ketiga area kemudian disambung menggunakan kondisi sambung asimtotik. Kondisi sambung nan pertama adalah kontinuitas daya radiasi. nan kedua adalah kontinuitas fluks radiasi. nan ketiga adalah konservasi total energi, artinya integral dari divergensi fluks radiasi di seluruh domain sama dengan nol setelah dikurangi dengan suku kopling ke materi. Ketiga kondisi ini menghasilkan satu nilai unik untuk fluks pada antarmuka. Nilai ini disebut fluks-konservatif singular. Secara matematis, fluks ini merupakan solusi dari persamaan integral Fredholm jenis kedua dengan kernel nan mengandung opasitas lokal. Tanpa menggunakan ekspansi asimtotik, nilai ini tidak dapat diperoleh dari skema beda hingga biasa lantaran skala panjang lapisan pemisah jauh di bawah resolusi grid. Metode ekspansi asimtotik memberikan rumus tertutup untuk fluks antarmuka nan hanya berjuntai pada nilai daya radiasi di kedua sisi dan opasitas rata-rata. Rumus ini menjadi dasar bagi semua skema numerik modern untuk Radiation-GRMHD. Metode Elemen Hingga Berbasis Geometri Simplektik untuk Sistem Hiperbolik Setelah formulasi asimtotik selesai, langkah berikutnya adalah diskretisasi numerik. Persamaan GRMHD sendiri adalah sistem hiperbolik nonlinear. Artinya, dia mempunyai gelombang karakter dengan kecepatan terbatas, tidak lebih dari kecepatan cahaya. Metode numerik untuk sistem hiperbolik kudu mempertahankan sifat kekekalan. Metode volume hingga biasa melakukannya dengan baik tetapi kandas pada pemisah singular lantaran memerlukan resolusi ekstrem. Solusinya adalah menggunakan metode komponen hingga dengan kegunaan pedoman polinomial orde tinggi, tetapi dengan modifikasi penting, bahwa pengetahuan ukur simplektik kudu dipertahankan. Geometri simplektik adalah struktur geometris pada ruang fase nan menjamin kekekalan volume ruang fase terhadap aliran Hamiltonian. Persamaan GRMHD tanpa radiasi dapat ditulis dalam corak Hamiltonian non-kanonik dengan bracket Poisson nan memuat divergensi medan magnet sebagai kendala. Dengan radiasi, sistem menjadi non-Hamiltonian lantaran ada disipasi, tetapi bagian konvektifnya tetap Hamiltonian. Metode simplektik mempertahankan bagian Hamiltonian secara eksak, sementara bagian disipatif ditambahkan sebagai suku terpisah nan diintegrasi dengan skema implisit. Dalam praktik numerik, ini dilakukan dengan memilih kegunaan pedoman komponen hingga nan ortogonal terhadap turunan simplektik dari Hamiltonian. Ortogonalitas ini menghasilkan matriks massa nan diagonal blok dengan struktur simplektik. Proses integrasi waktu menggunakan skema Runge-Kutta simplektik jenis Gauss-Legendre. Skema ini orde dua, empat, alias enam, tetapi selalu reversible dalam waktu untuk bagian Hamiltonian. Untuk bagian radiasi, digunakan skema implisit lantaran opasitas bisa sangat besar. Gabungan keduanya disebut IMEX-symplectic, di mana IMEX kependekan dari implicit-explicit. Hasil numerik nan sudah dilakukan oleh beberapa fisikawan menunjukkan bahwa metode ini mempertahankan konservasi daya total hingga ketelitian mesin untuk simulasi hingga satu juta langkah waktu, apalagi pada kasus dengan lapisan pemisah optik setajam apapun. Validasi dengan Masalah Uji Standar dan Aplikasi ke Piringan Akresi Metode di atas telah diuji pada tiga masalah standar dalam astrofisika komputasi. Masalah pertama adalah gelombang akustik linear dalam medium homogen dengan opasitas konstan. Solusi analitik untuk redaman gelombang akibat kopling radiasi tersedia dalam corak kegunaan Bessel. Metode komponen hingga simplektik memberikan galat relatif 10 pangkat minus lima untuk amplitudo dan 10 pangkat minus enam untuk frekuensi, pada grid dengan hanya sepuluh komponen per panjang gelombang. Metode volume hingga konvensional memerlukan seratus komponen untuk galat nan sama. Masalah kedua adalah kejutan radiatif, ialah diskontinuitas densitas dan suhu nan bergerak disertai emisi radiasi. Di depan kejutan, medium optik tipis. Di belakang kejutan, medium optik tebal. Lapisan transisi mempunyai ketebalan hanya beberapa kali jalur bebas. Metode simplektik sukses menangkap profil fluks radiasi kontinu tanpa osilasi, sementara metode konvensional menghasilkan loncatan semu pada fluks hingga 30 persen. Masalah ketiga adalah piringan akresi relativistik di sekitar lubang hitam Kerr dengan parameter spin 0,9. Simulasi dilakukan hingga seratus ribu waktu orbit. Hasilnya, fluks radiasi total nan mencapai tak terhingga sesuai dengan prediksi teori thin disk dari Shakura dan Sunyaev dengan aspek koreksi relativistik dari Novikov dan Thorne, tetapi dengan tambahan komponen non-termal dari korona nan sebelumnya tidak dapat dimodelkan secara bersamaan. Dengan demikian, metode asimtotik pemisah singular fluks-konservatif nan dipadukan dengan komponen hingga simplektik merupakan satu-satunya metode nan sampai saat ini bisa mensimulasikan Radiation-GRMHD dari horizon lubang hitam hingga skala jet dalam satu domain, tanpa dugaan perantara tentang corak spektrum radiasi. Semoga Bermanfaat.