SCF fails spinlock with POB-DZVP-REV2

jquertin

Hi,

I'm currently running some molecular calculation as tests before going to periodic one of porphyrin with Co2+.
I'm using spin DFT with a spinloc2 1 -5 setting (with HSEsol).
With POB-DZVP, the calculation finishes successfully while with rev2 it fails after the 4th SCF cycle (sometimes the 3rd) with the good old "NaN at SCF".
I'm not sure what to try beside giving it more memory (which gives the same results) as those are my first calculation with metals.

I'm not sure if it's even a bug but since POB-DZVP works and not POB-DZVP-REV2, I figure that there's something that I cannot control at play.

Best,
Jean

PS: I also tried with FMIXING 60, spinlock instead of spinloc2, tighter tolinteg but got the same problem.

por_co-opt-pob_dzvp_rev2.d12 por_co-opt-pob_dzvp_rev2.out

Aleks

Hi,

Have you tried running a calculation where you specify the full basis set via the general route? This way one could rule out whether the internally stored basis set has some intrinsic error associated with it.

You also have DIIS turned on, which can cause issues with spin locking and cause instabilities. You can try switching it off (NODIIS) or alternatively remove spin lock, set up the spins at the beginning (ATOMSPIN) and monitor whether the DIIS is keeping the right solution.

Still might be a bug, but worth ruling out the simpler possibilities.

Hope it helps.

Cheers,
Aleks

jquertin

Aleks Thanks for the suggestions. I'll try them and keep the post updated.

jquertin

So, in short, everything that I tried (ATOMSPIN, NODIIS, explicitly give the basis set and combining all three) failed. Either with a similar error ("NaN" during SCF) or SIGSEGV ("forrtl: severe (174): SIGSEGV, segmentation fault occurred").

GiacomoAmbrogio

Hi jquertin,
As Aleks correctly pointed out, it’s better not to use DIIS (which is active by default) when using SPINLOCK.
A good way to combine the two is like this:

SPINLOCK  
1  
-2  
THREDIIS  
0.005

This way, DIIS is only activated when SPINLOCK is turned off.

However, I tried running your input, and it’s likely that this is actually a bug. Unfortunately, I don’t have a better solution at the moment other than changing the functional.