Pendahuluan

Bayangkan Anda ingin memahami mengapa silikon dapat menjadi dasar industri semikonduktor, mengapa besi bersifat magnetik, mengapa suatu oksida dapat menghantarkan listrik hanya setelah diberi dopan, atau mengapa molekul tertentu menempel kuat pada permukaan katalis. Semua pertanyaan itu tampak berbeda, tetapi di tingkat mikroskopik semuanya kembali pada persoalan yang sama: bagaimana elektron dan inti atom tersusun, bergerak, dan saling berinteraksi di dalam material.

Di sinilah density functional theory—dalam bahasa Indonesia akan kita sebut teori fungsional kerapatan, atau singkatnya DFT—menjadi sangat penting. DFT adalah metode mekanika kuantum yang digunakan untuk menghitung struktur elektronik dan sifat fisis sistem atom, molekul, permukaan, dan kristal. Dalam fisika zat padat komputasi, DFT menjadi salah satu alat utama karena mampu menghubungkan struktur atomik dengan energi, ikatan, pita energi, momen magnetik, gaya pada atom, tekanan, dan banyak sifat fisis lain yang relevan untuk penelitian material. Fondasi modern DFT berasal dari teorema Hohenberg–Kohn dan formulasi Kohn–Sham pada tahun 1960-an (Hohenberg dan Kohn, 1964; Kohn dan Sham, 1965).

Buku ini ditulis untuk pembaca dewasa pemula: Anda mungkin sudah pernah mendengar mekanika kuantum, orbital, energi total, atau struktur pita, tetapi belum tentu merasa yakin bagaimana semua itu berubah menjadi sebuah penelitian DFT yang rapi dan siap dipublikasikan. Tujuan buku ini bukan membuat DFT tampak mudah secara palsu. Tujuan yang lebih sehat adalah membuat DFT terbaca, terstruktur, dapat diuji, dan dapat digunakan dengan tanggung jawab ilmiah.

Mengapa DFT perlu dipelajari dengan serius?

DFT sering terlihat praktis: kita menyiapkan struktur, memilih fungsional, memilih pseudopotential, mengatur cutoff energi dan k-point, menjalankan perangkat lunak, lalu memperoleh energi total, struktur pita, density of states, atau muatan atom. Karena banyak paket DFT modern menyediakan alur kerja yang relatif nyaman, seseorang dapat menghasilkan angka dan gambar cukup cepat.

Namun, menghasilkan keluaran bukan berarti menghasilkan pengetahuan.

Misalnya, sebuah perhitungan dapat memberi nilai celah pita sebesar 0,7 eV untuk suatu semikonduktor. Pertanyaan ilmiahnya tidak berhenti di sana. Kita harus bertanya:

• Apakah struktur kristalnya sudah direlaksasi dengan benar?
• Apakah cutoff energi dan k-point sudah konvergen?
• Apakah fungsional tukar-korelasi yang dipakai memang sesuai?
• Apakah nilai eigen Kohn–Sham boleh langsung ditafsirkan sebagai celah eksitasi eksperimen?
• Apakah efek spin–orbit coupling, korelasi kuat, atau interaksi van der Waals penting untuk sistem ini?
• Apakah hasil tersebut konsisten dengan literatur atau data eksperimen?

Pertanyaan-pertanyaan seperti ini membedakan pengguna perangkat lunak dari peneliti komputasi. Buku ini bertujuan membawa Anda ke arah yang kedua.

Dalam fisika zat padat, sifat material sering ditentukan oleh struktur elektroniknya: bagaimana keadaan elektron terisi, bagaimana pita energi terbentuk dalam kristal periodik, dan bagaimana elektron merespons perubahan posisi atom, medan eksternal, atau gangguan lain. Bahasa seperti kisi kristal, ruang resiprok, zona Brillouin, dan pita energi adalah bahasa dasar untuk menjelaskan elektron di padatan periodik (Kittel, 2005). DFT menyediakan cara komputasional untuk mengisi bahasa itu dengan angka, peta, dan prediksi.

Sebagai contoh, jika kita ingin membandingkan dua struktur kristal dari material yang sama, DFT dapat digunakan untuk menghitung energi total masing-masing struktur. Struktur dengan energi lebih rendah biasanya lebih stabil pada kondisi yang dimodelkan, selama asumsi perhitungannya tepat. Jika kita ingin mengetahui apakah sebuah adsorbat menempel kuat pada permukaan logam, DFT dapat membantu menghitung energi adsorpsi. Jika kita ingin melihat apakah material bersifat logam atau isolator, kita dapat menghitung struktur pita dan density of states. Semua contoh ini akan dibahas lebih rinci dalam bab-bab berikutnya.

