| 查看: 3359 | 回复: 20 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
cenwanglai荣誉版主 (知名作家)
老和山猥琐派九段
|
[交流]
【求助完毕】硫酸根离子的结构优化，偶极修正及势能修正。 已有2人参与
|
|||
|
SO4放在1nm边长的立方盒子中间,中性分子簇价电子数为30. 对于离子的计算，按照结构优化和自洽计算两个阶段是否带电可分为一下四种情况 带电（NELECT=32） 不带电(NELECT=30) 结构优化 静态自洽计算 A) 带电 带电 E1=-35.994eV B） 不带电 带电 E2=-35.765eV C） 带电 不带电 -- D） 不带电 不带电 E4=-26.539eV 以上结果说明，硫酸根是可以稳定单独存在的。而且，结构优化是否带电对体系基态是有显著影响的，二者能量差|E1-E2|=0.23eV. 但是这里没有加偶极修正。 在结构优化和自洽计算都带电且考虑偶极修正（NELECT=32, IDIPOL=4)后，在静态自洽文件夹中： grep "TOTEN" OUTCAR free energy TOTEN = -34.6451eV *********** adding dipol energy to TOTEN NOW ************* free energy TOTEN = -28.091eV 看起来考虑偶极修正后，差异还是很大的。 问题： （1）对于slab结构吸附带电离子，结构优化和静态计算的时候需要加电荷。在结构优化的同时就需要考虑偶极修正（slab IDIPOL =1~3,带电小分子 IDIPOL=4）？这样计算好像要慢很多，各位一般是怎么做的？ （2）还有一项势能修正（LDIPOL=.TRUE.），这个手册上说对于带电体系只对立方盒子有效。因为带电体系要做四极矩修正，这个只能在立方盒子超胞下才能做的。但是大部分超胞吸附模型都不是立方盒子，比如我的是六方的。吸附结构引入高介电常数的底材，势能修正项是否可以忽略呢？各位在计算的时候一般是怎么考虑势能修正这一项的呢？ （3）偶极修正（dipol correction）好理解，势能修正(potential correction)是什么意思呢？ （4）手册上讲：For charged cells or for calculations of molecules and surfaces with a large dipole moment, the energy converges very slowly with respect to the size L of the supercell. ...for charged systems with excess electrons more reliable results can be obtained if the energy after correction of the linear error (1/L) is plotted against 1/L^3 to extrapolates results manually for L->infinite. 这个多级矩修正与超胞尺寸L的关系怎么理解呢？ 完全蒙。请过路高手跟帖解惑或者提供参考资料。 [ Last edited by cenwanglai on 2010-10-22 at 19:45 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 缺陷稳定性分析 |
» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
感觉是下周放榜了
已经有8人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有81人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 进口色谱纯级的甲醇、乙腈有没有可能带有很多的杂离子，大家交流交流经验
已经有15人回复
- MS里怎么用原胞的坐标建结构
已经有6人回复
- “声子液体”优化热电性能
已经有9人回复
- gaussian 是如何进行部分优化的？
已经有5人回复
- 请教一个问题，优化结构键长过长怎么解决
已经有5人回复
- 结构-热分析时如何施加温度梯度载荷？
已经有6人回复
- 非水体系银/银离子电极的制作方法
已经有7人回复
- GaAs等异质结构太赫兹辐射材料
已经有6人回复
- 用Dmol优化收敛的结构，用vasp不收敛？
已经有9人回复
- 先进锂离子电池正极材料LiFePO4的研究进展
已经有237人回复
- 过渡元素离子都能形成配合物吗？
已经有1人回复
- 化学中可以通过什么作用来结合重金属离子？
已经有11人回复
- 请教用molpro打印氢分子离子的重叠积矩阵、单电子积分矩阵、偶极矩阵
已经有3人回复
- 【求助】平坦势能面下的结构优化
已经有4人回复
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
3楼2010-09-21 09:12:41
sunyang1988
至尊木虫 (著名写手)
- 1ST强帖: 1
- 应助: 24 (小学生)
- 贵宾: 0.984
- 金币: 9397.7
- 散金: 1593
- 红花: 23
- 沙发: 1
- 帖子: 1744
- 在线: 1255.8小时
- 虫号: 965316
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
5楼2010-09-25 17:59:11
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
6楼2010-09-25 18:58:12
sunyang1988
至尊木虫 (著名写手)
- 1ST强帖: 1
- 应助: 24 (小学生)
- 贵宾: 0.984
- 金币: 9397.7
- 散金: 1593
- 红花: 23
- 沙发: 1
- 帖子: 1744
- 在线: 1255.8小时
- 虫号: 965316
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
7楼2010-09-25 19:38:39
25