El muy corto respuesta es: cada vez que usted está en duda.
La comprobación de estabilidad es probablemente innecesaria cuando se trata de lidiar con el cierre de la shell de moléculas, donde usted puede estar bastante seguro de que un restringido metodología es razonable. Por ejemplo, al calcular el agua en la ICR nivel de la teoría, usted realmente no necesita buscar una mayor estabilidad de la función de onda, porque no sería físicamente razonable. Una de UHF enfoque en principio se llevan a un menor de energía, ya que tiene más flexibilidad.
Cada vez que calcular las moléculas en las que se podría suponer que algunos concha abierta personaje debe considerar la comprobación de la función de onda de la estabilidad. Normalmente no hace daño a realizar una verificación demasiados. El descubrimiento de inestabilidad posterior puede ser mucho más frustrante. Que no significa necesariamente que luchan por el correcto estado del suelo, como usted todavía está restringido en la vuelta. Voy a incluir un pequeño, pero bien conocido el ejemplo a continuación: $\ce{O2}$. Estoy usando Gaussian 09 Apo. D. 01 para el cálculo, pero se debe trabajar con algo más reciente. (No funciona de la misma manera con nada más.) Tenga en cuenta también, que no estoy realizando optimizaciones de geometría, pero el uso de $\mathbf{d}(\ce{O-O})=\pu{120pm}$ a mantenerlo simple.
Primero vamos a ejecutar una ordinaria, restringido cálculo. Si usted no está utilizando el SCF(xqc) de palabras clave, usted debe inmediatamente el aviso de que algo está mal.
%chk=bp86svp.la icr.cantar.chk
#p RBP86/Def2SVP/W06
DenFit ! Densidad de uso ajuste a la velocidad de cálculo
gfinput ! Imprime conjunto de base en un formato adecuado para su uso como base general de entrada
gfoldprint ! Imprime la base de establecer la información en el formato de Gauss
la pio(6/7=3) ! Los interruptores de impresión de todos los MO
simetría(suelto) ! El uso correcto de la simetría
scf(xqc,conver=10) ! Cuadrática convergencia después convencional (para casos difíciles)
INT(SuperFineGrid) ! Más grande de la cuadrícula, más precisa
tarjeta de título no utilizados
0 1
O 0.6 0.0 0.0
O -0.6 0.0 0.0
Usted encontrará un resumen de las energías en la final, si usted desea comprobar y ejecutar los cálculos.
A continuación, vamos a comprobar y optimizar la función de onda.
%chk=bp86svp.la icr.cantar.puñalada.opt.chk
%oldchk=bp86svp.la icr.cantar.chk
#p RBP86/Def2SVP/W06
DenFit
geom=comprobar adivinar=leer
gfinput gfoldprint de la pio(6/7=3)
simetría(suelta)
scf(xqc,conver=10) INT(SuperFineGrid)
estable=opt
tarjeta de título no utilizados
0 1
Usted debe encontrar lo siguiente en su salida:
La función de onda tiene un ICR -> UHF inestabilidad.
Seguido por otro SCF ejecutar. Preste atención a los siguientes:
SCF Hecho: E(UB-P86) = -150.201660211 una.u. después de 6 ciclos
Convg = 0.3033 D-11 52 Fock formaciones.
S**2 = 1.0016 -V/T = 2.0049
= 0.0000 = 0.0000 = 0.0000 = 1.0016 S= 0.6187
= 0.000000000000 E+00
La aniquilación de la primera tirada de contaminantes:
S**2 antes de la aniquilación 1.0016, después de 0.0126
Usted todavía encontrar un estado singlete, pero ahora como UHF de la función de onda. Usted debe buscar más tarde:
La función de onda es estable bajo las perturbaciones consideradas.
Lo que han encontrado es una ruptura de la simetría singlete. Sabemos que el estado debe ser un triplete. Por lo tanto, vamos a realizar una restringido concha abierta de cálculo. Ya sabemos que esto va a tener y RHF -> UHF inestabilidad, y ya que no podemos usar la convergencia cuadrática, se puede utilizar estable. Pero nos proporcionan un buen punto de partida para la UHF de cálculo, o al menos de energía para comparar otros.
%chk=bp86svp.rohf.viaje.chk
%oldchk=bp86svp.la icr.cantar.chk
#p ROBP86/Def2SVP/W06
DenFit
geom=comprobar adivinar=alfa
gfinput gfoldprint de la pio(6/7=3)
simetría(suelta)
scf(conver=10) ! Cuadrática convergencia no está disponible para RO
INT(SuperFineGrid)
tarjeta de título no utilizados
0 3
La parte importante aquí es, que genera un puro spin estado onda sunction:
= 0.0000 = 0.0000 = 1.0000 = 2.0000 S= 1.0000
Otro cálculo es mediante la mezcla de HOMO y LUMO para generar una ruptura de la simetría ansatz de la función de onda directamente. Usted encontrará el mismo resultado que con el método anterior.
%chk=bp86svp.uhf.cantar.chk
%oldchk=bp86svp.la icr.cantar.chk
#p UBP86/Def2SVP/W06
DenFit
gfinput gfoldprint de la pio(6/7=3)
simetría(suelta)
geom=comprobar adivinar=(leer,mezcla)
scf(xqc,conver=10)
INT(SuperFineGrid)
estable=opt
tarjeta de título no utilizados
0 1
Por lo tanto, encontrar de nuevo:
SCF Hecho: E(UB-P86) = -150.201660211 A. U. después de 11 ciclos
NFock= 11 Conv=0.88 D-12 -V/T= 2.0049
= 0.0000 = 0.0000 = 0.0000 = 1.0016 S= 0.6187
= 0.000000000000 E+00
KE= 1.494731832724 D+02 PE=-4.114759744301 D+02 EE= 8.357834648903 D+01
La aniquilación de la primera tirada de contaminantes:
S**2 antes de la aniquilación 1.0016, después de 0.0126
La función de onda es estable bajo las perturbaciones consideradas.
Ahora para la última parte, ejecute un verdadero UHF de cálculo:
%chk=bp86svp.uhf.viaje.puñalada.opt.chk
%nproc=8
%mem=8000MB
%oldchk=bp86svp.la icr.cantar.chk
#p UBP86/Def2SVP/W06
DenFit
geom=comprobar adivinar=leer
gfinput gfoldprint de la pio(6/7=3)
simetría(suelta)
scf(xqc,conver=10)
INT(SuperFineGrid)
estable=opt
tarjeta de título no utilizados
0 3
Usted encontrará como se esperaba
La aniquilación de la primera tirada de contaminantes:
S**2 antes de la aniquilación 2.0031, después 2.0000
La función de onda es estable bajo las perturbaciones consideradas.
[...]
La función de onda es ya estable.
De risa, ¿por qué no considerar un quinteto?
Archivo de comandos de Energía Funcional / Hartree ( ciclos )
bp86svp.la icr.cantar E(RB-P86) = -150.155030787 ( 4 )
bp86svp.la icr.cantar.puñalada.opt E(UB-P86) = -150.201660211 ( 6 )
bp86svp.rohf.viaje E(ROB-P86) = -150.213771982 ( 10 )
bp86svp.uhf.cantar E(UB-P86) = -150.201660211 ( 11 )
bp86svp.uhf.viaje.puñalada.opt E(UB-P86) = -150.215679340 ( 9 )
bp86svp.uhf.quin.estable E(UB-P86) = -149.651916836 ( 9 )
En conclusión, cuando se sospecha de algún concha abierta y de carácter, de verificación para la estabilidad. Si usted tiene buen comportamiento de las moléculas (singlete), como la mayoría de los productos orgánicos son, probablemente puede pasar por alto.
Debe ser consciente de que la estabilidad de optimización no necesariamente te da el estado del suelo.
