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

【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 1429  |  回复: 5

sealanlan

铁虫 (正式写手)

[求助] 两个结构类似的分子做柔性势能面扫描,结果差别很大,可否帮忙看下怎么回事?已有1人参与

两个类似的分子,固定开环两个碳原子之间的距离,用关键词opt=modredundant,做开环反应的柔性势能面扫描,结果差别很大,可否帮忙看下怎么回事?

第一个分子能够从闭合状态顺利地优化到开环状态,高斯view打开看中间结构,中间文件过渡态一切正常,如下:
势能面扫描结果
http://muchongimg.xmcimg.com/oss ... _1604036712_548.png
起始闭环结构
http://muchongimg.xmcimg.com/oss ... _1604036788_386.png
曲线最高点-过渡态结构
http://muchongimg.xmcimg.com/oss ... _1604036826_228.png
最终闭环结构
http://muchongimg.xmcimg.com/oss ... _1604036854_507.png

第二个分子能够从闭合状态顺利地优化到过渡态,从过渡态最高点之后,所有结构均发生扭曲变形,导致扫描不能进行到开环状态,详细结果如下:
势能面扫描结果
http://muchongimg.xmcimg.com/oss ... _1604037000_277.png
起始闭环结构
http://muchongimg.xmcimg.com/oss ... _1604037084_623.png
曲线最高点--过渡态正常结构
http://muchongimg.xmcimg.com/oss ... _1604037128_409.png
过渡态之后--扭曲变形结构


这是第二个分子热不稳定吗?还是计算方法出了问题?

两个分子的开环状态和闭环状态都是稳定的,我已经得到稳定结构构型,为啥势能面扫描会出这样的问题呢?
接下来我该如何寻找分子2的过渡态呢?opt=ts已经试过了,不太行,不得已做的势能面扫描,我该如何改进势能面扫描呢?

两个结构类似的分子做柔性势能面扫描,结果差别很大,可否帮忙看下怎么回事?
势能面扫描曲线.png


两个结构类似的分子做柔性势能面扫描,结果差别很大,可否帮忙看下怎么回事?-1
分子1起始闭环结构.png


两个结构类似的分子做柔性势能面扫描,结果差别很大,可否帮忙看下怎么回事?-2
分子1中间过渡态结构.png


两个结构类似的分子做柔性势能面扫描,结果差别很大,可否帮忙看下怎么回事?-3
分子1结束闭环结构.png


两个结构类似的分子做柔性势能面扫描,结果差别很大,可否帮忙看下怎么回事?-4
分子2势能面扫描曲线.png


两个结构类似的分子做柔性势能面扫描,结果差别很大,可否帮忙看下怎么回事?-5
分子2起始闭环结构.png


两个结构类似的分子做柔性势能面扫描,结果差别很大,可否帮忙看下怎么回事?-6
分子2中间曲线最高点过渡态结构.png


两个结构类似的分子做柔性势能面扫描,结果差别很大,可否帮忙看下怎么回事?-7
分子2过渡态之后的扭曲变形结构.png
回复此楼
努力
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
sealanlan: 金币+15, ★★★很有帮助, 非常感谢,我试试看 2020-10-30 15:31:54
这个扫描是沿着两个原子间距离还进行的。而开环 往往不单纯是就是严格沿着两个原子之间的距离来进行开环的。单一一个内坐标变量的扫描不能完全重现开环过程。但是既然扫描已经出现顶点,是潜在的可能的过渡态结构。不妨就把两个过渡态都优化出来, 跑IRC曲线。如果过渡态结构优化正确。IRC曲线应该就能更好的体现开环的过程了。
2楼2020-10-30 15:12:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sealanlan

铁虫 (正式写手)

引用回帖:
2楼: Originally posted by paramecium86 at 2020-10-30 15:12:38
这个扫描是沿着两个原子间距离还进行的。而开环 往往不单纯是就是严格沿着两个原子之间的距离来进行开环的。单一一个内坐标变量的扫描不能完全重现开环过程。但是既然扫描已经出现顶点,是潜在的可能的过渡态结构。 ...

扫描是沿着两个C-C原子间距离进行的
努力
3楼2020-10-30 15:32:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sealanlan

铁虫 (正式写手)

引用回帖:
2楼: Originally posted by paramecium86 at 2020-10-30 15:12:38
这个扫描是沿着两个原子间距离还进行的。而开环 往往不单纯是就是严格沿着两个原子之间的距离来进行开环的。单一一个内坐标变量的扫描不能完全重现开环过程。但是既然扫描已经出现顶点,是潜在的可能的过渡态结构。 ...

现在还有个问题就是opt=TS,优化不出来过渡态,
附件里分别是opt=TS优化的分子1和分子2的log文件,目前为止优化好多天了,一直没有结果
努力
4楼2020-10-30 15:40:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sealanlan

铁虫 (正式写手)

引用回帖:
4楼: Originally posted by sealanlan at 2020-10-30 15:40:51
现在还有个问题就是opt=TS,优化不出来过渡态,
附件里分别是opt=TS优化的分子1和分子2的log文件,目前为止优化好多天了,一直没有结果...

附件文件刚上传失败了,再传一次

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:libolin3@tal.com
  • 附件 1 : 01_NH_ts.log
  • 2020-10-30 15:42:32, 2.56 M
  • 附件 2 : 02_CH2_ts.log
  • 2020-10-30 15:43:44, 3.3 M
努力
5楼2020-10-30 15:44:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

引用回帖:
5楼: Originally posted by sealanlan at 2020-10-30 02:44:07
附件文件刚上传失败了,再传一次...

那看起来顶点的地方不是 过渡态的那种顶点,而是结构有很大变化时候出现能能量突跃。
6楼2020-11-01 17:12:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sealanlan 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[留学DIY] 10043以后咋出去啊 +3 RuiLove 2024-04-22 3/150 2024-04-27 12:16 by icm639
[教师之家] 博士论文被抄袭 +26 和尚敲小木鱼 2024-04-22 43/2150 2024-04-27 11:14 by liuyuzhou
[论文投稿] 投稿RSC旗下杂志突然间看不到投稿状态,里面投稿的文章信息不见了,有遇到过的吗? 50+3 dlying 2024-04-26 3/150 2024-04-27 09:58 by bobvan
[论文投稿] AC审稿 +3 马winner 2024-04-26 6/300 2024-04-27 09:08 by 马winner
[硕博家园] 聊天 +12 暮色恋伊人 2024-04-22 13/650 2024-04-27 08:40 by WASM
[硕博家园] 博士白读了 +47 Da_Meng_Zi 2024-04-21 52/2600 2024-04-27 08:25 by shl2112501
[论文投稿] LWT投 +3 AChen92 2024-04-26 3/150 2024-04-26 22:16 by hizifu
[考研] 没学上 +6 季向阳 2024-04-26 12/600 2024-04-26 21:06 by 季向阳
[有机交流] 环肽的合成 +3 徐来不惊 2024-04-25 5/250 2024-04-26 16:56 by 徐来不惊
[考博] 申博求助 +4 dskabdh 2024-04-24 11/550 2024-04-26 15:54 by dskabdh
[硕博家园] 考研,求职还是考编? +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+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 雪径踏青
[基金申请] 化学结清有情发出来了 +3 starboy7286 2024-04-20 3/150 2024-04-24 15:52 by 嘿,黑贝
[考博] 博士招生 +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
[考博] 申博成果界定是根据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
信息提示
请填处理意见