Articulo de referencia

pesos de Chirgwin-Coulson

En los cálculos modernos de la teoría del enlace de valencia (VB) , los pesos de Chirgwin-Coulson (también llamados pesos de Mulliken) son los pesos relativos de un conjunto de ...

En los cálculos modernos de la teoría del enlace de valencia (VB) , los pesos de Chirgwin-Coulson (también llamados pesos de Mulliken) son los pesos relativos de un conjunto de posibles estructuras de enlace de valencia de una molécula. Métodos relacionados para encontrar los pesos relativos de las estructuras de enlace de valencia son el método de Löwdin [ 1 ] y los pesos inversos. [ 2 ]

Las pesas llevan el nombre de BH Chirgwin y Charles Coulson , quienes las desarrollaron en 1950. [ 3 ]

Fondo

Para una función de ondaΨ=idoiΦi{\displaystyle \Psi =\sum \limits _{i}C_{i}\Phi _{i}}dóndeΦ1,Φ2,,Φnorte{\displaystyle \Phi _{1},\Phi _{2},\dots,\Phi _{n}}son un conjunto ortogonal y linealmente independiente de orbitales base , el peso de un orbital constituyenteΨi{\displaystyle \Psi _{i}}seríadoi2{\displaystyle C_{i}^{2}}dado que la integral de superposición,Sij{\displaystyle S_{ij}}, entre dos funciones de ondaΨi,Ψj{\displaystyle \Psi _{i},\Psi _{j}}sería 1 parai=j{\displaystyle i=j}y 0 paraij{\displaystyle i\neq j}En la teoría del enlace de valencia , sin embargo, las estructuras generadas no son necesariamente ortogonales entre sí y, a menudo, presentan una superposición sustancial entre las dos estructuras. Por lo tanto, al considerar orbitales constituyentes no ortogonales (es decir, orbitales con superposición distinta de cero), los términos no diagonales en la matriz de superposición serían distintos de cero y deben incluirse para determinar el peso de un orbital constituyente. Un método para calcular el peso de un orbital constituyente,Φi{\displaystyle \Phi _{i}}, propuesto por Chirgwin y Coulson sería: [ 3 ]

Fórmula de Chirgwin-Coulson

Wi=doiΦi|Ψ=doijdojΨi|Ψj=jdoidojSij{\displaystyle {\begin{aligned}W_{i}&=C_{i}\langle \Phi _{i}\vert \Psi \rangle =C_{i}\sum \limits _{j}C_{j}\langle \Psi _{i}\vert \Psi _{j}\rangle \\&=\sum \limits _{j}C_{i}C_{j}S_{ij}\end{aligned}}}

La aplicación de la fórmula de Chirgwin-Coulson a un orbital molecular produce la población de Mulliken del orbital molecular. [ 4 ]

Formulación rigurosa

Determinación de estructuras de enlace de valencia

El método de Rumer

Las 20 posibles estructuras de enlace de valencia para el butadieno y sus círculos de Rumer.

Yuri Rumer propuso un método para crear un conjunto completo e independiente de estructuras de enlace de valencia para una molécula . Para un sistema con n electrones y n orbitales, el método de Rumer consiste en disponer los orbitales en un círculo y conectarlos entre sí con líneas que no se intersecan. [ 5 ] Las estructuras covalentes, o sin carga, se pueden crear conectando todos los orbitales entre sí. Las estructuras iónicas, o cargadas, para un átomo dado se pueden determinar asignando una carga a una molécula y siguiendo el método de Rumer. Para el caso del butadieno , se muestran las 20 posibles estructuras de Rumer, donde 1 y 2 son las estructuras covalentes, 3-14 son las estructuras monoiónicas y 15-20 son las estructuras diiónicas. Las estructuras de enlace de valencia resultantes se pueden representar mediante una combinación lineal de determinantes.|ab¯dod¯|{\displaystyle |a{\overline {b}}c{\overline {d}}|}, donde una letra sin línea superior indica un electrón conα{\displaystyle \alpha }espín, mientras que una letra con línea superior indica un electrón conβ{\displaystyle \beta }espín . La estructura VB para 1, por ejemplo, sería una combinación lineal de los determinantes.|12¯34¯|{\displaystyle |1{\overline {2}}3{\overline {4}}|},|21¯34¯|{\displaystyle |2{\overline {1}}3{\overline {4}}|},|12¯43¯|{\displaystyle |1{\overline {2}}4{\overline {3}}|}, y|21¯43¯|{\displaystyle |2{\overline {1}}4{\overline {3}}|}. Para una especie monoaniónica, la estructura VB para 11 sería una combinación lineal de|12¯44¯|{\displaystyle |1{\overline {2}}4{\overline {4}}|}y|21¯44¯|{\displaystyle |2{\overline {1}}4{\overline {4}}|}, a saber:

ϕ11=12(|12¯44¯|+|21¯44¯|){\displaystyle \phi _{11}={\frac {1}{\sqrt {2}}}(|1{\overline {2}}4{\overline {4}}|+|2{\overline {1}}4{\overline {4}}|)}

Representación matricial de estructuras VB

Una estructura VB arbitraria|φ1φ2¯φ3φ4¯|{\displaystyle |\varphi _{1}{\overline {\varphi _{2}}}\varphi _{3}{\overline {\varphi _{4}}}\dots |}que contienenorte{\displaystyle n}electrones, representados por los índices de electrones1,2,,norte{\displaystyle 1,2,\dots ,n}, ynorte{\displaystyle n}orbitales, representados porφ1,φ2,,φnorte{\displaystyle \varphi _{1},\varphi _{2},\dots ,\varphi _{n}}, puede representarse mediante el siguiente determinante de Slater :

