计算化学公社

标题: 求助CASTEP做LEuH-Cl的结构优化的Optimize cell的设置问题 [打印本页]
作者
Author: yoyo00    时间: 2020-5-8 11:58
标题: 求助CASTEP做LEuH-Cl的结构优化的Optimize cell的设置问题
CASTEP中在做LEuH-Cl的结构优化(结构与层状双氢氧化物LDH相似），勾选Optimize cell,计算不收敛，不勾选Optimize cell 则可以收敛，但是计算出来的能量这样的（energy not corrected for finite basis set)，希望大神们指导指导！！！这是LEuH-Cl的结构，对称性为P21212
(, 下载次数 Times of downloads: 53)

(, 下载次数 Times of downloads: 620)


作者
Author: sobereva    时间: 2020-5-9 02:58
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并反映出帖子具体内容，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当、又啰嗦、又没体现帖子主要内容、又因为太长而没显示完整的标题“CASTEP中的结构优化，勾选Optimize cell,计算不收敛，不勾选Optimize cell,则可以s...”改了，以后务必注意恰当设置标题。

那个提示不是报错，看这里了解原理
https://www.tcm.phy.cam.ac.uk/castep/MD/node8.html
作者
Author: 卡开发发    时间: 2020-5-10 00:28
1、先确定结构是否正确。
2、考虑色散修正，并且做好动能截断的收敛性测试。
3、尽量使用对称性。
作者
Author: yoyo00    时间: 2020-5-11 20:28
sobereva 发表于 2020-5-9 02:58
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并反 ...

嗯嗯，好的
作者
Author: yoyo00    时间: 2020-5-11 20:30
卡开发发 发表于 2020-5-10 00:28
1、先确定结构是否正确。
2、考虑色散修正，并且做好动能截断的收敛性测试。
3、尽量使用对称性。

大神，CASTEP中的没有Eu的色散校正的参数怎么办啊
作者
Author: 卡开发发    时间: 2020-5-11 22:31
yoyo00 发表于 2020-5-11 20:30
大神，CASTEP中的没有Eu的色散校正的参数怎么办啊

Grimme的DFT-D3的程序包含DFT-D2的参数（https://www.chemie.uni-bonn.de/p ... t-version-of-dft-d3），Eu在其中有定义，不过单位是不同的，你需要自行换算一下。
  1. c the published radii in S.Grimme, J.Comput.Chem. 27, (2006), 1787-1799 (tab 1)
  2. c refer to the following values multiplied by 1.1 (rs6 in this code)
  3. c H, He
  4.          r0(1:86) = (/ 0.91d0,0.92d0,
  5. c Li-Ne
  6.      .      0.75d0,1.28d0,1.35d0,1.32d0,1.27d0,1.22d0,1.17d0,1.13d0,
  7. c Na-Ar
  8.      .      1.04d0,1.24d0,1.49d0,1.56d0,1.55d0,1.53d0,1.49d0,1.45d0,
  9. c K, Ca old
  10.      .      1.35d0,1.34d0,
  11. c Sc-Zn
  12.      .      1.42d0,1.42d0,1.42d0,1.42d0,1.42d0,
  13.      .      1.42d0,1.42d0,1.42d0,1.42d0,1.42d0,
  14. c Ga-Kr
  15.      .      1.50d0,1.57d0,1.60d0,1.61d0,1.59d0,1.57d0,
  16. c Rb, Sr
  17.      .      1.48d0,1.46d0,
  18. c Y-Cd
  19.      .      1.49d0,1.49d0,1.49d0,1.49d0,1.49d0,
  20.      .      1.49d0,1.49d0,1.49d0,1.49d0,1.49d0,
  21. c In, Sn, Sb, Te, I, Xe
  22.      .      1.52d0,1.64d0,1.71d0,1.72d0,1.72d0,1.71d0,
  23. c Cs,Ba,La,Ce-Lu
  24.      .      1.638d0,1.602d0,1.564d0,1.594d0,1.594d0,1.594d0,1.594d0,
  25.      .      1.594d0,1.594d0,1.594d0,1.594d0,1.594d0,1.594d0,1.594d0,
  26.      .      1.594d0,1.594d0,1.594d0,
  27. c Hf, Ta-Au
  28.      .      1.625d0,1.611d0,1.611d0,1.611d0,1.611d0,1.611d0,1.611d0,
  29.      .      1.611d0,
  30. c Hg,Tl,Pb,Bi,Po,At,Rn
  31.      .      1.598d0,1.805d0,1.767d0,1.725d0,1.823d0,1.810d0,1.749d0/)

  33.        c6(1:86) = (/0.14d0,0.08d0,
  34.      .   1.61d0,1.61d0,3.13d0,1.75d0,1.23d0,0.70d0,0.75d0,0.63d0,
  35.      .   5.71d0,5.71d0,10.79d0,9.23d0,7.84d0,5.57d0,5.07d0,4.61d0,
  36.      .   10.8d0,10.8d0,10.8d0,10.8d0,10.8d0,
  37.      .   10.8d0,10.8d0,10.8d0,10.8d0,10.8d0,10.8d0,10.8d0,16.99d0,
  38.      .   17.10d0,16.37d0,12.64d0,12.47d0,12.01d0,24.67d0,24.67d0,
  39.      .   24.67d0,24.67d0,24.67d0,24.67d0,24.67d0,24.67d0,24.67d0,
  40.      .   24.67d0,24.67d0,24.67d0,37.32d0,38.71d0,38.44d0,31.74d0,
  41.      .   31.50d0,29.99d0,315.275d0,226.994d0,176.252d0,
  42.      .  140.68d0,140.68d0,140.68d0,140.68d0,140.68d0,140.68d0,140.68d0,
  43.      .  140.68d0,140.68d0,140.68d0,140.68d0,140.68d0,140.68d0,140.68d0,
  44.      .  105.112d0,
  45.      .  81.24d0,81.24d0,81.24d0,81.24d0,81.24d0,81.24d0,81.24d0,
  46.      .  57.364d0,57.254d0,63.162d0,63.540d0,55.283d0,57.171d0,56.64d0 /)
复制代码
作者
Author: yoyo00    时间: 2020-5-12 08:11
卡开发发 发表于 2020-5-11 22:31
Grimme的DFT-D3的程序包含DFT-D2的参数（https://www.chemie.uni-bonn.de/pctc/mulliken-center/software ...

谢谢，大神



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