Dalton

[kbob]

Скачал сей фантастический пакет. Аналитические градиенты для CC-методов таки есть.
Не могу понят как расчитывается Hessian в CCSD.
Пытаюсь задать следующий инпут файл
[cut]**DALTON INPUT
.WALK
**INTEGRAL
.DIPLEN
.DEROVL
.DERHAM
**WAVE FUNCTION
.CC
*CC INPUT
.CCSD
**PROPERTIES
.NUMHES
*VIBANA
.HESPUN
*END OF INPUT[/cut]

получаю ошибку

Код: Выделить всё

At line 2906 of file cc_fop.F (unit = 12, file = 'CC_FCKH')
Fortran runtime error: End of file
Error in /opt/Dalton2011_release/bin/dalton.x, exit code 2

при этом вижу, что файл CC_FCKH есть, но он пустой.
Уважаемые коллеги подскажите что неправильно в моем инпуте?

[endlesslake]

вроде вот такой инпут считает численный hessian из аналитических градиентов (для уже оптимизированной геометрии)

[cut]**DALTON INPUT
.NMDDRV
**NMDDRV
.DORDR
1 1
.VIBANA
**INTEGRAL
.DIPLEN
.DEROVL
.DERHAM
**WAVE FUNCTIONS
.TITLE
test case for CCSD seminumerical hessian
.CC
*SCF INPUT
.THRESH
1.0d-12
*CC INPUT
.CCSD
.THRENR
1.0d-10
.THRLEQ
1.0d-10
.PRINT
2
**END OF DALTON INPUT[/cut]

[kbob]

Спасибо, директива DORDR действительно полезна.

Код: Выделить всё

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0  0x7F64968780F7
#1  0x7F64968786D4
#2  0x7F64961DA49F
#3  0x87BA6D in numdr1_
#4  0x881862 in numdrv_
#5  0x47BCCD in exedrv_
#6  0x404CA5 in MAIN__ at dalton.F:0
/opt/Dalton2011_release/bin/dalton: строка 294: 20101 Ошибка сегментирования                   (подготовлен дамп ядра) $DALTON

появился более осмысленный результат!

[kbob]

Опечатка - запятая между единицами обязательна, иначе программа безуспешно пытается найти 11-ю производную энергии по координатам!

Код: Выделить всё

**NMDDRV
.DORDR
1, 1

[endlesslake]

странно, в тестовых задачах - без запятой(с пробелами), да и у меня нормально тот инпут работал.

[kbob]

У меня заработала с запятой, причем вполне адекватный результат расчета.
Еще один вопрос.
Во всех DFT расчетах выходит сообщение об ошибке.

Number of electrons from numerical integration: 55.542987
Number of electrons from orbital occupations: 56.000000
Error in the number of electrons: -0.457013
Error larger than DFTELS (set input): 0.001000
Calculation aborted.

Откуда дробное количество электронов берется?

[endlesslake]

в прилагающихся тестах тоже ?

[kbob]

В тестах dft_hcth120, dft_qr_sym таже ошибка.

[kbob]

DALTON2013 is released. Note that users with a licence for Dalton2011 or previous versions must reregister to obtain a licence for DALTON2013.

Fourth, some of the new features:
• lsDalton:

o Dynamics using HF and DFT
o More response properties, including some magnetic properties like MCD
o DEC-MP2 energy, density and gradient
o Local orbitals
o Improved SCF optimization routines
o Massively parallel CCSD
o MPI parallelization of HF and DFT
o Improved integral code including MPI parallelism
o Matrix operation parallelization using PBLAS/SCALAPACK
o ADMM exchange (energy, gradients)

• Dalton:

o Subsystems CC using Cholesky decomposition
o Multiscale modeling using the Polarizable Embedding (PE) library.
o Static exchange (STEX) for X-ray spectroscopy
o Damped response via Complex Polarization Propagator (CPP)
o Quadratic response for open-shell DFT
o Relativistic corrections to nuclear shielding constants
o Empirical dispersion corrections DFT-D2, DFT-D3 and DFT-D3BJ
o Various performance improvements and a few bug fixes

We have enjoyed working out the theory for all the new features and implementing them. We hope that you will also find them useful.
On behalf of all the DALTON2013 authors,

Hans Jørgen Aa. Jensen