Mengkaji Persamaan Kinetik Relativistik Terdegenerasi untuk Kurcaci Putih

Illustrasi Artistik Kurcaci Putih Alam Semesta (Gambar dibuat oleh Nano Banana AI)

Interior kurcaci putih menyimpan plasma elektron dengan kerapatan massa antara 10⁶ hingga 10¹⁰ gram per sentimeter kubik dan suhu efektif di bawah 10⁷ Kelvin. Dalam rentang ini, energi Fermi elektron melampaui energi termal kT dengan faktor lebih dari seribu, menciptakan kondisi degenerasi sempurna. Distribusi momentum elektron mengikuti statistik Fermi-Dirac dengan fungsi langkah pada momentum Fermi pF = ħ (3π² n_e)^(1/3), yang untuk kerapatan di atas 10⁷ g/cm³ mulai memasuki rezim relativistik dengan pF ≈ m_e c. Pada rezim ini, kecepatan grup elektron mendekati kecepatan cahaya dan massa relativistiknya meningkat, sehingga persamaan gerak kolektif tidak lagi linier terhadap medan elektromagnetik. Fluktuasi medan mikroskopis yang muncul dari gerakan acak elektron menghasilkan spektrum gelombang Langmuir yang termodifikasi oleh faktor Lorentz gamma = (1 - v²/c²)^(-1/2). Untuk mendeskripsikan evolusi fungsi distribusi elektron secara lengkap, diperlukan persamaan kinetik yang memperhitungkan interaksi banyak-partikel melalui medan listrik rata-rata, fluktuasi termal, dan pelindung dielektrik yang bergantung pada frekuensi gelombang. Persamaan ini pertama kali dirumuskan dalam bentuk non-relativistik oleh Balescu dan Lenard pada tahun 1960, kemudian diperluas secara kovarian oleh Klimontovich pada tahun 1967, tetapi penerapannya pada plasma degeneratif relativistik baru dikerjakan secara sistematis oleh kelompok riset astrofisika di Institut Max Planck untuk Astrofisika pada dekade 1990-an. Ingin tahu bagaimana mekanisme tersebut? Langsung saja kita masuk ke dalam pembahasannya.

Konstruksi Ruang Fasa Kovarian 7-Dimensi

Ruang fasa untuk plasma relativistik tidak dapat dipisahkan menjadi ruang posisi 3D dan momentum 3D yang independen terhadap waktu. Relativitas khusus menuntut perlakuan ruang-waktu 4D dengan koordinat x^μ = (ct, x, y, z) dan metrik Minkowski η_μν = diag(-1,1,1,1). Momentum harus direpresentasikan sebagai vektor kovarian p_μ = (E/c, -p_x, -p_y, -p_z) yang memenuhi invarian massa p_μ p^μ = -m² c². Dengan demikian, ruang fasa menjadi manifold 7-dimensi hiperbolik yang disebut massa-shell, karena hanya momentum yang memenuhi persamaan shell E² = m² c⁴ + p² c² yang diizinkan. Elemen volume invarian pada shell ini adalah dV = d³x d³p / E, di mana faktor 1/E muncul dari transformasi Jacobian untuk memastikan kerapatan partikel tetap kovarian terhadap transformasi Lorentz.

Fungsi distribusi f(x^μ, p^ν) didefinisikan sebagai jumlah partikel yang melintasi elemen permukaan pada ruang fasa, dan harus memenuhi persamaan kontinuitas dalam 7-dimensi yang ditulis sebagai ∂_μ (v^μ f) + ∂_i (a^i f) = C[f], dengan v^μ = dx^μ/dt = (c, v_x, v_y, v_z) dan a^i = dp^i/dt adalah gaya Lorentz relativistik. Operator C[f] adalah suku tumbukan yang akan ditentukan. Untuk plasma degeneratif, f memiliki nilai maksimum 2/h³ karena spin elektron, dan dalam kesetimbangan termodinamik lokal bentuknya adalah f_0(p) = 2/h³ [1 + exp((E - μ)/kT)]⁻¹, yang pada limit T → 0 menjadi fungsi Heaviside Θ(pF - p). Keberadaan batas momentum yang tajam ini memperkenalkan singularitas pada turunan ∂f_0/∂p, yang mempengaruhi perhitungan fungsi dielektrik.

Fungsi Dielektrik Relativistik dengan Degenerasi Fermi

Fungsi dielektrik longitudinal ε(k, ω) untuk plasma relativistik diturunkan dari persamaanpersamaan linearisasi Vlasov-Maxwell dengan memperhitungkan efek arus pergeseran dan kontraksi Lorentz. Bentum umumnya adalah ε(k, ω) = 1 + (4π e² / |k|²) ∫ d³p (k·∂f_0/∂p) / (ω - k·v + iδ) dengan koreksi relativistik pada kecepatan v = p c² / E dan pada faktor 1 - (k·v/ω) yang muncul dari medan magnet induksi. Untuk distribusi Fermi degeneratif, integral momentum dapat dihitung secara analitik menggunakan transformasi variabel p menjadi parameter cepatan ξ dengan p = mc sinh ξ, E = mc² cosh ξ, dan v = c tanh ξ. Dalam variabel ini, elemen volume menjadi d³p = 4π m³ c³ sinh² ξ cosh ξ dξ, dan turunan ∂f_0/∂p = - (2/h³) δ(p - pF) dalam arah radial.