Dari fungsi gelombang ke kerapatan elektron

Untuk memahami mengapa DFT menarik, kita perlu mulai dari ide dasar mekanika kuantum. Dalam mekanika kuantum, keadaan suatu sistem dijelaskan oleh fungsi gelombang. Fungsi gelombang adalah objek matematis yang memuat informasi tentang kemungkinan hasil pengukuran. Untuk satu elektron, fungsi gelombang bergantung pada posisi elektron. Untuk banyak elektron, fungsi gelombang bergantung pada posisi semua elektron sekaligus.

Di sinilah kesulitannya muncul. Jika suatu sistem memiliki \(N\) elektron, fungsi gelombangnya tidak hanya bergantung pada satu koordinat ruang, tetapi pada koordinat semua elektron. Secara sederhana, untuk \(N\) elektron dalam tiga dimensi, fungsi gelombang bergantung pada \(3N\) koordinat ruang, belum termasuk spin. Artinya, kompleksitas masalah bertambah sangat cepat ketika jumlah elektron bertambah. Untuk kristal nyata dengan banyak atom, penyelesaian eksak masalah banyak-elektron hampir selalu tidak mungkin dilakukan secara langsung dengan fungsi gelombang penuh.

DFT menawarkan perubahan sudut pandang. Alih-alih menjadikan fungsi gelombang banyak-elektron sebagai variabel utama, DFT menjadikan kerapatan elektron sebagai variabel utama. Kerapatan elektron, biasanya ditulis \(\rho(\mathbf{r})\), menyatakan banyaknya elektron per satuan volume di sekitar posisi \(\mathbf{r}\). Jika \(\rho(\mathbf{r})\) besar di suatu daerah, daerah itu memiliki peluang lebih besar untuk ditempati elektron. Jika \(\rho(\mathbf{r})\) kecil, daerah itu relatif jarang ditempati elektron.

Contoh sederhana: pada atom, kerapatan elektron biasanya besar di sekitar inti karena elektron tertarik oleh muatan positif inti. Pada ikatan kovalen, kerapatan elektron dapat terkumpul di antara dua atom, menandakan adanya pembagian elektron yang membantu mengikat atom-atom tersebut. Pada material ionik, kerapatan dapat lebih terkonsentrasi di sekitar ion yang lebih elektronegatif.

Kerapatan elektron jauh lebih sederhana daripada fungsi gelombang banyak-elektron karena hanya bergantung pada tiga koordinat ruang, bukan pada semua koordinat elektron sekaligus. Keajaiban DFT adalah bahwa, untuk keadaan dasar, kerapatan elektron bukan sekadar ringkasan kasar. Teorema Hohenberg–Kohn menunjukkan bahwa kerapatan elektron keadaan dasar menentukan potensial eksternal, hingga konstanta aditif, dan karena itu menentukan sifat keadaan dasar sistem elektron dalam kerangka asumsi teorema tersebut (Hohenberg dan Kohn, 1964).

Pernyataan ini sangat dalam. Ia tidak mengatakan bahwa semua persoalan otomatis mudah. Ia mengatakan bahwa, secara prinsip, kerapatan elektron memuat informasi yang cukup untuk menentukan energi keadaan dasar. Karena itu energi dapat dipandang sebagai fungsional dari kerapatan.

Istilah fungsional berarti aturan yang menerima sebuah fungsi sebagai masukan dan menghasilkan sebuah bilangan sebagai keluaran. Fungsi biasa, misalnya \(f(x)\), menerima bilangan \(x\) dan menghasilkan bilangan. Fungsional, misalnya \(E[\rho]\), menerima fungsi kerapatan \(\rho(\mathbf{r})\) dan menghasilkan energi \(E\). Tanda kurung siku sering dipakai untuk menekankan bahwa masukan fungsional adalah sebuah fungsi, bukan sekadar bilangan.

Sebagai analogi, bayangkan Anda memiliki peta suhu \(T(\mathbf{r})\) di seluruh ruangan. Sebuah aturan yang menghitung suhu rata-rata ruangan dari seluruh peta suhu itu adalah fungsional: masukannya fungsi suhu di ruang, keluarannya satu bilangan. Dalam DFT, kita mencari aturan energi \(E[\rho]\) yang memberi energi total dari kerapatan elektron.

Peran Kohn–Sham: jembatan antara teori dan perhitungan

