ccenergy - Online na nuvem

PROGRAMA:

NOME

ccenergy - Programa de energia dupla e simples de cluster acoplado

DESCRIÇÃO

O programa ccenergia resolve o cluster acoplado simples e duplica as equações de amplitude.
Este código funcionará para restrição de spin, shell fechado Hartree-Fock (RHF), restrição de spin
Hartree-Fock (ROHF), Hartree-Fock (UHF) sem restrição de rotação e Brueckner
referências. O algoritmo faz uso de simetria espacial (D2h e seus subgrupos) e
técnicas de álgebra matricial para eficiência. Além disso, o programa pode ser usado para
simular métodos de correlação local.

REFERÊNCIAS

1. Uma introdução à teoria de agrupamento acoplado para químicos computacionais, TD Crawford
e HF Schaefer, Rev. Comp. Chem. 14 33-136 (2000).

2. Uma abordagem direta de decomposição de produto para exploração de simetria em muitos corpos
métodos. I. Cálculos de energia, JF Stanton, J. Gauss, JD Watts e RJ
Bartlett, J. Chem. Física 94 4334 (1991).

INPUT FORMATO

A entrada para este programa é lida do arquivo input.dat. As seguintes palavras-chave são válidas:

WFN = corda
Especifica o tipo de cálculo de cluster acoplado desejado. Entradas válidas em
presentes são CCSD, CCSD_T, BCCD e BCCD_T. Não há padrão.

REFERÊNCIA = corda
Especifica o tipo de orbitais usados ​​para a função de referência de determinante único.
As entradas avaliadas no momento são rhf, UHF e ROHF. Não há padrão.

MEMÓRIA = (reais MB)
A quantidade máxima de memória central a ser usada, em megabytes. Veja também o TIPO DE CACHE
e CACHELEV palavras-chave abaixo.

CONVERGÊNCIA = número inteiro
Convergência desejada nas amplitudes CC. A convergência é alcançada quando o RMS de
o erro na função de onda CC é menor que 10 ** (- n). O padrão é 7.

BRUECKNER_CONV = número inteiro
Convergência desejada para os orbitais nos cálculos de Brueckner-CC. Convergência é
alcançado quando o maior valor absoluto de uma única ampltiudes T1 é menor que
10 ** (- n). O padrão é 5.

MAXITER = número inteiro
Número máximo de iterações para resolver as equações de cluster acopladas. Padrões para
50.

FREEZE_CORE = corda
Especifica se os orbitais centrais devem ser congelados automaticamente na onda CC
função. Se esta opção for definida (de preferência na seção "padrão" de entrada),
então o cscf programa identificará os orbitais centrais. Códigos subsequentes começando
com transmissão irá congelar esses orbitais. Observe que essa palavra-chave geralmente torna o
especificação adicional de FROZEN_DOCC desnecessária, mas ainda pode ser necessária
em alguns casos especiais. Os valores permitidos são verdadeiro, falso, sim, não, 1, 0, grande e
pequeno.

FROZEN_DOCC = array_inteiro
O número de orbitais duplamente ocupados de menor energia em cada irredutível
representação da qual não haverá excitações. O pedido de algodão do
representações irredicíveis é usado. O padrão é o vetor zero. NB: Esta
palavra-chave é realmente lida pelo programa de transformação integral transmissão e os votos de
programa de classificação integral classificar. Veja também a palavra-chave FREEZE_CORE.

FROZEN_UOCC = vetor_inteiro
O número de orbitais desocupados de maior energia em cada representação irredutível
no qual não haverá excitações. O padrão é o vetor zero. NB: Esta
palavra-chave é realmente lida pelo programa de transformação integral transmissão e os votos de
programa de classificação integral classificar.

IMPRESSÃO = número inteiro
Determina o detalhamento da saída. Um valor de 0 (o padrão) especifica
impressão mínima e um valor de 2 fornecerá atualizações de saída conforme cada termo do CC
equações está completo. Nenhum outro valor foi usado ainda.

NUM_AMPS = número inteiro
Especifica o número das maiores amplitudes T1 e T2 para imprimir no final do CC
procedimento. O padrão é 10.

PRINT_MP2_AMPS = booleano
Indica se as amplitudes MP2 (que são estimativas iniciais para RHF e UHF
funções de referência) devem ser impressos. O padrão é falso.