|φ1φ2¯φ3φ4¯|=1norte¡|φ1(1)α(1)φ1(2)α(2)φ1(norte)α(norte)φ2(1)β(1)φ2(2)β(2)φ2(norte)β(norte)|{\displaystyle |\varphi _{1}{\overline {\varphi _{2}}}\varphi _{3}{\overline {\varphi _{4}}}\dots |={\frac {1}{\sqrt {n!}}}{\begin{vmatrix}\varphi _{1}(1)\alpha (1)&\varphi _{1}(2)\alpha (2)&\dots &\varphi _{1}(n)\alpha (n)\\\varphi _{2}(1)\beta (1)&\varphi _{2}(2)\beta (2)&\dots &\varphi _{2}(n)\beta (n)\\\vdots &\vdots &\ddots &\vdots \end{vmatrix}}}

Dóndeα(k){\displaystyle \alpha (k)}yβ(k){\displaystyle \beta (k)}representar unα{\displaystyle \alpha }oβ{\displaystyle \beta }giro en elkel{\displaystyle k^{\text{th}}}electrón, respectivamente. Para el caso de un sistema de dos electrones con orbitalesa{\displaystyle a}yb{\displaystyle b}, la estructura VB,|ab¯|{\displaystyle |a{\overline {b}}|}, puede representarse:|ab¯|=12|a(1)α(1)a(2)α(2)b(1)β(1)b(2)β(2)|{\displaystyle |a{\overline {b}}|={\frac {1}{\sqrt {2}}}{\begin{vmatrix}a(1)\alpha (1)&a(2)\alpha (2)\\b(1)\beta (1)&b(2)\beta (2)\end{vmatrix}}}

La evaluación del determinante arroja: [ 6 ]

|ab¯|=12(a(1)b(2)[α(1)β(2)]a(2)b(1)[α(2)β(1)]){\displaystyle |a{\overline {b}}|={\frac {1}{\sqrt {2}}}(a(1)b(2)[\alpha (1)\beta (2)]-a(2)b(1)[\alpha (2)\beta (1)])}

Definición de los pesos de Chirgwin-Coulson

Dada una función de ondaΨ=idoiΦi{\displaystyle \Psi =\sum \limits _{i}C_{i}\Phi _{i}}dóndeΦ1,Φ2,,Φnorte{\displaystyle \Phi _{1},\Phi _{2},\dots ,\Phi _{N}}es un conjunto completo y linealmente independiente de estructuras VB ydok{\displaystyle C_{k}}es el coeficiente de cada estructura, el peso de Chirgwin-CoulsonWK{\displaystyle W_{K}}de una estructura VBΦK{\displaystyle \Phi _{K}}se puede calcular de la siguiente manera: [ 3 ]

Wi=jdoidojΦi|Φj=jdoidojSij{\displaystyle W_{i}=\sum \limits _{j}C_{i}C_{j}\langle \Phi _{i}|\Phi _{j}\rangle =\sum \limits _{j}C_{i}C_{j}S_{ij}}

DóndeS{\displaystyle S}es la matriz de superposición que satisfaceΦi|Φj=Sij{\displaystyle \langle \Phi _{i}|\Phi _{j}\rangle =S_{ij}}.

Otros métodos para calcular los pesos de la estructura VB incluyen los pesos de Löwdin , dondeWiLowdin=j,kSij1/2dojSik1/2dok{\displaystyle W_{i}^{\text{Lowdin}}=\sum \limits _{j,k}S_{ij}^{1/2}C_{j}S_{ik}^{1/2}C_{k}}, [ 1 ] y pesos inversos, dondeWiinverso=1norte(doi2(S1)ii){\displaystyle W_{i}^{\text{inverse}}={\frac {1}{N}}{\bigg (}{\frac {C_{i}^{2}}{(S^{-1})_{ii}}}{\bigg )}}connorte{\displaystyle N}ser un factor de normalización definido pornorte=idoi2(S1)ii{\displaystyle N=\sum \limits _{i}{\frac {C_{i}^{2}}{(S^{-1})_{ii}}}}. [ 2 ] El uso de pesos de Löwdin e inversos es apropiado cuando los pesos de Chirgwin-Coulson superan 1 o son negativos. [ 2 ]

Descomposición de semideterminantes de orbitales moleculares

Dado un conjunto de orbitales moleculares,Ψ1,Ψ2,,Ψmetro{\displaystyle \Psi _{1},\Psi _{2},\dots ,\Psi _{m}}Para una molécula, consideremos el determinante de una población orbital dada, representada porDmes{\displaystyle D_{\text{MO}}}El determinante se puede escribir como el siguiente determinante de Slater:

Dmes=|Ψ1Ψ¯1Ψ2Ψ¯2|{\displaystyle D_{\text{MO}}=|\Psi _{1}{\overline {\Psi }}_{1}\Psi _{2}{\overline {\Psi }}_{2}\dots |}

Calcular el determinante explícitamente multiplicando esta expresión puede ser una tarea computacionalmente difícil, dado que cada orbital molecular está compuesto por una combinación de orbitales atómicos. Por otro lado, debido a que el determinante de un producto de matrices es igual al producto de determinantes, el determinante puede reagruparse en semideterminantes, uno de los cuales contiene solo electrones conα{\displaystyle \alpha }espín y el único con electrones deβ{\displaystyle \beta }giro, eso es:Dmes=hmesαhmesβ{\displaystyle D_{\text{MO}}=h_{\text{MO}}^{\alpha }h_{\text{MO}}^{\beta }}dóndehmesα=|ϕ1ϕ2|{\displaystyle h_{\text{MO}}^{\alpha }=|\phi _{1}\phi _{2}\dots |}yhmesβ=|ϕ¯1ϕ¯2|{\displaystyle h_{\text{MO}}^{\beta }=|{\overline {\phi }}_{1}{\overline {\phi }}_{2}\dots |}. [ 6 ] [ 7 ] [ 8 ]