Teorema Hohenberg–Kohn memberi dasar konseptual, tetapi belum langsung memberi resep komputasi yang praktis. Masalah utamanya adalah kita tidak mengetahui bentuk eksak fungsional energi untuk sistem banyak-elektron nyata. Kohn dan Sham kemudian memperkenalkan pendekatan yang sangat penting: mengganti masalah elektron berinteraksi dengan sistem bantu berupa elektron tak berinteraksi yang menghasilkan kerapatan elektron yang sama dengan sistem asli (Kohn dan Sham, 1965).

Sistem bantu ini bukan berarti elektron nyata tiba-tiba tidak saling berinteraksi. Interaksi elektron tetap ada, tetapi sebagian besar efeknya dimasukkan ke dalam beberapa suku energi, terutama energi tukar-korelasi. Energi tukar-korelasi mencakup bagian rumit dari efek mekanika kuantum banyak-elektron yang tidak tertangkap oleh energi kinetik elektron tak berinteraksi dan interaksi Coulomb klasik secara langsung.

Kata tukar berkaitan dengan sifat elektron sebagai fermion. Elektron memiliki spin setengah bilangan bulat dan fungsi gelombang banyak-elektronnya harus antisimetri ketika dua elektron dipertukarkan. Sifat ini menimbulkan efek tukar. Kata korelasi mengacu pada kenyataan bahwa gerak satu elektron tidak bebas dari gerak elektron lain karena adanya tolakan Coulomb dan efek kuantum. Dua elektron cenderung menghindari satu sama lain dengan cara yang lebih halus daripada yang dapat dijelaskan oleh medan rata-rata sederhana.

Dalam praktik, kita harus memilih pendekatan untuk energi tukar-korelasi. Contohnya adalah local density approximation, generalized gradient approximation, meta-GGA, fungsional hibrida, DFT+U, dan koreksi dispersi van der Waals. Pilihan ini bukan detail kecil. Pilihan fungsional dapat mengubah parameter kisi, energi formasi, celah pita, urutan kestabilan fase, atau prediksi magnetisme. Keterbatasan fungsional aproksimasi, termasuk kesalahan delokalisasi dan kesulitan menangani beberapa sistem berkorelasi kuat, telah dibahas secara luas dalam literatur DFT modern (Cohen, Mori-Sánchez, dan Yang, 2008).

Karena itu, salah satu pesan penting buku ini adalah: DFT bukan satu tombol, melainkan kerangka teori dan praktik numerik yang harus dikendalikan dengan pemahaman.

Apa yang dapat dihitung dengan DFT?

DFT paling sering digunakan untuk menghitung sifat keadaan dasar atau sifat yang berkaitan erat dengan keadaan dasar. Keadaan dasar adalah keadaan energi terendah suatu sistem pada kondisi tertentu. Banyak sifat material dapat diturunkan dari keadaan ini.

Beberapa contoh penting adalah sebagai berikut.

Pertama, DFT dapat menghitung energi total. Energi total memungkinkan kita membandingkan kestabilan relatif beberapa struktur. Misalnya, jika dua fase kristal memiliki komposisi sama, energi total per atom dapat membantu menentukan fase mana yang lebih stabil pada tekanan mendekati nol, selama efek suhu dan entropi tidak dominan.

Kedua, DFT dapat menghitung gaya pada atom. Gaya adalah turunan energi terhadap posisi atom. Dengan gaya ini, kita dapat melakukan relaksasi geometri: atom-atom digerakkan secara numerik sampai gaya menjadi kecil. Contohnya, ketika memodelkan molekul yang teradsorpsi pada permukaan, posisi awal adsorbat mungkin hanya tebakan. Relaksasi DFT membantu menemukan konfigurasi lokal yang lebih stabil.

Ketiga, DFT dapat menghitung tegangan dan tekanan dalam kristal. Kuantitas ini penting untuk optimasi parameter kisi, perhitungan modulus bulk, dan studi transisi fase akibat tekanan.

Keempat, DFT dapat menghasilkan struktur pita dan density of states. Struktur pita menunjukkan hubungan antara energi elektron dan vektor gelombang di ruang resiprok. Density of states menunjukkan banyaknya keadaan elektronik pada rentang energi tertentu. Keduanya sangat penting untuk menafsirkan apakah suatu material cenderung bersifat logam, semikonduktor, atau isolator.

Kelima, DFT dapat digunakan untuk menganalisis kerapatan muatan, ikatan, dan momen magnetik. Misalnya, pada material magnetik, perhitungan spin-polarized dapat menunjukkan perbedaan populasi elektron spin-up dan spin-down serta menghasilkan estimasi momen magnetik lokal.

