计算化学公社

标题: DFT-D4色散校正的简介与使用 [打印本页]
作者
Author: sobereva    时间: 2019-2-13 19:50
标题: DFT-D4色散校正的简介与使用
DFT-D4色散校正的简介与使用
Introduction and use of DFT-D4 dispersion correction

文/Sobereva@北京科音
First release: 2019-Feb-13  Last update: 2024-Sep-5


之前笔者在《DFT-D色散校正的使用》（http://sobereva.com/210）中详细介绍了DFT-D3的函数形式以及使用方法，还有其它相关的文章：《乱谈DFT-D》（http://sobereva.com/83）、《谈谈“计算时是否需要加DFT-D3色散校正？”》（http://sobereva.com/413）。大家如果对DFT-D3尚不了解的话请先阅读以上几篇。本文笔者简要介绍一下DFT-D3的后继者DFT-D4以及其实际使用方法。由于笔者精力有限，此文在原理方面就不写得很详细了。

1 关于DFT-D4方法

DFT-D3色散校正如今已极度流行，不光可以使那些计算弱相互作用差的泛函在计算弱相互作用精度方面有质的改进，对于计算很多中、大体系，特别是涉及过渡金属反应的热力学数据往往也有明显改进。DFT-D3最大的一个软肋就是没法考虑体系的实际电子结构。同种元素的原子，在不同化学体系中的电子结构显然是不同的，这从原子电荷上就可以直接体现出来。特别是过渡金属，在不同配合物里往往以不同氧化态形式存在，虽然其原子电荷的差异没有形式上的氧化态的差异那么夸张（参见《使用Multiwfn通过LOBA方法计算氧化态》http://sobereva.com/362里的一些论述），但是差异还是相当显著的。DFT-D方法计算色散校正能是基于原子的色散系数计算的，这从《DFT-D色散校正的使用》一文的公式里可以看得很清楚。实际上原子的色散系数是显著依赖于原子的电子结构的，特别是原子电荷这个直接表征原子在分子中带的净电荷的量对色散系数的影响是明显不可忽视的。但是由于DFT-D3没有考虑这点，而只纯粹依赖于体系的几何结构，当体系中存在显著带电的原子的时候，比如过渡金属配合物、离子化合物，DFT-D3在原理上就不是特别理想了（尽管从实际表现上来看，对于牵扯过渡金属的体系，用DFT-D3比不用一般还是明显有改进的）。

色散校正的方法很多，有的是像DFT-D3一样不考虑实际电子结构，如OBS、DFT-ulg；也有的是能体现实际电子结构的，如XDM、dDSC、TS、LRD、vdW-DF/VV10、MBD@rsSCS等。那些能体现实际电子结构的做法普遍比形式超级简单的DFT-D3更昂贵，编程实现往往更费劲，大多没有解析梯度，被主流程序支持得也比较少，比如笔者撰文时最新的G16 B.01一个都不直接支持，而ORCA 4.1.1版只支持VV10（也被叫做DFT-NL）。

为了解决DFT-D3的上述不足，Grimme在JCP, 147, 034112 (2017)中提出了DFT-D4，关键性的改进就是让原子色散系数与原子电荷挂钩，从而能更好地反应实际情况。2017年这篇文章实际上是DFT-D4的最初版，发表之后，Grimme一直也没公开提供能做DFT-D4计算的程序，令人觉得很莫名其妙，不禁猜测是不是DFT-D4还存在什么问题。后来在2018年7月，Grimme往ChemRxiv上上传了介绍DFT-D4最终版文章的v1版，在2019年1月25日又上传了这篇文章的v2版，参见https://doi.org/10.26434/chemrxiv.7430216.v2。我估计过不多久，这篇介绍DFT-D4最终版本的重量级文章就会正式在期刊上发表了。目前计算DFT-D4的程序已经公开了。
后记：DFT-D4最终版的正式文章已发表于J. Chem. Phys. 150, 154122 (2019)。

