Gamess и гибридные функционалы

[Shorku]

Боюсь показаться человеком, не умеющим читать мануалы, но сам почему-то найти не смог: можно ли в GamessUS подкручивать коэффициенты при DFT- и HF- обмене в гибридных функционалах? Например, выкрутить в PBE0 обмен до 40% или 50%.

[Ferom]

В двойных гибридных можно. Поищите в мануале подраздел
* * * B2x-PLYP double hybrid functionals * * *

Если же очень надо именно PBE0, то можно попробовать в файле dftxca.f исходного кода GAMESS исправить доли в строках

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

      ELSEIF(DFTYPE.EQ.PBE0) THEN
         NEXFG     = 4
         NCORFG    = 4
         EXCHR     = 'PBE&HFX'
         CORCHR    = 'PBE'
         CPBEX     = 0.75D+00
         DFTTYP(3) = 0.25D+00

[Shorku]

В двойных гибридных можно. Поищите в мануале подраздел
* * * B2x-PLYP double hybrid functionals * * *

Если же очень надо именно PBE0, то можно попробовать в файле dftxca.f исходного кода GAMESS исправить доли в строках

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

      ELSEIF(DFTYPE.EQ.PBE0) THEN
         NEXFG     = 4
         NCORFG    = 4
         EXCHR     = 'PBE&HFX'
         CORCHR    = 'PBE'
         CPBEX     = 0.75D+00
         DFTTYP(3) = 0.25D+00

Отлично, спасибо. А Вы не в курсе, в исходниках только одна (в смысле две, но в одном месте) цифирь и в одном месте зашита? Впрочем, поэкспериментирую)

[Гесс]

Мы накручивали в ADF, но он небесплатный. С гамессом не работал.

[Shorku]

Мы накручивали в ADF, но он небесплатный. С гамессом не работал.

В том и фишка, что много где это делается штатными средствами, а в гамессе я такого не нашел)

[sanya1024]

В FireFly можно:

HFX N1,N2 The weight of the exact (i.e. Hartree-Fock) exchange. This option is supported only for hybrid functionals. The default value is set as is appropriate for a particular value of DFTTYP. For non-self-consistent double-hybrid functionals an extended form of this keyword exists, namely HFX(2), which allows one to change the weight of the exact exchange in the parent DFT functional (see the description of PARENT keyword below).

PARENT For non-self-consistent double-hybrid functionals, this keyword defines the parent DFT functional that is used to generate the MOs and density. This keyword can be used to give the name of any supported non-local functional. The default value is PARENT=DEFAULT which means to use the parent functional used in the original formulation of the double-hybrid functional used. This option pertains to non-self-consistent double hybrid functionals only.

Сразу говорю: сама не пробовала (не видела необходимости).

[alxyppv]

подтверждаю, в FireFly очень просто делается. Вот, например, как сделать 15% в PBE0:

$CONTRL RUNTYP=energy MAXIT=240 SCFTYP=RHF DFTTYP=PBE0 $END
$DFT HFX=0.15 $END

[Shorku]

FireFly, к сожалению, не подходит

[sanya1024]

А чем не подходит? если нужны TDDFT градиенты, то в последней версии они уже есть.

[Shorku]

А чем не подходит? если нужны TDDFT градиенты, то в последней версии они уже есть.

SF TD-DFT, упоминаний о нем в firefly не нашел

[ Post made via Android ]

[sanya1024]

А, это да, только Гамесс, только хардкор (ну или Q-chem, но за деньги).

[Shorku]

(ну или Q-chem, но за деньги).

они дают демку на пару месяцев пощупать, по ощущениям работает гораздо быстрее гамесса)