Keenam, dengan metode lanjutan, DFT dapat menjadi dasar untuk mempelajari fonon, sifat termal, respons optik, permukaan, cacat, antarmuka, dan material berkorelasi. Namun, setiap perluasan membawa asumsi dan parameter tambahan yang harus diuji.

Apa yang tidak boleh dijanjikan DFT secara sembarangan?

DFT kuat, tetapi tidak mahakuasa. Kesalahan umum pemula adalah menganggap semua keluaran perangkat lunak sebagai kebenaran langsung. Sikap ini berbahaya.

Pertama, hasil DFT bergantung pada aproksimasi fungsional tukar-korelasi. Tidak ada satu fungsional yang paling baik untuk semua sistem dan semua sifat. Fungsional semilokal seperti LDA dan GGA sering berguna, tetapi dapat gagal untuk beberapa celah pita, interaksi dispersi, transfer muatan, atau sistem dengan korelasi elektron kuat.

Kedua, hasil DFT bergantung pada parameter numerik. Cutoff energi yang terlalu rendah, sampling k-point yang terlalu kasar, supercell yang terlalu kecil, atau kriteria konvergensi yang longgar dapat menghasilkan kesimpulan salah. Dua perhitungan dengan fungsional sama dapat berbeda hanya karena parameter numeriknya tidak setara.

Ketiga, beberapa besaran yang tampak mudah dibaca dari keluaran sebenarnya memiliki batas interpretasi. Misalnya, nilai eigen Kohn–Sham sering digunakan untuk menggambar struktur pita, tetapi tidak semua nilai eigen itu dapat disamakan begitu saja dengan energi eksitasi quasiparticle eksperimen. Dalam banyak kasus, celah pita dari DFT semilokal lebih kecil daripada celah eksperimen. Ini bukan sekadar “kesalahan perangkat lunak”, melainkan konsekuensi dari aproksimasi dan makna formal besaran yang dihitung.

Keempat, model atomik yang salah tidak dapat diselamatkan oleh perhitungan yang rapi. Jika struktur awal tidak relevan, keadaan spin tidak diuji, muatan cacat tidak dipertimbangkan, atau permukaan dibuat terlalu tipis, hasil akhir dapat tampak presisi tetapi tidak bermakna secara fisis.

Karena itu, buku ini akan berulang kali menekankan tiga kata: konvergensi, validasi, dan interpretasi.

• Konvergensi berarti hasil tidak berubah secara signifikan ketika parameter numerik diperketat.
• Validasi berarti hasil diperiksa terhadap eksperimen, literatur, model teori lain, atau uji konsistensi internal.
• Interpretasi berarti kita memahami apa yang benar-benar dihitung, bukan hanya apa yang ingin kita klaim.

Cara berpikir yang akan dipakai dalam buku ini

Buku ini dibangun dari dasar menuju praktik penelitian. Kita akan memulai dari mekanika kuantum karena DFT bukan sekadar teknik numerik; DFT adalah metode kuantum. Kita akan membahas fungsi gelombang, operator, persamaan Schrödinger, keadaan eigen, nilai harap, dan prinsip variasi. Setelah itu, kita masuk ke padatan periodik: kisi kristal, sel satuan, ruang resiprok, zona Brillouin, dan teorema Bloch.

Kemudian kita membahas masalah banyak-elektron, karena DFT lahir sebagai jawaban terhadap kesulitan tersebut. Dari sana, kita masuk ke teorema Hohenberg–Kohn, persamaan Kohn–Sham, fungsional tukar-korelasi, basis set, pseudopotential, metode PAW, sampling k-point, dan prosedur self-consistent field.

Setelah fondasi itu kuat, buku ini bergerak ke praktik penelitian: alur kerja DFT, uji konvergensi, optimasi struktur, persamaan keadaan, struktur pita, density of states, analisis muatan, magnetisme, cacat, permukaan, fonon, sifat optik, transport, perangkat lunak, dan desain proyek penelitian.

Urutannya sengaja dibuat bertahap. Jika Anda langsung memulai dari “cara menjalankan input file”, Anda mungkin bisa menghasilkan keluaran, tetapi sulit menilai apakah keluaran itu benar. Jika Anda hanya mempelajari teori tanpa praktik, Anda mungkin paham persamaan tetapi tidak siap melakukan penelitian. Buku ini mencoba menjaga keduanya tetap terhubung.

Contoh peta riset kecil

Untuk memberi gambaran, bayangkan sebuah proyek sederhana: “Apakah substitusi atom B pada material A dapat menurunkan energi formasi cacat dan mengubah sifat elektroniknya?”

