პარამეტრები
🔬 PySCF
ატომი x (Å) y (Å) z (Å)
მოლეკულის 3D სტრუქტურა
მოლეკულური ორბიტალების ენერგეტიკული დონეები
ელექტრონის სიმკვრივე ატომებზე (Mulliken)
ფიზიკური ინტუიცია — რას ითვლის DFT?

კვანტური მექანიკა გვეუბნება: ატომების სისტემის სრული აღწერა შეიძლება ტალღური ფუნქციით Ψ(r₁,r₂,...,rₙ). მაგრამ N ელექტრონის სისტემაში ეს 3N განზომილებიანი ობიექტია — პრაქტიკულად გამოუთვლელი N>10 შემთხვევაში. DFT-ის იდეა: ტალღური ფუნქციის ნაცვლად გამოვიყენოთ ელექტრონის სიმკვრივე ρ(r) — მხოლოდ 3 განზომილება! Walter Kohn-მა 1964 წელს დაამტკიცა, რომ სისტემის ყველა თვისება განისაზღვრება ρ(r)-ით. ეს არის DFT-ის საფუძველი.

ისტორიული კონტექსტი

1927 — Thomas-Fermi მოდელი: პირველი მცდელობა ელექტრონული სიმკვრივის გამოყენებისა — მარტივი, მაგრამ არაზუსტი. 1964 — Pierre Hohenberg და Walter Kohn: "Hohenberg-Kohn თეორემები" — DFT-ის მათემატიკური საფუძველი (Physical Review, 1964). 1965 — Walter Kohn და Lu Jeu Sham: "Kohn-Sham განტოლებები" — DFT-ის პრაქტიკული გამოთვლის სქემა. 1998 — Walter Kohn-მა მიიღო ნობელის პრემია ქიმიაში DFT-ის განვითარებისთვის. დღეს DFT ყველაზე გავრცელებული გამოთვლითი მეთოდია ქიმიასა და მასალათმცოდნეობაში.

Hohenberg-Kohn თეორემები

პირველი თეორემა: გარე პოტენციალი V_ext(r) ერთმნიშვნელოვნად განისაზღვრება ელექტრონის სიმკვრივით ρ(r). შედეგად — სრული ენერგია არის ρ(r)-ის ფუნქციონალი. მეორე თეორემა: ძირეული მდგომარეობის ენერგია მინიმალურია — ვარიაციული პრინციპი. სრული ენერგიის ფუნქციონალი:

$$ E[\rho] = T[\rho] + V_{ext}[\rho] + V_{ee}[\rho] $$

T[ρ] — კინეტიკური ენერგია, V_ext[ρ] — გარე პოტენციალი (ბირთვები), V_ee[ρ] — ელექტრონ-ელექტრონური ურთიერთქმედება (კულონი + გაცვლა-კორელაცია).

Kohn-Sham განტოლებები

Kohn-Sham (1965): რეალური ურთიერთქმედი ელექტრონების სისტემა ჩანაცვლდება ეფექტური პოტენციალის V_eff(r) მქონე არაურთიერთქმედი ელექტრონების სისტემით. ეს იძლევა ერთ-ელექტრონის შრედინგერის განტოლებებს:

