:orphan:


Nb atom
=======

DIRAC presently does not perform an open-shell Hartree--Fock in a strict sense,
rather an average-of-configuration calculation, which amounts to the
optimization of the average energy of a set of configurations (or determinants)
generated from the specification of a given number of open shells and their
electron occupations. As an example of such a calculation, showing in
particular tricks to help convergence, we consider the Nb atom.

The ground state configuration of Nb is
[Kr].4d\ :sup:`4`\ .5s\ :sup:`1`\  and so a reasonable input for an
average Hartree--Fock calculation would be::

  **DIRAC
  .WAVE FUNCTION
  .ANALYZE
  **HAMILTONIAN
  .X2C
  **INTEGRALS
  *READIN
  .UNCONTRACT
  **WAVE FUNCTION
  .SCF
  *SCF
  .CLOSED SHELL
   18 18
  .OPEN SHELL
   2
   4/10,0
   1/2,0
  **ANALYZE
  .MULPOP
  *MULPOP
  .VECPOP
   1..oo
   1..oo
  *END OF

Here we create separate open shells for the four *4d* electrons and the
single *5s* electron to reduce the number of states in the average.

We use the following mol file::

    DIRAC


    C   1
           41.    1
    Nb      0.0000000000        0.0000000000        0.0000000000
    LARGE BASIS dyall.v2z
    FINISH

Unfortunately this calculation does not converge::

  ERROR, SCF not converged, Dirac stops!

  --- SEVERE ERROR, PROGRAM WILL BE ABORTED ---

The reason is that the open-shell *4d* and *5s* orbitals mix.  We can see this
in the Mulliken population analysis::

  * Electronic eigenvalue no. 10: -0.2308050163326       (Occupation : f = 0.4000)  m_j=  1/2
  ==========================================================================================

  * Gross populations greater than 0.00010

  Gross     Total   |    L Ag Nb s      L Ag Nb dxx    L Ag Nb dyy    L Ag Nb dzz
  -----------------------------------------------------------------------------------
   alpha    1.0000  |      0.8647        -0.0540        -0.0540         0.2432
   beta     0.0000  |      0.0000         0.0000         0.0000         0.0000

  * Electronic eigenvalue no. 11: -0.2230586181555       (Occupation : f = 0.4000)  m_j= -3/2
  ==========================================================================================

  * Gross populations greater than 0.00010

  Gross     Total   |    L Ag Nb dxx    L Ag Nb dyy    L B1gNb dxy    L B2gNb dxz    L B3gNb dyz
  --------------------------------------------------------------------------------------------------
   alpha    0.2650  |      0.0663         0.0663         0.1325         0.0000         0.0000
   beta     0.7350  |      0.0000         0.0000         0.0000         0.3675         0.3675

  * Electronic eigenvalue no. 12: -0.2205588496606       (Occupation : f = 0.4000)  m_j=  1/2
  ==========================================================================================

  * Gross populations greater than 0.00010

  Gross     Total   |    L Ag Nb s      L Ag Nb dxx    L Ag Nb dyy    L Ag Nb dzz    L B2gNb dxz    L B3gNb dyz
  -----------------------------------------------------------------------------------------------------------------
   alpha    0.2178  |      0.0285         0.0481         0.0481         0.0931         0.0000         0.0000
   beta     0.7822  |      0.0000         0.0000         0.0000         0.0000         0.3911         0.3911

  * Electronic eigenvalue no. 13: -0.2186062869816       (Occupation : f = 0.4000)  m_j= -3/2
  ==========================================================================================

  * Gross populations greater than 0.00010

  Gross     Total   |    L Ag Nb dxx    L Ag Nb dyy    L B1gNb dxy    L B2gNb dxz    L B3gNb dyz
  --------------------------------------------------------------------------------------------------
   alpha    0.7350  |      0.1837         0.1837         0.3675         0.0000         0.0000
   beta     0.2650  |      0.0000         0.0000         0.0000         0.1325         0.1325

  * Electronic eigenvalue no. 14: -0.2178694489641       (Occupation : f = 0.4000)  m_j=  1/2
  ==========================================================================================

  * Gross populations greater than 0.00010

  Gross     Total   |    L Ag Nb s      L Ag Nb dxx    L Ag Nb dyy    L Ag Nb dzz    L B2gNb dxz    L B3gNb dyz
  -----------------------------------------------------------------------------------------------------------------
   alpha    0.7822  |      0.1068         0.1728         0.1728         0.3298         0.0000         0.0000
   beta     0.2178  |      0.0000         0.0000         0.0000         0.0000         0.1089         0.1089

  * Electronic eigenvalue no. 15: -0.2161776019252       (Occupation : f = 0.5000)  m_j=  5/2
  ==========================================================================================

  * Gross populations greater than 0.00010

  Gross     Total   |    L Ag Nb dxx    L Ag Nb dyy    L B1gNb dxy
  --------------------------------------------------------------------
   alpha    1.0000  |      0.2500         0.2500         0.5000
   beta     0.0000  |      0.0000         0.0000         0.0000

We see significant mixing of *4d* and *5s* orbitals.  This mixing is forbidden
from atomic symmetry. The problem is that we run the atom only in linear
supersymmetry and so the mixing is introduced by numerical noise due to the
energetic proximity of these orbitals.

