|
·
Evaluar
la exactitud del método CPHF aplicado al C4H6
en función de la utilización de diferentes bases, según
referente experimental.
·
Evaluar
la base más eficiente en términos de exactitud según referente
experimental contra tamaño para el tipo de moléculas del caso.
·
Obtener
una relación entre el largo de la molécula y las propiedades del
caso de modo de realizar extrapolaciones a moléculas mayores y en
particular al poliacetileno infinito (CH)¥.
En este trabajo se encara el estudio de polienos, que conforman una familia
con destacables propiedades NLO, propiedades que se incrementan con
el largo de su cadena molecular.
En
el momento de su publicación se contaba con datos experimentales sólo
para el C4H6, por lo que el cálculo de
propiedades para polienos más largos, y la deducción de una
posible función que permita extrapolaciones directas sobre las
propiedades de polienos de cualquier longitud son un aporte del todo
original.
Nótese que dada la simetría central de estas moléculas el
dipolo y la primera hiperpolarizabilidad son nulos, y dada las
particularidades de su geometría sólo hay cuatro valores
independientes no nulos para el tensor polarizabilidad y 7 para el
tensor segunda hiperpolarizabilidad.
Por ser también el C4H6 el único
polieno con propiedades obtenidas experimentalmente, esta molécula
será entonces la que dirija la selección de la base de cálculo más
eficiente.
Cartilla
de parámetros de cálculo
·
Software:
HONDO
·
Hardware:
Sistemas paralelos de computadoras ICAP-1 e ICAP-2.
·
Método
de cálculo de las propiedades NLO: CPHF.
·
Geometría:
ver Ilustración
V-2. Las coordenadas de los átomos fue tomada de una
optimización de geometría hecha en base 6-31G y simetría C2h.
·
Base:
Ver Tabla
V-6 y comentario.
·
Presición
de las integrales: 10-9 u.a.
·
Criterio
de dependencia lineal de funciones base: 10-5 u.a.
·
Criterio
de convergencia: 10-6 u.a. en la etapa SCF, 10-4
u.a. en la etapa CPHF.
Ilustración
V-2
Tabla
V-6
|
Etiqueta de
referencia
|
Base fundamental
|
Orbitales difusos y
exponentes según átomo
|
|
|
|
C
|
H
|
|
|
|
s
|
p
|
d
|
s
|
p
|
|
6-31G
|
6-31G
|
|
|
|
|
|
|
6-31G*
|
“
|
|
|
0.8
|
|
|
|
A
|
“
|
0.05
|
|
|
|
|
|
B
|
“
|
|
0.05
|
|
|
|
|
C
|
“
|
0.05
|
0.05
|
|
|
|
|
D
|
“
|
0.05
|
0.05
|
|
0.04
|
|
|
E
|
“
|
0.05
|
0.05
|
|
|
|
|
|
|
0.015
|
|
|
|
|
|
F
|
“
|
0.05
|
0.05
|
|
|
|
|
|
|
|
0.015
|
|
|
|
|
G
|
“
|
0.05
|
0.05
|
|
0.04
|
|
|
|
|
|
0.015
|
|
|
|
|
H
|
“
|
0.05
|
0.05
|
|
0.04
|
|
|
|
|
0.015
|
0.015
|
|
|
|
|
I
|
“
|
0.05
|
0.05
|
|
0.04
|
|
|
|
|
0.015
|
0.015
|
|
0.01
|
|
|
J
|
“
|
0.05
|
0.05
|
0.05
|
0.04
|
|
|
K
|
“
|
0.05
|
0.05
|
0.8
|
0.04
|
|
|
|
|
|
0.015
|
0.2
|
|
|
|
|
|
|
|
0.05
|
|
|
|
|
|
|
|
0.0125
|
|
|
|
L
|
“
|
0.05
|
0.05
|
0.8
|
0.04
|
0.9
|
|
|
|
|
0.015
|
0.2
|
|
0.1
|
|
|
|
|
|
0.05
|
|
0.01
|
|
|
|
|
|
0.0125
|
|
|
|
6-31G-PD
|
“
|
|
0.05
|
0.05
|
|
|
Como
ya expresamos en el Caso 1, el enriquecimiento de las bases estándar
con orbitales difusos es de relevante influencia en el cálculo de
las propiedades del caso. Es por ello que en su trabajo los autores
consideran una amplia g
de configuraciones de orbitales difusos aplicados sobre la base
6-31G estándar de modo de discernir entre ellas la más eficiente.
Tabla
V-7
|
Método
|
<a>
|
<g>
|
|
|
C4H6
|
C4H6
|
C6H8
|
|
Experimental dcSHG
|
58.3
|
27400
|
89700
|
En la Tabla V-7 se observan los
valores experimentales disponibles a tomar como referentes en los
casos correspondientes, los valores se expresan en u.a. Las
expresiones <a>
y <g>
corresponden a la suerte de promedios arrojados por los experimentos definidos por la Ecuación V-1.
Tabla
V-8
|
Etiqueta de base
|
<a>
|
Desviación del
referente experimental
|
<g>
|
Desviación del
referente experimental
|
|
6-31G
|
42.82
|
-27%
|
1098
|
-96%
|
|
6-31G*
|
43.24
|
-26%
|
1055
|
-96%
|
|
A
|
45.85
|
-21%
|
1908
|
| |