частоты в Gamess

[Amora]

Здравствуйте. Я столкнулась со следующей проблемой. У меня есть комплекс метан-азот. Нужно с помощью Gamess посчитать частоты колебательные. Я делаю оптимизацию, и одна частота всегда мнимая. (Оптимизировала в Gaussiane и считала частоты- все ок). Подскажите как нужно оптимизировать правильно геометрию?
я делала так:


$CONTRL SCFTYP=RHF RUNTYP=OPTIMIZE COORD=UNIQUE MPLEVL=2 UNITS=ANGS $END
$SYSTEM TIMLIM=600 MEMDDI=100 MWORDS=1000 $END
$BASIS GBASIS=ACCT $END
$DATA
CH4...N2  opt
C1
C 6.0 0.0000000000 0.0000000000 0.0000000000
H 1.0 -1.0860100000 0.0000000000 0.0000000000
H 1.0 0.3631400000 1.0236800000 0.0000000000
H 1.0 0.3616700000 -0.5123400000 -0.8874900000
H 1.0 0.3616700000 -0.5123400000 0.8874900000
N 7.0 3.6197000000 0.0003800000 -0.5571200000
N 7.0 3.6197000000 0.0003800000 0.5571200000
$END

[VTur]

Добавьте
$STATPT OPTTOL =1.0d-5 HESS =CALC IHREP=1 HSSEND=.t. $END

[Reizen]

IHREP=1

это жестоко

Мне как раз то и нравится в GAMESS-е, что можно указать, допустим, IHREP=4 или 7, и будет тратиться намного меньше времени на расчёт гессиана (тем более, что в данном случае он будет считаться "трёхточечно"!).

Amora, главная проблема здесь, как кажется, в том, что из-за пологости ППЭ оптимизация всё время останавливается не в самом минимуме, а на каком-то расстоянии от него.

Поэтому главное лекарство в предложенной VTur-ом строке (имхо) - это OPTTOL =1.0d-5. Можете даже поставить OPTTOL =1.0d-7. А HESS=CALC и IHREP=1 оставьте на крайний случай, если первое средство не поможет.

[Amora]

!!!!!!!!!!!!!!!