| 查看: 1315 | 回复: 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 ] |
回复此楼» 猜你喜欢
寻博导
已经有13人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有9人回复
-
物理化学论文润色/翻译怎么收费?
已经有227人回复
-
哈尔滨工业大学（深圳）绿色化学团队招收2026年秋季入学博士生
已经有11人回复
-
澳大利亚昆士兰科技大学PhD博士奖学金
已经有0人回复
-
【211博士招生】环境化学、地学、毒理方向，擅长污染物迁移转化降解、计算模拟者优先
已经有16人回复
-
求助如何利用vasp计算磁光克尔效应
已经有1人回复
-
070300总分298求调剂
已经有1人回复
-
一志愿厦门大学化学学硕307求调剂
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+4,VIP+0):感谢ganlh老师的回答。其实，我一直是“手动改变原子位置”这一流派的人 10-10 09:12
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+4,VIP+0):感谢ganlh老师的回答。其实，我一直是“手动改变原子位置”这一流派的人 10-10 09:12
 本帖仅楼主可见 |
11楼2009-10-10 08:40:04
hairan
木虫 (著名写手)
瞌睡木虫
- QC强帖: 3
- 应助: 61 (初中生)
- 贵宾: 1.06
- 金币: 4974.3
- 散金: 8072
- 红花: 44
- 帖子: 2901
- 在线: 733.4小时
- 虫号: 432069
- 注册: 2007-08-12
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
loovfnd(金币+2,VIP+0):感谢参与！ 10-8 18:21
小木虫(金币+0.5):给个红包，谢谢回帖交流
loovfnd(金币+2,VIP+0):感谢参与！ 10-8 18:21
|
版主误解 我说用这个关键词没有成功过是针对那天争论的出现特殊错误的体系说的，意思是用这个关键词对那个错误不管用。 其他体系，比如D4h对称性分子降低到D2h对称性计算，我算过，可以进行 不过也有个别体系存在使用了这个关键词，程序还是判断为原有对称性的情况。 我感觉程序选择对称轴肯定是选那条对称性最高的轴，比如把D4h降低对称性的时候，肯定选择的是那条D4的轴，而不会去选择其他的D2或者C2轴。 这可能也是为什么很多分子无法用这个关键词降低对称性计算的原因，比如把D4h降为C2，沿D4轴有很多种可能的C2模式，程序不知道该选哪一个，然后就干脆还用D4h来计算。 至于如何让程序锁定这个轴，没有专门研究过。不敢妄论。 [ Last edited by hairan on 2009-10-8 at 17:21 ] |
2楼2009-10-08 17:08:24
yytsnake
至尊木虫 (职业作家)
- 应助: 357 (硕士)
- 贵宾: 0.01
- 金币: 18204.3
- 散金: 5
- 红花: 36
- 帖子: 3189
- 在线: 444.6小时
- 虫号: 92066
- 注册: 2005-09-08
- 性别: GG
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
loovfnd(金币+2,VIP+0):感谢参与！ 10-8 18:21
yjcmwgk(金币+2,VIP+0):有用有用！你说的太好了。我经常作为斑竹的身份去教训别人“去看高斯说明书中关于……的部分”，这次你也该狠狠教训我“去看GaussianIOP中关于……的部分” 10-9 10:55
小木虫(金币+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
loovfnd
至尊木虫 (著名写手)
- 应助: 35 (小学生)
- 贵宾: 1.111
- 金币: 14361.4
- 散金: 236
- 红花: 4
- 帖子: 2790
- 在线: 187小时
- 虫号: 219674
- 注册: 2006-03-14
- 专业: 物理有机化学
4楼2009-10-08 18:22:37