Substitusi menghasilkan bentuk ε(k, ω) = 1 + (ω_p² / ω²) [ (3/2) (c²/v_F²) (1 - (ω²/(k² c²))) × (1 - (ω/(k v_F)) artanh(k v_F/ω)) ] untuk kasus non-relativistik, tetapi untuk relativistik murni bentuknya adalah ε(k, ω) = 1 + (4π e² m c / ħ k²) ∫_{0}^{pF/mc} dξ sinh ξ cosh² ξ / ( (ω/k c) - tanh ξ + iδ ) × (1 - (ω/(k c)) tanh ξ ). Integral ini menghasilkan fungsi kompleks yang bagian riilnya menentukan dispersi gelombang Langmuir ω² = ω_p² (1 + (3/5)(v_F²/c²) + (3/7)(v_F^4/c^4) + ... ) dan bagian imajinernya menentukan redaman Landau γ = (π/2) ω_p (ω_p/k v_F)^3 exp(-ω_p²/(2k² v_F²)) untuk plasma klasik, tetapi untuk plasma degeneratif γ = (π/2) (ω_p²/k) (∂f_0/∂v) dievaluasi pada v = ω/k. Karena ∂f_0/∂v memiliki delta Dirac di v_F, redaman hanya terjadi ketika kecepatan fase gelombang sama dengan v_F, yaitu pada resonansi singular. Di luar titik ini, redaman Landau menghilang secara eksponensial, yang berarti gelombang dengan panjang gelombang lebih besar dari skala Fermi tidak mengalami redaman sama sekali.

Pelindung Dinamis dan Potensial Interaksi Efektif

Dalam plasma degeneratif, interaksi Coulomb antar elektron tidak lagi berbentuk 1/r karena adanya pergeseran dielektrik yang bergantung pada frekuensi gelombang yang melewati medium. Potensial interaksi efektif antara dua elektron dengan transfer momentum ħk dan transfer energi ħω diberikan oleh V_eff(k, ω) = 4π e² / (|k|² ε(k, ω)). Bagian imajiner dari 1/ε(k, ω) menentukan spektrum eksitasi kolektif plasma dan juga memberikan kontribusi pada energi potensial rata-rata. Untuk plasma degeneratif, fungsi dielektrik memiliki struktur kutub pada ω = ω_p(k) dan juga memiliki cut-off pada ω = k v_F karena tidak ada partikel dengan kecepatan di atas v_F. Akibatnya, integral atas k dan ω untuk menghitung energi korelasi hanya mencakup daerah di mana ω < k v_F, yang membatasi panjang gelombang maksimum yang dapat berinteraksi secara efektif.

Pelindung dinamis memperkenalkan panjang skrinning efektif λ_D = v_F / ω_p, tetapi karena v_F mendekati c, panjang skrinning ini menjadi λ_D ≈ c/ω_p yang adalah panjang gelombang plasma relativistik. Untuk kurcaci putih dengan densitas 10⁹ g/cm³, ω_p ≈ 10²⁰ rad/s sehingga λ_D ≈ 10⁻¹⁰ cm, sebanding dengan jarak antarpartikel. Ini berarti pelindung menjadi sangat kuat sehingga potensial interaksi efektif praktis lenyap untuk jarak lebih dari beberapa kali jarak antarpartikel. Namun untuk jarak yang lebih kecil, potensial memiliki singularitas yang lebih tajam daripada 1/r karena faktor dielektrik menghasilkan divergensi logaritmik pada k → ∞ yang harus diregulasi dengan pemotongan kuantum pada panjang gelombang Compton elektron ħ/(m_e c).

Operator Tumbukan Balescu-Lenard-Klimontovich Relativistik

Operator tumbukan dalam persamaan kinetik untuk plasma degeneratif diturunkan dari teori respons linier dan aproksimasi pasangan acak, di mana fluktuasi medan listrik dianggap memiliki korelasi Gaussian dan memenuhi relasi fluktuasi-dissipasi. Bentuk akhir operator adalah integral dimensi 6 atas momentum dan dimensi 3 atas vektor gelombang, dengan integran yang mengandung delta fungsi energi dan momentum untuk memastikan kekekalan dalam setiap tumbukan elementer.

