Using escf I am trying to calculate dynamic polarisabilities, setting frequencies with $scfinstab dynpol.

I find if I try to do "too many" frequencies at once, then the calculation crashes with, for instance (for chloride):

..........................................................................................

Iteration IRREP Converged Max. Euclidean

roots residual norm

0 2 e - integrals written in 0 blocks requiring 0 k-byte

1 a 17 1.007987072749249D-02

Warning! RPA metric is ill-conditioned: cond = 1.618844632179714E+016

2 a 20 1.639881319565660D-03

MODTRACE: no modules on stack

singular metric in <ritz2>

escf ended abnormally

..........................................................................................

If I reduce the calculation down to say 1 frequency, then the calculation succeeds. In some cases I can successfully push through 15 frequencies, in other cases (for the same system) even 10 frequencies crashes. There's no obvious pattern explaining the crash.

I'm not sure that it's related to the amount of memory available. Increasing $rpacor or $maxcor makes no difference, and the turbomole is not even using that (uses only 35MB of the 500MB I gave it).

Is there any sensible fix or input parameter for handling many dynpol frequencies at once or are single frequency calculations the only way to work around the problem?

Thanks,

Drew