Teori fungsi rapatan
Teori fungsi kerapatan ("DFT, Density functional theory") merupakan salah satu dari beberapa pendekatan populer untuk perhitungan struktur elektron banyak-partikel secara mekanika kuantum untuk sistem molekul dan bahan rapat. Teori Fungsi Kerapatan (DFT) adalah teori mekanika kuantum yang digunakan dalam fisika dan kimia untuk mengamati keadaan dasar dari sistem banyak partikel.
Metode tradisional dalam perhitungan struktur elektron, seperti teori Hartree-Fock didasarkan pada fungsi gelombang banyak-elektron yang rumit. Sasaran utama dari teori fungsi kerapatan adalah menggantikan fungsi gelombang elektron banyak-partikel dengan kerapatan elektron sebagai besaran dasarnya. Fungsi gelombang partikel-banyak bergantung pada 3N variabel, yaitu tiga variabel ruang untuk masing-masing N elektron, sedangkan kerapatan hanya merupakan fungsi dari 3 variabel, jadi merupakan suatu besaran yang sederhana untuk ditangani, baik secara konsep maupun secara praktis.
Walaupun teori fungsi kerapatan memiliki dasar konseptualnya dalam model Thomas-Fermi, DFT tidak berlandaskan pijakan teoretis yang kuat sampai munculnya teorema Hohenberg-Kohn (HK) yang menunjukkan adanya pemetaan satu-satu antara kerapatan elektron keadaan dasar dengan fungsi gelombang keadaan dasar dari sistem banyak-partikel. Selain itu, teorema HK membuktikan bahwa kerapatan keadaan dasar meminimalkan energi elektron total sistem tersebut. Karena teorema HK berlaku hanya untuk keadaan dasar, DFT juga merupakan sebuah teorema keadaan dasar.
Teorema Hohenberg-Kohn hanya suatu teorema keberadaan, yang menyatakan bahwa penggambaran itu ada, tetapi tidak menghasilkan penggambaran apapun yang tepat seperti itu. Teorema tersebut dalam penggambaran ini dibuat pendekatan. Penggambaran yang paling terkenal adalah pendekatan kerapatan lokal (LDA) yang memberikan pendekatan penggambaran dari kerapatan sistem terhadap energi total. LDA digunakan untuk gas elektron yang seragam, dikenal juga sebagai jellium.
Pada kenyataannya, teorema HK jarang digunakan secara langsung untuk membuat perhitungan. Sebagai gantinya, implementasi teori fungsi kerapatan yang paling umum digunakan saat ini adalah metode Kohn-Sham. Dalam kerangka DFT Kohn-Sham, masalah interaksi elektron banyak partikel, potensial statis eksternal direduksi menjadi sebuah masalah yang mudah dikerjakan dengan penggantian elektron yang tidak berinteraksi menjadi sebuah potensial efektif. Potensial efektif meliputi potensial eksternal dan pengaruh interaksi “Colomb” antar elektron.
Dalam banyak kasus, DFT dengan pendekatan kerapatan lokal memberikan hasil yang memuaskan jika dibandingkan dengan data eksperimen pada daya komputasi yang relatif rendah, ketika dibandingkan dengan cara-cara penyelesaian masalah mekanika kuantum banyak-partikel yang lain.
DFT menjadi sangat terkenal untuk perhitungan dalam fisika keadaan padat sejak tahun 1970. Akan tetapi, DFT tersebut tidak dapat dipertimbangkan cukup akurat untuk perhitungan kimia kuantum sampai tahun 1990, ketika pendekatan digunakan dalam teori dihasilkan perbaikan yang lebih baik. DFT kini merupakan suatu metode yang mengarahkan pada perhitungan struktur elektron dalam berbagai bidang.
Akan tetapi, masih ada sistem yang tidak dapat dijelaskan dengan baik dengan LDA. LDA tidak dapat menjelaskan dengan baik interaksi antar molekul, terutama gaya van der Waals (dispersi). Hasil lain yang terkenal adalah perhitungan celah pita dalam semikonduktor, tetapi larangan ini tidak dapat memperlihatkan kegagalan, karena DFT adalah teori keadaan dasar dan celah pita adalah sifat keadaan tereksitasi.
Teori fungsi kerapatan pertama kali dikembangkan oleh Thomas dan Fermi pada tahun 1920. Mereka menghitung energi sebuah atom dengan energi kinetiknya diwakili sebagai fungsi dari kerapatan elektron, ini dikombinasikan dengan ungkapan klasik untuk interaksi inti-elektron dan elektron-elektron (yang dapat juga diwakili dalam hubungan kerapatan elektron).
Walaupun ini adalah tahap awal yang penting, ketelitian persamaan Thomas-Fermi terbatas karena persamaan tersebut tidak dapat memperlihatkan pertukaran energi dari sebuah atom yang diramalkan dengan teori Hartree-Fock. Fungsi pertukaran energi ditambahkan oleh Dirac pada tahun 1928.
Akan tetapi, teori Thomas-Fermi-Dirac tetap tidak akurat untuk beberapa aplikasi, karena teori tersebut sulit untuk memperlihatkan energi kinetik dengan sebuah fungsi kerapatan, dan teori tersebut mengabaikan hubungan antara elektron keseluruhan.
Pada umumunya dalam perhitungan struktur elektron banyak partikel, inti yang dimiliki molekul atau 'cluster' terlihat tetap (pendekatan Born-Oppenheimer), menghasilkan sebuah potensial eksternal statis formula_1 di mana elektron berpindah. Keadaan elektron yang mantap 'stationer' dijelaskan dalam fungsi gelombang formula_2 penyelesaian Persamaan Schrodinger banyak elektron.
formula_3
dengan formula_4 adalah jumlah elektron dan formula_5
adalah interaksi elektron-elektron. Operator formula_6
dan formula_5 biasa disebut sebagai operator umum 'universal' karena operator tersebut berlaku untuk semua sistem, dengan formula_1 adalah sistem yang tidak dapat berdiri sendiri atau tidak umum/'universal'. Dapat terlihat perbedaan yang nyata antara permasalahan partikel tunggal dengan permasalahan banyak partikel yang rumit, yaitu pada persamaan interaksi formula_5. Sekarang, banyak metode yang pintar untuk menyelesaikan persamaan Schrodinger, misalnya teori gangguan diagramatis dalam fisika, sedangkan dalam kimia kuantum sering menggunakan metode interaksi konfigurasi (CI), yang didasarkan pada pengembangan sistematis fungsi gelombang dalam determinan Slater. Akan tetapi, permasalahan dengan metode ini adalah membutuhkan kemampuan komputasi yang sangat besar, yang membuatnya tidak mungkin untuk menerapkannya pada sistem kompleks yang lebih besar.
Dalam hal ini DFT memberikan suatu alternatif yang menarik, yang lebih bermanfaat karena DFT memberikan suatu cara sistematis pemetaan permasalahan banyak-partikel, dengan formula_5, menjadi permasalahan partikel-tunggal tanpa formula_5. Dalam DFT, variabel kunci adalah kerapatan partikel formula_12 yang diberikan oleh persamaan:
formula_13.
Hohenberg dan Kohn membuktikan pada tahun 1964 bahwa hubungan yang dinyatakan di atas dapat dibalikkan, yaitu menjadi kerapatan keadaan dasar formula_14. Prinsip ini memungkinkan untuk menghitung fungsi gelombang keadaan dasar yang bersesuaian formula_15. Yang berarti, formula_16 merupakan fungsi yang unik dari formula_17, yang diberikan oleh persamaan:
formula_18
Dan sebagai konsekuensinya, semua keadaan dasar yang dapat teramati formula_19 juga merupakan fungsi dari formula_17
formula_21
Dari sini dengan mengikuti fakta-fakta yang ada, bahwa energi keadaan dasar juga merupakan fungsi dari formula_17
formula_23,
Dengan kontribusi potensial eksternal formula_24 dapat ditulis secara eksplisit dalam persamaan kerapatan.
formula_25
Fungsi-fungsi formula_26 dan formula_27 disebut fungsi umum ‘universal’, sedangkan formula_28 merupakan fungsi yang tidak umum ‘tidak universal’, karena fungsi tersebut bergantung pada sistem yang sedang dipelajari. Setelah sebuah sistem ditetapkan, yaitu ... diketahui, langkah selanjutnya meminimalkan fungsi tersebut.
formula_29
Dengan mengacu pada formula_12, dengan asumsi memiliki ungkapan yang dapat dipercaya untuk formula_26 dan formula_27. Keberhasilan meminimalkan fungsi energi akan menghasilkan kerapatan keadaan dasar formula_17 dan semua keadaan dasar yang dapat teramati.
Permasalahan variasi dalam meminimalkan fungsi energi formula_34 dapat diselesaikan dengan menerapkan metode Lagrangian dari pengali yang tidak dapat ditentukan, yang telah dilakukan oleh Kohn dan Sham pada tahun 1965. Dengan metode ini, salah satu kegunaan nyata fungsi tersebut dalam persamaan di atas dapat ditulis sebagai fungsi kerapatan fiktif dari fungsi yang tidak berinteraksi.
formula_35,
Dengan formula_36 notasi untuk energi kinetik yang tidak berinteraksi dan formula_37 adalah potensial efektif eksternal di mana partikel berpindah. Secara nyata, formula_38 jika formula_37 dipilih untuk menjadi:
formula_40.
Sehingga, dapat menyelesaikan persamaan Kohn-Sham dari sistem pelengkap yang tidak saling berinteraksi.
formula_41,
Yang menghasilkan orbital formula_42 itu dihasilkan kembali kerapatan
formula_12 dari sistem banyak-partikel yang sebenarnya.
formula_44.
Potensial partikel-tunggal efektif formula_37 dapat ditulis lebih detail sebagai :
formula_46,
Dengan notasi persamaan kedua disebut persamaan Hartree yang menjelaskan tolakan ‘Coulomd’ antar elektron-elektron, sedangkan persamaan terakhir formula_47 disebut potensial hubungan pertukaran. Dengan formula_47 meliputi semua interaksi banyak partikel. Karena persamaan Hartree dan formula_47 bergantung pada formula_50, yang bergantung pada formula_42, yang bergantung pula pada formula_37, masalah penyelesaian persamaan Kohn-Sham dapat dilakukan dengan sebuah cara ""self-consistent"". Biasanya dimulai dengan perkiraan awal untuk formula_12, lalu menghitung kesesuaian formula_37 dan penyelesaian persamaan Kohn-Sham untuk formula_42. Dari perhitungan ini sebuah kerapatan baru dan mulai kembali. Prosedur ini diulang terus menerus sampai konvergensi dicapai.
Masalah utama dengan DFT adalah bahwa fungsi yang tepat untuk hubungan dan pertukaran tidak diketahui kecuali untuk gas elektron bebas. Akan tetapi, pendekatan yang ada membolehkan perhitungan kuantitas fisik tertentu secara akurat. Dalam fisika, pendekatan yang digunakan secara luas adalah pendekatan kerapatan lokal (LDA), dengan fungsi yang hanya bergantung pada kerapatan pada koordinat dengan fungsi yang telah dievaluasi.
formula_56
Pendekatan kerapatan spin lokal (LSDA) adalah penyetaraan secara langsung dari LDA yang meliputi spin elektron:
formula_57
Tingginya akurasi rumusan untuk kerapatan energi hubungan-pertukaran formula_58 telah dibangun dari simulasi gas elektron bebas.
Pendekatan penyetaraan gradien (GGA) masih tetap lokal, tetapi mempertimbangkan gradien dari kerapatan pada koordinat yang sama:
formula_59
Menggunakan (GGA) tersebut memberikan hasil yang sangat baik untuk geometri molekul dan energi keadaan dasar yang telah dicapai. Banyak tambahan peningkatan perhitungan telah dibuat untuk DFT dengan pengembangan yang lebih representatif dari fungsi.
Penyetaraan relativistik untuk formalisasi DFT mengarahkan pada teori fungsi kerapatan arus.
Dalam kenyataannya, teori Koh-Sham dapat diaplikasikan dalam 2 cara yang berbeda yang bergantung pada apa yang diteliti. Dalam keadaan padat, basis set gelombang datar digunakan pada syarat batas berkala. Selain itu, tekanan besar diletakkan pada konsistensi yang tepat dengan model idealisasi dari ‘gas elektron seragam’, yang memperlihatkan perilaku yang sama untuk keadaan padat yang tanpa batas. Dalam fase gas dan cair, tekanan ini sedikit mengalami pengurangan, karena gas elektron yang seragam tersebut merupakan model yang lemah untuk perilaku atom atau molekul yang diskrit. Karena batasan yang dikurangi, sebuah variasi besar dari fungsi hubungan-pertukaran telah dikembangkan untuk aplikasi-aplikasi kimia. Yang paling populer dari fungsi tersebut dikenal sebagai B3LYP. Parameter-parameter yang dapat terukur dari fungsi-fungsi tersebut umumnya dicobakan untuk ‘training set’ molekul. Walaupun hasil yang diperoleh dengan fungsi tersebut biasanya relatif akurat untuk beberapa aplikasi, tidak ada cara sistematis untuk meningkatkannya (berlawanan dengan beberapa metode yang didasarkan pada fungsi gelombang tradisional seperti interaksi konfigurasi atau penggabungan metode ""cluster""). Oleh karena itu, dalam pendekatan DFT saat ini tidak mungkin untuk menghitung kesalahan perhitungan tanpa membandingkannya dengan metode lain ataupun dengan eksperimen.