24小时热门版块排行榜    

查看: 2817  |  回复: 6

forkn

银虫 (小有名气)

[交流] 【求助】如何利用分子动力学加速DFT计算的收敛 已有5人参与

各位高手,利用castep对过渡金属体系加自旋进行结构优化,由于是自己建立的模型,没有相关的文献可供参考,因此初始结构可能不够合理,SCF很难收敛。大多数时间做电子最小化时第一步就跳出了。
我已经参考过别人的意见,先不加自旋几何优化,在此结构上再算加自旋优化。现在有3个问题:

1 好像听人说过,第一性原理计算指出,加不加自旋对对体系的结构变化不大,是这样吗?如果是这样的话,那理应在不加自旋优化出的结构上加自旋进一步优化是比较容易收敛的啊?因为平衡结构很接近了啊,但是我这样算的时候,加上自旋还是很难收敛?

2 加自旋几何优化时,自旋的设置对收敛有没有大的影响?是不是自旋设置的不合适导致了不收敛?我在modify-spin里设置的自旋,formal charge 选择默认,只是对formal spin 选项进行设置,分别选了higt spin ,low spin ,其他默认。计算时也对 use formal spin as initial 打钩了,这样设置正确吗?

3 如果是初始结构不合理,那么不加自旋几何优化也会很难收敛,听一些专家和版主说过最好利用分子动力学先对初始结构进行优化,再拿到castep中优化。像过渡金属体系,还是磁性的(Fe/Ni体系),用ms中的哪个模块比较合适?如何用分子动力学的相关模块,设置参数都有哪些啊?好像需要经验参数,力场什么的?偶对分子动力学一点不懂啊,如何上手?
回复此楼
山阻石拦,大江毕竟东流去;雪辱霜欺,梅花依旧向阳开!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yb2006

铁杆木虫 (正式写手)

★ ★
forkn(金币+2):谢谢你的建议 2010-06-23 19:55:58
qasd(金币+2):谢谢 2010-06-23 20:27:22
forkn(金币+3):谢谢 2010-07-05 21:25:34
我是用VASP的,不太知道castep中参数怎么设置,但一般我在算比较复杂的体系的时候都是先结构优化跑几十步,然后用MD跑个几千步,这样的话比较容易收敛,而且能使结构跳出local minimum,找到真正的稳定结构~~你可以用siesta这种比较快的程序包来跑MD,然后把结构导入castep中进行优化~~我也不太了解castep,帮你顶了!
2楼2010-06-23 18:51:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forkn

银虫 (小有名气)

引用回帖:
Originally posted by yb2006 at 2010-06-23 18:51:36:
我是用VASP的,不太知道castep中参数怎么设置,但一般我在算比较复杂的体系的时候都是先结构优化跑几十步,然后用MD跑个几千步,这样的话比较容易收敛,而且能使结构跳出local minimum,找到真正的稳定结构~~你可 ...

谢谢你的经验,那为什么先优化跑几十步再MD跑而不直接MD跑啊?呵呵,不懂啊!另外,siesta这个软件好用吗?我没有任何基础,自己当干,不知容易上手不?
山阻石拦,大江毕竟东流去;雪辱霜欺,梅花依旧向阳开!
3楼2010-06-23 20:02:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher


hedaors(金币+1):感谢交流 2010-06-24 07:53:25
forkn(金币+1):你指的是直接用dmol加自旋优化吗?还是指先用dmol对初始结构优化 2010-06-24 15:14:42
引用回帖:
Originally posted by forkn at 2010-06-23 20:02:33:

谢谢你的经验,那为什么先优化跑几十步再MD跑而不直接MD跑啊?呵呵,不懂啊!另外,siesta这个软件好用吗?我没有任何基础,自己当干,不知容易上手不?

结构不大用dmol就行了。
4楼2010-06-24 07:22:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

宏hong

木虫 (正式写手)


gavinliu7390(金币+1):谢谢交流! 2010-06-24 10:05:18
castep也可以算的
5楼2010-06-24 09:37:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

condensed

木虫 (著名写手)

★ ★
gavinliu7390(金币+2):谢谢交流! 2010-06-24 11:05:35
forkn(金币+2):谢谢,可惜我没有vasp的程序啊,也不会用,还是想要castep,学校有这个软件的版权 2010-06-24 15:16:21
没有对自旋的经验,不过有对复杂体系优化的经验.如果你用vasp优化可以分为几步来优化.比如首先固定晶格,光优化原子位置.然后固定原子位置优化晶格.最后在一起优化.这样比一开始就开放所有的自由度要快得多.
6楼2010-06-24 10:54:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forkn

银虫 (小有名气)

我的意思是想知道过渡金属磁性体系的初始结构用MS中的哪个分子动力学模块能给先优化下,然后在castep中跑,这样比较容易收敛,也同时易于找到真正能量最低的结构,而不是local minimum。或是其他能够快速几何优化的程序包也行。
山阻石拦,大江毕竟东流去;雪辱霜欺,梅花依旧向阳开!
7楼2010-06-24 15:22:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 forkn 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +31 医学老男孩 2026-08-13 69/3450 2026-08-16 00:16 by mo_mo_mo_mo
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见