New HF code convergence problem

[thierry....@gmail.com]

Dear all

we hit a serious convergence problem with the HF code in Molpro2025.

We tried to redo some old calculations which converged without any problem using the old Molpro2009 HF code. The systems studied are carbon nanotubes and nanocones.

Unfortunately, the only way to make these converge with the new HF is to remove random (it depends on the system) pairs of electrons to obtain starting orbitals.

We tried several option of the new code (SO-SCI, level shifts, algorithm,...) nothing helped.

I join the input and the corresponding Molpro2009 and Molpro2025 outputs.

Thank you for any hint you may bring.

Best wishes,

 Thierry

[tibo...@gmail.com]

Dear  Thierry ,

We have also been seeing strangely worse HF/KS convergence ever since the default HF program was replaced in version 2021.1.

But for us SO-SCI, or SO usually works, albeit slower than the old code.

My recommendation is to try using the old SCF code, accessible by the old option, eg. {hf,old}.

Best regards,

Tibor Győri