24小时热门版块排行榜    

查看: 2409  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 306求0703调剂一志愿华中师范 +10 纸鱼ly 2026-03-21 11/550 2026-03-24 17:22 by qingfeng258
[考研] 080500求调剂 +3 zzzzfan 2026-03-24 3/150 2026-03-24 16:38 by barlinike
[考研] 321求调剂 +4 Ymlll 2026-03-24 4/200 2026-03-24 14:44 by sprinining
[考研] 277分求调剂,跨调材料 +3 考研调剂lxh 2026-03-24 3/150 2026-03-24 13:52 by JourneyLucky
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[考研] 材料专硕英一数二306 +8 z1z2z3879 2026-03-18 8/400 2026-03-23 20:49 by baobaoye
[考研] 306求调剂 +9 chuanzhu川烛 2026-03-18 9/450 2026-03-23 13:17 by luoyongfeng
[考研] 317求调剂 +12 申子申申 2026-03-19 18/900 2026-03-22 22:23 by luoyongfeng
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
信息提示
请填处理意见