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 ondadóndeson un conjunto ortogonal y linealmente independiente de orbitales base , el peso de un orbital constituyenteseríadado que la integral de superposición,, entre dos funciones de ondasería 1 paray 0 paraEn 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,, propuesto por Chirgwin y Coulson sería: [ 3 ]
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

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., donde una letra sin línea superior indica un electrón conespín, mientras que una letra con línea superior indica un electrón conespín . La estructura VB para 1, por ejemplo, sería una combinación lineal de los determinantes.,,, y. Para una especie monoaniónica, la estructura VB para 11 sería una combinación lineal dey, a saber:
Representación matricial de estructuras VB
Una estructura VB arbitrariaque contieneelectrones, representados por los índices de electrones, yorbitales, representados por, puede representarse mediante el siguiente determinante de Slater :
Dóndeyrepresentar unogiro en elelectrón, respectivamente. Para el caso de un sistema de dos electrones con orbitalesy, la estructura VB,, puede representarse:
La evaluación del determinante arroja: [ 6 ]
Definición de los pesos de Chirgwin-Coulson
Dada una función de ondadóndees un conjunto completo y linealmente independiente de estructuras VB yes el coeficiente de cada estructura, el peso de Chirgwin-Coulsonde una estructura VBse puede calcular de la siguiente manera: [ 3 ]
Dóndees la matriz de superposición que satisface.
Otros métodos para calcular los pesos de la estructura VB incluyen los pesos de Löwdin , donde, [ 1 ] y pesos inversos, dondeconser un factor de normalización definido por. [ 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,Para una molécula, consideremos el determinante de una población orbital dada, representada porEl determinante se puede escribir como el siguiente determinante de Slater:
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 conespín y el único con electrones degiro, eso es:dóndey. [ 6 ] [ 7 ] [ 8 ]
Tenga en cuenta que cualquier orbital molecular dadopuede escribirse como una combinación lineal de orbitales atómicos, es decir, para cada, existende tal manera que. Por lo tanto, el semideterminantepuede descomponerse aún más en los semideterminantes para un ordenamiento de orbitales atómicoscorrespondiente a una estructura VB. Como tal, el orbital molecularpuede representarse como una combinación de los semideterminantes de los orbitales atómicos,. El coeficientese puede determinar evaluando la siguiente matriz: [ 6 ] [ 7 ] [ 8 ]
El mismo método puede utilizarse para evaluar el semideterminante para elelectrones,. Por lo tanto, el determinantepuede expresarse como, dóndeí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 dosorbitales, indicados comoy. Las posibles estructuras VB parason las dos estructuras covalentes,yindicadas como 1 y 2 respectivamente, así como las estructuras iónicas.yindicadas como 3 y 4 respectivamente, como se muestra a continuación.

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 en,, se puede utilizar en lugar de las estructuras VBy: [ 9 ]
Donde el signo negativo surge de la antisimetría del intercambio electrónico. Como tal, la función de onda para lamolécula,, puede considerarse como una combinación lineal de la estructura de Heitler-London y las dos estructuras de enlace de valencia iónica.
La matriz de solapamiento entre los orbitales atómicos entre las tres configuraciones de enlace de valencia,, yse proporciona en la salida para los cálculos de enlace de valencia. A continuación se muestra un ejemplo de salida: [ 6 ]
Hallar los autovectores de la matriz, dóndees el hamiltoniano yes energía debida al solapamiento orbital , produce el vector VB, que satisface: [ 10 ]
Resolviendo para el vector VBEl uso de la teoría del funcional de la densidad produce los coeficientesy. Por lo tanto, se pueden calcular los pesos de Coulson-Chrigwin: [ 6 ]
Para comprobar la consistencia, los pesos inversos se pueden calcular determinando primero la inversa de la matriz de superposición:
A continuación, la constante de normalización.se puede determinar:
Los pesos finales son:, y.
De manera informal, los pesos calculados indican que la función de onda para laLa 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

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 para. Considerando únicamente elorbitales de oxígeno y etiquetando elorbital en eloxígeno,tiene 6 posibles estructuras VB según el método de Rumer. Suponiendo que no hay solapamiento de orbitales atómicos, elLa estructura puede representarse mediante los determinantes.: [ 6 ]
tiene los siguientes tres orbitales moleculares , uno donde todo el oxígenoLos orbitales están en fase, uno donde hay un nodo en el oxígeno central y otro donde todo el oxígenoLos orbitales están desfasados, como se muestra a continuación:

Las funciones de onda para cada uno de los orbitales molecularesse puede escribir como una combinación lineal de cada uno de los oxígenosorbitales de la siguiente manera: [ 6 ]
Dóndeindica el coeficiente deen un orbital molecular. Consideremos las contribuciones VB para el estado fundamental de,Utilizando el método de semideterminantes, los semideterminantes para el estado fundamental son:
Mediante el método de expansión de semideterminante, el coeficiente,, para una estructuraes:
Lo cual implica que el estado fundamental tiene los siguientes coeficientes:
Dada la siguiente matriz de solapamiento para los semideterminantes: [ 6 ]
La superposición entre dos estructuras VB representadas por el producto de dos determinantes VBse puede evaluar hallando el producto de la superposición entre los dos semideterminantes, es decir:
Por ejemplo, la superposición entre los orbitalesysería:
Los pesos de las estructuras de Lewis estándar paraseríayLos pesos se pueden obtener calculando primero los pesos de Chirgwin-Coulson para sus determinantes constituyentes:
Los pesos para las estructuras de Lewis estándar serían la suma de los pesos de los determinantes constituyentes. Por lo tanto: [ 3 ]
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,, el peso es:
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)) 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:

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 ]

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. sistema. A continuación se muestran las estructuras de resonancia dirradicales primarias (1 y 2) y una estructura zwitteriónica secundaria (3):

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 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 .
- 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 .
- 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 .
- ↑ 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 .
- ^ Rumer, Georg (1932). "Zur Theorie der Spinvalenz" . Göttinger Nachr. : 337.
- 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.
- 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.
- 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 .
- ^ 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 .
- ^ J. Verbeek, JH Langenberg, CP Byrman, F. Dijkstra, JH van Lenthe, TORTUGA: Un programa Ab Initio VB/VBSCF (1998-2000)
- 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 .
- ↑ Engelberts, Jeroen Johan. “Análisis del enlace químico mediante la teoría del enlace de valencia ab initio”. Universidad de Utrecht , 2017.
- ↑ 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 .
- ^ 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 .
- Química