24小时热门版块排行榜    

查看: 2518  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

sealanlan

铁虫 (正式写手)

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

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

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

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

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的回帖
查看全部 6 个回答

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
这个扫描是沿着两个原子间距离还进行的。而开环 往往不单纯是就是严格沿着两个原子之间的距离来进行开环的。单一一个内坐标变量的扫描不能完全重现开环过程。但是既然扫描已经出现顶点,是潜在的可能的过渡态结构。 ...

现在还有个问题就是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文件,目前为止优化好多天了,一直没有结果...

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

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿江南大学085701环境工程专硕总分287求调剂 +5 18266118446 2026-03-09 5/250 2026-03-11 16:51 by 2020015
[基金申请] 提交后的基金本子,已让学校撤回了,可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
[考研] 265求调剂 +8 小木虫085600 2026-03-09 11/550 2026-03-10 19:44 by 绿色低碳化学
[考研] 085602化工求调剂 +7 董boxing 2026-03-10 7/350 2026-03-10 17:07 by BruceLiu320
[考研] 求调剂! +3 朔朔话 2026-03-09 3/150 2026-03-10 15:48 by houyaoxu
[考研] 材料与化工求调剂 +11 与冬清宁 2026-03-07 14/700 2026-03-10 14:23 by Linda Hu
[考研] 0817一志愿苏州大学280复试调剂 +9 kk扛 2026-03-07 9/450 2026-03-10 14:20 by peike
[考研] 研究生招生 +3 徐海涛11 2026-03-10 6/300 2026-03-10 12:43 by 徐海涛11
[考研] 085600材料与化工,一志愿广州985,求调剂 +15 qqyyaill 2026-03-05 15/750 2026-03-10 08:23 by Linda Hu
[考研] 家人们 调剂不迷路 看这里 +8 likeihood 2026-03-09 13/650 2026-03-10 08:09 by likeihood
[考研] 考研调剂,一志愿山东大学材料与化工,328分,政治51 +5 关你西红柿929 2026-03-08 6/300 2026-03-09 13:50 by 新篇章DFSS
[考研] 081700学硕一志愿北京化工大学数二英一过六级有竞赛求调剂 +5 galaxary 2026-03-07 7/350 2026-03-08 04:37 by wang_dand
[考研] 一志愿211 化学305分求调剂 +3 0703杨悦305分 2026-03-05 3/150 2026-03-07 21:15 by yinhuanshun
[考研] 一志愿武汉理工085601,初试301分,请问能调剂到湖北吗 +3 肖yang 2026-03-06 3/150 2026-03-07 19:08 by ms629
[考研] 诚求调剂,323分 有专利/科研/实习经历 +6 橙子cyx 2026-03-06 8/400 2026-03-07 19:06 by 新篇章DFSS
[考研] 求调剂 +4 呼呼?~+123456 2026-03-06 4/200 2026-03-06 23:11 by L135790
[考研] 不限学校专业的调剂同学看过来 +5 啊摆啊摆 2026-03-05 9/450 2026-03-06 12:06 by 啊摆啊摆
[考研] 求材料调剂 +4 berdmond 2026-03-05 4/200 2026-03-05 19:45 by 黑衣馒头人
[考研] 复试调剂 +5 呼呼?~+123456 2026-03-05 5/250 2026-03-05 19:43 by caszguilin
[考研] 纺织、生物、化学、材料等专业 +3 Eember. 2026-03-05 7/350 2026-03-05 16:35 by Eember.
信息提示
请填处理意见