Tenga en cuenta que cualquier orbital molecular dadoΨmes{\displaystyle \Psi _{\text{MO}}}puede escribirse como una combinación lineal de orbitales atómicosϕ1,ϕ2,,ϕnorte{\displaystyle \phi _{1},\phi _{2},\dots ,\phi _{n}}, es decir, para cadaΨi{\displaystyle \Psi _{i}}, existendoij{\displaystyle C_{ij}}de tal manera queΨi=jdoijϕj{\displaystyle \Psi _{i}=\sum \limits _{j}C_{ij}\phi _{j}}. Por lo tanto, el semideterminantehmesα{\displaystyle h_{\text{MO}}^{\alpha }}puede descomponerse aún más en los semideterminantes para un ordenamiento de orbitales atómicoshrα=|ϕ1,ϕ2,,ϕnorte|{\displaystyle h_{r}^{\alpha }=|\phi _{1},\phi _{2},\dots ,\phi _{n}|}correspondiente a una estructura VBr{\displaystyle r}. Como tal, el orbital molecularΨi{\displaystyle \Psi _{i}}puede representarse como una combinación de los semideterminantes de los orbitales atómicos,hmesα=rdorαhrα{\displaystyle h_{\text{MO}}^{\alpha }=\sum \limits _{r}C_{r}^{\alpha }h_{r}^{\alpha }}. El coeficientedorα{\displaystyle C_{r}^{\alpha }}se puede determinar evaluando la siguiente matriz: [ 6 ] [ 7 ] [ 8 ]

dorα=|do11do21donorte1do12do22donorte2do1nortedo2nortedonortenorte|{\displaystyle C_{r}^{\alpha }={\begin{vmatrix}C_{11}&C_{21}&\dots C_{n1}\\C_{12}&C_{22}&\dots C_{n2}\\\vdots &\vdots &\ddots \\C_{1n}&C_{2n}&\dots C_{nn}\\\end{vmatrix}}}

El mismo método puede utilizarse para evaluar el semideterminante para elβ{\displaystyle \beta }electrones,hmesβ{\displaystyle h_{\text{MO}}^{\beta }}. Por lo tanto, el determinanteDmes{\displaystyle D_{\text{MO}}}puede expresarse comoDmes=r,sdorαdorβhrαhsβ{\displaystyle D_{\text{MO}}=\sum \limits _{r,s}C_{r}^{\alpha }C_{r}^{\beta }h_{r}^{\alpha }h_{s}^{\beta }}, dónder,s{\displaystyle r,s}índice en todas las posibles estructuras VB. [ 6 ] [ 7 ] [ 8 ]

Ejemplos de cálculos para moléculas simples

Cálculos para la molécula de hidrógeno

La molécula de hidrógeno puede considerarse una combinación lineal de dosH{\displaystyle {\ce {H}}}1s{\displaystyle 1s}orbitales, indicados comoφ1{\displaystyle \varphi _{1}}yφ2{\displaystyle \varphi _{2}}. Las posibles estructuras VB paraH2{\displaystyle {\ce {H_2}}}son las dos estructuras covalentes,|φ1φ2¯|{\displaystyle |\varphi _{1}{\overline {\varphi _{2}}}|}y|φ2φ1¯|{\displaystyle |\varphi _{2}{\overline {\varphi _{1}}}|}indicadas como 1 y 2 respectivamente, así como las estructuras iónicas.|φ1φ1¯|{\displaystyle |\varphi _{1}{\overline {\varphi _{1}}}|}y|φ2φ2¯|{\displaystyle |\varphi _{2}{\overline {\varphi _{2}}}|}indicadas como 3 y 4 respectivamente, como se muestra a continuación.

Posibles configuraciones de espín de la molécula de hidrógeno

Debido a que las estructuras 1 y 2 representan enlaces covalentes en la molécula de hidrógeno y el intercambio de electrones de la estructura 1 produce la estructura 2, las dos estructuras covalentes se pueden combinar en una función de onda. Como tal, el modelo de Heitler-London para el enlace enH2{\displaystyle {\ce {H_2}}},ΦHL{\displaystyle \Phi _{HL}}, se puede utilizar en lugar de las estructuras VB|φ1φ2¯|{\displaystyle |\varphi _{1}{\overline {\varphi _{2}}}|}y|φ1¯φ2|{\displaystyle |{\overline {\varphi _{1}}}\varphi _{2}|}: [ 9 ]

ΦHL=|φ1φ2¯||φ1¯φ2|{\displaystyle \Phi _{HL}=|\varphi _{1}{\overline {\varphi _{2}}}|-|{\overline {\varphi _{1}}}\varphi _{2}|}

Donde el signo negativo surge de la antisimetría del intercambio electrónico. Como tal, la función de onda para laH2{\displaystyle {\ce {H_2}}}molécula,ΨH2{\displaystyle \Psi _{{\text{H}}_{2}}}, puede considerarse como una combinación lineal de la estructura de Heitler-London y las dos estructuras de enlace de valencia iónica.

ΨH2=do1ΦHL+do2|φ1φ1¯|+do3|φ2φ2¯|{\displaystyle \Psi _{{\text{H}}_{2}}=C_{1}\Phi _{HL}+C_{2}|\varphi _{1}{\overline {\varphi _{1}}}|+C_{3}|\varphi _{2}{\overline {\varphi _{2}}}|}

La matriz de solapamiento entre los orbitales atómicos entre las tres configuraciones de enlace de valenciaΦHL{\displaystyle \Phi _{HL}},|φ1φ1¯|{\displaystyle |\varphi _{1}{\overline {\varphi _{1}}}|}, y|φ2φ2¯|{\displaystyle |\varphi _{2}{\overline {\varphi _{2}}}|}se proporciona en la salida para los cálculos de enlace de valencia. A continuación se muestra un ejemplo de salida: [ 6 ]

