应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 如何让计算离子从晶体内部扩散到表面的能量？
152/2返回列表上一页12
查看: 2810  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

【答案】应助回帖

★ ★ ★ ★ ★
franch: 金币+5, 应助指数+1, 谢谢回帖交流。。。 2014-03-31 22:26:03
引用回帖:
10楼: Originally posted by 1164028932 at 2014-03-30 14:16:14
麻烦你看一下，我根据你的意思建立的反应物与生成物的结构如下图所示，你看一下合理吗?对不对？还有就是优化的时候要固定晶格优化吗?优化时固定其他原子只优化表面原子还是都优化？

reaction.png

production.png
...

OB轴向取的太少了。建议把OB轴向取得再大一些，否则相当于OB向的一条线上的Ba都脱出，这不合理。另外真空层建议至少10A，如果严格的话建议做测试，CASTEP的计算量和晶格大小有关，所以要看能承受多大的计算量了，如果是Dmol3就无所谓，真空层取个三五十A问题都不大。

优化应当作不opt cell（但是因此也有可能造成不收敛）的全几何优化，不要固定，这样拖出前后晶格弛豫的能量会包含在脱出的能量中，也很容易理解，离子脱出后晶格需要重组、再构。不过也有可能会不容易收敛。

不过我估计CASTEP计算这个体系有点危险，速度相当慢（我估计单节点每个structure都要几个小时），因为大致上从尺寸上看K-point基本是gamma了，这样CASTEP就算跨节点也不会有好的效率（CASTEP计算划分是依据K-point的）。速度如果比较慢就不要用Ultrafine，可以降低到fine。另外可以尝试Dmol3（也有可能收敛造成问题）甚至开源的SIESTA或Openmx。差不多的精确程度来看，VASP速度比CASTEP好得多，这个你就要好好请教楼上前辈了。

老菜鸟的一点建议，你先试试吧。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
11楼2014-03-30 16:02:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1164028932

金虫 (初入文坛)
引用回帖:
11楼: Originally posted by 卡开发发 at 2014-03-30 16:02:13
OB轴向取的太少了。建议把OB轴向取得再大一些，否则相当于OB向的一条线上的Ba都脱出，这不合理。另外真空层建议至少10A，如果严格的话建议做测试，CASTEP的计算量和晶格大小有关，所以要看能承受多大的计算量了，如 ...

恩好的，不要谦虚，在我眼里你已经是大神了，我还是请教一下，我这个能用Dmol算过渡态吗?我前几天就是这个结构分别用CASTEP和Dmol计算了同一个路径的扩散能，参数设置基本一样，为什么结果差距很大呀？我看文献没发现有用Dmol算晶体内部离子扩散的?麻烦大神指导一下。谢谢。
赞一下 回复此楼
12楼2014-03-31 21:28:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
12楼: Originally posted by 1164028932 at 2014-03-31 21:28:52
恩好的，不要谦虚，在我眼里你已经是大神了，我还是请教一下，我这个能用Dmol算过渡态吗?我前几天就是这个结构分别用CASTEP和Dmol计算了同一个路径的扩散能，参数设置基本一样，为什么结果差距很大呀？我看文献没发 ...

原子轨道的方法计算固体是不太准确，优化晶格结构都有点困难，castep对于PBC体系肯定好些的。你指的结果差距指的是扩散能垒？除了平面波方法中的pseudopot不一样，另外k-point和Ecut一般需要测试我们是知道的，Dmol3中的k-point以及原子轨道的截断也是要测试的，这些都不一样结果不一样也很正常。不过我觉得你这个体系偏大，CASTEP算起来有点危险，完全准确不可能，趋势正确，计算参数有力有据即可。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
13楼2014-04-01 09:37:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1164028932

金虫 (初入文坛)
引用回帖:
13楼: Originally posted by 卡开发发 at 2014-04-01 09:37:04
原子轨道的方法计算固体是不太准确，优化晶格结构都有点困难，castep对于PBC体系肯定好些的。你指的结果差距指的是扩散能垒？除了平面波方法中的pseudopot不一样，另外k-point和Ecut一般需要测试我们是知道的，Dmo ...

