VALBOND - VALBOND

En mecánica molecular , VALBOND es un método para calcular la energía de flexión del ángulo que se basa en la teoría del enlace de valencia . Se basa en funciones de fuerza orbital , que se maximizan cuando los orbitales híbridos del átomo son ortogonales . La hibridación de los orbitales de enlace se obtiene a partir de fórmulas empíricas basadas en la regla de Bent , que relaciona la preferencia por el carácter p con la electronegatividad.

Las funciones VALBOND son adecuadas para describir la energía de la distorsión del ángulo de enlace no solo alrededor de los ángulos de equilibrio, sino también en distorsiones muy grandes. Esto representa una ventaja sobre la aproximación de oscilador armónico más simple utilizada por muchos campos de fuerza, y permite que el método VALBOND maneje moléculas hipervalentes y complejos de metales de transición . El término de energía VALBOND se ha combinado con campos de fuerza como CHARMM y UFF para proporcionar una forma funcional completa que incluye también estiramiento de enlaces, torsiones e interacciones no unidas.

Forma funcional

Moléculas no hipervalentes

Para un ángulo α entre enlaces normales (no hipervalentes) que involucran un orbital híbrido sp ^m d ⁿ , la contribución de energía es

${\ Displaystyle E (\ alpha) = k (S ^ {max} -S (\ alpha))}$ ,

donde k es un factor de escala empírico que depende de los elementos involucrados en el enlace, S ^max , la función de resistencia máxima , es

${\ Displaystyle S ^ {max} = {\ sqrt {\ frac {1} {1 + m + n}}} (1 + {\ sqrt {3m}} + {\ sqrt {5n}})}$

y S (α) es la función de fuerza

${\ Displaystyle S (\ alpha) = S ^ {max} {\ sqrt {1 - {\ frac {1 - {\ sqrt {1- \ Delta ^ {2}}}} {2}}}}}$

que depende de la integral de no ortogonalidad Δ:

${\ Displaystyle \ Delta = {\ frac {1} {1 + m + n}} \ left [1 + m \ cos \ alpha + {\ frac {n} {2}} (3 \ cos ^ {2} \ alfa -1) \ derecha]}$

La contribución de energía se agrega dos veces, una por cada uno de los orbitales de enlace involucrados en el ángulo (que pueden tener diferentes hibridaciones y diferentes valores de k ).

Para átomos de bloque p no hipervalentes, el valor de hibridación n es cero (sin contribución del orbital d), y m se obtiene como% p (1-% p), donde% p es el carácter p del orbital obtenido de

${\ Displaystyle \% p_ {i} = {\ frac {n_ {p} wt_ {i}} {\ sum _ {j} wt_ {j}}}}$

donde la suma de j incluye todos los ligandos, pares solitarios y radicales en el átomo, n _p es la "hibridación bruta" (por ejemplo, para un átomo "sp ² ", n _p = 2). El peso wt _i depende de los dos elementos involucrados en el enlace (o solo uno para un par solitario o radicales) y representa la preferencia por el carácter p de diferentes elementos. Los valores de los pesos son empíricos, pero pueden racionalizarse en términos de la regla de Bent.

Moléculas hipervalentes

Para las moléculas hipervalentes, la energía se representa como una combinación de configuraciones de VALBOND, que son similares a las estructuras de resonancia que colocan enlaces de cuatro electrones de tres centros (3c4e) de diferentes maneras. Por ejemplo, ClF ₃ se representa con un enlace de dos centros "normal" y un enlace 3c4e. Hay tres configuraciones diferentes para ClF ₃ , cada una con un enlace Cl-F diferente como enlace de dos centros. Para sistemas más complicados, el número de combinaciones aumenta rápidamente; SF ₆ tiene 45 configuraciones.

${\ Displaystyle E_ {tot} = \ sum _ {j} c_ {j} E_ {j}}$

donde la suma es sobre todas las configuraciones j , y el coeficiente c _j está definido por la función

${\ Displaystyle c_ {j} = {\ frac {\ Displaystyle \ prod _ {i = 1} ^ {bombo} \ cos ^ {2} \ alpha _ {i}} {\ Displaystyle \ sum _ {j = 1} ^ {config} \ prod _ {i = 1} ^ {bombo} \ cos ^ {2} \ alpha _ {i}}}}$

donde "hype" se refiere a los enlaces 3c4e. Esta función asegura que se favorezcan las configuraciones donde los enlaces 3c4e son lineales.

Los términos de energía se modifican multiplicándolos por un factor de orden de enlace, BOF, que es el producto de los órdenes de enlace formales de los dos enlaces involucrados en el ángulo (para enlaces 3c4e, el orden de enlace es 0.5). Para los enlaces 3c4e, la energía se calcula como

${\ Displaystyle E (\ alpha) = BOF \ times k _ {\ alpha} [1- \ Delta (\ alpha + \ pi) ^ {2}]}$

donde Δ es nuevamente la función de no ortogonalidad, pero aquí el ángulo α está compensado por 180 grados (π radianes).

Finalmente, para asegurar que se reproduce la preferencia axial frente a ecuatorial de diferentes ligandos en compuestos hipervalentes, se resta un término de "energía de compensación". Tiene la forma

${\ Displaystyle E_ {offset} = \ sum _ {i = 1} ^ {config} c_ {i} \ sum _ {j = 1} ^ {hype} {\ frac {EN_ {ija} + EN_ {ijb}} {2}}}$

donde los términos EN dependen de la diferencia de electronegatividad entre el ligando y el átomo central de la siguiente manera:

${\ displaystyle EN_ {ija} = 30 \ times (en_ {lig} -en_ {ca}) \ times ss}$

donde ss es 1 si la diferencia de electronegatividad es positiva y 2 si es negativa.

Para las moléculas hipervalentes de bloque p, no se utilizan orbitales d, por lo que n = 0. La contribución p m se estima a partir de métodos de química cuántica ab initio y un análisis de orbitales de enlace natural (NBO).

Extensión

Las extensiones más recientes, disponibles en el conjunto de códigos CHARMM , incluyen la influencia trans (o efecto trans ) dentro de VALBOND-TRANS y la posibilidad de ejecutar dinámicas moleculares reactivas con "VALBOND multiestado".

Referencias