根据这篇预印版的DFT-D4最终形式的文章，可以看到最终形式的DFT-D4会有以下特征：
·和DFT-D3(BJ)一样使用BJ阻尼函数。目前DFT-D4已经对所有主流泛函都拟合了阻尼参数
·带有和DFT-D3一样的Axilrod–Teller–Muto (ATM)形式的三体校正项。DFT-D4也可以结合MBD形式的三体校正项，但更为昂贵，而且实测的精度并没有ATM那么好
·色散系数依赖于electronegativity equilibration (EEQ)方法计算的原子电荷。EEQ这个方法和已经很流行的EEM方法差不多，都是基于电负性均衡原理的思想，依赖于和元素有关的经验参数，可以瞬间计算很大体系的原子电荷。EEM的一些介绍可以看Multiwfn程序手册的3.9.15节
·有解析梯度（Hessian目前只有半数值的）。这是VV10之类其它依赖于电子结构的色散校正方法所不具备的
·最高支持到86号元素（Rn）
·和DFT-D3一样基本上都是零耗时
·根据Grimme的测试，DFT-D4的精度能达到VV10这样更复杂的色散校正方法的精度，甚至有时更好

对于那些不存在显著带电原子的情况，D4和D3精度半斤八两，至少不输于D3。而对于有显著带电原子的情况，D4比D3的改进还是挺明显的，除了上述提到的Grimme的D4相关文章里有体现外，还可以看他的Acc. Chem. Res., 52, 258 (2019)，里面专门展示了DFT-D3、D4在金属有机体系上的效果。

DFT-D4原理上也可以结合其它原子电荷。2017年那篇JCP里初代DFT-D4结合的是GFN-xTB方法（一种半经验形式的DFT方法，类似DFTB）计算时顺带产生的Mulliken电荷，结合这种电荷时被称为DFT-D4(TB)，但测试发现这种形式不仅算原子电荷的过程更麻烦（需要依赖于专做GFN-xTB的xtb程序），增加了耗时，结果还不如结合EEQ电荷的最终版DFT-D4好。

Grimme的文章也表示，DFT-D这种色散校正形式到了DFT-D4这一代就基本没有油水进一步可榨了。由于最终形式的DFT-D4计算又快效果又好，应该在未来会很快普及、被大量量化程序直接支持。而且笔者也认为，有了DFT-D4，VV10等绝大部分以往提出的能够展现实际电子结构的色散校正方法就可以彻底退出历史舞台了（但有个别这些方法由于仗着更高级的计算形式，精度可能做到比DFT-D4更高一些，所以还不是完全没存在意义）。

当然，未来肯定有很多DFT-D4相关的测试文章接踵而至，也不排除DFT-D4方法的定义又出现一些改变，以上说法只是笔者撰文时的情况和看法。


2 Grimme的dftd4程序

2.1 简介

Grimme写的独立的做DFT-D4的程序dftd4已经发布在了https://www.chemie.uni-bonn.de/pctc/mulliken-center/software/dftd4，可以直接下载源代码包。下文说的是笔者2019年2月13日在此处下载到的dftd4.2.0.tar.xz的情况，以后版本可能有些地方会有变动。

这个DFT-D4程序用法、用处和DFT-D3很相似。由于此程序里面没有包含xtb程序，因此做的DFT-D4只能是基于EEQ电荷的（如上所述，这也正是标准的DFT-D4的情况）。此程序的三体校正项支持前述的ATM和MBD形式，因此确切来说，此程序可以算的DFT-D4具体形式包括DFT-D4(EEQ)-ATM和DFT-D4(EEQ)-MBD。

dftd4得先编译才能运行。为便于读者使用，按照本节做法编译好的Linux版和Windows版我直接提供在此：http://sobereva.com/soft/dftd4.2.0_binary.rar。其中有.exe的是Win32版（libiomp5md.dll是其运行时需要的动态库文件，应放到与之相同的目录），没后缀的是Linux版，都是可执行文件，直接就能运行。

2.2 最简单粗暴的编译方法

