24小时热门版块排行榜     石溪大学接受考研调剂申请>

【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 1877  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

quantumor

金虫 (著名写手)

快乐兔子

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

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

mujie2319

新虫 (小有名气)

请问一下楼主,我想用qe计算声子谱  该怎么固定原子呢
6楼2019-11-18 17:03:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

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的回帖

wzhsun

新虫 (小有名气)

引用回帖:
3楼: Originally posted by quantumor at 2019-09-12 08:28:08
感谢版主详细的回复!谢谢!具体该如何操作才能将30金币交给你?请明示,我很久不来木虫,不记得该如何操作了,谢谢!

人家解惑不是为了你的30金币 是乐于助人 不要那么物质 人性是好的

发自小木虫IOS客户端
4楼2019-10-17 12:07:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] LWT投 +3 AChen92 2024-04-26 3/150 2024-04-26 22:16 by hizifu
[基金申请] 基金开始函评了吗? +16 wych1103 2024-04-25 16/800 2024-04-26 21:32 by 淀粉搬运工
[考研] 0854-0855调剂 +8 shangannum1 2024-04-21 12/600 2024-04-26 16:42 by yz仔
[论文投稿] with editor 两个月了,什么原因? +7 yiersan9 2024-04-24 15/750 2024-04-26 16:13 by jonewore
[硕博家园] 考研,求职还是考编? +15 xizj 2024-04-21 24/1200 2024-04-26 11:49 by Kan客
[论文投稿] Chemical Engineering Journal投稿3周了,一直显示With editor状态。这是送审了吗? 10+4 yifeng11 2024-04-20 13/650 2024-04-26 09:48 by yifeng11
[基金申请] "颜宁:基础研究应顶天立地"能做到基础研究同时顶天立地的才是牛人 +5 zju2000 2024-04-24 5/250 2024-04-26 09:36 by LittleBush
[考博] 取博导收留 5+4 zzb777888 2024-04-20 10/500 2024-04-26 08:52 by polymerfriend
[论文投稿] Nature一直在编辑手里,考虑好几天了,是悬了吗 +12 彩虹初见 2024-04-24 12/600 2024-04-25 19:21 by 雪径踏青
[教师之家] 期末给学生划重点都是什么话术啊 +16 luokereng 2024-04-20 18/900 2024-04-25 15:46 by BusyGer
[考博] 求博导 +6 好okjh 2024-04-21 10/500 2024-04-25 14:04 by 好okjh
[基金申请] 国社科项目,你们学校都限额申报吗? +7 屡战屡败 2024-04-21 10/500 2024-04-25 12:10 by 屡战屡败
[访问学者] CSC的访问学者申请,没有个评审意见,也不知道怎么改,还有必要申请吗 +4 flyingship 2024-04-20 4/200 2024-04-24 21:09 by 59038mute
[考博] 博士招生 +4 zx179 2024-04-24 7/350 2024-04-24 15:01 by H考研成功
[教师之家] 大家访学都是怎么找的啊? +3 luokereng 2024-04-22 3/150 2024-04-24 11:40 by xuechenli
[论文投稿] 期刊推荐 20+4 木颜尘ip 2024-04-22 7/350 2024-04-24 10:06 by bobvan
[考博] 申博成果界定是根据Jcr分区还是中科院分区 +4 我属驴核动力驴 2024-04-22 5/250 2024-04-24 08:47 by 晓目崇
[教师之家] 找工作过程中碰到好多人问南京某新换校长211的政策,规定如下 +15 ZHONGWU_U 2024-04-21 17/850 2024-04-24 05:12 by juta256
[考博] 研二光催化6月底4篇2区 +7 wjtab 2024-04-22 11/550 2024-04-23 06:59 by byron2012
[论文投稿] 编辑是选国外的好还是国内的好。 +8 lizhengke06 2024-04-20 8/400 2024-04-22 08:58 by cuiyunjian
信息提示
请填处理意见