detcas - Online sa Cloud

PROGRAMA:

NAME

detcas - Determinant Configuration Interaction / Complete Active
Space Self-Consistent-Field Program

DESCRIPTION

Ang programa detcas gumaganap ng complete-active-space self-consistent-field (CASSCF)
pag-optimize ng mga molecular orbital sa pamamagitan ng dalawang hakbang na pamamaraan kasabay ng
determinant configuration program sa pakikipag-ugnayan DETCI. Ang programa ay medyo simple at
kasalukuyang gumagamit ng Newton-Raphson approach para i-update ang mga orbital, na gumagamit ng simple
tinatayang orbital na Hessian. Ang convergence ay pinabilis gamit ang direktang pagbabaligtad ng Pulay ng
ang iterative subspace (DIIS) na pamamaraan. Ang code ay isinulat upang payagan ang mas pangkalahatan
mga wavefunction na hindi kinakailangang nagtatampok ng buong CI na paggamot sa aktibong espasyo. Sa
partikular, ang anumang restricted active space (RAS) CI wavefunction na sinusuportahan ng DETCI ay maaaring
ginamit, na nagpapahintulot sa mga wavefunction ng MCSCF ng uri ng RASSCF.

Mga sanggunian

Tinatayang Orbital Hessian:

1. G. Chaban, MW Schmidt, at MS Gordon, Theor. Chim. Acta 97, 88-95 (1997).

Restricted Active Space CI:

1. Determinant Based Configuration Interaction Algorithm para sa Kumpleto at Pinaghihigpitan
Configuration Interaction Spaces, J. Olsen, BO Roos, P. Jorgensen, at HJ Aa.
Jensen, J. Chem. Phys. 89, 2185 (1988).

Restricted Active Space SCF:

1. P.-A. Malmqvist, A. Rendell, at BO Roos, J. Phys. Chem. 94, 5477 (1990).

Programa ng DETCI:

1. CD Sherrill, Computational Algorithms para sa Large-Scale Full at Multi-Reference
Configuration Interaction Wavefunctions, PhD thesis, University of Georgia, Athens,
GA, 1996.

2. CD Sherrill at HF ​​Schaefer, The Configuration Interaction Method: Advances
sa Highly Correlated Approaches, Adv. Quantum Chem. 34, 143-269 (1999).

MGA FILE KAILANGAN

input.dat - Input na file
file78 - MO one-electron integrals (fzc operator)
file72 - MO dalawang-electron integral
file73 - MO one-particle density matrix
file74 - MO dalawang-particle density matrix
file75 - MO Lagrangian

MGA FILE Nai-update

output.dat - Output file
file14.dat - Record ng mga energies at orbital gradients

INPUT FORMAT

Ang mga sumusunod na argumento sa command-line ay magagamit:

-tahimik Nagbibigay ito ng parehong resulta bilang I-PRINT=0.

-o fname
Nagbibigay ng filename para sa output file. Default sa output.dat.

Ang karagdagang input para sa program na ito ay binabasa mula sa file input.dat. Ang mas karaniwang ginagamit
ang mga keyword ay:

WFN = pisi
Ang mga sinusuportahang uri ng function ng wave ay CASSCF at RASSCF.

CONVERGENSIYA = kabuuan
Convergence ninanais sa orbital gradient. Nakakamit ang convergence kapag ang RMS
ng error sa orbital gradient ay mas mababa sa 10**(-n). Ang default ay 4 para sa
mga kalkulasyon ng enerhiya at 7 para sa mga gradient.

DOCC = integer_array
Ibinibigay ng vector na ito ang bilang ng mga dobleng orbital na inookupahan sa bawat irrep. meron
walang default.

SOCC = integer_array
Ang vector na ito ay nagbibigay ng bilang ng mga single occupied orbital sa bawat irrep. meron
walang default.

ENERGY_CONVERGENCE = kabuuan
Convergence ninanais sa kabuuang MCSCF enerhiya. Ang default ay 7.

FROZEN_DOCC = integer_array
Ang bilang ng pinakamababang enerhiya na dobleng inookupahan ng mga orbital sa bawat hindi mababawasan
representasyon na literal na mapi-freeze (hindi na-update sa MCSCF). Ang
Ang pag-order ng cotton ng mga hindi mapag-aalinlanganang representasyon ay ginagamit. Ang default ay ang
zero vector.

FROZEN_UOCC = integer_vector
Ang bilang ng pinakamataas na enerhiya na walang tao na mga orbital sa bawat hindi mababawasan na representasyon
na literal na mapi-freeze (hindi na-update sa MCSCF). Ang default ay ang zero
vector

RESTRICTED_DOCC = integer_array
Ang bilang ng pinakamababang enerhiya na dobleng inookupahan ng mga orbital sa bawat hindi mababawasan
representasyon na i-optimize ngunit pananatiling dobleng inookupahan sa MCSCF.
Ang mga orbital na ito ay kasunod ng FROZEN_DOCC mga orbital. Ang default ay ang zero
vector

RESTRICTED_UOCC = integer_array
Ang bilang ng pinakamataas na enerhiya na walang tao na mga orbital sa bawat hindi mababawasan na representasyon
na i-optimize ngunit pananatiling walang tao sa MCSCF. Dumating ang mga orbital na ito
bago ang FROZEN_UOCC mga orbital. Ang default ay ang zero vector.