下面说一下如何以最简单的方式编译。

Windows 32bit版：笔者机子里装了Visual Studio 2017以及Intel Parallel Studio XE 2019 Update 1里面的Intel Fortran compiler和Intel MKL库。启动Visual Studio 2017，新建一个项目，把dftd4压缩包里的include目录下的.f文件和source目录下的.f90文件都放到项目的文件夹里，并且加入项目里。把编译的配置从默认Debug改为Release。之后在项目属性里，在Fortran标签页下把Libaries里的Runtime Library设成Multithreaded，Use Intel Math Kernel Library设Sequential。Language标签页里Process OpenMP Directives设Generate Parallel Code。之后如果你愿意的话还可以去设对CPU和指令集的优化选项，这里就不用了。然后对项目直接编译，项目目录下的Release目录里就产生了dftd4.exe可执行文件，直接就能用了，而且传给别人别人也能用。

Linux版：笔者的机子里已经装了Intel Fortran compiler 19.0.1.144和MKL数学库（是按照《VASP最简单的安装方法》http://sobereva.com/455里面的做法装的Intel Parallel Studio XE 2019 Update 1里面的）。将dftd4压缩包解压后，将include目录下的.f文件拷入source目录，再进入source目录，直接运行以下命令即在当前目录下产生了dftd4可执行文件
ifort *.f90 -qopenmp -qopenmp-link=static -mkl -static-intel -stand f08 -check bounds -o dftd4
其中那些带static的选项是为了让编译出的可执行文件不依赖于动态库，从而在别人的机子里也确保能运行。

2.3 标准但麻烦的编译方法

dftd4基于meson+ninja的代码生成系统组合，如果你想形式上最“优雅”地编译（吃饱了撑的。其实没任何好处），需要系统里有meson和ninja。这俩玩意儿比较新，可能很多人不熟悉，这俩程序其实就相当于常见的cmake和make，只不过有一些额外的特点而已。笔者认为DFT-D4这么简单的程序用mason+ninja纯粹是莫名其妙，也给用户添麻烦，一般系统里都不自带它们，还得现装。下面我们来通过meson+ninja编译DFT-D4，机子里已经装了上述Intel Fortran编译器和MKL库。笔者也尝试了用此系统自带的gfortran 4.8.5，发现编译不过去。

先运行以下命令以安装meson、lapack静态库（顺带会安装blas静态库）、ninja。由于yum这么安装之后ninja的可执行文件名是ninja-build，所以得做个符号链接成为ninja
yum install lapack-devel
yum install meson
yum install ninja-build
ln -s /usr/bin/ninja-build /usr/bin/ninja

将下载到的DFT-D4程序的压缩包解压，进入此目录，直接将以下内容输入命令行窗口，就会通过ifort来编译
FC=ifort meson setup build && ninja -C build
然后在当前目录下的build目录下就出现了名为dftd4的可执行程序，可以直接用了。编译过程中终端里会提示什么不识别-pipe之类选项的警告，不用管它，只要程序能用就行了。如果想运行得更快，而且你的机子里有MKL库的话，可以把编译时依赖的lapack给替换为MKL库，把dftd4目录下自带的meson.build替换为我提供的这个即可：http://sobereva.com/attach/464/meson.build。

2.4 基本用法

在含有dftd4的目录下直接用./dftd4就可以显示出此程序可以接的选项，一看就懂，就不多说了。比如下面的命令是对H2O.xyz这个文件里的体系基于BLYP泛函的参数计算标准形式的DFT-D4，即DFT-D4(EEQ)-ATM形式的色散校正能
./dftd4 H2O.xyz --func blyp

瞬间看到以下信息
           -------------------------------------------------
          |                Calculation Setup                |
           -------------------------------------------------
   coordinate file      : H2O.xyz
   number of atoms      :      3
   charge               :      0
   non-additivity corr. : ATM
   charge model         : EEQ
   functional           : b-lyp
   omp threads          :      1
   memory needed (est.) :    0.0 Mb


