The full CI program
This module is the determinant full CI program, as described in
P.J. Knowles and N.C. Handy, Chem. Phys. Lett. 111, 315 (1984),
P.J. Knowles and N.C. Handy, Comp. Phys. Commun. 54 (1989) 75.
Published work resulting from the use of this program should cite these references.
The program in normal use finds the lowest eigenvector of the complete CI hamiltonian matrix; more sophisticated use is possible, but not documented here. The program is interfaced to free standing versions such as supplied in the CPC program library by use of the DUMP
option.
The program is called with the command FCI
.
Defining the orbitals
ORBIT
,name.file;
name.file specifies the record from which orbitals are read. The default is the set of orbitals from the last SCF, MCSCF or CI calculation.
Occupied orbitals
OCC
,$n_1,n_2,\ldots,n_8$;
$n_i$ specifies numbers of occupied orbitals (including CORE
) in irreducible representation number $i$. If not given, the default is the complete basis set.
Frozen-core orbitals
CORE
,$n_1,n_2,\ldots,n_8$;
$n_i$ is the number of frozen-core orbitals in irrep number $i$. These orbitals are doubly occupied in all configurations, i.e., not correlated. If no CORE
card is given, the program uses the same core orbitals as the last CI calculation; if there was none, then the atomic inner shells are taken as core. To avoid this behaviour and correlate all electrons, specify
CORE
Defining the state symmetry
The number of electrons and the total symmetry of the wavefunction are specified on the WF
card:
WF
,elec,sym,spin
where
- elec: is the number of electrons
- sym: is the number of the irreducible representation
- spin: defines the spin symmetry, spin$=2S$ (singlet=0, doublet=1, triplet=2, etc.)
Density matrix
The 1-electron density matrix can be computed with
DM
,record.file
and then subsequently used to calculate one-electron properties using the PROP
program.
Excited states
STATE
,value;
Specifies the number of states to be sought (default 1).
Convergence threshold
THR
,value;
Convergence threshold (default $10^{-5}$). If changes in the normalised wavefunction become less than this value, then the iterative eigenvector solver terminates.
Printing options
PRINT
,code,value;
Print options. Generally, the value determines how much intermediate information is printed. value=-1 means no print (default for all codes). if value is omitted, it is taken as zero, which is usually appropriate. Specification of higher values will generate more output. The following codes are allowed:
ORBITAL
Print molecular orbitalsINTEGRAL
Print integralsTIMING
Print extra timing informationDIAGONAL
Print diagonal elements of HamiltonianHAMILTONIAN
Print much intermediate information
Interface to other programs
FCI,DUMP=
filename;
causes the FCI diagonalization to be bypassed, with input information and transformed integrals being written to a formatted file filename. The format is as described in Comp. Phys. Commun. 54 (1989) 75, and in an external program should be accessed using the FCIdump class or module available from https:%%//%%bitbucket.org/pjknowles/fcidump.
Example
- examples/h2o_fci.inp
***,FCI calculation for water in dz basis set basis=dz geometry={ o h,,0.0, 1.1,1.4 h,,0.0,-1.1,1.4} hf {fci wf,10,1;core,1} ! these options could be omitted ! FCI energy = -76.08088039 ---;