$$ \left[-\frac{\hbar^2}{2m}\nabla^2 + V_{eff}(\mathbf{r})\right]\psi_i(\mathbf{r}) = \varepsilon_i\psi_i(\mathbf{r}) $$
$$ V_{eff}(\mathbf{r}) = V_{ext}(\mathbf{r}) + \int\frac{\rho(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}d\mathbf{r}' + V_{xc}[\rho] $$

V_eff(r) = V_ext(r) + V_Hartree(r) + V_xc(r). V_xc — გაცვლა-კორელაციის პოტენციალი (exchange-correlation) — DFT-ის ერთადერთი "უცნობი", სხვადასხვა ფუნქციონალებით ახლოვდება (LDA, GGA, hybrid).

ელექტრონის სიმკვრივე ρ(r)

ელექტრონის სიმკვრივე გამოითვლება დაკავებული ორბიტალებიდან:

$$ \rho(\mathbf{r}) = \sum_{i=1}^{N} |\psi_i(\mathbf{r})|^2 $$

ρ(r) — ელექტრონის პოვნის ალბათობის სიმკვრივე. ∫ρ(r)dr = N (ელექტრონთა სრული რაოდენობა). Mulliken population analysis — ρ ატომებზე განაწილება.

HOMO და LUMO — მოლეკულური ფრონტალური ორბიტალები

HOMO (Highest Occupied Molecular Orbital) — ყველაზე მაღალი დაკავებული ორბიტალი. LUMO (Lowest Unoccupied Molecular Orbital) — ყველაზე დაბალი თავისუფალი ორბიტალი. HOMO-LUMO სხვაობა = მოლეკულის "ელექტრონული Band Gap":

$$ E_g = E_{LUMO} - E_{HOMO} $$
ორბიტალი განმარტება როლი
HOMO ყველაზე მაღალი დაკავებული მოლეკულური ორბიტალი ელექტრონის დონორი, ნუკლეოფილობა
LUMO ყველაზე დაბალი თავისუფალი მოლეკულური ორბიტალი ელექტრონის აქცეპტორი, ელექტროფილობა
სიდიდეთა ცხრილი
სიმბოლო სახელი ერთეული
ρ(r)ელექტრონის სიმკვრივეe/ų
E[ρ]სრული ენერგიის ფუნქციონალიeV, Ha
ψᵢ(r)Kohn-Sham ორბიტალიÅ⁻³/²
εᵢორბიტალური ენერგიაeV
Vxcგაცვლა-კორელაციის პოტენციალიeV
Haჰარტრი — ენერგიის ატომური ერთეული27.2114 eV
a₀ბორის რადიუსი — სიგრძის ატომური ერთეული0.529 Å
კავშირი სხვა მოდულებთან

DFT → კრისტალური სტრუქტურა: DFT გამოთვლა ყოველთვის იწყება ატომების კოორდინატებით — ელემენტარული უჯრედის განსაზღვრით. DFT → ზონური სტრუქტურა: DFT იძლევა ზუსტ Band Structure-ს — tight-binding-ზე გაცილებით ზუსტს. DFT → NEGF: NEGF კვანტური მოძრაობის გამოთვლა იყენებს DFT-ის ჰამილტონიანს — ეს არის DFT+NEGF მიდგომა ნანომოწყობილობებისთვის.

რეალური გამოყენება

ფარმაცევტიკა: ახალი წამლების მოლეკულების დიზაინი — HOMO-LUMO სხვაობა განსაზღვრავს ქიმიურ რეაქტიულობას. კატალიზი: ადსორბციის ენერგიების გამოთვლა კატალიზატორის ზედაპირზე. ბატარეები: Li-ion ბატარეების ელექტროდული მასალების DFT მოდელირება. ნახევარგამტარები: Si, GaAs Band Gap-ის ზუსტი გამოთვლა მოწყობილობების დიზაინისთვის.

ნაბიჯ-ნაბიჯ გზამკვლევი
1
მოლეკულის არჩევა
ჩამოსაშლელი მენიუდან აირჩიეთ მოლეკულა. ატომების კოორდინატები ანგსტრომებში ავტომატურად ჩაიტვირთება.
2
კოორდინატების რედაქტირება
ცხრილში შეგიძლიათ x, y, z კოორდინატების შეცვლა. 3D მოლეკულა ავტომატურად განახლდება.
3
გამოთვლა (Celery)
"გამოთვლა" ღილაკი Celery task queue-ში გაუშვებს გამოთვლას. სტატუსი ყოველ წამს განახლდება.
4
ენერგეტიკული დონეები
მწვანე ხაზები — დაკავებული ორბიტალები (HOMO ჩათვლით). ცისფერი — თავისუფალი (LUMO). ყვითელი წყვეტილი ხაზი — HOMO-LUMO სხვაობა.
5
ელექტრონის სიმკვრივე
ბარ-გრაფიკი გვიჩვენებს Mulliken-ის ანალიზს — რამდენი ელექტრონი "ეკუთვნის" თითოეულ ატომს. მეტი ელექტრონი = უფრო ელექტრონეგატიური ატომი.
შედეგის ინტერპრეტაცია

HOMO-LUMO მცირე სხვაობა → მოლეკულა ქიმიურად რეაქტიულია. დიდი სხვაობა → სტაბილური, ნაკლებ რეაქტიული. H₂-ისთვის დიდი gap — ძალიან სტაბილური. ბენზოლისთვის მცირე gap — ელექტრონები დელოკალიზებულია.

თვითშეფასება

Q1. DFT = Density Functional Theory. "Density" რას ნიშნავს?

A მასალის სიმკვრივე (g/cm³)
B ელექტრონის სიმკვრივე ρ(r)
C ატომთა რაოდენობა
D ენერგიის სიმკვრივე

Q2. ვინ მიიღო ნობელის პრემია DFT-ისთვის?

A Albert Einstein
B Erwin Schrödinger
C Walter Kohn
D Niels Bohr

Q3. HOMO ნიშნავს:

A ყველაზე მაღალი დაკავებული მოლეკულური ორბიტალი
B ყველაზე დაბალი დაკავებული მოლეკულური ორბიტალი
C ყველაზე მაღალი თავისუფალი ორბიტალი
D ჰამილტონიანის ოპერატორი

Q4. DFT-ის მთავარი უპირატესობა კლასიკური კვანტური მექანიკის მეთოდებთან შედარებით?

A უფრო ზუსტია
B მარტივი ფორმულებია
C მხოლოდ ბირთვებს ითვლის
D 3N-ის ნაცვლად 3 განზომილება — გაცილებით სწრაფი გამოთვლა