Page **1** of **1**

### Problems with vibrational analysis with DFT and PCM

Posted: **08 Oct 2019, 14:51**

by **sauer**

Hi everybody,

we are trying to do run a geometry optimization and vibrational analysis of a molecule in PCM at the DFT level, but have problems. Apparently one has to use numerical derivatives for the Hessian. But asking for this we get a

"MEMCHK ERROR, not a valid memget id in work(kalloc-1)". See the attached output file, which also includes in the input files.

But maybe we are using the wrong keywords to get this.

Any help would be appreciated.

Stephan

### Re: Problems with vibrational analysis with DFT and PCM

Posted: **10 Oct 2019, 13:43**

by **magnus**

Looks like a bug which is possibly connected with switches between symmetries during numerical differentiation. Unfortunately, turning off symmetry doesn't work because of another bug.

In any case, I would run this in two steps. First do a geometry optimization and then perform vibrational analysis using the optimized geometry (and numerical Hessian).

### Re: Problems with vibrational analysis with DFT and PCM

Posted: **10 Oct 2019, 19:45**

by **sauer**

Hi Magnus,

thanks for looking at this. Unfortunately, running it in two steps does not help either. It runs in the same bugs. It is a bit strange, because it must mean, that nobody does DFT/PCM geometry optimizations with Dalton.

Stephan

### Re: Problems with vibrational analysis with DFT and PCM

Posted: **10 Oct 2019, 19:57**

by **magnus**

Sorry I didn't mean to imply that it works if you run out in two steps. That was to avoid the calculation of a numerical Hessian in the geometry optimization step where it is not needed.

The geometry optimization works if you do not ask for the numerical Hessian. However, it is done with a fixed cavity according to Kenneth and I'm not exactly sure what the implication of that is.

The error you get is because of a bug in the calculation of the numerical Hessian when you include PCM.

### Re: Problems with vibrational analysis with DFT and PCM

Posted: **10 Oct 2019, 20:00**

by **magnus**

If the bug is related to switching symmetry as I suspect then perhaps it can be avoided if symmetry detection is turned off. Unfortunately specifying "NoSymmetry" in the molecule file is not enough.

### Re: Problems with vibrational analysis with DFT and PCM

Posted: **10 Oct 2019, 20:29**

by **sauer**

And you do not happen to know where I could turn it off hard-code in the program?

### Re: Problems with vibrational analysis with DFT and PCM

Posted: **10 Oct 2019, 23:21**

by **xiongyan21**

GAMESS can easily incorporate PCM into DFT and TDDFT calculations, and I just have successfully tried a UHF spin-flip TDDFT geometry optimization of a pesticide molecule for its triplet state in water, obtaining 5 states preset.

I will look into the s-squared.

Perhaps I also can use the geometry obtained for adiabatically using pcm for its phosphorescence calculation, vertical or adiabatic, in water with Dalton.

### Re: Problems with vibrational analysis with DFT and PCM

Posted: **11 Oct 2019, 08:21**

by **magnus**

sauer wrote: ↑10 Oct 2019, 20:29

And you do not happen to know where I could turn it off hard-code in the program?

Not very specific but somewhere in DALTON/abacus/abander.F.

### Re: Problems with vibrational analysis with DFT and PCM

Posted: **11 Oct 2019, 09:30**

by **sauer**

Ok, thanks.