2- and 4-component relativistic KR-MCSCF module
Original implementation written by Joern Thyssen, Timo Fleig and Hans Joergen Aagaard Jensen, parallelization and continuous improvement by Stefan Knecht and Hans Joergen Aagaard Jensen. Linear symmetry adaptation by Stefan Knecht and Hans Joergen Aagaard Jensen.
For more details on the theory and actual implementation see [Jensen1996], [Thyssen2004] , [Thyssen2008], [Knecht2009].
specifies the CI module behind KRMCSCF, choices are LUCIAREL and GASCIP. Note: GASCIP does not work with linear symmetry and for large-scale MCSCF (> \(1.0 \times 10^7\) determinants) only LUCIAREL can be used efficiently. default:
.CI PROGRAM
GASCIP
Specification of CI calculation with electron distribution in orbital (GAS) spaces. The first line contains the number of GA spaces to be used (1-7), followed by one line per GAS with a separation by a “/” of the min/max number of electrons in each GAS and the number of orbitals per fermion correp (either one (no inversion symmetry) or two (inversion symmetry: gerade ungerade) entries per line).
The first entry before the “/” gives the minimal number of accumulated (!) electrons after consideration of this GAS, the second the corresponding maximum number, separated by blanks. The minimum and maximum accumulated occupations allow for a very flexible parameterization of the wave function. All determinants fulfilling the occupation constraints will be constructed. The second entry after the “/” gives the number of orbitals per fermion correp. See also the open-shell input in SCF which is similar to the syntax used here. The design of a GAS scheme is non-trivial and should be motivated by the electronic structure of the system (e.g. inner core, outer core, valence, virtual space). Sometimes it is useful to subdivide the valence space, for scientific reasons, or/and the virtual space, for technical reasons (save core memory). See reference [Fleig2006a], pp. 27 for more details. example for a CASSCF run with 10 electrons in 10 orbitals also known as CAS(10,10):
.GAS SHELLS
1
10 10 / 5 5
if all orbitals (fullMCSCF, feasible only for systems with 2-3 active electrons) should be included, use:
.GAS SHELLS
1
2 2 / all
Inactive orbitals per fermion correp. default: all orbitals active. example for 4 gerade and 2 ungerade orbitals frozen in a molecule with inversion center:
.INACTIVE
4 2
convergence threshold in the MCSCF gradient given as double precision value. The present default value is quite tight and might be adapted if one is only interested in MCSCF energies, e.g. \(1.d-03\). default:
0.5d-05
integer value giving the (boson/fermion) symmetry of the state to optimize on. default:
1
If you run the calculation in linear symmetry you have to specify 2 \(\times \Omega\) value of the state to optimize on. The doubling stems from the fact that we want to avoid non-integer input, e.g. in case of an odd number of electrons we might have \(\Omega\) = 1/2, 3/2, 5/2, etc. values and the corresponding input for the \(\Omega\) = 1/2 would then read as
1
If we have a system with an even number of electrons and inversion symmetry the input for an \(\Omega\) = 2g state would read as
4g
integer value giving the maximum number of MACRO iterations in the MCSCF optimization. default:
25
integer value giving the maximum number of MICRO iterations in the MCSCF optimization. default:
50
delete active-secondary e-e rotations (rotations between electronic-electronic spinors) in the gradient and Hessian calculation specified by an orbital string of virtual (secondary) orbitals for each fermion correp. default: include all active-secondary e-e rotations. example: (deleting all e-e rotations for secondary orbitals 20,21 and 22 in fermion correp 1 (gerade) and 24,25,26,27 and 28 in fermion correp 2 (ungerade)):
.DELETE
20,21,22
24..28
skip e-e rotations (rotations between electronic-electronic spinors) in the gradient and Hessian calculation. default: include e-e rotations.
skip e-p rotations (rotations between electronic-positronic spinors) in the gradient and Hessian calculation. default: include e-p rotations.
raise default print level to the given integer value. Please use with care as you may get millions of output lines if you choose a too high value. default:
.PRINT
0
NOTE: the following keywords must be placed under the input deck
*OPTIMI
compute natural orbital occupation numbers for the final optimized electronic state. default: do not compute natural orbital occupation numbers.
analyze the final CI wave function printing the coefficients for each determinant above a given threshold \(10^{-2}\). default: do not analyze the final CI wave function.
maximum size of Davidson subspace. default: 3 times the number of eigenstates (see CIROOTS) to optimize on. example:
.MXCIVE
3