S=|S11S21S22S31S32S33|=|10,7789042310,778904230,435432581|{\displaystyle S={\begin{vmatrix}S_{11}\\S_{21}&S_{22}\\S_{31}&S_{32}&S_{33}\\\end{vmatrix}}={\begin{vmatrix}1\\0.77890423&1\\0.77890423&0.43543258&1\\\end{vmatrix}}}

Hallar los autovectores de la matrizHmiS=0{\displaystyle H-ES=0}, dóndeH{\displaystyle H}es el hamiltoniano ymi{\displaystyle E}es energía debida al solapamiento orbital , produce el vector VBdo{\displaystyle {\vec {c}}}, que satisface: [ 10 ]

ΨH=do{ΦHL,|φ1φ1¯|,|φ2φ2¯|}=do1ΦHL+do2|φ1φ1¯|+do3|φ2φ2¯|{\displaystyle \Psi _{H}={\vec {c}}\{\Phi _{HL},|\varphi _{1}{\overline {\varphi _{1}}}|,|\varphi _{2}{\overline {\varphi _{2}}}|\}=C_{1}\Phi _{HL}+C_{2}|\varphi _{1}{\overline {\varphi _{1}}}|+C_{3}|\varphi _{2}{\overline {\varphi _{2}}}|}

Resolviendo para el vector VBdo{\displaystyle {\vec {c}}}El uso de la teoría del funcional de la densidad produce los coeficientesdo1=0,787469{\displaystyle C_{1}=0.787469}ydo2=do3=0,133870{\displaystyle C_{2}=C_{3}=0.133870}. Por lo tanto, se pueden calcular los pesos de Coulson-Chrigwin: [ 6 ]

W1=do12S11+do1do2S12+do1do3S13=0,784{\displaystyle W_{1}=C_{1}^{2}S_{11}+C_{1}C_{2}S_{12}+C_{1}C_{3}S_{13}=0.784}
W2=W3=0,108{\displaystyle W_{2}=W_{3}=0.108}

Para comprobar la consistencia, los pesos inversos se pueden calcular determinando primero la inversa de la matriz de superposición:

S1=|6.464493.50783.137393.50781.366123.13739|{\displaystyle S^{-1}={\begin{vmatrix}6.46449\\-3.5078&3.13739\\-3.5078&1.36612&3.13739\\\end{vmatrix}}}

A continuación, la constante de normalización.norte{\displaystyle N}se puede determinar:

norte=KdoK2(S1)KK=0,0185{\displaystyle N=\sum \limits _{K}{\frac {C_{K}^{2}}{(S^{-1})_{KK}}}=0.0185}

Los pesos finales son:W1=1norte(do12(S1)11)=0,803{\displaystyle W_{1}={\frac {1}{N}}{\bigg (}{\frac {C_{1}^{2}}{(S^{-1})_{11}}}{\bigg )}=0.803}, yW2=W3=0,098{\displaystyle W_{2}=W_{3}=0.098}.

De manera informal, los pesos calculados indican que la función de onda para laH2{\displaystyle {\ce {H_2}}}La molécula tiene una contribución menor de una especie iónica no predicha por un modelo estrictamente de orbitales moleculares para el enlace.

Cálculos para el ozono

Las 6 posibles estructuras VB y los correspondientes círculos de Rumer para el ozono.

Para determinar los pesos relativos de cada estructura de resonancia del ozono se requiere, en primer lugar, la determinación de las posibles estructuras VB paraO3{\displaystyle {\ce {O_3}}}. Considerando únicamente elpag{\displaystyle p}orbitales de oxígeno y etiquetando elpag{\displaystyle p}orbital en eliel{\displaystyle i^{\text{th}}}oxígenoϕi{\displaystyle \phi _{i}},O3{\displaystyle {\ce {O_3}}}tiene 6 posibles estructuras VB según el método de Rumer. Suponiendo que no hay solapamiento de orbitales atómicos, elkel{\displaystyle k^{\text{th}}}La estructura puede representarse mediante los determinantes.Φk{\displaystyle \Phi _{k}}: [ 6 ]

Φ1=12(|ϕ2ϕ2¯ϕ1ϕ3¯|+|ϕ2ϕ2¯ϕ3ϕ1¯|){\displaystyle \Phi _{1}={\frac {1}{\sqrt {2}}}(|\phi _{2}{\overline {\phi _{2}}}\phi _{1}{\overline {\phi _{3}}}|+|\phi _{2}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{1}}}|)}
Φ2=12(|ϕ1ϕ1¯ϕ2ϕ3¯|+|ϕ1ϕ1¯ϕ3ϕ2¯|){\displaystyle \Phi _{2}={\frac {1}{\sqrt {2}}}(|\phi _{1}{\overline {\phi _{1}}}\phi _{2}{\overline {\phi _{3}}}|+|\phi _{1}{\overline {\phi _{1}}}\phi _{3}{\overline {\phi _{2}}}|)}
Φ3=12(|ϕ1ϕ2¯ϕ3ϕ3¯|+|ϕ2ϕ1¯ϕ3ϕ3¯|){\displaystyle \Phi _{3}={\frac {1}{\sqrt {2}}}(|\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|+|\phi _{2}{\overline {\phi _{1}}}\phi _{3}{\overline {\phi _{3}}}|)}
Φ4=|ϕ1ϕ1¯ϕ2ϕ2¯|{\displaystyle \Phi _{4}=|\phi _{1}{\overline {\phi _{1}}}\phi _{2}{\overline {\phi _{2}}}|}
Φ5=|ϕ2ϕ2¯ϕ3ϕ3¯|{\displaystyle \Phi _{5}=|\phi _{2}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|}
Φ6=|ϕ1ϕ1¯ϕ3ϕ3¯|{\displaystyle \Phi _{6}=|\phi _{1}{\overline {\phi _{1}}}\phi _{3}{\overline {\phi _{3}}}|}

