# *KRMCSCF¶

*2- and 4-component relativistic KR-MCSCF module*¶

The 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].

## Mandatory keywords under *KRMCSCF¶

### .GAS SHELLS¶

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 corep (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 corep. 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¶

Inactive orbitals per fermion corep. *default*: all orbitals are active. The example below for a molecule with inversion
center marks the lowest 4 *gerade* and 2 *ungerade* orbitals as inactive, i.e. they are always kept doubly occupied in all determinants
of the CI expansion:

```
.INACTIVE
4 2
```

All remaining orbitals above the inactive + active space are considered as secondary orbitals, i.e. they are kept empty in all determinants of the CI expansion. By default all orbital rotations between the inactive-active and active-secondary space are included.

## Optional keywords under *KRMCSCF¶

### .CI PROGRAM¶

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
```

### .THRESH¶

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
```

### .SYMMETRY¶

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 and use LUCIAREL as CI PROGRAM. 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
```

### .MAX MACRO¶

integer value giving the maximum number of MACRO iterations in the MCSCF optimization. *default*:

```
25
```

### .MAX MICRO¶

integer value giving the maximum number of MICRO iterations in the MCSCF optimization. *default*:

```
50
```

### .FROZEN¶

specifies how many of the inactive orbitals for each fermion corep should be *frozen*, that is, not change during the KRMCSCF optimization.
The frozen orbitals are not included in 2-electron integral transformation, which saves some CPU time.
*default*: include all inactive orbital rotations.
*example:* Freeze the first two inactive orbitals in fermion corep 1 (*gerade*) and the first inactive orbital in fermion corep 2 (*ungerade*):

```
.FROZEN
2 1
```

### .DELETE¶

delete occupied-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 corep.
*default*: include all active-secondary e-e rotations.
*example*: Delete all e-e rotations for secondary orbitals 20,21 and
22 in fermion corep 1 (*gerade*) and 24,25,26,27 and 28 in fermion corep 2 (*ungerade*):

```
.DELETE
20,21,22
24..28
```

### .SKIPEE¶

skip e-e rotations (rotations between electronic-electronic spinors) in the gradient and Hessian calculation. *default*: include e-e rotations.

### .WITHEP¶

include e-p rotations (rotations between electronic-positronic spinors) in the gradient and Hessian calculation. *default*: skip e-p rotations.

### .SKIPEP¶

skip e-p rotations (rotations between electronic-positronic spinors) in the gradient and Hessian calculation. *default*: skip e-p rotations.

### .MVOFAC¶

generate modified inactive and virtual orbitals based on block-diagonalization of FC + “MVOFAC”*FV.
The value does not change the MCSCF wave functions,
but the non-active orbitals for post-MCSCF calculations (e.g. *LUCITA or *KRCI) are modified.
*default*: MVOFAC = 1.0
*example*: Mimic Bauchschlicher’s MVOs for HF-CI for the MCSCF wave function:

```
.MVOFAC
0.0
```

### .PRINT¶

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
```

# *OPTIMI¶

## optional keywords under *OPTIMI¶

**NOTE: the following keywords must be placed under the input deck**

```
*OPTIMI
```

### .NOOCCN¶

compute natural orbital occupation numbers for the final optimized electronic state. *default*: do not compute natural orbital occupation numbers.

### .ANALYZ¶

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.

### .MAX CI¶

maximum number of initial CI iterations. *default*:

```
.MAX CI
5
```

### .MXCIVE¶

maximum size of Davidson subspace. *default*: 3 times the number of eigenstates (see .CIROOTS) to optimize on.
example:

```
.MXCIVE
3
```