24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2149  |  回复: 13

flowers123

银虫 (小有名气)

[求助] 求助:关于二聚体的结合能和nosymm关键字的问题

目前计算遇到困难,不知如何展开下去,问了实验室的师兄也都没有给出实用的建议,毕竟他们都没做过类似的题目。希望有遇到过的虫友给分析分析。目前在做二聚体的结合能问题如下;
1、二聚体具有严格的分子对称性,在计算二聚体的总能量时,发现在输入文件中加入nosymm比不加计算得到的能量大0.001hatree左右,计算扣除bsse的结合能时,结合能差0.3kcalmol左右。查了资料,nosymm是去掉分子轨道的对称性,不明白为什么对高对称性的二聚体总能量影响那么大?
2、单体的能量计算时,因为要加入别的东西,破坏了对称性,所以一直是加入nosymm的,不加就算不动。由加了nosymm的总能量减去单体的能量得到的结合能,跟文献资料相差0.002kcal\mol,是不是说明,为了保持一致,必须两个文件都要加入nosymm?
老板很忙,自己年纪一大把学点东西很难,希望各位虫友帮忙指导一下,我看了一下午小木虫上nosymm相关的内容,还是没有搞懂。希望有人能帮我解决,在此拜谢!!

[ Last edited by flowers123 on 2011-7-3 at 17:56 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

盼望长成一只小菜鸟......
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ieee18

捐助贵宾 (著名写手)

桥东会务


【答案】应助回帖

★ ★
uuv2010(金币-2): 请针对问题回帖 2011-07-03 20:04:07
ben_ladeng: 建议去学术会议发布栏发帖。希望您不要这样发布消息。 2011-07-03 20:04:59
内容已删除
2楼2011-07-03 19:28:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flowers123

银虫 (小有名气)

怎么没有人回答我呢
求助帖,非常焦急,无关的内容请不要贴进来,谢谢啦
盼望长成一只小菜鸟......
3楼2011-07-04 08:43:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤狼

金虫 (正式写手)

独行者

【答案】应助回帖

★ ★
flowers123(金币+5): 感谢! 2011-07-04 14:28:59
gmy1990(金币+2): 非常感谢你对论坛的贡献,奖励一下! 2011-07-05 11:37:34
加上利用对称性会极大的减少计算量,对于高分子有对称性还是利用对称性比较好(个人见解)~~但是对于某些体系的最稳定状态并不一定是对称性最高的状态,所以对称性的选择要注意下。另外这个能量差很大么?感觉用g03和g09算的能量误差有时候就比这个大了~~看看文献上用的方法和软件版本与你的是否有出入
轮回的过往,注定的未来;独孤的行者,孤独的路途……
4楼2011-07-04 10:24:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

★ ★
flowers123(金币+5): 感谢 2011-07-04 14:29:19
flowers123(金币+2): 追加,感谢回帖内容如此长 2011-07-04 15:01:45
gmy1990(金币+2): 非常感谢你对论坛的贡献,奖励一下! 2011-07-05 11:38:17
我的体会是,程序计算时默认进行对称性操作,其目的仅仅是利用对称性操作,可以对对称性相同的轨道之间的积分避免重复计算,从而减少计算量。

nosymm仅只是说程序对分子不进行对称性操作,不是说不承认其固有的对称性,无非是多一些重复计算而已。所以加不加nosymm影响不大。

当然由于加nosymm,没有强调先找出对称性,重复计算对称性相同的轨道之间的积分会多少有点计算误差。

nosymm只是在需要保持分子前后坐标保持不变时加,如作电子密度差。因为不加nosymm,程序会进行对称性操作,将分子改变坐标轴取向、挪动分子电子密度重心于坐标原点。

对于单个原子,要慎加nosymm。因为默认进行对称性操作,计算结果会将某一个单占或双占(孤对电子)轨道放在某一个坐标轴方向,有利于我们观察。
但加nosymm,因不计较对称性,这种单占或双占(孤对电子)轨道可能用多个坐标轴分量来表达,不利于我们观察。

[ Last edited by zhou2009 on 2011-7-4 at 10:46 ]
5楼2011-07-04 10:33:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flowers123

银虫 (小有名气)

引用回帖:
Originally posted by 独孤狼 at 2011-07-04 10:24:32:
加上利用对称性会极大的减少计算量,对于高分子有对称性还是利用对称性比较好(个人见解)~~但是对于某些体系的最稳定状态并不一定是对称性最高的状态,所以对称性的选择要注意下。另外这个能量差很大么?感觉用g ...

分子体系并不是很大,总共有三十个原子左右,对称性非常好。我认为这个差值很大的,因为结合能只有3kcal\mol ,差0.3kcal\mo就是10%,只是无论如何不能允许的。
我知道加上利用对称性会极大的减少计算量,但是加与不加相差比较大,不知道怎么选择。我原来一直认为,对称性这么好,应该加不加都一样的。请问你的意思是说,我的体系对称性低的状态更稳定么?
我为了重复文献上的计算结果,用的它结构数据,方法和基组都一样。只有计算包它是g03,我用g09,目前实验室计算系统没有g03。没办法直接重复计算。

另外阁下能不能指点下,g03和g09记算能量一般都差别这么大么?怎么调整g09的参数是更接近g03的结果?我目前调整了scf收敛精度,发现基本没有影响。
盼望长成一只小菜鸟......
6楼2011-07-04 14:57:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flowers123

银虫 (小有名气)

引用回帖:
Originally posted by zhou2009 at 2011-07-04 10:33:04:
我的体会是,程序计算时默认进行对称性操作,其目的仅仅是利用对称性操作,可以对对称性相同的轨道之间的积分避免重复计算,从而减少计算量。

nosymm仅只是说程序对分子不进行对称性操作,不是说不承认其固有的 ...

感谢斑竹大人回帖指导。
我也本来认为加了 nosymm,最多只是更昂贵些。使我的结构问题么?也不是啊,结构是用的文献上的数据,自己手写的z矩阵坐标,构建的模型,对称性是不容置疑的。可不可以认为,我的构型对称性是肯定的确定没有问题的,不用程序对它进行对称性操作,将分子改变坐标轴取向、挪动分子电子密度重心于坐标原点等等,所以要使用nosymm。那这样的话,加不加nosymm应该差别不大的!为啥我的计算结果差别如此巨大!!是在不能理解。

弱弱的问一句版主,我这样的数据差距也属于计算误差范围么?(⊙﹏⊙b汗,我觉得数值差距有点太大了)因为加了nosymm的结果更接近文献值,我是否可以认为,不加nosymm得到的总能量是有问题的?也就是说,我因为计算单体能量的时候使用了nosymm,为了保持一致性计算总能量是也要加呢?
在此拜谢
盼望长成一只小菜鸟......
7楼2011-07-04 15:35:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flowers123

银虫 (小有名气)

发出去就不能修改了,恳请懂行的各位指导我,金币不断追加中.....钱不是问题,虽然钱不多。
盼望长成一只小菜鸟......
8楼2011-07-04 15:41:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤狼

金虫 (正式写手)

独行者

引用回帖:
Originally posted by flowers123 at 2011-07-04 14:57:46:
分子体系并不是很大,总共有三十个原子左右,对称性非常好。我认为这个差值很大的,因为结合能只有3kcal\mol ,差0.3kcal\mo就是10%,只是无论如何不能允许的。
我知道加上利用对称性会极大的减少计算量,但是 ...

能把你的头文件贴一下么
轮回的过往,注定的未来;独孤的行者,孤独的路途……
9楼2011-07-04 18:22:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flowers123

银虫 (小有名气)

%mem=256mb
%nprocl=4
%nproc=2
#gfinput iop(6/7=3)
#p mp2/gen nosymm
我不是很明白 阁下说的 头文件是不是坐标前面的这部分。
今天下午我个饭gfinput输出了一下,dz基组的详细内容,发现g09里面的都在和g94差好多,这么说是不是g03中用dz也和g09的dz也会有差别,是个是不是也会引起计算结果不同呢?
盼望长成一只小菜鸟......
10楼2011-07-04 20:35:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 flowers123 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 295分求调剂 +6 ?要上岸? 2026-04-17 6/300 2026-04-21 08:18 by Equinoxhua
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[论文投稿] 期刊推荐 +3 材料研究生 2026-04-15 5/250 2026-04-20 16:02 by 豆豆7758
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +4 yexuqing 2026-04-19 4/200 2026-04-20 14:47 by brantleo
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 5/250 2026-04-20 10:47 by YuY66
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 297,工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 291求调剂 +11 关忆北. 2026-04-14 11/550 2026-04-19 17:16 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
信息提示
请填处理意见