[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Convergence of SCF Loop in FEFF8 for TiO2



Dear Users,

 I'm sending again - previous msg had incorrect copy of Moreno's list:

 Sergio Moreno reported a lack of convergence of the SCF loop
for TiO2 with the following list of POTENTIALS:

POTENTIALS
*    ipot   Z  element            l_scmt  l_fms   stoichiometry
        0   22   Ti                2       2       0.01
        1   22   Ti                2       2       2
        2    8   O                 1       1       4

 I was able to reproduce the problem on our SGI with his input file.
It seems that with this input, the SCF loop is almost converging, but
somehow going into a limit cycle.

 Then I changed the POTENTIALS as follows, and it converged in 11 iterations.
 The -1 for l_scmt means the FEFF default, which is usually recommended.
For this system it is 2 ( i.e., s,p,d basis) both for Ti and O.
Looks like the potentials for the O in TiO2 want/need a little d character
to converge.

POTENTIALS
*    ipot   Z  element            l_scmt  l_fms   stoichiometry
        0   22   Ti                -1       2       0.01
        1   22   Ti                -1       2       2
        2    8   O                 -1       1       4

 I also changed RPATH to 0.1 - you don't need or should use
the path expansion with FMS calculations, so set RPATH less than
distance to first shell.

*RPATH     4.6
RPATH     0.1

  J. Rehr

----- End of forwarded message from John J. Rehr -----