NumForce module with direct SCF DFT methods

[Dempsey]

Hi,

I'm a new Turbomole user and I've been having trouble getting the aoforce module to run with COSMO. I switched over to NumForce and found that NumForce calculations only run with RI-DFT. This could be a problem for me as I want to use the PBE0 functional which doesn't work well with RI-DFT and thus I'm in this situation where I need direct SCF + COSMO + NumForce/aoforce but I cannot have it all.

Is it intentional that the NumForce module only works for RI-DFT i.e. "NumForce -ri > force.out" or should it work with direct SCF methods too?

[antti_karttunen]

Hi,

It is true that NumForce+COSMO runs only with ridft (not dscf). But you can use PBE0 functional with ridft just fine (unless you are using some rather old version of TM). Just turn on ri in define and ridft will run PBE0 calculations fine (control file will have the $rij keyword). Some more details in this white paper: https://arxiv.org/pdf/1610.07779.pdf

Best,
Antti

[Dempsey]

Hi Antti,

Thank you for your rapid response and apologies for mine being so delayed, I somehow didn't see that you replied despite checking many times. Part of my issue was resolved, I realised that I have to call NumForce -ri > force.out to ensure NumForce runs with the ridft module and thus works with COSMO calculations. I also have these calculations working with the PBE0 functional so thank you.

Should the aoforce module work with COSMO? It seems to fail for me saying "CPKS iterations did not converge". I have tried increasing $forceiterlimit and changing $forceconv. I find it quite useful to be able to view my vibrational modes using aoforce2g98 however I haven't been able to find a method like this if I perform the vibrational frequency analysis using NumForce. The NumForce module works well on my system but aoforce does not. Could you provide a reason for this?

Kind Regards,
Dempsey

[antti_karttunen]

Hi,

aoforce module does work with COSMO. But if you encounter convergence issues, you may need to tune some keywords as discussed here: https://forum.turbomole.org/index.php?topic=1105.0. I normally use NumForce with COSMO. You will still be able to use aoforce2g98 with COSMO. After running the NumForce job, go to the subdirectory numforce and execute
aoforce2g98 aoforce.out > freq.log
This will give you the "synthetic" Gaussian output as usual.

Best,
Antti

[Dempsey]

Hi,

Is it possible to visualize vibrations produced using NumForce? There's no numforce2g98 or equivalent is there?

Best,
Dempsey

[antti_karttunen]

Hi,

sorry, I was bit unclear. You can always find aoforce.out also after NumForce run, fromthe numforce directory. The final step of NumForce is always run with aoforce. So you can use aoforce2g98 also with NumForce.

Best,
Antti

[Dempsey]

Hi,

Ah, I see! Much appreciated!

Best,
Dempsey