O3{\displaystyle {\ce {O_3}}}tiene los siguientes tres orbitales moleculares , uno donde todo el oxígenopag{\displaystyle p}Los orbitales están en fase, uno donde hay un nodo en el oxígeno central y otro donde todo el oxígenopag{\displaystyle p}Los orbitales están desfasados, como se muestra a continuación:

La molecularπ{\displaystyle \pi }órbitas del ozono, de izquierda a derecha,π1{\displaystyle \pi _{1}},π2{\displaystyle \pi _{2}}, yπ3{\displaystyle \pi _{3}}

Las funciones de onda para cada uno de los orbitales molecularesπi{\displaystyle \pi _{i}}se puede escribir como una combinación lineal de cada uno de los oxígenospag{\displaystyle p}orbitales de la siguiente manera: [ 6 ]

|π1π2π3|=|do11do12do13do21do22do23do31do32do33||ϕ1ϕ2ϕ3|=|0,3680,7640,3680,71000,7100,6140,6710,614||ϕ1ϕ2ϕ3|{\displaystyle {\begin{vmatrix}\pi _{1}\\\pi _{2}\\\pi _{3}\\\end{vmatrix}}={\begin{vmatrix}C_{11}&C_{12}&C_{13}\\C_{21}&C_{22}&C_{23}\\C_{31}&C_{32}&C_{33}\\\end{vmatrix}}{\begin{vmatrix}\phi _{1}\\\phi _{2}\\\phi _{3}\\\end{vmatrix}}={\begin{vmatrix}0.368&0.764&0.368\\0.710&0&-0.710\\0.614&-0.671&0.614\\\end{vmatrix}}{\begin{vmatrix}\phi _{1}\\\phi _{2}\\\phi _{3}\\\end{vmatrix}}}

Dóndedoij{\displaystyle C_{ij}}indica el coeficiente deϕj{\displaystyle \phi _{j}}en un orbital molecularπi{\displaystyle \pi _{i}}. Consideremos las contribuciones VB para el estado fundamental deO3{\displaystyle {\ce {O_3}}},|π1π1¯π2π2¯|{\displaystyle |\pi _{1}{\overline {\pi _{1}}}\pi _{2}{\overline {\pi _{2}}}|}Utilizando el método de semideterminantes, los semideterminantes para el estado fundamental son:

|ϕ1ϕ2|gramo=do11do12do21do22=0,542{\displaystyle |\phi _{1}\phi _{2}|_{g}={\begin{Vmatrix}C_{11}&C_{12}\\C_{21}&C_{22}\\\end{Vmatrix}}=-0.542}
|ϕ2ϕ3|gramo=do12do13do22do23=0,542{\displaystyle |\phi _{2}\phi _{3}|_{g}={\begin{Vmatrix}C_{12}&C_{13}\\C_{22}&C_{23}\\\end{Vmatrix}}=-0.542}
|ϕ1ϕ3|gramo=do11do13do21do23=0,523{\displaystyle |\phi _{1}\phi _{3}|_{g}={\begin{Vmatrix}C_{11}&C_{13}\\C_{21}&C_{23}\\\end{Vmatrix}}=-0.523}

Mediante el método de expansión de semideterminante, el coeficiente,doi{\displaystyle C_{i}}, para una estructura|ϕiϕj¯ϕkϕl¯|{\displaystyle |\phi _{i}{\overline {\phi _{j}}}\phi _{k}{\overline {\phi _{l}}}|}es:

|ϕiϕj¯ϕkϕl¯|=|ϕiϕk||ϕjϕl|{\displaystyle |\phi _{i}{\overline {\phi _{j}}}\phi _{k}{\overline {\phi _{l}}}|=|\phi _{i}\phi _{k}||\phi _{j}\phi _{l}|}

Lo cual implica que el estado fundamental tiene los siguientes coeficientes:

Ψgramo=0,416Φ1+0,400Φ2+0,400Φ3+0,294Φ4+0,294Φ5+0,274Φ6=0,294(|ϕ2ϕ2¯ϕ1ϕ3¯|+|ϕ2ϕ2¯ϕ3ϕ1¯|)+0,283(|ϕ1ϕ1¯ϕ2ϕ3¯|+|ϕ1ϕ1¯ϕ3ϕ2¯|)+0,283(|ϕ1ϕ2¯ϕ3ϕ3¯|+|ϕ2ϕ1¯ϕ3ϕ3¯|)+0,294|ϕ1ϕ1¯ϕ2ϕ2¯|+0,294|ϕ2ϕ2¯ϕ3ϕ3¯|+0,274|ϕ1ϕ1¯ϕ3ϕ3¯|{\displaystyle {\begin{aligned}\Psi _{g}&=-0.416\Phi _{1}+0.400\Phi _{2}+0.400\Phi _{3}+0.294\Phi _{4}+0.294\Phi _{5}+0.274\Phi _{6}\\&=-0.294(|\phi _{2}{\overline {\phi _{2}}}\phi _{1}{\overline {\phi _{3}}}|+|\phi _{2}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{1}}}|)+0.283(|\phi _{1}{\overline {\phi _{1}}}\phi _{2}{\overline {\phi _{3}}}|+|\phi _{1}{\overline {\phi _{1}}}\phi _{3}{\overline {\phi _{2}}}|)+0.283(|\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|+|\phi _{2}{\overline {\phi _{1}}}\phi _{3}{\overline {\phi _{3}}}|)+\\&\quad \quad 0.294|\phi _{1}{\overline {\phi _{1}}}\phi _{2}{\overline {\phi _{2}}}|+0.294|\phi _{2}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|+0.274|\phi _{1}{\overline {\phi _{1}}}\phi _{3}{\overline {\phi _{3}}}|\end{aligned}}}