Dengan DFT, proyek seperti ini dapat dipecah menjadi langkah-langkah ilmiah:

1. Menentukan struktur kristal material A.
2. Menguji konvergensi cutoff energi dan k-point.
3. Merelaksasi struktur murni.
4. Membuat supercell dan mengganti satu atom dengan atom B.
5. Merelaksasi struktur cacat.
6. Menghitung energi formasi dengan rumus yang sesuai.
7. Membandingkan density of states sebelum dan sesudah substitusi.
8. Memeriksa apakah muncul keadaan elektronik baru di dekat tingkat Fermi.
9. Membandingkan hasil dengan literatur atau eksperimen jika tersedia.
10. Menyatakan batasan model, misalnya ukuran supercell, pilihan fungsional, dan kemungkinan keadaan muatan cacat.

Contoh ini menunjukkan bahwa penelitian DFT bukan hanya “menghitung satu angka”. Penelitian DFT adalah rangkaian keputusan. Setiap keputusan harus dapat dijelaskan.

Sikap ilmiah yang diperlukan

DFT mengajarkan kerendahan hati. Hasil komputasi sering terlihat sangat presisi: banyak digit, gambar berwarna, kurva halus, dan file keluaran panjang. Tetapi presisi numerik tidak sama dengan kebenaran fisis. Angka 3,742819 eV tidak otomatis lebih benar daripada 3,7 eV jika model, fungsional, atau parameter numeriknya belum diuji.

Sikap ilmiah yang baik dalam DFT mencakup beberapa kebiasaan:

• tidak mempercayai hasil pertama tanpa uji konvergensi;
• mencatat semua parameter penting;
• membedakan hasil numerik, asumsi model, dan interpretasi fisis;
• membaca makalah pembanding secara kritis;
• melaporkan metode dengan cukup rinci agar dapat direproduksi;
• berani mengatakan bahwa suatu hasil masih terbatas.

Dalam penelitian modern, reproduktibilitas menjadi semakin penting. Basis data material dan alur kerja komputasi skala besar menunjukkan bahwa DFT dapat digunakan secara sistematis untuk menjelajahi banyak material, tetapi keberhasilan pendekatan semacam itu tetap bergantung pada standar metodologi, dokumentasi parameter, dan validasi yang baik (Jain et al., 2013).

Apa yang diharapkan setelah membaca buku ini?

Setelah menyelesaikan buku ini, Anda diharapkan mampu membaca makalah DFT dengan lebih kritis. Ketika melihat kalimat seperti “perhitungan dilakukan dengan GGA-PBE, cutoff 500 eV, k-point 6×6×6, dan struktur direlaksasi hingga gaya lebih kecil dari 0,01 eV/Å”, Anda tidak hanya mengenali istilahnya, tetapi juga memahami mengapa parameter itu penting dan apa yang perlu ditanyakan.

Anda juga diharapkan mampu merancang perhitungan DFT dasar untuk material kristalin: memilih struktur, menjalankan SCF, melakukan relaksasi, menguji konvergensi, menghitung energi relatif, membuat struktur pita dan density of states, serta menulis bagian metode komputasi secara jujur dan lengkap.

Lebih jauh lagi, Anda diharapkan memahami batas DFT. Peneliti yang baik bukan hanya tahu kapan DFT dapat digunakan, tetapi juga tahu kapan DFT perlu diperluas, dibandingkan dengan metode lain, atau ditafsirkan dengan hati-hati.

DFT adalah jembatan antara persamaan kuantum dan sifat material nyata. Jembatan itu kuat, tetapi harus dilalui dengan peta yang benar. Buku ini adalah peta awal Anda.

References

Cohen, A. J., Mori-Sánchez, P., dan Yang, W. (2008). Insights into current limitations of density functional theory. Science, 321(5890), 792–794.

Hohenberg, P., dan Kohn, W. (1964). Inhomogeneous electron gas. Physical Review, 136(3B), B864–B871.

Jain, A., Ong, S. P., Hautier, G., Chen, W., Richards, W. D., Dacek, S., Cholia, S., Gunter, D., Skinner, D., Ceder, G., dan Persson, K. A. (2013). Commentary: The Materials Project: A materials genome approach to accelerating materials innovation. APL Materials, 1, 011002.

Kittel, C. (2005). Introduction to Solid State Physics (8th ed.). John Wiley & Sons.

Kohn, W., dan Sham, L. J. (1965). Self-consistent equations including exchange and correlation effects. Physical Review, 140(4A), A1133–A1138.