O      q(ref)        CN(ref)   cov. CN(ref)     α(AIM,ref)
     0.0000000      0.0000000      0.0000000      5.1967092
    -0.3547147      0.9924594      0.8041679      6.1223924
    -0.5988153      1.9886928      1.6112394      6.7491941
     0.0000000      0.9985672      0.9798440      5.1987934
H      q(ref)        CN(ref)   cov. CN(ref)     α(AIM,ref)
     0.0000000      0.0000000      0.0000000      5.0540161
     0.0000000      0.9117922      0.8942243      2.7207580

           -------------------------------------------------
          |               Damping Parameters                |
           -------------------------------------------------
   s6                   :     1.0000
   s8                   :     2.3408
   a1                   :     0.4449
   a2                   :     4.0933

           -------------------------------------------------
          |                     Results                     |
           -------------------------------------------------
Edisp  /kcal,au:    -0.2980  -0.00047485

dispersion energy written to file .EDISP

即色散校正能为-0.00047485 Hartree。如果你用cat .EDISP显示刚产生的.EDISP隐藏文件，会看到更高精度的数据：-0.00047485257046358403。

dftd4另外常用的选项是--chrg，用于设定体系净电荷；--grad，用于计算色散校正能的梯度；--hess，用于计算色散校正能的半数值的Hessian。

DFT-D4程序支持的泛函在dftd4源代码包里dfuncpar.f90里可以看到。在这个文件里搜get_d4eeqbjatm，凡是下面有名字的都可以用。值得一提的是，截止到2019-May-8，目前官网上的dftd4的这个参数库文件里仍没有常用的M06-2X泛函的参数，因此DFT-D4没法结合M06-2X使用。


3 在量子化学程序中使用DFT-D4

后记：从ORCA 4.2开始，支持且只支持DFT-D4(EEQ)-ATM形式的DFT-D4，直接写个D4关键词即可使用。下面的关于ORCA的文字可以忽略了。

ORCA从4.1开始支持了DFT-D4，实际上是调用程序目录里otool_dftd4可执行文件算的，这个文件其实就是专供ORCA的dftd4。只要在ORCA输入文件里简简单单地写上D4关键词，计算时就会自动计算出D4校正能并加到总能量上。ORCA也支持DFT-D4下的优化任务。

对于笔者撰文时最新的ORCA 4.1.1，大家会发现ORCA输出的DFT-D4校正能和直接用Grimme的dftd4程序给出的不同（选的都是相同泛函的情况下）。这是因为起码对于ORCA 4.1.1来说，默认算的是DFT-D4(TB)形式的校正能，因此计算时会自动调用ORCA自带的otool_xtb（专供ORCA的xtb程序）做GFN-xTB计算。由于DFT-D4(EEQ)-ATM形式更便宜且结果更好，我预感ORCA以后早晚会把D4关键词的默认情况改成DFT-D4(EEQ)-ATM形式。


Molpro从2018.1开始支持了DFT-D4。Turbomole从7.3也支持了DFT-D4。

CP2K从2024.2开始支持了DFT-D4。

不知Gaussian什么时候会正式支持DFT-D4，估计是早晚的事。对于尚不支持DFT-D4的Gaussian，如果想在Gaussian计算中使用的话，可以利用《将Gaussian与ORCA联用搜索过渡态、产生IRC、做振动分析》（http://sobereva.com/422）和《将Gaussian与Grimme的xtb程序联用搜索过渡态、产生IRC、做振动分析》（http://sobereva.com/421）里介绍的Gaussian的external关键词自己写个接口让Gaussian可以在计算能量、受力、Hessian的时候调用DFT-D4，实现起来很简单。
作者
Author: ene    时间: 2019-2-13 23:37
本帖最后由 ene 于 2019-2-13 23:51 编辑

有个非常zz但是管用的编译方法，可以不用依赖meson和ninja。
首先进入 source/ 目录，执行
  1. cp ../include/param_ref.f ./