REINICIAR = booleano
Permite que o programa use os vetores T1 e T2 antigos de cálculos anteriores como o
suposições iniciais para um novo cálculo. Isso é particularmente útil para geometria
otimizações. O padrão é VERDADEIRO. NB: A capacidade de reiniciar um cálculo é
também controlado pelas fases dos orbitais moleculares. Se o parâmetro FASE
do arquivo chkpt (cf. libchkpt) é definido como FALSE, então o usuário REINICIAR entrada
será ignorado. Este comportamento pode ser substituído, porém com o FORCE_RESTART
opção.

FORCE_RESTART = booleano
Força possível reinício das equações CC das amplitudes T1 e T2 anteriores.
Esta opção deve ser usada apenas por especialistas.

AO_BASIS = corda
Calcule as contribuições de quatro integrais de índice virtual, , do atômico
integrais de base orbital. Esta opção retarda o cálculo um pouco, mas muito
reduz os requisitos de espaço em disco. Os valores permitidos são NONE (padrão), DISK (use o
Integrais de base AO armazenados no disco), ou DIRETOS (use centavos para calcular o AO-
integrais dinâmicos em cada iteração CC). A opção DIRETA ainda é
experimental e deve ser usado apenas por especialistas.

CACHETYPE = corda
Seleciona o tipo de prioridade para manter o cache de memória automático usado pelo
Códigos DPD. (Consulte libdpd.html para obter mais detalhes.) Um valor de LOW (o padrão)
seleciona um esquema de "baixa prioridade" em que a exclusão de itens do cache é
com base em prioridades pré-programadas. Um valor de LRU seleciona um "menos usado recentemente"
esquema em que o item mais antigo do cache será o primeiro a ser excluído.

CACHELEV = número inteiro
Seleciona o nível de cache automático desejado no armazenamento de vários
amplitudes, integrais e intermediários no procedimento de agrupamento acoplado. Um valor
de 0 não retém nenhuma quantidade no cache, enquanto um nível de 6 tenta armazenar todos
quantidades em cache. Para cálculos particularmente grandes, um valor de 0 pode ajudar
com certos tipos de problemas de memória. O padrão é 2, o que significa que todos
quantites de quatro índices com até dois índices orbitais virtuais (por exemplo,
integrais) podem ser mantidos no cache.

DIIS = booleano
Seleciona o uso da técnica de inversão direta no subespaço iterativo (DIIS)
de Pulay para aceleração de convergência. O padrão é VERDADEIRO.

LOCAL = booleano
Simule os efeitos das técnicas de correlação local. O padrão é FALSE.

LOCAL_CUTOFF = reais
O valor de tolerância (sempre entre um e zero) para o Broughton-Pulay
verificação de integridade usada para construir domínios orbitais para cálculos de CC local. o
o padrão é 0.02. Ver J. Broughton e P. Pulay, J. Comp. Chem. 14, 736-740 (1993)
e C. Hampel e H.-J. Werner, J. Chem. Phys. 1046286-6297 (1996).

LOCAL_METHOD = corda
Seleciona o tipo de esquema CCSD local a ser simulado. As entradas válidas são WERNER
para o método desenvolvido por H.-J. Werner e colegas de trabalho ou AOBASE para o método
desenvolvido pela GE Scuseria e colaboradores (atualmente inoperante). O padrão é
WERNER.

LOCAL_WEAKP = corda
Seleciona o tratamento desejado de "pares fracos" no método local-CCSD. Um valor de
NEGLIGÊNCIA ignora pares fracos inteiramente. Um valor de NENHUM trata pares fracos no mesmo
maneira como pares fortes. Um valor de MP2 (o padrão) usa segunda ordem
teoria de perturbação para corrigir a energia CCSD local calculada com pares fracos
ignorado.

PRINT_PAIR_ENERGIES = booleano
Indica se as energias do par MP2 e CCSD devem ser impressas. Este é apenas
possível para referências RHF. O padrão é falso.

SPINADAPT_ENERGIES = booleano
Indica se as energias do par adaptado ao spin devem ser impressas. O padrão é falso.

18 Junho, 2002 ccenergia(1)

Use ccenergy online usando serviços onworks.net