恩就是的，扩散能垒差距很大，你说的PBC体系是什么体系？
赞一下 回复此楼
14楼2014-04-02 22:35:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
14楼: Originally posted by 1164028932 at 2014-04-02 22:35:31
恩就是的，扩散能垒差距很大，你说的PBC体系是什么体系？...

周期边界条件（periodic boundary condition）。我觉得你先看看趋势（和试验方面），比较纯粹的绝对值意义不大，这个本身也依赖泛函的选取，很显然的就是即便同一种泛函采用不同基函数展开，特性也很不同。PBC下Dmol3计算某些性质我听说不是很靠谱。

castep大家很注意Ecut的选取，但是Dmol3大家对Rcut的选取似乎并不是很重视，两者的可比性我觉得都只有这些相关的参数极限外推才能够作比较。

另外，我不确定的是你现在过渡态构型上两者做出来的定性是否正确？Dmol3可以在TS search后计算Freq，然后进行TS opt，经过优化后的构型是否比较过？

靠谱点的话还是VASP。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
15楼2014-04-02 23:00:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1164028932 的主题更新
152/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 299求调剂 +5 ZVVZ13 2026-04-08 5/250 2026-04-09 10:13 by 逆水乘风
[考研] 0703调剂 +17 拾玖壹 2026-04-04 19/950 2026-04-09 09:18 by may_新宇
[基金申请] 有爆料，一个青年教师卖房得400万，然后换了一个四青帽子 +6 babu2015 2026-04-08 6/300 2026-04-09 08:14 by kudofaye
[考研] 308求调剂 +16 倘若起风了呢 2026-04-05 16/800 2026-04-08 23:03 by 猪会飞
[考研] 专硕310求调剂 +7 捞捞我…. 2026-04-04 8/400 2026-04-08 20:24 by yutian743
[考研] 085600材料与化工301分求调剂院校 +27 刺痛jk 2026-04-06 28/1400 2026-04-08 16:16 by luoyongfeng
[考研] 312 化工或制药调剂 +9 小小墨123 2026-04-02 10/500 2026-04-08 10:30 by screening
[考研] 315求调剂 +17 小羊小羊_ 2026-04-02 18/900 2026-04-07 22:01 by lijunpoly
[考研] 22408 318分求调剂 +4 勤奋的小笼包 2026-04-06 6/300 2026-04-07 15:05 by 纸鹤555
[论文投稿] Decision: Revise for Editor还会送审吗 100+3 CccccccccFD 2026-04-04 5/250 2026-04-07 10:58 by 北京莱茵润色
[考研] 一志愿华中农业大学0710（A）初试329分 求调剂 +5 一名26考研生 2026-04-04 5/250 2026-04-07 08:54 by 18828373951
[考研] 生物学调剂 可调剂到生物与医药 +3 李政莹 2026-04-06 3/150 2026-04-06 19:02 by macy2011
[考研] （调剂）一志愿报考哈尔滨工业大学0857资源与环境专业378分考生 +7 狠狠加油 2026-04-05 8/400 2026-04-06 16:52 by momo皓
[考研] 0857大类环境工程B区求调剂 +3 龚禹铭 2026-04-05 3/150 2026-04-06 10:22 by 蓝云思雨
[考研] 301求调剂 +6 细胞相关蛋白 2026-04-02 10/500 2026-04-06 08:34 by jp9609
[考研] 考研生物学考A区211，初试322，科目生化和生物综合，求调剂 +6 。。。54 2026-04-03 6/300 2026-04-05 14:54 by JOKER0401
[考研] 286求调剂 +3 草木不言 2026-04-04 3/150 2026-04-04 22:40 by lbsjt
[考研] 323求调剂 +8 李佳乐1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
[考研] [调剂信息]085408光电信息 求调剂 总分291分数一英一 +3 iz11az 2026-04-02 3/150 2026-04-04 19:09 by 蓝云思雨
[考研] 279求调剂 +5 傅文秋 2026-04-02 5/250 2026-04-02 18:10 by 笔落锦州
信息提示
关闭
请填处理意见
关闭