24小时热门版块排行榜    

查看: 2182  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lzn_t

木虫 (小有名气)

[求助] 怎样在DMol3 结构优化中使构型基本保持不变?

现已知同一物质存在两种不同结构a和b,已知晶胞(supercell)在不带电荷时a结构更为稳定;带有+1电荷时b更稳定。问题是:b结构的supercell在电中性下优化所得结果会变为a结构的supercell,但我需要得到b构型在电中性时的最优supercell,该如何实现?

[ Last edited by lzn_t on 2012-2-20 at 17:37 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yadong80

铜虫 (小有名气)

WDD880227(金币+1): 感谢交流 2012-02-21 08:34:25
结构优化的原则应该是尊重实验结构,也就是优化后的结构与优化前的结构(实验结构)越接近越好,如果出入太大,就得检查优化的设置,或则程序本身了。如果真如楼主所说,结构优化能够产生这样大的结构调整,楼主应该检查一下设置了。
PS: 个人觉得,楼主的描述是在考察我们的思辨能力,a和b是同一物质,已经说明a的supercell在电中性时最稳定,b的supercell在带一个单位电荷时最稳定,那么在电中性条件下优化b的supercell,得到a的supercell有什么好奇怪的?我想程序想告诉你a的supercell就是电中性条件下b的最优supercell
楼主最好说的具体点,a和b的差异是什么?
4楼2012-02-20 23:16:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

1006094470

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
WDD880227(金币+1): 感谢提示 2012-02-21 08:33:54
把所有原子选上,然后modify-constraints可不可以?
2楼2012-02-20 21:51:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
WDD880227(金币+1): 感谢提示 2012-02-21 08:34:11
不知道LZ对于“基本构型”的理解是怎样的,这个词汇很暧昧,很难量化表示。
如果说是比例尺一样放大缩小,但各原子之间的空间对称关系保持不变,可以用分数坐标固定;如果说要各原子的绝对距离都不变,就要迪卡尔坐标固定。
相应的,还有在某轴向上固定其他轴向上变动之类的。
如果LZ想要的是,可以稍微构造变化一点,但是不要变化太多。比如正方体变成底比高为100:99的四方也可以接受,但是如果变成了不一定什么样子的结构就不能接受的话,那LZ就对电脑的智能期待太高了,或者,LZ的行为已经类似于“向着自己想要的结果上取优化”,从而与造结果没什么区别了……
思想重于技巧,内涵重于表象
3楼2012-02-20 22:07:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lzn_t

木虫 (小有名气)

是这样的,在电中性时两种结构均存在,虽然a结构的能量更低,也更稳定,但由于b结构向a结构转变时存在能量势垒,所以并非b结构在得到一个电子后一定会转变为a结构。本计算的初衷也是在建立电中性时两种结构的合理模型,然后寻找结构转变的能量路径(energy path of chemical reaction)。由于两种结构的supercell在建立时在某些细节之处没有可靠的数据可供使用,故而需先对所建立的模型结构优化,再计算其反应势垒。但现在问题是电中性时b优化结果即为a结构,无法将优化后的b作为reactant来对待。
    另外,已验证在带有+1电荷时,a结构优化后不会变为b结构。因此,估计在电中性时两者之间势垒很小,而带有+1电荷时则势垒稍高。尽管如此,依然想得到其具体的势垒值,但怎样得到电中性时合理的b结构呢?
5楼2012-02-21 11:09:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 293求调剂 +4 世界首富 2026-03-11 4/200 2026-03-17 09:23 by 雾散后相遇lc
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 材料专硕326求调剂 +5 墨煜姒莘 2026-03-15 5/250 2026-03-16 21:30 by 木瓜膏
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +5 Liwangman 2026-03-15 5/250 2026-03-16 17:10 by 我的船我的海
[考研] 0703化学调剂 ,六级已过,有科研经历 +7 曦熙兮 2026-03-15 7/350 2026-03-16 16:34 by houyaoxu
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 326求调剂 +4 上岸的小葡 2026-03-15 5/250 2026-03-16 08:39 by Linda Hu
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[考研] 255求调剂 +3 李嘉慧, 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考研] 279求调剂 +3 莫xiao 2026-03-10 4/200 2026-03-11 08:06 by 斩魂滴兔子!
信息提示
请填处理意见