Поиск ПС реакции циклизации

[sorcal]

Добрый день, уважаемые коллеги.
Недавно столкнулся со следующей проблемой.

Есть реакция галогенциклизации. Интермедиатами могут быть бромониевый или тиираниевые ионы.

pic1.gif

Методами PM3 и MNDO ПС реакции перехода тиираниевого иона в тиазиновый цикл локализуется.

pic2.gif

А вот при попытке отыскания ПС этой же реакции в рамках HF, DFT (с разными функционалами и базисами) и MP2,
процедура sadpoint в firefly сваливает структуру в минимум, соответствующий циклическому продукту.
Есть у кого-нибудь какие-нибудь идеи как разрешить эту проблему и оценить E_A ab initio методами?

Использую такую шапку:
$contrl aimpac=.t. pltorb=.t. scftyp=rhf runtyp=sadpoint
DFTTYP=B3LYP COORD=cart icharg=1 mult=1 MAXIT=999 $end
$statpt ststep=0.05 opttol=0.0001 nprt=-2 nstep=999 $end
$SCF diis=.true. damp=.t. soscf=.false. dirscf=.t. $END
$system timlim=50000 mwords=50 $end
$basis gbasis=n31 ngauss=6 npfunc=1 ndfunc=1 diffsp=.t.
diffs=.t. $end

[VTur]

ПС этой же реакции в рамках HF, DFT (с разными функционалами и базисами) и MP2

и полуэмпирические очень часто (обычно) не совпадают. Как раз сейчас мы с одним уважаемым участником форума бьёмся над такой задачей.
Метод один - в полуэмпирике получаете строение ПС и начинаете снова искать ПС уже HF, DFT или MP2. Лучше начинать с очень маленьких базисов, тогда строение ПС будет ближе к полуэмпирике. А потом постепенно наращиваете базис (чаще именно так приходится делать). Если резко нарастить, можно пролететь мимо мнимой частоты.
Хотя, если представлете себе нужную структуру, то можно сразу и без полуэмпирики и в больших базисах.

[sorcal]

Метод один - в полуэмпирике получаете строение ПС и начинаете снова искать ПС уже HF, DFT или MP2.

Да, я обычно так и поступаю.
А вот постепенно наращивать базис не пробовал. Попробую. Спасибо.