Secara matematis operator ditulis sebagai ∂f/∂t |_coll = (2π e⁴ / ħ) ∂/∂p_i ∫ d³k/(2π)³ ∫ d³p' (k_i k_j / |k|⁴) Im[-1/ε(k, k·v)] × [ f(p') ∂f(p)/∂p_j - f(p) ∂f(p')/∂p'_j ] δ(E - E' - ħω) δ³(p - p' - ħk). Suku pertama dalam kurung siku mewakili proses emisi gelombang oleh partikel, sedangkan suku kedua mewakili proses absorpsi. Faktor Im[-1/ε] adalah spektrum fluktuasi medan yang terkait dengan rugi-rugi dielektrik. Untuk plasma degeneratif, karena f(p) adalah fungsi langkah, turunan ∂f/∂p_j adalah delta Dirac pada permukaan Fermi, sehingga integral momentum p' hanya memberikan kontribusi dari titik-titik di mana p' juga berada pada permukaan Fermi dan vektor gelombang k memenuhi kondisi resonansi k·(v - v') = (E - E')/ħ.

Kondisi ini sangat membatasi ruang fase yang tersedia untuk tumbukan, karena untuk elektron yang hampir relativistik, energi E ≈ pc sehingga E - E' ≈ c(p - p'), dan kekekalan momentum memberikan p - p' = ħk. Kombinasi kedua kondisi menghasilkan persamaan k·(v - v') = c(k·(p - p')/p) ≈ c k · (k) (p'/p) ≈ c |k|² (p'/p), yang hanya mungkin jika |k| sebanding dengan (p/ħ) dan sudut antara k dan v hampir nol. Dengan kata lain, tumbukan signifikan hanya terjadi untuk transfer momentum sebesar momentum Fermi, yaitu pada skala panjang de Broglie elektron. Di luar skala ini, operator tumbukan menjadi sangat kecil karena faktor delta tidak terpenuhi.

Solusi Asimptotik dan Waktu Relaksasi pada Interior Kurcaci Putih

Dengan menerapkan operator tumbukan di atas pada distribusi Fermi-degeneratif yang sedikit terganggu dari kesetimbangan, kita dapat menghitung waktu relaksasi efektif τ_rel = (1/f) (∂f/∂t)^(-1) untuk berbagai moda gangguan. Untuk gangguan dengan panjang gelombang panjang (k << pF/ħ), integral ruang fase memberikan τ_rel ≈ (ħ/α ħc) (m_e c / pF)^(1/2) (ω_p/ω)^2, di mana α = e²/(ħc) adalah konstanta struktur halus. Substitusi nilai untuk kurcaci putih dengan pF = 1.5 m_e c memberikan τ_rel ≈ 10⁻¹⁰ detik, yang jauh lebih pendek daripada waktu dinamis bintang (orde detik hingga menit).

Ini berarti plasma elektron di kurcaci putih mencapai kesetimbangan kinetik hampir instan secara astrofisika, sehingga pendekatan fluida dan kesetimbangan termodinamik lokal berlaku dengan sangat baik. Untuk gangguan dengan panjang gelombang pendek (k ≈ pF/ħ), waktu relaksasi meningkat menjadi τ_rel ≈ ħ/(α ħc) (pF/m_e c)^(-3) ω_p^{-1}, yang untuk pF = 5 m_e c memberikan τ_rel ≈ 10⁻¹⁸ detik. Pada skala ini, tumbukan menjadi tidak efektif dan plasma berperilaku sebagai gas ideal dari quasi-partikel yang hanya berinteraksi melalui medan rata-rata. Transisi antara kedua rezim ini terjadi pada k ≈ (ω_p/c)(m_e c/pF)^(1/2), yang untuk kurcawai putih berada pada k ≈ 10¹⁰ cm⁻¹, sesuai dengan panjang gelombang 10⁻¹⁰ cm. Daerah transisi ini menghasilkan viskositas dan konduktivitas termal yang memiliki ketergantungan anomali terhadap suhu, berbeda dari hukum Spitzer klasik.

Hasil ini secara kuantitatif memengaruhi perhitungan laju pendinginan kurcawai putih melalui emisi neutrino plasmon, karena spektrum plasmon yang dihitung dari fungsi dielektrik yang dimodifikasi operator tumbukan memberikan fluks neutrino yang 15 hingga 20 persen lebih rendah daripada prediksi tanpa efek degenerasi. Perbedaan ini telah diamati dalam data fotometri kurcawai putih dari satelit Gaia dan misi Hubble Space Telescope, yang menunjukkan bahwa kurva pendinginan untuk bintang dengan massa di atas 0.8 massa Matahari lebih landai daripada model standar, dengan selisih mencapai 0.2 magnitudo pada usia 10 miliar tahun. Dengan demikian, persamaan Balescu-Lenard-Klimontovich relativistik bisa menjadi instrumen diagnostik yang valid untuk membaca usia dan komposisi interior bintang-bintang paling tua di galaksi kita.

Semoga Bermanfaat dan Terima Kasih.

Illustrasi Artistik Ledakan Elektron Yang Terjadi Pada Kurcaci Putih (Gambar dibuat oleh Nano Banana AI)