24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1342  |  回复: 15
当前主题已经存档。

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

[交流] 【讨论】关于symmetry(pg=?)命令的一个疑惑,向各位专家讨教

向各位大侠请教一个问题:
以前很少用symmetry(pg=?)这个命令。前天看到hairan专家说他用symmetry(pg=?)从来没有成功过。于是我的兴趣就被激发出来了。

用了这个命令我做了一些测试:
(1)我的师弟magic6做了一个测试。他的一个分子的输入文件是D2对称性的。然后我要求他加上symmetry(pg=c2)进行优化(加“%kjob L301 1”控制字段和“#p”命令字段以达到测试目的,后来发现#p根本没派上用场,呵呵)。
(2)我自己以前做的一个C3h对称性的分子的输入文件,加nosymm的,我自己再调出来。对比着看了看。
(3)我自己以前做的一个C2v对称性的分子的输入文件,没加nosymm,我自己再调出来。对比着看了看。

结果如下:

(1)D2对称性分子,加symmetry(pg=c2),在Framework这一栏中,显示如下:
Stoichiometry     C74H64N4O8   
Framework group   D2[C2(NCNCC.CCNCN),X(C68H64O8)]
Deg. of freedom   110
Full point group  D2


(2)C3h对称性分子,加nosymm,在Framework这一栏中,显示如下:
Stoichiometry     C60H84N30Pt6
Framework group   C3H[X(C60H84N30Pt6)]
Deg. of freedom   89
Full point group  C3H


(3)C2v对称性分子,没加nosymm,在Framework这一栏中,显示如下:
Stoichiometry     C92H16CoN8(2)
Framework group   C2V[C2(Co),SGV(C4N2),SGV'(C4N2),X(C84H16N4)]
Deg. of freedom   90
Full point group  C2V


三个文件对比看,有个问题就想问了:如果一个分子对称性较高,它会选择哪个轴来定低对称性的轴呢?举例来讲:如果一个D4d对称性的分子,我强制定为C2对称性进行优化,那么Gaussian将选择哪个C2轴来建立它的density matrix呢?——毕竟D4d分子有多个C2轴
还有个问题:在第一个测试中,“NCNCC.CCNCN”代表的是什么意思呢?
第三个问题:如果一个D4d对称性的分子,我强制定为C2对称性进行优化,能不能控制Gaussian用特定的对称轴来构建density matrix呢?
请诸位大侠不吝赐教


[ Last edited by zeoliters on 2009-11-13 at 18:29 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)

瞌睡木虫

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
loovfnd(金币+2,VIP+0):感谢参与! 10-8 18:21
引用回帖:
Originally posted by yjcmwgk at 2009-10-8 14:25:
向各位大侠请教一个问题:
以前很少用symmetry(pg=?)这个命令。前天看到hairan专家说他用symmetry(pg=?)从来没有成功过。于是我的兴趣就被激发出来了。

用了这个命令我做了一些测试:
(1 ...

版主误解
我说用这个关键词没有成功过是针对那天争论的出现特殊错误的体系说的,意思是用这个关键词对那个错误不管用。
其他体系,比如D4h对称性分子降低到D2h对称性计算,我算过,可以进行
不过也有个别体系存在使用了这个关键词,程序还是判断为原有对称性的情况。
我感觉程序选择对称轴肯定是选那条对称性最高的轴,比如把D4h降低对称性的时候,肯定选择的是那条D4的轴,而不会去选择其他的D2或者C2轴。
这可能也是为什么很多分子无法用这个关键词降低对称性计算的原因,比如把D4h降为C2,沿D4轴有很多种可能的C2模式,程序不知道该选哪一个,然后就干脆还用D4h来计算。
至于如何让程序锁定这个轴,没有专门研究过。不敢妄论。

[ Last edited by hairan on 2009-10-8 at 17:21 ]
据权威报纸调查,84%的青少年犯罪者有玩电子游戏的经历,故应该禁止经营电子游戏。据我们所知,100%的青少年犯罪者有穿鞋的经历,故制鞋厂应该停业。
2楼2009-10-08 17:08:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yytsnake

至尊木虫 (职业作家)

★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
loovfnd(金币+2,VIP+0):感谢参与! 10-8 18:21
yjcmwgk(金币+2,VIP+0):有用有用!你说的太好了。我经常作为斑竹的身份去教训别人“去看高斯说明书中关于……的部分”,这次你也该狠狠教训我“去看GaussianIOP中关于……的部分” 10-9 10:55
对点群不在行,对高斯的对称性计算更是一头雾水(正在看资料学习呢)。
当体系有多个子群时,应该用一个选项来指定要哪一个吧。资料显示:
By using a lower symmetry than the molecule has there might be several equal
symmetry planes possible to use as main symmetry plane.
IOP(2/20) is needed to define which one is the main symmetry plane.

IOp(2/19)
Largest allowed point group, as Hollerith string.

IOp(2/20)
Number (1-3 for X-Z) of axis to help specify which subgroup of the type specified in IOp(19) to use.

[ Last edited by yytsnake on 2009-10-8 at 17:45 ]
3楼2009-10-08 17:30:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

loovfnd

至尊木虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
这个symmetry(pg=?),我一般使用在,高对称性时有虚频,然后稍微降低一点对称性时优化。成功过
莫强求!
4楼2009-10-08 18:22:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yyx19840628

木虫 (著名写手)

鹰的眼睛,狼的耳朵,豹的速 ...

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
loovfnd(金币+1,VIP+0):感谢参与! 10-9 08:08
我一般优化都不用这个的,要不要紧啊?对结果有没有影响
人生的最大遗憾莫过于错误地坚持了不该坚持的,轻易地放弃了不该放弃的……
5楼2009-10-08 19:26:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

引用回帖:
Originally posted by yyx19840628 at 2009-10-8 19:26:
我一般优化都不用这个的,要不要紧啊?对结果有没有影响

没影响
6楼2009-10-08 19:42:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
loovfnd(金币+1,VIP+0):感谢参与! 10-9 08:08
yjcmwgk(金币+0,VIP+0):不是事事求助于黑箱,而是想探讨一下:如果让黑箱自己选择对称轴,它会怎么选 10-9 10:55
引用回帖:
Originally posted by loovfnd at 2009-10-8 18:22:
这个symmetry(pg=?),我一般使用在,高对称性时有虚频,然后稍微降低一点对称性时优化。成功过

同意。感觉这个功能用处不大。手动起来又可靠又不麻烦,何必事事都求助于黑箱呢。
Simpler, stronger, more tolerant
7楼2009-10-08 19:51:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)

瞌睡木虫

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
loovfnd(金币+1,VIP+0):感谢参与! 10-9 08:08
引用回帖:
Originally posted by yytsnake at 2009-10-8 17:30:
对点群不在行,对高斯的对称性计算更是一头雾水(正在看资料学习呢)。
当体系有多个子群时,应该用一个选项来指定要哪一个吧。资料显示:
By using a lower symmetry than the molecule has there might be several equal
symmetry planes possible to use as main symmetry plane.
IOP(2/20) is needed to define which one is the main symmetry plane.

IOp(2/19)
Largest allowed point group, as Hollerith string.

IOp(2/20)
Number (1-3 for X-Z) of axis to help specify which subgroup of the type specified in IOp(19) to use.

这个讲到点子上了,估计能解决版主的问题,应该多奖金币
据权威报纸调查,84%的青少年犯罪者有玩电子游戏的经历,故应该禁止经营电子游戏。据我们所知,100%的青少年犯罪者有穿鞋的经历,故制鞋厂应该停业。
8楼2009-10-08 20:48:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
loovfnd(金币+1,VIP+0):感谢参与! 10-9 08:08
yjcmwgk(金币+0,VIP+0):我知道。我只是想探讨一下:如果让黑箱自己选择对称轴,它会怎么选 10-10 09:13
如果非要用低点群,可以加上一个小位移降低对称性,或者换程序——如pc-gamess/gamess(us)
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
9楼2009-10-08 22:46:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

引用回帖:
Originally posted by neweroica at 2009-10-8 19:51:


同意。感觉这个功能用处不大。手动起来又可靠又不麻烦,何必事事都求助于黑箱呢。

不是事事求助于黑箱,而是想探讨一下:如果让黑箱自己选择对称轴,它会怎么选
10楼2009-10-09 08:59:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yjcmwgk 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程085601,270求调剂 +41 @ASDF1234 2026-04-08 44/2200 2026-04-14 09:23 by barlinike
[考研] 085404 298分求调剂 +11 呼啦呼啦呼呼呼 2026-04-10 12/600 2026-04-14 08:38 by wfj257
[基金申请] 2026 WR青拔 +3 liuchb715 2026-04-09 6/300 2026-04-13 18:40 by liuchb715
[考研] 291求调剂 +8 关忆北. 2026-04-11 8/400 2026-04-12 09:32 by 逆水乘风
[考研] 一志愿郑州大学 22408 305分求调剂 +5 安小满zzz 2026-04-08 5/250 2026-04-12 00:41 by 蓝云思雨
[考研] 药学305求调剂 +8 玛卡巴卡boom 2026-04-10 8/400 2026-04-12 00:07 by zhouwenxian
[考研] 材料工程日语考生求调剂 +7 0856?调剂 2026-04-10 7/350 2026-04-11 21:33 by 蓝云思雨
[考研] 291分调剂 +5 上岸小莹加油 2026-04-09 6/300 2026-04-11 21:06 by 逆水乘风
[考研] 调剂求助 +6 果然有我 2026-04-11 7/350 2026-04-11 16:22 by 明月此时有
[考研] 调剂 +5 文道星台 2026-04-11 5/250 2026-04-11 15:01 by 凯凯要变帅
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 368求调剂 +3 17385968172 2026-04-10 3/150 2026-04-10 20:12 by 电子奥德彪
[考研] 266求调剂 +29 阳阳哇塞 2026-04-07 29/1450 2026-04-10 16:20 by 高维春
[考研] 一志愿沪9,326生物学求相关专业调剂 +4 刘墨墨 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
[考研] 江苏大学 工科调剂 捡漏 +3 Evan_Liu 2026-04-09 5/250 2026-04-10 10:22 by Evan_Liu
[考研] 292求调剂 +9 笑笑袁 2026-04-09 9/450 2026-04-10 10:05 by LHGeng
[考研] 初试分332,一志愿报考西北工业大学, +11 故人?? 2026-04-09 11/550 2026-04-09 21:54 by JineShine
[考研] 调剂 +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
[考研] 275 求调剂 +8 Lei812514 2026-04-07 8/400 2026-04-08 12:46 by chemisry
信息提示
请填处理意见