Пытаюсь оптимизировать геометрию одной стильбеноподобной молекулы из 26 атомов в US GAMESS. Нужно для публикации.
Планировал сначала полуэмпирикой (AM1, PM3), а потом DFT оптимизировать, но ни в одном из расчетов, что я пробовал, не получилось достигнуть равновесной геометрии (***** FAILURE TO LOCATE STATIONARY POINT, TOO MANY STEPS TAKEN *****).
Пробовал в разных вариациях: увеличить NSTEP, MAXIT, уменьшить TRMIN, пробовал Ньютона-Рафсона, включал NOSYMM (хотя тут и так С1), дописывал UNITS=ANGS и COORD=UNIQUE, хотя это и так значения по умолчанию. Короче, много чего попробовал.
Предположил, что, может быть, полуэмпирика чересчур старается выплощить молекулу и не достигает равновесной геометрии из-за невыгодных взаимодействий этиленовых водородов с арильными орто-водородами, и решил попробовать DFT в нехитром базисе.
Инпут был следующий:
Код: Выделить всё
! File created by the GAMESS Input Deck Generator Plugin for Avogadro
$BASIS GBASIS=N31 NGAUSS=6 DIFFSP=.TRUE. $END
$CONTRL SCFTYP=RHF RUNTYP=OPTIMIZE DFTTYP=B3LYP MAXIT=50 $END
$SCF DIRSCF=.TRUE. $END
$SYSTEM MWORDS=100 MEMDDI=100 $END
$DATA
С1
// декартовы координаты атомов; все значения меньше 10 по абсолютной величине, ароматические кольца (транс-конфигурация) молекулы немного выведены из плоскости этиленового фрагмента при построении //
$END
Код: Выделить всё
MAXIMUM GRADIENT = 0.0009005 RMS GRADIENT = 0.0003724
NSERCH: 50 E= -743.4487139282 GRAD. MAX= 0.0009005 R.M.S.= 0.0003724
HESSIAN UPDATED USING THE BFGS FORMULA
ACTUAL ENERGY CHANGE WAS -0.0001455143
PREDICTED ENERGY CHANGE WAS -0.0001662654 RATIO= 0.875
MIN SEARCH, CORRECT HESSIAN, TRYING PURE NR STEP
NR STEP HAS LENGTH = 74.274857
TRIM/QA LAMBDA FOR NON-TS MODES = -0.01210643
TRIM/QA STEP HAS LENGTH = 0.070711
RADIUS OF STEP TAKEN= 0.07071
CURRENT TRUST RADIUS= 0.07071
Подскажите, пожалуйста, как мне выполнить расчет корректно.