24小时热门版块排行榜    

查看: 2400  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 335求调剂 +3 yuyu宇 2026-03-23 4/200 2026-03-23 19:03 by macy2011
[考研] 生物学一志愿985,分数349求调剂 +6 zxts12 2026-03-21 9/450 2026-03-23 18:37 by macy2011
[考研] 284求调剂 +6 Zhao anqi 2026-03-22 6/300 2026-03-23 09:23 by king123!
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 求调剂一志愿海大,0703化学学硕304分,有大创项目,四级已过 +6 幸运哩哩 2026-03-22 10/500 2026-03-22 20:10 by edmund7
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[考研] 318求调剂 +4 plum李子 2026-03-21 7/350 2026-03-22 14:17 by ColorlessPI
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[考研] 268求调剂 +9 简单点0 2026-03-17 9/450 2026-03-21 15:37 by lature00
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
信息提示
请填处理意见