24小时热门版块排行榜    

查看: 2521  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
信息提示
请填处理意见