复制代码
确保系统内有lapack.a blas.a，把这两个库拷贝到当前目录下
最后编译（P.S. 如果要使用 -static 命令，我这边不能同时使用 -fopenmp 命令，否则编译出来的文件没法运行。不清楚具体是WSL的bug还是代码本身的问题。）

编译的前几次会失败，反复执行编译命令到没有报错为止，在当前目录下就会产生dftd4这个可执行文件。

  1. gfortran -O3 -fbounds-check -fopenmp -std=f2008 -o dftd4 program_dftd.f90 printout.f90 mctc_global.f90 mctc_systools.f90 mctc_readin.f90 mctc_econv.f90 mctc_timings.f90 error.f90 class_set.f90 class_molecule.f90 class_param.f90 argparser.f90 geometry_reader.f90 disp_output.f90 eeq_model.f90 coordination_number.f90 dftd4.f90 dfuncpar.f90 lapack.a blas.a
复制代码






作者
Author: Accelerator    时间: 2019-2-14 01:19
“日落西山的turbomole”xswl
和turbomole什么仇什么怨系列(
作者
Author: sobereva    时间: 2019-2-14 06:24
Accelerator 发表于 2019-2-14 01:19
“日落西山的turbomole”xswl
和turbomole什么仇什么怨系列(

都快被ORCA挤下悬崖了，还在收费
母公司还被收购了
作者
Author: sobereva    时间: 2019-2-14 07:53
今日对本文进行了更新，直接添加了已经编译好的Linux和Windows下的dftd4可执行文件的下载链接，并且加入了如何在Windows和Linux下极其简单地编译dftd4的介绍。在文末提了一嘴在Gaussian计算中利用DFT-D4的思路。
作者
Author: granvia    时间: 2019-2-14 12:33
采用EEQ估算电荷就是为了能支持解析梯度吧？
作者
Author: Accelerator    时间: 2019-2-14 14:57
sobereva 发表于 2019-2-14 06:24
都快被ORCA挤下悬崖了，还在收费
母公司还被收购了

不过很奇怪我在德国这边看到turbomole的用户非常多 相反他们不太用ORCA算DFT 不知道为什么
作者
Author: sobereva    时间: 2019-2-14 19:38
Accelerator 发表于 2019-2-14 14:57
不过很奇怪我在德国这边看到turbomole的用户非常多 相反他们不太用ORCA算DFT 不知道为什么

毕竟是德国开发的，历史较久，因此必然有很多非学术层面的因素以及历史性因素
作者
Author: michaelm    时间: 2019-2-27 00:48
这个版本只能算总校正能？能像dft-d3一样，算片段之间的校正么？
作者
Author: sobereva    时间: 2019-2-27 00:50
michaelm 发表于 2019-2-27 00:48
这个版本只能算总校正能？能像dft-d3一样，算片段之间的校正么？

目前版本没有现成选项
作者
Author: michaelm    时间: 2019-3-2 14:06
问个土问题：这个色散校正能，理论上应该怎么正确地得到？是比较dft和更高级的计算方法例如ccsd之间的能差吗？
试着玩了一下orca里的dft-4和vv10，得到的校正能天差地别，一个正一个负。怎么知道哪个更正确？
作者
Author: sobereva    时间: 2019-3-2 15:06
michaelm 发表于 2019-3-2 14:06
问个土问题：这个色散校正能，理论上应该怎么正确地得到？是比较dft和更高级的计算方法例如ccsd之间的能差 ...

理论上DFT-D3/4是以近似方式计算原子的动态极化率，然后通过Casimir–Polder公式计算出原子间的C6系数（然后再获得C8等系数），带入进形式上略类似于分子力场的DFT-D计算公式得到。虽然这个过程每一步都牵扯近似，但是思想是严格的。

DFT-D4是没问题的。
vv10用的是明显不同的思路。化学上感兴趣的是相对能量，因此比较两种校正得到的绝对能量是没太大意义的。
作者
Author: michaelm    时间: 2019-3-5 14:39
今天试了一下算双甲烷分子之间的色散修正。

sapt2/aug-cc-pvtz大概是-0.93kcal/mol。
b3lyp-d3bj给出的修正是-0.90，vv10是-0.71，都还算靠谱。orca里的b3lyp-d4给出的修正是+2.86，完全不靠谱，可能是现在orca里用的GFN2-xTB模型的电荷不适合用在d4？建议大家暂时不要用orca目前版本里的d4。

作者
Author: sobereva    时间: 2019-3-5 19:26
michaelm 发表于 2019-3-5 14:39
今天试了一下算双甲烷分子之间的色散修正。

sapt2/aug-cc-pvtz大概是-0.93kcal/mol。

独立的dftd4算的结果是多少？
作者
Author: michaelm    时间: 2019-3-5 22:20
本帖最后由 michaelm 于 2019-3-6 04:13 编辑

独立的dft-d4没有片段之间的选项。用super molecule的方法算的话，是-0.81。我还没有仔细读d4的paper。不知道d4有没有类似BSSE的问题需要考虑。例如EEQ方法，同一个原子的电荷在单体和双体时候，是不是会有不同。如果有的话，用super molecule方法会引入额外的误差，这种误差在独立的dft-d4不知道怎么搞掉。
作者
Author: sobereva    时间: 2019-3-6 07:20
michaelm 发表于 2019-3-5 22:20
独立的dft-d4没有片段之间的选项。用super molecule的方法算的话，是-0.81。我还没有仔细读d4的paper。不知 ...

supramolecular方式结果看来合理。预感早晚会给dftd4加上计算片段间项的功能，毕竟dftd4刚出来。原理上讲supramolecular算片段间D4校正能是不妥的。
这和BSSE问题没有直接联系
作者
Author: hzamber    时间: 2019-5-11 20:59
使用orca进行DFT-D4计算，orca版本：Program Version 4.1.1  - RELEASE  -
输入文件：
!b3lyp d4 def2-tzvp  pal2
*xyz 0 1
O                 -1.29527542   -0.95244604   -0.08439711
H                 -1.96684629   -1.20450472    0.55440905
H                 -1.21891600    0.01667075   -0.01273795
O                  1.47586351   -0.63755538   -0.07513260
H                  2.03414372   -1.07560919    0.57227133
H                  0.61081072   -1.08144807   -0.01749750
O                 -0.17970784    1.59357918    0.12071627
H                 -0.09171077    2.24303363   -0.58191728
H                  0.62465016    1.04596562    0.06646441
*

计算结果：
-------------------------------------------------------------------------------
                          DFT DISPERSION CORRECTION                           
                                                                              
                                 DFTD4 V1.0                                    
-------------------------------------------------------------------------------
The B3LYP functional is recognized
Active option DFTDOPT                   ...         5   

-------------------------   ----------------
Dispersion correction            0.000000000
-------------------------   ----------------


-------------------------   --------------------
FINAL SINGLE POINT ENERGY      -229.305319165147
-------------------------   --------------------
计算结果并没有进行D4校正计算，计算设置有问题？

作者
Author: yflchx    时间: 2019-5-11 21:36
本帖最后由 yflchx 于 2019-5-11 21:42 编辑
michaelm 发表于 2019-3-2 14:06
问个土问题：这个色散校正能，理论上应该怎么正确地得到？是比较dft和更高级的计算方法例如ccsd之间的能差 ...

Ec(VV10) = Ec(nl) + βN          nl即为nonlocal

β是与b有关的参数，N是电子数；上述两项之和为正。Ec(nl)为负，可以与-D相类比。其实，单纯作为色散校正Ec(nl)即可。
作者
Author: sobereva    时间: 2019-5-12 10:27
hzamber 发表于 2019-5-11 20:59
使用orca进行DFT-D4计算，orca版本：Program Version 4.1.1  - RELEASE  -
输入文件：
!b3lyp d4 def2-tz ...

你得确保dftd4模块正常运行了才行，比如没出现找不到库文件之类。如果没正常运行，返回的数值就为0。

-------------------------------------------------------------------------------
                          DFT DISPERSION CORRECTION                           
                                                                              
                                 DFTD4 V1.0                                    
-------------------------------------------------------------------------------
The B3LYP functional is recognized
Active option DFTDOPT                   ...         5   

OMP: Warning #181: OMP_STACKSIZE: ignored because KMP_STACKSIZE has been defined
normal termination of dftd4
-------------------------   ----------------
Dispersion correction           -0.004195560
-------------------------   ----------------


-------------------------   --------------------
FINAL SINGLE POINT ENERGY      -229.309514679398
-------------------------   --------------------

作者
Author: hzamber    时间: 2019-5-12 14:48
sobereva 发表于 2019-5-12 10:27
你得确保dftd4模块正常运行了才行，比如没出现找不到库文件之类。如果没正常运行，返回的数值就为0。

...

使用批处理脚本进行的计算，没有注意到提示，确实有提示没有找到库文件，如下：
error while loading shared libraries: libmkl_intel_lp64.so
: cannot open shared object file: No such file or directory
需要安装intel编译器?
作者
Author: sobereva    时间: 2019-5-12 17:01
hzamber 发表于 2019-5-12 14:48
使用批处理脚本进行的计算，没有注意到提示，确实有提示没有找到库文件，如下：
error while loading s ...

是。或者从别人那里把这个文件拷过来
作者
Author: hzamber    时间: 2019-5-12 18:43
sobereva 发表于 2019-5-12 17:01
是。或者从别人那里把这个文件拷过来

谢谢社长回复。
这个文件是否可以跨平台使用，如Centos的编译文件复制到Suse linux中使用？
作者
Author: sobereva    时间: 2019-5-13 00:12
hzamber 发表于 2019-5-12 18:43
谢谢社长回复。
这个文件是否可以跨平台使用，如Centos的编译文件复制到Suse linux中使用？

可能可以
作者
Author: Novice    时间: 2019-8-10 18:06
社长果然厉害，ORCA4.2版The charge dependent atom-pairwise dispersion correction (keyword ! D4) is using the D4(EEQ)-ATM dispersion model [123], other D4 versions, using tight-binding partial charges, are currently only available with the standalone DFT-D4 program.
作者
Author: sobereva    时间: 2019-8-11 06:41
Novice 发表于 2019-8-10 18:06
社长果然厉害，ORCA4.2版The charge dependent atom-pairwise dispersion correction (keyword ! D4) is us ...

实测，4.2版的D4对应的确实是D4(EEQ)-ATM结果。但程序输出有个bug，一开头显示的仍是
Your calculation utilizes the atom-pairwise dispersion correction
based on tight binding partial charges (D4)
Cite in your paper:
E.Caldeweyher, C.Bannwarth, S.Grimme, J.Chem.Phys., (2017), 147, 034112.
E.Caldeweyher, S.Ehlert, A.Hansen, H.Neugebauer, S.Grimme, J.Chem.Phys., (2018), in progress.
这很有误导性
作者
Author: Reminder    时间: 2020-11-27 19:26
Sob老师您好，请问dftd4是否不支持86号元素以后的，即包含锕系元素的部分？
作者
Author: sobereva    时间: 2020-11-30 10:08
Reminder 发表于 2020-11-27 19:26
Sob老师您好，请问dftd4是否不支持86号元素以后的，即包含锕系元素的部分？

不支持。原文摘要部分就已经写明了：J. Chem. Phys. 150, 154122 (2019)
作者
Author: lao7    时间: 2020-11-30 12:00
orca能管理好自己的内存使用就好了，有时候设置内存不合理，会导致算了一段时间，到振动分析时中断报错。这点确实不如Gaussian。
作者
Author: scottwu    时间: 2020-12-22 17:37
想问一下各位大神，DFT-D4算出来的校正能，和总能”直接相加“就可以了吧？（是相加关系，不用改变符号了吧？）
作者
Author: wzkchem5    时间: 2020-12-22 17:46
scottwu 发表于 2020-12-22 17:37
想问一下各位大神，DFT-D4算出来的校正能，和总能”直接相加“就可以了吧？（是相加关系，不用改变符号了吧 ...

对，一般如果一个能量项为负，这个负号在程序输出的数值里就会带着
作者
Author: scottwu    时间: 2020-12-22 18:18
wzkchem5 发表于 2020-12-22 17:46
对，一般如果一个能量项为负，这个负号在程序输出的数值里就会带着

多谢！！！
作者
Author: scottwu    时间: 2021-4-14 11:55
对于涉及结构优化，如果不能在优化过程中把D4考虑进去，那么对优化完成后的结构才进行D4校正是不是就没意义了？
作者
Author: sobereva    时间: 2021-4-14 12:48
scottwu 发表于 2021-4-14 11:55
对于涉及结构优化，如果不能在优化过程中把D4考虑进去，那么对优化完成后的结构才进行D4校正是不是就没意义 ...

否
得把前提明确清楚，用的什么泛函，算的什么问题
作者
Author: scottwu    时间: 2021-5-4 20:03
sobereva 发表于 2021-4-14 12:48
否
得把前提明确清楚，用的什么泛函，算的什么问题

就对于一个用PBE优化得到稳定结构的例子来说，因为我发现OUTCAR中对于能量的收敛标准，参考的是范德华校正后的能量
作者
Author: sobereva    时间: 2021-5-5 00:08
scottwu 发表于 2021-5-4 20:03
就对于一个用PBE优化得到稳定结构的例子来说，因为我发现OUTCAR中对于能量的收敛标准，参考的是范德华校 ...

如果你的体系用PBE优化结果就合理，算能量再考虑DFT-D4校正原理上是没有问题的
作者
Author: scottwu    时间: 2021-5-6 09:12
sobereva 发表于 2021-5-5 00:08
如果你的体系用PBE优化结果就合理，算能量再考虑DFT-D4校正原理上是没有问题的

受教受教，多谢sob老师
作者
Author: 吕安琪    时间: 2021-9-6 12:34
sobereva 发表于 2019-2-14 07:53
今日对本文进行了更新，直接添加了已经编译好的Linux和Windows下的dftd4可执行文件的下载链接，并且加入了 ...

sob老师 请问现在高斯可以支持DFT-D4计算吗
作者
Author: sobereva    时间: 2021-9-6 15:11
吕安琪 发表于 2021-9-6 12:34
sob老师 请问现在高斯可以支持DFT-D4计算吗

不
能不能这种问题自行看官网上的手册便知
作者
Author: 七尺贱    时间: 2023-2-20 14:07
sobereva 发表于 2021-9-6 15:11
不
能不能这种问题自行看官网上的手册便知

老师我想问一下M06-2X里面描述色散的部分和DFT-D3这种有什么区别吗，计算弱相互作用加上D3是否更为稳妥
作者
Author: sobereva    时间: 2023-2-20 14:09
七尺贱 发表于 2023-2-20 14:07
老师我想问一下M06-2X里面描述色散的部分和DFT-D3这种有什么区别吗，计算弱相互作用加上D3是否更为稳妥

M06-2X对色散作用的描述能力是体现在交换-相关泛函的定义的层面上的，不是后校正
此文说了
谈谈“计算时是否需要加DFT-D3色散校正？”
http://sobereva.com/413（http://bbs.keinsci.com/thread-9772-1-1.html）
作者
Author: 七尺贱    时间: 2023-2-20 14:23
sobereva 发表于 2023-2-20 14:09
M06-2X对色散作用的描述能力是体现在交换-相关泛函的定义的层面上的，不是后校正
此文说了
谈谈“计算 ...

谢谢老师
作者
Author: zhangbo11011    时间: 2024-7-23 10:38
Reminder 发表于 2020-11-27 19:26
Sob老师您好，请问dftd4是否不支持86号元素以后的，即包含锕系元素的部分？

您好，请问您对于锕系元素怎么进行色散矫正的



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3