An efficient way to handle such a case is to strip off the open-shell electrons
and simply calculate first Nb\ :sup:`5+`\  with the electron configuration of
Kr and to save the converged DFCOEF file::

  **DIRAC
  .WAVE FUNCTION
  .ANALYZE
  **HAMILTONIAN
  .X2C
  **INTEGRALS
  *READIN
  .UNCONTRACT
  **WAVE FUNCTION
  .SCF
  *SCF
  .CLOSED SHELL
   18 18
  #.OPEN SHELL
  # 2
  # 4/10,0
  # 1/2,0
  **ANALYZE
  .MULPOP
  *MULPOP
  .VECPOP
   1..oo
   1..oo
  *END OF

This calculation converges smoothly and
now the Mulliken population analysis tells us that
the gerade orbitals 10--14 are purely *d*, whereas
gerade orbital 15 is pure *s*. Please verify this.

We will now restart from the converged DFCOEF
and reinsert the open-shell electrons and impose
vector selection by overlap::

  **DIRAC
  .WAVE FUNCTION
  .ANALYZE
  **HAMILTONIAN
  .X2C
  **INTEGRALS
  *READIN
  .UNCONTRACT
  **WAVE FUNCTION
  .SCF
  *SCF
  .CLOSED SHELL
   18 18
  .OPEN SHELL
   2
   4/10,0
   1/2,0
  .OVLSEL
  .NODYNSEL
  **ANALYZE
  .MULPOP
  *MULPOP
  .VECPOP
   1..oo
   1..oo
  *END OF

This makes sure that the orbitals obtained in the calculation on
Nb\ :sup:`5+`\  stay in place. Again, we keep the DFCOEF file
for further restart.

The final step is again to restart from DFCOEF but this time
we relax the SCF by removing overlap selection::

  **DIRAC
  .WAVE FUNCTION
  .ANALYZE
  **HAMILTONIAN
  .X2C
  **INTEGRALS
  *READIN
  .UNCONTRACT
  **WAVE FUNCTION
  .SCF
  *SCF
  .CLOSED SHELL
   18 18
  .OPEN SHELL
   2
   4/10,0
   1/2,0
  #.OVLSEL
  #.NODYNSEL
  **ANALYZE
  .MULPOP
  *MULPOP
  .VECPOP
   1..oo
   1..oo
  *END OF

We obtain convergence and our open-shell orbitals are pure *d* and *s*. This is
what we wanted::

  * Electronic eigenvalue no. 10: -0.2557501903224       (Occupation : f = 0.4000)  m_j= -3/2
  ==========================================================================================

  * Gross populations greater than 0.00010

  Gross     Total   |    L Ag Nb dxx    L Ag Nb dyy    L B1gNb dxy    L B2gNb dxz    L B3gNb dyz
  --------------------------------------------------------------------------------------------------
   alpha    0.8000  |      0.2000         0.2000         0.4000         0.0000         0.0000
   beta     0.2000  |      0.0000         0.0000         0.0000         0.1000         0.1000

  * Electronic eigenvalue no. 11: -0.2557501876526       (Occupation : f = 0.4000)  m_j=  1/2
  ==========================================================================================

  * Gross populations greater than 0.00010

  Gross     Total   |    L Ag Nb dxx    L Ag Nb dyy    L Ag Nb dzz    L B2gNb dxz    L B3gNb dyz
  --------------------------------------------------------------------------------------------------
   alpha    0.4000  |      0.0667         0.0667         0.2667         0.0000         0.0000
   beta     0.6000  |      0.0000         0.0000         0.0000         0.3000         0.3000

  * Electronic eigenvalue no. 12: -0.2505426815019       (Occupation : f = 0.4000)  m_j=  5/2
  ==========================================================================================

  * Gross populations greater than 0.00010

  Gross     Total   |    L Ag Nb dxx    L Ag Nb dyy    L B1gNb dxy
  --------------------------------------------------------------------
   alpha    1.0000  |      0.2500         0.2500         0.5000
   beta     0.0000  |      0.0000         0.0000         0.0000

  * Electronic eigenvalue no. 13: -0.2505426814916       (Occupation : f = 0.4000)  m_j= -3/2
  ==========================================================================================

  * Gross populations greater than 0.00010

  Gross     Total   |    L Ag Nb dxx    L Ag Nb dyy    L B1gNb dxy    L B2gNb dxz    L B3gNb dyz
  --------------------------------------------------------------------------------------------------
   alpha    0.2000  |      0.0500         0.0500         0.1000         0.0000         0.0000
   beta     0.8000  |      0.0000         0.0000         0.0000         0.4000         0.4000

  * Electronic eigenvalue no. 14: -0.2505426791823       (Occupation : f = 0.4000)  m_j=  1/2
  ==========================================================================================

  * Gross populations greater than 0.00010

  Gross     Total   |    L Ag Nb dxx    L Ag Nb dyy    L Ag Nb dzz    L B2gNb dxz    L B3gNb dyz
  --------------------------------------------------------------------------------------------------
   alpha    0.6000  |      0.1000         0.1000         0.4000         0.0000         0.0000
   beta     0.4000  |      0.0000         0.0000         0.0000         0.2000         0.2000

  * Electronic eigenvalue no. 15: -0.2044678705701       (Occupation : f = 0.5000)  m_j=  1/2
  ==========================================================================================

  * Gross populations greater than 0.00010

  Gross     Total   |    L Ag Nb s
  --------------------------------------
   alpha    1.0000  |      1.0000
   beta     0.0000  |      0.0000