24小时热门版块排行榜    

查看: 2795  |  回复: 9

vigaryang

木虫 (正式写手)

[求助] 关于结构优化中的对称性问题

如果要优化一个像NH3这种明确知道对称性的分子,是不是在优化时需要事先指定优化过程中对称性不变?如果是的话,是通过symmetry(pg=?)实现呢,还是在分子内坐标输入部分设置变量实现?网上查了一下,发现很难见到与symmetry(pg=?)这个命令有关的信息,不知是为什么。难道对于这种分子也不需要指定对称性吗,而任由其自由优化吗。谢谢。
回复此楼

» 收录本帖的淘帖专辑推荐

量化程序使用方法——Caussian,NBO,Mult,S

» 猜你喜欢

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

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhficcas

银虫 (正式写手)

【答案】应助回帖


gmy1990(金币+1): 非常感谢你对论坛的贡献,奖励一下! 2011-05-23 08:43:08
vigaryang(金币+2): 2011-05-23 21:48:17
通过内坐标来控制变量, 这个NH3可能还要引入虚原子来精确定义二面角.
2楼2011-05-23 06:20:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vigaryang

木虫 (正式写手)

引用回帖:
Originally posted by zhficcas at 2011-05-23 06:20:29:
通过内坐标来控制变量, 这个NH3可能还要引入虚原子来精确定义二面角.

为什么不用symmetry(pg=?)呢?
3楼2011-05-23 06:47:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小白如水

铁杆木虫 (小有名气)

关于对称性方面也是我迷惑的 关注一下。
为了消除虚频有种尝试是忽略对称性(nosymm)这是不是以降低计算精度为代价呢?
4楼2011-05-23 09:59:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

【答案】应助回帖

★ ★ ★
gmy1990(金币+3): 非常感谢你对论坛的贡献,奖励一下! 2011-05-23 12:26:39
vigaryang(金币+5): 2011-05-23 21:48:42
vigaryang(金币+3): 非常感谢 2011-05-25 14:49:08
Gaussian在优化的过程中,默认保持对称性不变。
也就是说,你的输入结构在gaussian读入之后,会首先做一个对称性判断,如果它判定为C3v了,那么,优化过程中,他就一直按照C3v来优化,如果按C3v优化不出来,他就按出错退出——一直到g03都是如此。。。但是最近发现g09似乎有所改变了,他会尝试改变对称性的,这一点我还没弄清楚他的机制。
不一定是按照内坐标输入,即便你按照笛卡尔坐标输入也是可以的,当然,这个时候gaussian在判断对称性的时候,是会因为你的笛卡尔坐标的精度不够而容易产生误判的,就好象gview在判断对称性的时候,也会有Default,tight,very tight, loose,very loose等选项一样,gaussian在判读对称性的时候也有tight.loose等标准,你可以自己指定,当然,它的标准与gview的标准是不同的。你的分子在gaussian读入后,开始跑202的时候就会给出它的判断结果——不过我猜测在跑103的时候就应该已经算出来了。如果你对gaussian的判读没有把握,或者对你输入的笛卡尔坐标的精度没有把握,那么在此时kill掉gaussian,看看输出文件就知道它怎么判定的了。
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
5楼2011-05-23 12:12:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vigaryang

木虫 (正式写手)

引用回帖:
Originally posted by zhangmt at 2011-05-23 12:12:42:
Gaussian在优化的过程中,默认保持对称性不变。
也就是说,你的输入结构在gaussian读入之后,会首先做一个对称性判断,如果它判定为C3v了,那么,优化过程中,他就一直按照C3v来优化,如果按C3v优化不出来,他就 ...

学习了。我试了下添加symmetry(pg=?),发现对称性在计算前后反而比不添加时发生了更大的变化,诧异。

还有,如果一开始被gaussian判定为C3V的话,是不是对于优化过程中的每一步构型,其对称都是C3V?在我的某一个计算中,初始构型和结束构型都是Oh对称(在gaussview中搭建初始分子时通过point group菜单控制了对称性),但我在.log文件中搜索point group关键词时,发现中间这些构型的point group都是C1。不知这正常不。问题有点多,麻烦您了。
6楼2011-05-25 11:25:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

引用回帖:
在我的某一个计算中,初始构型和结束构型都是Oh对称(在gaussview中搭建初始分子时通过point group菜单控制了对称性),但我在.log文件中搜索point group关键词时,发现中间这些构型的point group都是C1。不知这正常不。

如果你在log文件中搜索point group,看到gaussian判定为C1,那么它就按C1优化,也就是没有对称性。原因有二:
1、gaussian对于对称性的判定是比较笨的,尤其Oh它经常是判定不出来的,这个应该是个bug,不知道g09改善了没有;
2、就是判定标准的问题,如果不用内坐标的形式,那么即便gview给出正确的对称性,但gaussian的标准与gview不同,也可能给出不同的结果。通常是由于笛卡尔坐标的精度有限导致的。
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
7楼2011-05-25 14:32:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sangym

木虫 (著名写手)

【答案】应助回帖

在优化之前应该先调好对称性,调到能达到最高的对称性,这样优化起来会节省时间,而且出来的结果也更合理,呵呵!!!!!
8楼2011-05-25 14:48:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vigaryang

木虫 (正式写手)

引用回帖:
Originally posted by sangym at 2011-05-25 14:48:56:
在优化之前应该先调好对称性,调到能达到最高的对称性,这样优化起来会节省时间,而且出来的结果也更合理,呵呵!!!!!

那为什么很多时候在计算中会在命令行中特意加入nosymmetry呢?
9楼2011-05-25 15:38:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sangym

木虫 (著名写手)

引用回帖:
Originally posted by vigaryang at 2011-05-25 15:38:22:
那为什么很多时候在计算中会在命令行中特意加入nosymmetry呢?

是这样的,我的理解不管你的构型有没有对称性,你如果不加对称性的话,程序会自己找对称性,如果你的体系没有对称性你直接加nonsymmetry,程序就不用自己找,直接会默认没有对称性,这样节省时间!!!!这是我理解,如果有更好更详细的答案欢迎交流,呵呵!!祝好运
10楼2011-05-25 19:07:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 vigaryang 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 0703化学调剂 ,六级已过,有科研经历 +4 曦熙兮 2026-03-15 4/200 2026-03-15 18:01 by JourneyLucky
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 311求调剂 +3 26研0 2026-03-15 3/150 2026-03-15 09:12 by JourneyLucky
[考研] 一志愿天大化工(085600)调剂总分338 +6 蔡大美女 2026-03-09 6/300 2026-03-14 02:46 by JourneyLucky
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-09 6/300 2026-03-14 02:14 by JourneyLucky
[考研] 一志愿天津大学,英一数二305分求调剂,四六级已过 +8 小小番的茄 2026-03-09 8/400 2026-03-14 01:53 by JourneyLucky
[考研] 材料工程,326分,求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 材料调剂,307分 +13 张泳铭1 2026-03-09 17/850 2026-03-13 11:09 by 薛云鹏
[考研] 一志愿华中师范071000,325求调剂 +5 RuitingC 2026-03-12 5/250 2026-03-13 10:43 by hyswxzs
[考研] 296求调剂 +3 大口吃饭 身体健 2026-03-13 3/150 2026-03-13 10:31 by 学员8dgXkO
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
信息提示
请填处理意见