求助CASTEP做LEuH-Cl的结构优化的Optimize cell的设置问题

yoyo00

CASTEP中在做LEuH-Cl的结构优化(结构与层状双氢氧化物LDH相似），勾选Optimize cell,计算不收敛，不勾选Optimize cell 则可以收敛，但是计算出来的能量这样的（energy not corrected for finite basis set)，希望大神们指导指导！！！这是LEuH-Cl的结构，对称性为P21212

yoyo00

卡开发发 发表于 2020-5-11 22:31
Grimme的DFT-D3的程序包含DFT-D2的参数（https://www.chemie.uni-bonn.de/pctc/mulliken-center/software ...

谢谢，大神

卡开发发

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/)
    
32. 
    
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 /)
    

复制代码

yoyo00

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

大神，CASTEP中的没有Eu的色散校正的参数怎么办啊

yoyo00

sobereva 发表于 2020-5-9 02:58
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并反 ...

嗯嗯，好的

卡开发发

1、先确定结构是否正确。
2、考虑色散修正，并且做好动能截断的收敛性测试。
3、尽量使用对称性。

sobereva

如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并反映出帖子具体内容，仔细看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