Kimia komputasi menggunakan model teoretis dan simulasi komputer untuk memprediksi struktur molekul, sifat, dan reaktivitas. Bidang ini mencakup hierarki metode yang menyeimbangkan akurasi dengan biaya komputasi. Pada tingkat tertinggi, metode ab initio memecahkan persamaan Schrödinger hanya menggunakan konstanta fisik fundamental. Teori Hartree-Fock (HF) menyediakan aproksimasi medan-rata-rata di mana setiap elektron bergerak dalam medan rata-rata semua lainnya. Metode pasca-HF (MP2, CCSD, CCSD(T)) memulihkan energi korelasi elektron, dengan CCSD(T) dianggap sebagai standar emas untuk molekul kecil. Teori fungsional kerapatan (DFT) menawarkan rasio akurasi-terhadap-biaya yang menguntungkan dengan menyatakan energi sebagai fungsi kerapatan elektron daripada fungsi gelombang banyak-elektron.
Fungsional DFT dan Himpunan Basis yang Populer
Fungsional DFT umum meliputi B3LYP (fungsional hibrida yang menggabungkan pertukaran HF eksak), M06 (fungsional Minnesota yang diparameterisasi untuk termokimia golongan-utama dan logam transisi), dan ωB97X-D (hibrida pemisah-jarak dengan dispersi empiris). Setiap fungsional memiliki kekuatan dan kelemahan: B3LYP bekerja baik untuk molekul organik tetapi buruk untuk kompleks logam transisi dan sistem terikat dispersi; fungsional keluarga M06 menangani logam lebih baik; ωB97X-D mencakup koreksi dispersi yang kritis untuk interaksi non-kovalen. Himpunan basis mengaproksimasi orbital molekul sebagai kombinasi linier fungsi basis. Himpunan basis minimal (STO-3G) memiliki sedikit fungsi tetapi akurasi buruk. Himpunan valensi-terpisah gaya Pople (6-31G*, 6-311+G**) menambahkan fungsi polarisasi dan difus. Himpunan konsisten korelasi Dunning (cc-pVDZ, cc-pVTZ, cc-pVQZ) memungkinkan konvergensi sistematis menuju batas himpunan basis lengkap, penting untuk perhitungan akurasi-tinggi.
Permukaan Energi Potensial dan Optimasi Geometri
Permukaan energi potensial (PES) menggambarkan bagaimana energi sistem molekul bervariasi dengan koordinat nuklir. Titik stasioner pada PES sesuai dengan geometri kesetimbangan (minimum) dan keadaan transisi (titik pelana orde-pertama). Algoritma optimasi geometri mengikuti gradien negatif energi (gaya) untuk menemukan minimum terdekat. Perhitungan Hessian (turunan kedua) mengkonfirmasi sifat titik stasioner: minimum memiliki semua frekuensi vibrasi positif, sedangkan keadaan transisi memiliki tepat satu frekuensi imajiner yang sesuai dengan koordinat reaksi. Jalur reaksi yang menghubungkan minimum melalui keadaan transisi ditemukan menggunakan metode seperti pita elastis terdorong (NEB) atau perhitungan koordinat reaksi intrinsik (IRC).
Analisis Frekuensi Vibrasi dan Termokimia
Analisis frekuensi vibrasi memberikan validasi struktural dan data termokimia. Frekuensi yang dihitung dapat dibandingkan dengan spektrum IR dan Raman eksperimental setelah menerapkan faktor skala (biasanya 0,96-0,98 untuk DFT) untuk mengoreksi anharmonisitas dan ketidaklengkapan himpunan basis. Termodinamika statistik dari perhitungan yang sama menghasilkan energi titik nol, entalpi, dan energi bebas Gibbs pada suhu yang ditentukan. Tidak adanya frekuensi imajiner mengkonfirmasi minimum sejati, sedangkan perhitungan IRC memverifikasi bahwa keadaan transisi menghubungkan reaktan dan produk yang benar. Koreksi termokimia ini penting untuk menghitung laju reaksi menggunakan teori keadaan transisi: k = (k_BT/h) exp(-ΔG‡/RT).
Model Solvasi
Sebagian besar kimia terjadi dalam larutan, membuat pemodelan solvasi menjadi penting. Model solvasi kontinum memperlakukan pelarut sebagai介质 dielektrik terpolarisasi yang mengelilingi zat terlarut dalam rongga. Model kontinum polarizable (PCM) membangun rongga dari bola atom yang tumpang tindih. Model SMD menambahkan suku jarak-pendek untuk akurasi yang lebih baik dengan pelarut tertentu. Model ini memperhitungkan polarisasi elektrostatik pelarut (respons cepat) dan, dalam beberapa implementasi, kontribusi kavitasi dan dispersi. Namun, model kontinum tidak dapat menangkap interaksi spesifik pelarut-zat terlarut seperti ikatan hidrogen dengan kulit solvasi pertama. Untuk kasus seperti itu, solvasi eksplisit menggunakan molekul pelarut diskrit dalam skema QM/MM atau pendekatan klaster-kontinum diperlukan.
Metode QM/MM dan Dinamika Molekul
Metode QM/MM (mekanika kuantum/mekanika molekul) membagi sistem menjadi daerah kuantum (diperlakukan dengan DFT atau ab initio) dan daerah klasik (diperlakukan dengan medan gaya). Pendekatan hibrida ini memungkinkan studi akurat tentang reaksi enzimatik, di mana situs aktif (substrat, residu katalitik, kofaktor logam) diperlakukan secara mekanika kuantum sementara protein sekitarnya diperlakukan dengan MM. Simulasi dinamika molekul (MD) merambatkan posisi atom dari waktu ke waktu dengan mengintegrasikan persamaan Newton menggunakan medan gaya klasik (AMBER, CHARMM, OPLS). MD mengeksplorasi ruang konformasi dan memberikan informasi termodinamika dan kinetika. Ensemble umum adalah NVT (jumlah, volume, suhu konstan) dan NPT (tekanan konstan). Teknik pengambilan sampel yang ditingkatkan (metadinamika, umbrella sampling) mengatasi hambatan energi yang tidak akan dilewati dalam skala waktu MD konvensional.
Aplikasi dalam Desain Obat dan Ilmu Material
Kimia komputasi sangat diperlukan dalam pengembangan farmasi. Penambatan molekul memprediksi bagaimana molekul kecil mengikat target protein dengan memberi skor jutaan pose terhadap situs pengikatan. Penapisan virtual menyaring pustaka senyawa besar secara komputasi sebelum pengujian eksperimental, secara dramatis mengurangi biaya. Gangguan energi bebas (FEP) dan integrasi termodinamika memberikan prediksi afinitas pengikatan yang sangat akurat untuk optimasi pemimpin. Dalam ilmu material, metode komputasi memprediksi celah pita, sifat mekanik, dan mekanisme reaksi. Penapisan komputasi throughput-tinggi mempercepat penemuan katalis baru, material baterai, dan senyawa fotovoltaik. Potensial antara atom yang dipelajari mesin, dilatih pada data DFT, sekarang memungkinkan akurasi mendekati-DFT pada kecepatan medan gaya, membuka perbatasan baru dalam kimia komputasi.