Dada la siguiente matriz de solapamiento para los semideterminantes: [ 6 ]

S=||ϕ1ϕ2|||ϕ1ϕ2||ϕ1ϕ2|||ϕ1ϕ3||ϕ1ϕ3|||ϕ1ϕ3||ϕ1ϕ2|||ϕ2ϕ3||ϕ1ϕ3|||ϕ2ϕ3||ϕ2ϕ3|||ϕ2ϕ3||=|0,983770,126340,999930,008100,126340,98377|{\displaystyle S={\begin{vmatrix}\langle |\phi _{1}\phi _{2}|||\phi _{1}\phi _{2}|\rangle \\\langle |\phi _{1}\phi _{2}|||\phi _{1}\phi _{3}|\rangle &\langle |\phi _{1}\phi _{3}|||\phi _{1}\phi _{3}|\rangle \\\langle |\phi _{1}\phi _{2}|||\phi _{2}\phi _{3}|\rangle &\langle |\phi _{1}\phi _{3}|||\phi _{2}\phi _{3}|\rangle &\langle |\phi _{2}\phi _{3}|||\phi _{2}\phi _{3}|\rangle \end{vmatrix}}={\begin{vmatrix}0.98377\\0.12634&0.99993\\0.00810&0.12634&0.98377\end{vmatrix}}}

La superposición entre dos estructuras VB representadas por el producto de dos determinantes VB|ϕaϕb¯ϕdoϕd¯|||ϕwϕincógnita¯ϕyϕz¯|{\displaystyle \langle |\phi _{a}{\overline {\phi _{b}}}\phi _{c}{\overline {\phi _{d}}}|||\phi _{w}{\overline {\phi _{x}}}\phi _{y}{\overline {\phi _{z}}}|\rangle }se puede evaluar hallando el producto de la superposición entre los dos semideterminantes, es decir:

|ϕaϕb¯ϕdoϕd¯|||ϕwϕincógnita¯ϕyϕz¯|=(|ϕaϕdo|||ϕwϕy|)(|ϕbϕd|||ϕincógnitaϕz|){\displaystyle \langle |\phi _{a}{\overline {\phi _{b}}}\phi _{c}{\overline {\phi _{d}}}|||\phi _{w}{\overline {\phi _{x}}}\phi _{y}{\overline {\phi _{z}}}|\rangle =(\langle |\phi _{a}\phi _{c}|||\phi _{w}\phi _{y}|\rangle )(\langle |\phi _{b}\phi _{d}|||\phi _{x}\phi _{z}|\rangle )}

Por ejemplo, la superposición entre los orbitales|ϕ1ϕ2¯ϕ3ϕ3¯|{\displaystyle |\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|}y|ϕ1ϕ2¯ϕ2ϕ3¯|{\displaystyle |\phi _{1}{\overline {\phi _{2}}}\phi _{2}{\overline {\phi _{3}}}|}sería:

|ϕ1ϕ2¯ϕ3ϕ3¯|||ϕ1ϕ2¯ϕ2ϕ3¯|=(|ϕ1ϕ3|||ϕ1ϕ2|)(|ϕ2ϕ3|||ϕ2ϕ3|)=(0,12634)(0,98377)=0,12429{\displaystyle \langle |\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|||\phi _{1}{\overline {\phi _{2}}}\phi _{2}{\overline {\phi _{3}}}|\rangle =(\langle |\phi _{1}\phi _{3}|||\phi _{1}\phi _{2}|\rangle )(\langle |\phi _{2}\phi _{3}|||\phi _{2}\phi _{3}|\rangle )=(0.12634)(0.98377)=0.12429}

Los pesos de las estructuras de Lewis estándar paraO3{\displaystyle {\ce {O_3}}}seríaW(Ψ2){\displaystyle W(\Psi _{2})}yW(Ψ3){\displaystyle W(\Psi _{3})}Los pesos se pueden obtener calculando primero los pesos de Chirgwin-Coulson para sus determinantes constituyentes:

W(|ϕ1ϕ2¯ϕ3ϕ3¯|)=k0,283dok|ϕ1ϕ2¯ϕ3ϕ3¯|||Φk|=0,283[0,294(|ϕ1ϕ2¯ϕ3ϕ3¯|||ϕ2ϕ2¯ϕ1ϕ3¯|+|ϕ1ϕ2¯ϕ3ϕ3¯|||ϕ2ϕ2¯ϕ3ϕ1¯|)+0,283(|ϕ1ϕ2¯ϕ3ϕ3¯|||ϕ1ϕ1¯ϕ2ϕ3¯|+|ϕ1ϕ2¯ϕ3ϕ3¯|||ϕ1ϕ1¯ϕ3ϕ2¯|)+0,283(|ϕ1ϕ2¯ϕ3ϕ3¯|||ϕ1ϕ2¯ϕ3ϕ3¯|+|ϕ1ϕ2¯ϕ3ϕ3¯|||ϕ2ϕ1¯ϕ3ϕ3¯|)+0,294|ϕ1ϕ2¯ϕ3ϕ3¯|||ϕ1ϕ1¯ϕ2ϕ2¯|+0,294|ϕ1ϕ2¯ϕ3ϕ3¯|||ϕ2ϕ2¯ϕ3ϕ3¯|+0,274|ϕ1ϕ2¯ϕ3ϕ3¯|||ϕ1ϕ1¯ϕ3ϕ3¯|]=0,111{\displaystyle {\begin{aligned}W(|\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|)&=\sum \limits _{k}0.283C_{k}\langle |\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|||\Phi _{k}|\rangle \\&=0.283[-0.294(\langle |\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|||\phi _{2}{\overline {\phi _{2}}}\phi _{1}{\overline {\phi _{3}}}|\rangle +\langle |\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|||\phi _{2}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{1}}}|\rangle )+0.283(\langle |\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|||\phi _{1}{\overline {\phi _{1}}}\phi _{2}{\overline {\phi _{3}}}|\rangle +\langle |\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|||\phi _{1}{\overline {\phi _{1}}}\phi _{3}{\overline {\phi _{2}}}|\rangle )\\&\quad \quad +0.283(\langle |\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|||\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|\rangle +\langle |\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|||\phi _{2}{\overline {\phi _{1}}}\phi _{3}{\overline {\phi _{3}}}|\rangle )+0.294\langle |\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|||\phi _{1}{\overline {\phi _{1}}}\phi _{2}{\overline {\phi _{2}}}|\rangle +0.294\langle |\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|||\phi _{2}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|\rangle \\&\quad \quad +0.274\langle |\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|||\phi _{1}{\overline {\phi _{1}}}\phi _{3}{\overline {\phi _{3}}}|\rangle ]\\&=0.111\end{aligned}}}
W(|ϕ2ϕ1¯ϕ3ϕ3¯|)=W(|ϕ1ϕ1¯ϕ2ϕ3¯|)=W(|ϕ1ϕ1¯ϕ3ϕ2¯|)=0,111{\displaystyle W(|\phi _{2}{\overline {\phi _{1}}}\phi _{3}{\overline {\phi _{3}}}|)=W(|\phi _{1}{\overline {\phi _{1}}}\phi _{2}{\overline {\phi _{3}}}|)=W(|\phi _{1}{\overline {\phi _{1}}}\phi _{3}{\overline {\phi _{2}}}|)=0.111}

Los pesos para las estructuras de Lewis estándar serían la suma de los pesos de los determinantes constituyentes. Por lo tanto: [ 3 ]

W(Ψ2)=W(|ϕ1ϕ1¯ϕ2ϕ3¯|)+W(|ϕ1ϕ1¯ϕ3ϕ2¯|)=0,222{\displaystyle W(\Psi _{2})=W(|\phi _{1}{\overline {\phi _{1}}}\phi _{2}{\overline {\phi _{3}}}|)+W(|\phi _{1}{\overline {\phi _{1}}}\phi _{3}{\overline {\phi _{2}}}|)=0.222}
W(Ψ3)=W(|ϕ1ϕ2¯ϕ3ϕ3¯|)+W(|ϕ2ϕ1¯ϕ3ϕ3¯|)=0,222{\displaystyle W(\Psi _{3})=W(|\phi _{1}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{3}}}|)+W(|\phi _{2}{\overline {\phi _{1}}}\phi _{3}{\overline {\phi _{3}}}|)=0.222}

Esto se compara bien con los pesos de Chirgwin-Coulson reportados de 0,226 para la estructura de Lewis estándar del ozono en el estado fundamental. [ 8 ]

Para el estado dirradical,Ψ1{\displaystyle \Psi _{1}}, el peso es:

W(|ϕ2ϕ2¯ϕ1ϕ3¯|)=k0,294dok|ϕ2ϕ2¯ϕ1ϕ3¯||Φk|=0,106{\displaystyle W(|\phi _{2}{\overline {\phi _{2}}}\phi _{1}{\overline {\phi _{3}}}|)=\sum \limits _{k}-0.294C_{k}|\phi _{2}{\overline {\phi _{2}}}\phi _{1}{\overline {\phi _{3}}}||\Phi _{k}|=0.106}
W(|ϕ2ϕ2¯ϕ3ϕ1¯|)=0,106{\displaystyle W(|\phi _{2}{\overline {\phi _{2}}}\phi _{3}{\overline {\phi _{1}}}|)=0.106}
W(Ψ1)=W(|ϕ2ϕ¯2ϕ1ϕ¯3|)+W(|ϕ2ϕ¯2ϕ1ϕ¯3|)=0,106+0,106=0,212{\displaystyle W(\Psi _{1})=W(|\phi _{2}{\overline {\phi }}_{2}\phi _{1}{\overline {\phi }}_{3}|)+W(|\phi _{2}{\overline {\phi }}_{2}\phi _{1}{\overline {\phi }}_{3}|)=0.106+0.106=0.212}

Esto también se compara favorablemente con los pesos de Chirgwin-Coulson reportados de 0,213 para el estado dirradical del ozono en el estado fundamental. [ 8 ]

Aplicaciones a compuestos del grupo principal

Borazina

Borazina , ( fórmula química)B3norte3H6{\displaystyle {\ce {B_3N_3H_6}}}) es un compuesto cíclico y planar que es isoelectrónico con el benceno . Dado el par de electrones no enlazantes en el orbital p del nitrógeno fuera del plano y el orbital p vacío del boro, es posible la siguiente estructura de resonancia:

Mesómeros de borazina

Sin embargo, los cálculos de VB utilizando un conjunto de bases D95 de doble zeta indican que las estructuras de resonancia predominantes son la estructura con los tres pares de electrones no enlazantes en el nitrógeno (etiquetada como 1 a continuación) y las seis estructuras de resonancia con un doble enlace entre el boro y el nitrógeno (etiquetadas como 2 a continuación). Los pesos relativos de las dos estructuras son 0,17 y 0,08 respectivamente. [ 11 ] [ 12 ]

Estructuras de resonancia dominantes de la borazina

Por el contrario, las estructuras de resonancia dominantes del benceno son las dos estructuras de Kekulé, con un peso de 0,15, y 12 estructuras monozwitteriónicas con un peso de 0,03. Los datos, en conjunto, indican que, a pesar de la similitud en apariencia y estructura, los electrones en la borazina están menos deslocalizados que los del benceno. [ 11 ]

