PM3
PM3
PM3
Aprobado el 07-06-2014
INTRODUCCIN
Una de las tareas de la qumica computacional es la determinacin de estructuras y
propiedades moleculares mediante el uso de mecnica molecular, mtodos semiempricos y
teoras de orbitales moleculares a primeros principios (ab initio)2
a*
Facultad de Qumica e Ingeniera Qumica, Universidad Nacional Mayor de San Marcos, Ciudad
Universitaria, Av. Venezuela s/n Pab. B, Lima 100, Per. joracperu@yahoo.es
La gran dificultad para realizar clculos ms exactos como el mtodo ab initio y DFT (Teora
del funcional de densidad) para el estudio de molculas de mayor tamao han sido los mtodos
aproximados llamados semiempricos, esto debido a las restricciones tecnolgicas existentes3.
Entre los mtodos semiempricos desarrollados se encuentra el mtodo PM3 (Parameterized
Model number 3)4 . Este mtodo surgi por la reparametrizacin del mtodo AM1 (Modelo
Austin 1, llamado as por la Universidad de Texas en Austin), y est parametrado para los
elementos: H, C, Si, Ge, Sn, Pb, N, P, As, Sb, Bi, O, S, Se, Te, F, Cl, Br, I, Al, Ga, In, Tl, Be, Mg,
Zn, Cd, Hg5.
El mtodo PM3 toma la energa de valencia molecular al, que incluye la repulsin
Los mtodos tipo Dewar (MINDO, AM1, PM3) tratan a la molcula como una coleccin de
electrones de valencia y centros atmicos, donde cada centro (core) consiste de un tomo y sus
electrones internos. La aproximacin ms simple se puede escribir de la siguiente forma:
Donde CACB son las cargas del core A y B, adems ( SASA | SBSB ) es la integral de repulsin
electrnica que involucra los orbitales de valencia s de los tomos A y B; el ltimo trmino fAB
corresponde a un valor que vara segn el mtodo y que concuerda con los datos
experimentales.
El clculo de la entalpa para los mtodos tipo Dewar, entre ellos el PM3, se obtienen
mediante:
Donde HOMO es la energa del orbital molecular ocupado ms alto y LUMO es la energa del
orbital molecular desocupado ms bajo.
La presencia de una buena separacin del HOMO y LUMO en las molculas ha sido asociada
con una estabilidad estructural y cintica14, mientras que cuando la diferencia de energa
HOMO-LUMO es pequea es de esperar que el sistema sea poco estable y altamente reactivo15.
Pearson postul que las molculas que tienen una separacin entre HOMO y LUMO grande
tienden a ser por regla natural lo ms duras posibles e incrementan su estabilidad16.
Detalles computacionales
Se ha determinado la dureza (), la energa HOMO, energa LUMO, energa total del sistema a
0 K y a 298 K y el calor de formacin, de las geometras optimizadas de cada una de las tres
posibles estructuras de la figura 1. Los estudios computacionales fueron realizados por medio
de clculos semiempricos a nivel PM3 que se incluyen en el paquete de programas
Hyperchem 8.0 Profesional17, implementado en una computadora personal Intel (R) Atom
N450 1.66Ghz, 1Gb de RAM.
RESULTADOS Y DISCUSIN
Los datos experimentales1 indican que el complejo de plomo (II) tiene una estereoqumica 1:2
(metal: ligando).
El espectro IR del ligando muestra las bandas de absorcin a 1250 y 2507 cm-1 que
corresponde a la tensin v (C= S) y v (S-H); asimismo, se puede notar un hombro ancho que
corresponde a la seal del grupo NH a 3250 cm-1, lo que demuestra la coexistencia de dos
formas tautomricas con intercambio rpido del protn. En el espectro IR del complejo de
plomo, se puede distinguir la desaparicin de la seal del grupo NH y el grupo SH. En cambio
la seal del C=S se observa desplazado ligeramente a 1237 cm-1, confirmando de esta manera
la coordinacin del plomo con el ligando a travs de los tomos de S y del nitrgeno
hidraznico que se desprotona. En la figura 2, en el rango entre 3071 y 2729 cm-1 del espectro
del ligando se observa una serie de bandas de absorcin de estiramiento que corresponde a
CPh-H18,19,20,21 debido a la resonancia de los enlaces conjugados dentro del anillo aromtico y el
anillo heterocclico. En el espectro del complejo de plomo (II), slo se observa una seal
definida a 3068 cm-1 debido a que la formacin del nuevo enlace entre el plomo y los tomos
dadores S y N que genera un anillo quelato de cuatro miembros provoca la desaparicin de la
resonancia entre los enlaces conjugados y el anillo triazol.
En los espectros de RMN-H1 a diferentes temperaturas (-70C - 60C) en CDCl3 del ligando, se
observa la presencia de dos formas tautomricas. Se observa por ejemplo a 60C, la seal del
protn NH a 12.45 ppm y mientras que el protn CSH es detectado a 1.55 ppm.
En el espectro C13(H1) se nota un desplazamiento a campo bajo del grupo C=N (de 149,7 a
158,6 ppm) tpico de estructuras en donde el grupo C=N no est en conjugacin con otros
dobles enlaces. La transicin del grupo C=S se desplaza a campo alto (de 166,4 a 161,7 ppm),
ganando cierta proteccin debido a que la densidad electrnica, en ausencia de la capacidad de
resonancia, se concentra en los tomos de N-C-S y Pb consolidndose el enlace C=S. Por otra
parte, las seales correspondientes a los hidrgenos aromticos no han sufrido un
desplazamiento importante. La seal del hidrgeno amnico NH en el ligando (14,32 ppm) no
es observada en el complejo de plomo (II), lo que corrobora que se produjo la desprotonacin
del grupo NH hidraznico.
Estudio computacional
Se ha estudiado a nivel computacional por el mtodo semiemprico PM3 debido a que el
mtodo presenta parametrizacin para el tomo de plomo5 . Este mtodo es muy popular para
evaluar sistemas orgnicos (caso del ligando) y por ser muy preciso en clculos de energa y
longitud de enlace22.
En cada una de las estructuras de la figura 1: 1a), 1b) y 1c) se determin la energa HOMO-
LUMO, la dureza () energa a 0K y 298K adems del calor de formacin. Para realizar un
correcto modelado de las estructuras se tuvo que aadir el par libre presente en el tomo de
plomo, ya que existe influencia de este par libre en la estereoqumica de los complejos de
plomo (II)23,24 que presentan tomos dadores como l oxgeno25, azufre y selenio26.
En la figura 4 se muestra las geometras de coordinacin utilizadas para el modelar las tres
posibilidades de coordinacin vistas en la figura 1. El modelo I fue utilizado para modelar la
estructura figura 1a y el modelo II para las estructuras 1b y 1c. Estas geometras fueron
introducidas al programa y optimizadas hasta obtener la mnima energa para cada uno de los
casos.
En base a estas geometras optimizadas de mnima energa se obtuvieron las siguientes
propiedades mencionadas en la tabla 1.
Tabla 1. Energa total (a 0K, 298K), Calor de Formacin, Dureza (n),Energa HOMO,
Energa LUMO, Separacin HOMO-LUMO para las tres posibilidades de
coordinacin del complejos de Pb (II).
Se puede observar que de las tres estructuras propuestas, la estructura (1b) es la que menor
estabilidad termodinmica y dureza qumica posee, por lo tanto se podra considerar una
estructura bastante reactiva y poco estable, mientras que la estructura (1c) es la que posee
mayor estabilidad termodinmica, buena separacin HOMO-LUMO y mayor dureza de las
tres opciones; cumpliendo con el Principio de Mxima Dureza de Pearson en el cual se asocia
a la mayor separacin HOMO-LUMO con la mayor estabilidad estructural y cintica14.
Los datos experimentales muestran por medio de los espectros IR, C13 (H1)- APT y H1-RMN
que el ligando est totalmente desprotonado en el complejo, adems de que sus puntos de
coordinacin son por medio del tomo de azufre y del nitrgeno hidraznico (figura 1c).
En estudios estructurales con difraccin de rayos X de complejos de plomo (II)
tetracoordinados con tiosemicarbazona18 muestran la formacin de enlace Pb-S (2.730A) y Pb-
N (2.482 ), tambin los complejos de Pb (II) con 3-hexametileneiminiltiosemicarbazona19
exhiben puntos de coordinacin en los tomos de azufre Pb-S (2.7064 ) y nitrgeno
hidraznico (2.628 ). Esto verificara que la estructura (1c) es la que probablemente se
formara, ya que posee los mismos puntos de coordinacin antes mencionados, adems, los
clculos obtenidos muestran que dicha estructura posee mayor estabilidad y alta dureza global.
El principio de cido y base dura y blanda (HSAB)27 tambin nos puede ayudar a dilucidar cul
de las estructuras propuestas es la ms probable. Tomando en cuenta que el catin de plomo
Pb (II) es considerado como cido intermedio y que el ligando posee dos zonas de
coordinacin en los tomos de azufre (blando) y del nitrgeno hidraznico (duro) este catin
puede interactuar con ambas zonas de coordinacin o con las dos al mismo tiempo, dando un
complejo bastante estable28.
Otro aspecto a tomar en cuenta es la influencia del par libre en los complejos de Pb (II), ya que
el par libre provoca que aparezca una repulsin entre el par libre y los ligandos, mucho ms
fuerte, que entre los tomos del ligando.
Se propone entonces que debido a la influencia del par libre en la estereoqumica del complejo
de plomo (II) la geometra tetradrica no es la adecuada y que la geometra ms probable ser la
geometra disfenoidal (balancn) o bipirmide trigonal26.
CONCLUSIONES
Las estructuras propuestas para el complejo de plomo (II) con el ligando 1,5 difenil-1,2,4
triazol-3-tiona han sido evaluadas a nivel semiemprico utilizando el mtodo PM3 obteniendo
que la estructura (3c) propuesta con el par libre adicionado al tomo de plomo, present mayor
dureza y estabilidad termodinmica a comparacin de las dos estructuras restantes (3a y 3b),
cumpliendo as con el principio de mxima dureza de Pearson (PMH) y corroborando con lo
encontrado experimentalmente.
La propuesta de una geometra tetradrica para el complejo de plomo (II) parece no ser
correcta, ya que la influencia del par libre en los complejos de plomo (II) con tomos dadores
como el oxgeno, el azufre y el selenio, es importante y est avalada por la teora de repulsin
de pares de electrones de valencia (VSEPR) ya que el plomo no es un metal de transicin, este
estudio propone la geometra de balancn o bipirmide trigonal para el complejo Bis (1,5-
difenil-1,2, 4-triazol-3-tionato) plomo (II) (Pb(DTT)2).
REFERENCIAS
1. Angulo-Cornejo J.; Serrano-Ruiz M.; Henning L. Sntesis y caracterizacin de 1,5
difenil-1, 2, 4-triazol-3-tiona (HDTT) y Bis (1,5-difenil-1, 2, 4-triazol-3-tionato) Pb (II)
((Pb(DTT)2)). Rev. Soc. Quim 2003; 69(3): 131-142.
2. Clark, T.; Handbook of Computational Chemistry. New York: Wiley-Interscience; 1985.
3. Levine I.; Qumica Cuntica. 5ta edicin. Madrid: Pearson Prentice Hall; 2001
4. Stewar J.J.; Optimization of parameters for semiempirical methods I, II Method and
Applications. J. Compt. Chem 1989; 10 (2): 209-221.
5. Stewart J.J. Optimization of parameters for semiempirical methods. III. Extension of
PM3 to Be, Mg, Zn, Ga, Ge, As, Se, Cd, In, Sn, Sb, Te, Hg, Tl, Pb and Bi. J. Compt. Chem
1991; 12 (3): 320-341.