Ошибка при анализе hessiane

[AlexKsen]

Уважаемые форумчане, здравствуйте. Пытаюсь посчитать гессиан, но все попытки приводят к одному результату - к ошибке. Пытаюсь в CCSD. Нужно посмотреть термохимию. Что я делаю не так? И может более низкого по качеству расчета будет достаточно? Может MP2 попробовать?

[amge]

Уважаемые форумчане, здравствуйте. Пытаюсь посчитать гессиан, но все попытки приводят к одному результату - к ошибке. Пытаюсь в CCSD.

$FORCE METHOD=ANALYTIC $END
Ишь че хотите!^*) Аналитического гессиана для CCSD, наверное, ни в одной программе нет. Аналитические градиенты кое-где есть, но не в гамессе. Так что только полностью численный, и будет изнуряюще долго.
___________
^*) Кстати, примерно то же самое сказано в конце приложенного Вами out-файла

[Гесс]

Ишь че хотите!^*) Аналитического гессиана для CCSD, наверное, ни в одной программе нет. Аналитические градиенты кое-где есть, но не в гамессе. Так что только полностью численный, и будет изнуряюще долго.

А вот это ТС не подойдет? http://slater.chemie.uni-mainz.de/cfour ... Properties

Пытаюсь в CCSD.

YOUR INPUT MISTAKENLY REQUESTED METHOD=ANALYTIC
ERROR IN $FORCE INPUT, PLEASE FIX AND TRY AGAIN.

ANALYTIC COMPUTATION OF THE HESSIAN IS IMPLEMENTED
FOR S,P,D BASIS SETS (INCLUDING ECP USAGE) FOR AB INITIO RHF, ROHF,
GVB (OPEN SHELL OR TCSCF) AND MCSCF (CISTEP=ALDET/ORMAS) WAVEFUNCTIONS.

SEMI-NUMERIC HESSIANS CAN BE COMPUTED WHEN ANALYTIC GRADIENTS ARE AVAILABLE:
IF F,G FUNCTIONS ARE IN THE BASIS SET,
IF USING GRID-BASED OR GRID-FREE DFT,
IF UHF, OR IF GVB WITH MORE THAN ONE PERFECT PAIR,
IF MCSCF USES SOME CISTEP OTHER THAN ALDET OR ORMAS,
IF MPLEVL=2 FOR RHF OR UHF WAVEFUNCTIONS,
IF EXCITED STATE CIS FOR RHF WAVEFUNCTIONS,
IF GROUND OR EXCITED STATE GUGA CI FOR RHF WAVEFUNCTIONS,
IF MODEL CORE POTENTIAL INTEGRALS ARE USED,
IF SEMI-EMPIRICAL MODELS AM1, PM3, OR MNDO ARE USED,
IF EFP, PCM, ONSAGER, OR COSMO SOLVENT MODELS ARE USED,
IF THERE IS AN APPLIED EXTERNAL ELECTRIC FIELD,
IF SCALAR RELATIVISTIC TRANSFORMATIONS ARE IN USE, OR
IF RUNNING SCF-MI COMPUTATIONS.

FULLY NUMERIC HESSIANS CAN BE COMPUTED WHEN ONLY THE ENERGY IS AVAILABLE:
WHEN MPLEVL=2 FOR SCFTYP=ROHF/OSPT=RMP OR SCFTYP=MCSCF,
WHEN CITYP IS SELECTED AND THE SCFTYP IS NOT RHF,
WHEN THE CITYP SELECTED IS NOT CIS OR GUGA, OR
WHEN ANY CCTYP IS SELECTED

[AlexKsen]

Спасибо. А какими методами Вы бы посоветовали попробовать посчитать.

[Гесс]

Гм. Ну согласно написанному выше - численным в CCSD, получисленным в MP2, либо в других программах...
Только не сравнивайте результаты из разных программ даже в одном методе.

[AlexKsen]

Спасибо.

[ Post made via Windows Smartphone ]