|
|
本帖最后由 ChemiAndy 于 2018-3-22 06:29 编辑
有些初学DFTB的朋友表示对DFTB+的输入文件设置有些困惑，我这里分享一些经验。
首先说一下DFTB，SCC-DFTB这几个概念的不同。DFTB是指基于密度泛函的紧束缚模型，这个模型是50年代提出的用于周期固体体系一种近似模型，使用“轨道能+排斥势”表示总能并通过拟合高级计算结果获取经验参数的一种方法。后来发现这个模型用于分子体系计算非常节约时间，可惜键长、生成热等性质误差太大。这可以通过进一步增加含参数的势函数以提高精确性。Elsner等人提出了基于电荷自洽的DFTB：SCC-DFTB。这个模型尝试在原DFTB的能量式后面增加一个多项式，这个多项式是关于原子电荷的泰勒展开，用于修正DFTB对交换、相关作用的误差。然后用高水平DFT计算结果去训练（拟合）这个多项式的参数。这个多项式展开到2级叫SCC-DFTB2，展开到3级叫SCC-DFTB3。其中还分加不加氢键校正（加氢键应该使用DampXH=Yes并设置gamma参数DampXHExponent）SCC-DFTB3还分ThirdOrder和ThirdOderfull两种，其中ThirdOrder只包含对角元项(on-site)，而ThirdOrderFull还包括交叉项。正常的SCC-DFTB3应该使用ThirdOrderFull，并且正确给出HubbardDerivs。
下面是各种不同方法的Hamiltonian部分的写法：（注意尽量使用最新方法即SCC-DFTB3-h）
SCC-DFTB2
- Hamiltonian = DFTB {
- SCC = Yes
- SlaterKosterFiles = Type2FileNames {
- Prefix = "/home/user/bin/dftb/1.2/slako/mio-0-1/"
- Separator = "-"
- Suffix = ".skf"
- }
- charge = 1.0 // not neccessary your case, default 0
- MaxAngularMomentum {
- C = "p"
- N = "p"
- H = "s"
- O = "p"
- }
- Filling = Fermi {
- Temperature [Kelvin] = 50
- }
- DampXH = No
- }
SCC-DFTB2-h (带氢键校正)
- Hamiltonian = DFTB {
- SCC = Yes
- SlaterKosterFiles = Type2FileNames {
- Prefix = "/home/user/bin/dftb/1.2/slako/mio-0-1/"
- Separator = "-"
- Suffix = ".skf"
- }
- charge = 1.0 // not neccessary your case, default 0
- MaxAngularMomentum {
- C = "p"
- N = "p"
- H = "s"
- O = "p"
- }
- Filling = Fermi {
- Temperature [Kelvin] = 50
- }
- DampXH = Yes
- DampXHExponent = 4.53
- }
- Hamiltonian = DFTB {
- SCC = Yes
- SlaterKosterFiles = Type2FileNames {
- Prefix = "/home/user/bin/dftb/1.2/slako/mio-0-1/"
- Separator = "-"
- Suffix = ".skf"
- }
- charge = 1.0
- MaxAngularMomentum {
- C = "p"
- N = "p"
- H = "s"
- O = "p"
- }
- Filling = Fermi {
- Temperature [Kelvin] = 50
- }
- Dispersion = SlaterKirkwood {
- PolarRadiusCharge = HybridDependentPol {
- C = {
- CovalentRadius [Angstrom] = 0.8
- HybridPolarisations [Angstrom^3, Angstrom,] = {
- 1.382 1382 1.382 1.064 1.064 1.064 3.8 3.8 3.8 3.8 3.8 3.8 2.50
- }
- }
- H = {
- CovalentRadius [Angstrom] = 0.4
- HybridPolarisations [Angstrom^3, Angstrom,] = {
- 0.386 0.386 0.000 0.000 0.000 0.000 3.5 3.5 3.5 3.5 3.5 3.5 0.80
- }
- }
- N = {
- CovalentRadius [Angstrom] = 0.8
- HybridPolarisations [Angstrom^3, Angstrom,] = {
- 1.030 1.030 1.090 1.090 1.090 1.090 3.8 3.8 3.8 3.8 3.8 3.8 2.82
- }
- }
- O = {
- CovalentRadius [Angstrom] = 0.75
- HybridPolarisations [Angstrom^3, Angstrom,] = {
- 0.560 0.560 0.000 0.000 0.000 0.000 3.8 3.8 3.8 3.8 3.8 3.8 3.15
- }
- }
- }
- }
- DampXH = Yes
- DampXHExponent = 4.53
- }
（SCC-DFTB3-h 见二楼）
|
评分 Rate
-
查看全部评分 View all ratings
|
[复制链接 Copy URL]
发表于 Post on 2018-3-9 21:44:40
收藏 Add to favorites
楼主 Author