S 2 N 2

El dinitruro de disulfuro es un compuesto plano cuadrado que contiene un sistema conjugado de 6 electrones. π{\displaystyle \pi }sistema. A continuación se muestran las estructuras de resonancia dirradicales primarias (1 y 2) y una estructura zwitteriónica secundaria (3):

Estructuras de resonancia dominantes de S2N2

Los cálculos de enlace de valencia utilizando el conjunto de bases doble zeta completo D95 de Dunning indican que la estructura de resonancia dominante es el dirradical singlete con un enlace nitrógeno-nitrógeno largo (estructura 1), con un peso de Chirgwin-Coulson de 0,47. Este valor es sustancialmente mayor que el peso del dirradical singlete centrado en los azufres (estructura 2), que tiene un peso de Chirgwin-Coulson de 0,06. [ 13 ] Este resultado concuerda bien con las reglas generales sobre estructuras de Lewis, a saber, que las cargas formales deben minimizarse, y contrasta con resultados computacionales anteriores que indicaban que 1 es la estructura dominante. [ 14 ]

Referencias

  1. 1 2 Löwdin, Per-Olov (marzo de 1953). "Sobre la teoría de orbitales moleculares de compuestos orgánicos conjugados con aplicación al anillo de benceno perturbado". The Journal of Chemical Physics . 21 (3): 496– 515. Bibcode : 1953JChPh..21..496L . doi : 10.1063/1.1698934 . ISSN 0021-9606 . 
  2. 1 2 3 Gallup, GA; Norbeck, JM (septiembre de 1973). "Análisis de población de funciones de onda de enlace de valencia y BeH2". Chemical Physics Letters . 21 (3): 495– 500. Bibcode : 1973CPL....21..495G . doi : 10.1016/0009-2614(73)80292-1 . ISSN 0009-2614 . 
  3. 1 2 3 4 Chirgwin, BH; Coulson, CA (22 de marzo de 1950). "La estructura electrónica de los sistemas conjugados. VI". Actas de la Royal Society A: Ciencias matemáticas, físicas y de ingeniería . 201 (1065): 196– 209. Bibcode : 1950RSPSA.201..196C . doi : 10.1098/rspa.1950.0053 . ISSN 1364-5021 . 
  4. Mulliken, RS (octubre de 1955). "Análisis de población electrónica en funciones de onda moleculares LCAO–MO. I". The Journal of Chemical Physics . 23 (10): 1833– 1840. Bibcode : 1955JChPh..23.1833M . doi : 10.1063/1.1740588 . ISSN 0021-9606 . 
  5. ^ Rumer, Georg (1932). "Zur Theorie der Spinvalenz" . Göttinger Nachr. : 337.
  6. 1 2 3 4 5 6 7 8 9 Shaik, Sason; Hiberty, Philippe C. (16 de noviembre de 2007). Guía de un químico para la teoría del enlace de valencia . Hoboken, NJ, EE. UU.: John Wiley & Sons, Inc. doi : 10.1002/9780470192597 . ISBN 9780470192597.
  7. 1 2 3 Shaik, Sason; Hiberty, Philippe C. (17 de agosto de 2004). "Teoría del enlace de valencia, su historia, fundamentos y aplicaciones: una introducción". Reviews in Computational Chemistry . John Wiley & Sons, Inc. págs. 1–100 . doi : 10.1002/0471678856.ch1 . ISBN  9780471445258.
  8. 1 2 3 4 5 Hiberty, PC; Leforestier, C. (marzo de 1978). "Expansión de las funciones de onda de orbitales moleculares en funciones de onda de enlace de valencia. Un procedimiento simplificado". Journal of the American Chemical Society . 100 (7): 2012– 2017. Bibcode : 1978JAChS.100.2012H . doi : 10.1021/ja00475a007 . ISSN 0002-7863 . 
  9. ^ Heitler, W.; Londres, F. (junio de 1927). "Wechselwirkung neutraler Atome und homöopolare Bindung nach der Quantenmechanik". Zeitschrift für Physik . 44 ( 6– 7): 455– 472. Bibcode : 1927ZPhy...44..455H . doi : 10.1007/bf01397394 . ISSN 1434-6001 . 
  10. ^ J. Verbeek, JH Langenberg, CP Byrman, F. Dijkstra, JH van Lenthe, TORTUGA: Un programa Ab Initio VB/VBSCF (1998-2000)
  11. 1 2 Benker, Daniel; Klapötke, Thomas M.; Kuhn, Gerhard; Li, Jiabo; Miller, Christian (2005). "Un cálculo ab initio de enlace de valencia (VB) de la energía de deslocalización π en borazina, B3N3H6". Heteroatom Chemistry . 16 (5): 311– 315. doi : 10.1002/hc.20095 . ISSN 1042-7163 . 
  12. Engelberts, Jeroen Johan. “Análisis del enlace químico mediante la teoría del enlace de valencia ab initio”. Universidad de Utrecht , 2017.
  13. Klapoetke, Thomas M.; Li, Jiabo; Harcourt, Richard D. (2004-10-12). "Cálculos de enlace de valencia ab initio de doble zeta (D95) para los estados fundamentales de S2N2 y S42+". ChemInform . 35 (41). doi : 10.1002/chin.200441002 . ISSN 0931-7597 . 
  14. ^ Gerratt, J.; McNicholas, SJ; Karadakov, PB; Sironi, M.; Raimondi, M.; Cooper, DL (enero de 1996). "La extraordinaria estructura electrónica del N 2 S 2 ". Revista de la Sociedad Química Estadounidense . 118 (27): 6472– 6476. Bibcode : 1996JAChS.118.6472G . doi : 10.1021/ja953994f . ISSN 0002-7863 .