cscf: en línea en la nube

PROGRAMA:

NOMBRE

cscf: resuelve las ecuaciones de Hartree-Fock

DESCRIPCIÓN

El programa cscf realiza el procedimiento iterativo para resolver las ecuaciones de Hartree-Fock.

Este programa está restringido a la simetría D2h y sus subgrupos y las ocupaciones orbitales
deben ser números enteros. Por lo tanto, ciertos estados de momento angular puro derivados de
La ocupación parcial de orbitales degenerados no se puede obtener con los códigos actuales. Para
ejemplo, el estado 2PIu (doblete PI u) de ONO lineal derivado de la energía más baja
La configuración lineal (pi u) 1 solo se puede calcular como 2B2u (doblete B2u) o 2B3u
(doblete B 3u) componente del estado 2PIu (doblete PI u), y el resultado espacial
La función de onda no tendrá simetría PI. En cierto sentido, sin embargo, esto es deseable,
ya que la energía será una función continua del ángulo de flexión. Calculando la energía
de configuraciones dobladas como 2B2u (doblete B 2u) o 2B3u (doblete B 3u) y haciendo un 2PIu puro
El estado (doblete PI u) en geometrías lineales da como resultado una discontinuidad pronunciada.

En su mayor parte, estados tripletes resultantes de la doble ocupación de un doble degenerado
orbital, como el estado 3A2 (triplete A 2) resultante de (e ') 2 o (e ") 2
configuraciones en simetría D3h, o el estado 3SIGMAg (triplete SIGMA g) de a (pi g) 2 o (pi
u) 2 configuraciones en simetría Dinfh (D infinito h), tendrán la simetría espacial adecuada.
Los estados singlete que resultan de estas mismas configuraciones electrónicas son inherentemente
multiconfiguración y, como tales, no están bien representados por una configuración única
funciones de onda.

Referencias

Método de archivo PK:

1. RC Raffenetti, Chem. Phys. Letón. 20 (1973) 335.

Cálculos de simetría molecular y HF de capa cerrada:

1. M.Dupuis y HFKing, Int. J. Quant. Chem. 11 (1977) 613.

DIIS para carcasa cerrada:

1. P. Pulay, Chem. Phys. Letón. 73 (1980) 393.

2. P. Pulay, J. Comp. química 3 (1982) 556.

Coeficientes de acoplamiento (alfa y beta) para carcasa abierta:

1. CCJ Roothaan, Rev. Mod. Phys. 32 (1960) 179.

Mojadura:

1. DR Hartree, "El cálculo de estructuras atómicas" (Wiley: Nueva York) 1957.

2. MC Zerner y M. Hehenberger, Chem. Phys. Letón. 62 (1979) 550.

Cambio de nivel:

1. VR Saunders e IH Hillier, Int. J. Quant. Chem. 7 (1973) 699.

CONVERGENTE CSCF

Para casos de cáscara abierta difíciles, se recomienda que una cáscara cerrada adecuada
El cálculo se ejecute primero (agregue o elimine un electrón extra) y que este vector SCF luego
utilizarse como suposición para la función de onda de capa abierta deseada. Para los casos de TCSCF, siempre es
Es prudente ejecutar un shell cerrado (o quizás el triplete apropiado) SCF primero y luego usar
esto como una suposición para el TCSCF.

Para sistemas de carcasa abierta, se recomienda un valor de cambio de nivel de 0.5 a 3.0. Empiece con un
valor alto (2.0 - 3.0) para el primer cálculo de SCF y luego redúzcalo (a 0.5 - 1.0) para
ejecuciones posteriores que utilizan un vector SCF convergente como punto de partida.

Es extremadamente importante tener en cuenta que esta versión del código ya no admite
DE TIPO ABIERTO. Un deben use la nuevas las palabras claves REFERENCIA y MULTIP a especificar la tipo of SCF
necesario.

ENTRADA FORMATO

La pestaña cscf El programa busca a través de la ruta de palabra clave predeterminada (primero SCF y luego en DEFAULT)
para las siguientes palabras clave:

ETIQUETA = cadena
Esta es una cadena de caracteres que se incluirá en la salida. Esta cadena no se usa
por el programa. No hay ningún valor predeterminado.

WFNMás = cadena
Este es el tipo de función de onda que se desea en última instancia. El valor predeterminado es SCF.

DE TIPO ABIERTO is no longer apoyadas

REFERENCIA = cadena
Esto especifica el tipo de cálculo SCF que se desea hacer. Puede ser uno de FHD
(para una camiseta de concha cerrada), ROHF (para un cálculo de concha abierta restringida), UHF
(para un cálculo de capa abierta sin restricciones), DOSCON (para una configuración de dos
singlete), o SPECIAL. Si SPECIAL se da, luego el acoplamiento alfa y beta
los coeficientes se deben dar con el ALPHA y BETA palabras clave. El valor predeterminado es FHD.

MULTP = entero
Especifica la multiplicidad de la molécula. El valor predeterminado es singlete.

CARGA = entero
Especifica la carga de la molécula. Defauly es 0.

DOCC = vector_entero
Esto da el número de orbitales doblemente ocupados en cada irreducible
representación. No hay ningún valor predeterminado. Si no se da, CSCF intentará
adivinar las ocupaciones utilizando el núcleo hamiltoniano.

SOCC = vector_entero
Esto da el número de orbitales ocupados individualmente en cada irreducible
representación. No hay ningún valor predeterminado.

DERTIPO = cadena
Esto especifica el orden de la derivada que finalmente se realizará. Esta usado
según el scf programa para determinar si ciertos archivos se van a escribir y también es
utilizado para determinar la convergencia predeterminada de la función de onda. El valor predeterminado es
PRIMERO.

MAXITER = entero
Esto da el número máximo de iteraciones. El valor predeterminado es 40.

CONVERGENCIA = entero
Esto especifica qué tan estrechamente convergerá la función de onda. La convergencia es
determinado comparando el cambio RMS en la matriz de densidad ("delta P") con el
valor dado. El criterio de convergencia es 10 ** (-entero). El valor predeterminado es 7 si
ambas DERTIPO = NINGUNO y WFNMás = SCF se dan y 10 en caso contrario.

CAMBIO DE NIVEL = real
Esto especifica el cambio de nivel. El valor predeterminado es 1.

DIRECTO = booleano
Especifica si se debe realizar el cálculo de SCF con una técnica directa integral. los
el valor predeterminado es falso.

IMPRIMIR_MOS = booleano
Especifica si imprimir los orbitales moleculares o no. El valor predeterminado es falso.

También hay una gran cantidad de parámetros de entrada que se utilizan con menos frecuencia. Si no lo hace
comprender lo que significan las siguientes opciones, luego asegúrese de que no aparezcan en su
aporte. Los valores predeterminados funcionarán en la inmensa mayoría de los casos. Estos se especifican
con las siguientes palabras clave:

BORRAR_INTS = booleano
Los archivos integrales se borrarán si WFNMás = SCF y DERTIPO = PRIMERO or DERTIPO = NINGUNO.
Si desea mantener archivos integrales, configure BORRAR_INTS = falso. El valor predeterminado es
verdadera.

REORDENAR = cadena
El parámetro controla el reordenamiento de los orbitales moleculares. Si se establece en ANTES entonces el
adivinar los orbitales del archivo de punto de control se reordenarán. Si se establece en DESPUÉS DE , convergió
Los orbitales se reordenarán antes de escribirse en el archivo de punto de control. En cualquiera
case MOORDER Se debe dar un parámetro para especificar el mapa de reordenamiento. El valor predeterminado es
no reordenar los orbitales.

MOORDER = vector_entero
Esto especifica un vector de reordenamiento de orbitales moleculares. Solo se usará si
REORDENAR Está establecido. Este vector asigna cada orbital a su nuevo índice, p. Ej. MOORDER = (0
2 1) especifica que después de reordenar los orbitales 1 y 2 se intercambiarán. El rango de
este vector es el mismo que el número de OM. Los índices están en orden Pitzer
(ordenado por simetría, luego por energía dentro de cada bloque de simetría), base-0. CSCF
probablemente fallará si el MOORDER dado mezcla orbitales de diferentes irreps. Allí
no es predeterminado.

ALPHA = vector_real
If DE TIPO ABIERTO = SPECIAL, entonces este parámetro da los coeficientes de acoplamiento alfa.
El número de elementos en este vector es MM (MM + 1) / 2, donde MM es el número de
representaciones irreducibles que contienen orbitales moleculares ocupados individualmente. Allí
no es predeterminado.

BETA = vector_real
If DE TIPO ABIERTO = SPECIAL, entonces este parámetro da los coeficientes de acoplamiento beta.
El número de elementos en este vector es MM (MM + 1) / 2, donde MM es el número de
representaciones irreducibles que contienen orbitales moleculares ocupados individualmente. Allí
no es predeterminado.

GUESS = cadena
Esta opción determina el tipo de conjetura inicial en el vector propio que utilizará CSCF.
Las únicas opciones válidas en este momento son: (1) GUESS = NÚCLEO, lo que hace que use
vector propio hamiltoniano central para iniciar el cálculo; (2) GUESS = AUTO lo cual
da como resultado un intento de utilizar el vector MO en el archivo de punto de control, o recurre a
núcleo adivina si no hay un vector propio en ese archivo. El valor predeterminado si AUTO.

YO IMPRIMO = entero
Esta es una opción de impresión. El valor predeterminado es 0.

MO_SALIDA = booleano
Imprime los orbitales con simetría y ocupaciones al final de la
cálculo. El valor predeterminado es verdadero.

GIRAR = booleano
Los orbitales moleculares no se rotarán si esto es falso. La rotación solamente
afecta a los orbitales virtuales de los sistemas de capa abierta. Este parámetro debe ser verdadero
para gradientes correlacionados y debe ser falso para derivadas segundas y superiores.
El valor predeterminado es falso si WFNMás = SCF y cierto de lo contrario.

VERIFICAR_ROT = booleano
Compruebe la rotación del orbital molecular descrita anteriormente para asegurarse de que no haya columnas de
la matriz de vectores propios de SCF se intercambia mediante la rotación. No tiene ningún efecto si GIRAR =
false. El defecto es cierto.

CHECK_MO_ORTOGONALIDAD = booleano
Compruebe si los orbitales moleculares son ortonormales. Útil solo para depurar. los
el valor predeterminado es falso.

DIIS = booleano
Esto determina si se utilizará diis. El defecto es cierto.

DIISINICIO = entero
Esto da la primera iteración para la que se utilizará DIIS. El valor predeterminado es 0.

NDIIS = entero
Esto da el número de matrices de error a usar en el procedimiento diis. El valor por defecto
es 6 para carcasa cerrada, 4 para carcasa abierta y 3 para tcscf.

DIISDAMP = real
Esto da el factor de amortiguación para el procedimiento de diis. El valor predeterminado es 0.0 para
carcasa cerrada, 0.02 para carcasa abierta y 0.01 para tcscf.

INC = real
Esto se usa en tcscf para determinar la frecuencia con la que se recalculan los coeficientes ci.
Un número pequeño (~ 0.25) hará que se vuelvan a calcular casi cada scf
iteración. El valor predeterminado es 0.25.

DYN_ACC = booleano
Al realizar scf directo, esto especifica si los cortes de precisión integral dinámica
se utilizará. El valor predeterminado es verdadero (use cortes dinámicos). Las iteraciones iniciales son
realizado con integrales con precisión de seis dígitos. Después de que la densidad converja a
Se completan 10 ^ -5 o 30 iteraciones, se utiliza la precisión integral completa. Si scf
Los problemas de convergencia se experimentan deshabilitando los cortes dinámicos al establecer este
variable a falso podría ayudar.

ORTOG_SOLO = booleano
A veces, en CASSCF u otros esquemas no HF / KS para optimización orbital, es
útil para reortogonalizar los MO de otras geometrías para la geometría actual, por lo que
se pueden utilizar como una suposición inicial para los nuevos MO. Esto puede ser realizado por
ejecutando CSCF con ORTOG_SOLO = su verdadero. Después de que los orbitales estén ortogonalizados, el
El programa se cerrará sin realizar un cálculo SCF. Esta palabra clave será
se ignora si no hay orbitales anteriores en el archivo de punto de control. Predeterminado a su verdadero
if WFNMás = DETCAS.

30 mayo, 1991 cscf(1)

Utilice cscf en línea utilizando los servicios de onworks.net