24小时热门版块排行榜    

查看: 2179  |  回复: 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的回帖

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的回帖

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的回帖

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的回帖

lzn_t

木虫 (小有名气)

引用回帖:
: Originally posted by souledge at 2012-02-20 22:07:37:
不知道LZ对于“基本构型”的理解是怎样的,这个词汇很暧昧,很难量化表示。
如果说是比例尺一样放大缩小,但各原子之间的空间对称关系保持不变,可以用分数坐标固定;如果说要各原子的绝对距离都不变,就要迪卡尔 ...

见4L
6楼2012-02-21 12:55:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lzn_t

木虫 (小有名气)

引用回帖:
: Originally posted by yadong80 at 2012-02-20 23:16:44:
结构优化的原则应该是尊重实验结构,也就是优化后的结构与优化前的结构(实验结构)越接近越好,如果出入太大,就得检查优化的设置,或则程序本身了。如果真如楼主所说,结构优化能够产生这样大的结构调整,楼主应 ...

见5L
7楼2012-02-21 12:56:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yadong80

铜虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
zzy870720z(金币+1): 谢谢指教 2012-02-21 15:47:08
不太清楚你所说的a和b之间到底什么关系,相差多少,如果体现在晶胞参数上,楼主或许可以手动调节晶胞参数,计算势能曲线,得到一个相对能量较低的b的supercell
8楼2012-02-21 13:09:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lzn_t 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 255求调剂 +3 李嘉慧, 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 288求调剂 +14 王晓阳- 2026-03-09 19/950 2026-03-14 02:05 by JourneyLucky
[考研] 265求调剂 +9 小木虫085600 2026-03-09 12/600 2026-03-14 01:11 by JourneyLucky
[考研] 一志愿湖师大化学289求调剂 +6 XMCMM3.14159 2026-03-10 6/300 2026-03-14 00:28 by JourneyLucky
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 279求调剂 +3 抓着星星的女孩 2026-03-10 3/150 2026-03-13 23:47 by userper
[考研] 材料工程,326分,求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 328化工专硕求调剂 +4 。,。,。,。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 一志愿华中师范071000,325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考研] 08食品或轻工求调剂,本科发表3篇sci一区top论文,一志愿南师大食品科学与工程 +3 我是一个兵, 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考研] 纺织、生物、化学、材料相关专业招生了 +4 耶耶业 2026-03-09 7/350 2026-03-12 19:05 by Equinoxhua
[考研] 一志愿江南大学085701环境工程专硕总分287求调剂 +5 18266118446 2026-03-09 5/250 2026-03-11 16:51 by 2020015
信息提示
请填处理意见