24小时热门版块排行榜     意得论文润色满3千送无线充电宝

查看: 100  |  回复: 2

quantumor

金虫 (著名写手)

快乐兔子

[求助] 部分优化的原子固定问题已有1人参与

Gaussian的部分优化(利用readfreeze)中,发现有时候被固定的原子在优化过程中也在跑,请问有何解决办法?
据说是用GaussView的Atom List Editor中的Freeze设置数值也可以固定原子,但我没用过,在Freeze下有0,-1,-2,-3四个值,请问这四个值的含义是什么?该死的说明书被我翻遍了也没找数值的含义。
回复此楼
愿好运与快乐伴随你!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zhou2009: 金币+5 2019-09-12 15:58:35
quantumor(小红豆代发): 金币+30, 回帖应助 2019-09-12 17:14:56
在写分子的坐标部分的时候 在原子的名字和原子的坐标之间加上-1 就可以冻结该原子。而标记为 0 则为普通优化的原子。 而 -2 -3代表了当你的优化有不同的片段时候,片段内部也就是具有同样负数数值的原子是彼此保持相对位置不变的。但是标记为-2 -3的不同片段之间的相对位置是可以在优化中变化的。但是这个用法只能在做ONIOM计算时候用。一般时候不能用。

而且需要注意的是,高斯一开始会默认根据你分子的质心位置到原点 并且根据对称性 转动你的分子到一个合适的角度(standard orientation)。这时候你的坐标会发生变化。所以只看数值的话 坐标有时候会变动。但是优化的时候相对位置还是冻结的。 如果需要冻结的原子连初始坐标一点都不可以变 那么就要加入 nosymm关键词。 这样高斯就不会帮你 reorientation。

举个例子 比如下面是一个水分子 想要冻结O原子

%chk=E:\G16W1\h2o_freeze2.chk
# b3lyp/6-31g opt nosymm

Title Card Required

0 1
O         -1        -2.38259685   -1.24309390    0.00000000
H          0       -1.42259685   -1.24309390    0.00000000
H          0       -2.70305143   -0.33815807    0.00000000

结果优化算完之后的坐标是

O               -1   -2.38259700   -1.24309400    0.00000000
H                0   -1.40687400   -1.22715400    0.00000000
H                0   -2.70390200   -0.32165300    0.00000000

而如果不加nosymm关键词

算好之后的坐标是
O               -1    0.00000000    0.00000000    0.11424600
H                0    0.00000000    0.79120400   -0.45698500
H                0    0.00000000   -0.79120400   -0.45698500

这俩优化好的结构如果检查键长健角什么的其实是一样的。 只是后者在优化开始的时候被平移和转动了。
2楼2019-09-12 04:18:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quantumor

金虫 (著名写手)

快乐兔子

感谢版主详细的回复!谢谢!具体该如何操作才能将30金币交给你?请明示,我很久不来木虫,不记得该如何操作了,谢谢!
愿好运与快乐伴随你!
3楼2019-09-12 08:28:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 quantumor 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 大学当辅导员好还是高中老师好 +44 妹妹kln 2019-09-13 46/2300 2019-09-19 12:48 by 寒夜微曙
[教师之家] 中年学术危机不约而至~ +12 小纵队长 2019-09-17 24/1200 2019-09-19 12:32 by 向宽处行
[硕博家园] 第一次去招聘会,感受了社会的残酷,优秀的人那么多,找工作好难。 +11 Actor~1 2019-09-18 12/600 2019-09-19 11:37 by liuhao_FOERC
[论文投稿] 论文今天下午终于能够被检索了,哈哈,毕业有望,开帖纪念一下!!! +19 孤鹰大漠 2019-09-18 23/1150 2019-09-19 10:28 by 孤鹰大漠
[教师之家] 评职称老人老办法新人新办法是否不合理 +21 长江学者1 2019-09-16 25/1250 2019-09-19 10:26 by Zed-jie
[找工作] 奉劝没有家底的虫友不要来杭州高校,真实工资很低 +29 gx木木重 2019-09-14 32/1600 2019-09-19 09:58 by 大白狼
[材料综合] 推荐几个检测机构 +6 艾尚雪 2019-09-18 6/300 2019-09-19 09:06 by happyever83
[考博] 清华水利 +6 箕山 2019-09-15 10/500 2019-09-18 13:20 by 箕山
[有机交流] 求助下面的一个反应 10+5 muxiaomuku 2019-09-15 12/600 2019-09-18 11:17 by 88817753
[考博] 啊,终于又等到了一个消息 +4 城口的木匠 2019-09-15 8/400 2019-09-18 08:32 by 城口的木匠
[有机交流] 噻吩环上醛基 10+4 北国?? 2019-09-13 14/700 2019-09-17 20:37 by xinxin91730
[考博] 想读博士,求推荐导师 +8 党叮当 2019-09-17 11/550 2019-09-17 18:38 by BruceChum
[考博] 四级没过,化学,一篇一区,一篇二区sci,老师同意报考,学校审核能过么? +34 杨宇丫头 2019-09-14 62/3100 2019-09-17 09:49 by 希望大哥
[基金申请] 2020年山东省自然科学基金 +3 suger舒贝特 2019-09-16 4/200 2019-09-17 06:19 by dkljl
[分子生物] 求GFP 质粒 +4 这也行ww 2019-09-12 5/250 2019-09-17 00:05 by donkey2017
[论文投稿] 论文二次被拒 +27 若初见~ 2019-09-12 31/1550 2019-09-16 13:24 by 相视朦胧中
[论文投稿] SCI修改稿和答复意见里出现了几处明显错误,是否会影响最终论文录用 10+4 lijin2014 2019-09-15 17/850 2019-09-15 23:00 by lijin2014
[考博] 你们打算报几个学校? +5 WJSza 2019-09-15 27/1350 2019-09-15 15:41 by WJSza
[硕博家园] 坐标成都 +11 爱运动的女孩 2019-09-15 11/550 2019-09-15 10:11 by 踏雪飞花
[论文投稿] CEJ这个期刊现在什么情况?有人投过吗 +10 肿瘤靶向材料 2019-09-14 12/600 2019-09-15 09:44 by 肿瘤靶向材料
信息提示
请填处理意见