The \(MnO_6\) system¶
The \(MnO_6\) compound is good example of the oxide crystal structure.
In the following we show how one can get converged molecular spinors (or molecular orbitals) of its ground state using single determinant SCF method.
Here the noAcAc abbreviation (“no active-active”) means value of OPENFAC=0, what is no active-active interation among open-shell spinors. The wAcAc (“with active-active”) means OPENFAC=1, that is full active-active interation (or full open-shell contribution to the Fock matrix).
For simplicity, we employ the two-component Hamiltonian, see the keyword .X2C.
The way to the convergence¶
With the aim to achieve the SCF convergence we investigated these three ways:
pam --inp=MnO6.x2c.bare_noacac.inp --mol=MnO6.crystal.mol --get "DFCOEF=DFCOEF.MnO6.noAcAc"
2. Using AcAc with molecular spinors (the DFCOEF file) from the previous SCF run the \(MnO_6\) system does not converge with the default DIIS accelerator. The new input file to download is MnO6.x2c.wacac.inp:
pam --inp=MnO6.x2c.wacac.inp --mol=MnO6.crystal.mol --put "DFCOEF.MnO6.noAcAc=DFCOEF"
3. Using the AcAc interaction with the parameter of OLEVEL=0.2, and, with the DFCOEF file from the point 1, we finally obtain the convergence. The total SCF energy of the system is -1617.5839185788545 a.u. (see the corresponding output file, MnO6.x2c.wacac_olvlshift_MnO6.crystal.out) The new DIRAC input file to download is MnO6.x2c.wacac.inp:
pam --inp=MnO6.x2c.wacac_olvlshift.inp --mol=MnO6.crystal.mol --put "DFCOEF.MnO6.noAcAc=DFCOEF" --get "DFCOEF=DFCOEF.MnO6.wAcAc"
This tutorial shows that it is important to try all possible means to achieve the SCF convergence if the default setting is not working.