NCASITER = kabuuan
Pinakamataas na bilang ng mga pag-ulit upang ma-optimize ang mga orbital. Ang pagpipiliang ito dapat be
tinukoy in ang DEFAULT seksyon of input, dahil kailangan itong makita ng
control program PSI. Default sa 1.

PRINT = kabuuan
Tinutukoy ng opsyong ito ang verbosity ng output. Ang isang halaga ng 1 o 2 ay tumutukoy
minimal na pag-print, ang halaga na 3 ay tumutukoy sa verbose printing. Ang mga halaga ng 4 o 5 ay
ginagamit para sa pag-debug. Huwag gumamit ng level 5 maliban kung ang test case ay napakaliit (hal
STO H2O CISD).

Ang hindi gaanong ginagamit na mga keyword ay:

DIIS_FREQ = kabuuan
Ang pamamaraan ng extrapolation ng DIIS ay susubukan bawat n mga pag-ulit. Ang default
ay 1.

DIIS_MAX_VECS = kabuuan
Ang maximum na bilang ng mga subspace vectors para sa DIIS procedure. Pagkatapos ng numerong ito
ng mga vector ay naabot, ang mga mas lumang vector ay ihuhulog mula sa subspace bilang
kailangan. Ang default ay 8.

DIIS_MIN_VECS = kabuuan
Ang pinakamababang bilang ng mga subspace vectors bago ang isang interpolation ng DIIS ay maaaring
gumanap. Ang default ay 2. Hindi makatuwirang magkaroon ng mga value na mas mababa sa 2.

DIIS_START = kabuuan
Ang numero ng pag-ulit kapag Pulay's Direct Inversion of the Iterative Subspace (DIIS)
pamamaraan para sa acceleration ng convergence ay dapat na naka-on. Ang default ay 3.
Bago ang pag-ulit na ito, hindi idinaragdag ang mga vector sa subspace ng DIIS. Ang una
Ang hakbang ng DIIS ay hindi gagawin hanggang sa DIIS_MIN_VECS Ang mga vector ay nasa subspace ng DIIS.

LAG_FILE = kabuuan
File (unit number) para sa pagbabasa ng lagrangian matrix. Ang default na halaga ay
kasalukuyang 75.

OPDM_FILE = kabuuan
File (unit number) para sa pagbabasa ng one-particle density matrix. Ang default na halaga
ay kasalukuyang 73.

SCALE_GRAD = boolean
Sinasabi kung isusukat ang orbital gradient sa pamamagitan ng tinatayang (diagonal) na orbital
Hessian. Ang default ay TRUE.

SCALE_STEP = tunay
Scale factor para sa orbital rotation step. Ang default ay 1.0.

TPDM_FILE = kabuuan
File (unit number) para sa pagbabasa ng two-particle density matrix. Ang default na halaga
ay kasalukuyang 74.

LEVEL_SHIFT = boolean
Sinasabi kung papayagan ang paglipat ng antas ng hessian matrix. Ito ay maaaring gamitin
upang matiyak na ang hessian ay positibong tiyak para sa mga simulang pag-ulit. Kung
ang hessian ay hindi positive definite malapit sa convergence, kung gayon ito ay maaaring isang senyales ng
mga hindi katatagan ng numero sa MCSCF. Ang default ay TRUE.

SHIFT = tunay
Sinasabi kung gaano kalaki ang dapat ilipat sa mga diagonal na elemento ng hessian para sa antas
paglilipat. Ang default ay 0.01.

DETERM_MIN = tunay
Minimum na pinahihintulutang halaga para sa determinant ng hessian matrix kung level shifting
ay sa. Hindi ito ginagamit kung LEVEL_SHIFT=FALSE. Ang default ay 0.00001.

STEP_MAX = tunay
Ito ang maximum na pinapayagang solong pag-ikot ng orbital. Ang default ay 0.30.

FORCE_STEP = boolean
Binibigyang-daan nito ang user na lampasan ang kinakalkula na hakbang at puwersahin ang isang hakbang sa a
partikular na direksyon. Ito ay maaaring maging kapaki-pakinabang kung sinusubukang pilitin ang pagkalkula
mula sa isang saddle point. Ang default ay FALSE.

FORCE_PAIR = kabuuan
Ito ang index para sa independiyenteng pares na dapat paikutin kung
FORCE_STEP=TRUE. Ito ay hindi pinapansin kung hindi man. Walang default.

FORCE_VALUE = tunay
Ito ang orbital rotation value para sa pag-ikot ng mga orbital na tinukoy ni
FORCE_PAIR. Ito ay ginagamit lamang kung ang FORCE_STEP ay nakatakda sa TRUE. Ang default na halaga ay
0.0.

CHECK_HESSIAN = boolean
Kalkulahin ang eigenvalues ​​ng orbital hessian. Ito ay mabuti para sa pagsuri
mga solusyon sa saddle-point. Dapat gamitin nang walang level shifting. Ang default ay
MALI.

EIGEN_VECTORS = boolean
Gusto mo bang kalkulahin ang eigenvectors ng orbital hessian pati na rin. Ito ay
hindi pinansin kung CHECK_HESSIAN=FALSE. Ang default ay FALSE.

8 Mayo, 1998 detcas(1)

Gumamit ng mga detcas online gamit ang mga serbisyo ng onworks.net