a variacional del estado de transición a la - Páxinas persoais - USC ...
a variacional del estado de transición a la - Páxinas persoais - USC ...
a variacional del estado de transición a la - Páxinas persoais - USC ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
54 2.3. Métodos <strong>de</strong> <strong>la</strong> estructura electrónica<strong>la</strong> aproximación orbital que conduce a <strong>la</strong>s ecuaciones Hartree-Fock proporcionael resultado correcto. El hamiltoniano <strong>de</strong> este sistema sólo contendría lostérminos monoelectrónicos:Ĥ s =n∑ĥ(i) =i=1n∑− 1 2 ∇2 (i) +i=1n∑v s (i) (2.103)Para este sistema <strong>la</strong> función <strong>de</strong> onda <strong><strong>de</strong>l</strong> <strong>estado</strong> fundamental, Ψ s,0 ), se pue<strong>de</strong>escribir como un <strong>de</strong>terminante <strong>de</strong> S<strong>la</strong>ter <strong>de</strong> los espín-orbitales <strong>de</strong> Kohn-Shamψ i = φ KSi σ cuya parte espacial φ i KS resulta ser función propia <strong><strong>de</strong>l</strong> hamiltonianopropio <strong>de</strong> un sistema <strong>de</strong> partícu<strong>la</strong>s no interactuantes, esto es:(− 1 )2 ∇2 i + v s φ KSi = ǫ i φ KSi (2.104)A partir <strong>de</strong> dichos orbitales Khon-Sham se pue<strong>de</strong> <strong>de</strong>finir <strong>la</strong> <strong>de</strong>nsidad electrónicacomo,ρ = ∑ ∣∣φ KS∣ 2 i(2.105)ii=1obteniéndose una expresión para <strong>la</strong> energíaE v [ρ 0 ] =n∑i=1∫ǫ i = T s [ρ 0 ] +ρ(⃗r)v s (⃗r)d⃗r (2.106)Para un sistema real en el que los n electrones interaccionan, <strong>la</strong> ecuación <strong><strong>de</strong>l</strong>a energía tiene <strong>la</strong> forma:∫E v [ρ 0 ] = T [ρ 0 ] + ρ(⃗r)v n (⃗r)d⃗r + V ee [ρ 0 ] (2.107)Haciendo uso <strong>de</strong> este sistema <strong>de</strong> referencia y <strong>de</strong> <strong>la</strong>s siguientes igualda<strong>de</strong>s,δ ¯T [ρ 0 ] − ¯T s [ρ 0 ] (2.108)δ ¯V ee [ρ 0 ] = V ee [ρ] − 1 ∫ ∫ ρ(⃗r1 )ρ(⃗r 2 )d⃗r 1 d⃗r 22 r 12(2.109)así como <strong>de</strong> <strong>la</strong> función <strong>de</strong> <strong>la</strong> energía <strong>de</strong> corre<strong>la</strong>ción-intercambio (en inglésexchange-corre<strong>la</strong>tion):E xc = δ ¯T [ρ 0 ] + δ ¯V [ρ 0 ] (2.110)Kohn y Sham fueron capaces <strong>de</strong> reescribir <strong>la</strong> ecuación que <strong>de</strong>fine el funcional<strong>de</strong> <strong>la</strong> energía como sigue:∫E 0 = E v [ρ] = ρ(⃗r)v(⃗r)d(⃗r) + T s [ρ] + 1 ∫ ∫ ρ(⃗r1 )ρ(⃗r 2 )d⃗r 1 d⃗r 2 + E xc [ρ]2 r 12(2.111)Los tres primeros términos <strong>de</strong> esta ecuación son conocidos si po<strong>de</strong>mos obtenerρ a partir <strong>de</strong> los orbitales <strong>de</strong> Khom-Sham. Por lo tanto si conocemos cuál es