“第25届北京科音初级量子化学培训班” 将于6月11至14日于北京举办,是初学者真正从头一次性完整、系统性学习量子化学计算,从而能很快上手量子化学研究的重要机会,比起自己摸索、鼓捣能少走无数弯路。报名正在进行中,请点击此链接查看,欢迎参加!

计算化学公社

 找回密码 Forget password
 注册 Register

请问用alphafold预测的肽结构,改成线性结构后做几何优化无法收敛如何解决?

查看数: 29 | 评论数: 2 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2026-6-3 00:36

正文摘要:

老师我预测得结构具有二级结构,修饰非标准氨基酸残基,做几何优化后直接跑MD容易崩溃,所以我应gaussview修饰了主链的N-CA-C-N二面角为180度,再用avogadro的MMFF94s预优化,再用orca做的opt,但是一直无法收敛,请 ...

回复 Reply

sobereva 发表于 Post on 1 hour ago
“无法收敛” 指代不明。如果单纯是几何优化不收敛,常规该怎么解决怎么解决,本论坛上有无数讨论,论坛首页google框搜索,下文也说了基本解决思想
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

动力学崩溃的最大可能就是拓扑文件有bug,这是真正的重中之重,然而你怎么准备的拓扑文件却一丁点都没说(无数人在网上提问时都犯这个致命问题)。怎么排查和解决在http://sobereva.com/soft/Sobtop#FAQ8里明确说了。这事和量子化学做几何优化没有丝毫直接联系,当前做这么大体系的量子化学优化纯粹白浪费时间。
虽然初始结构里有不合理接触也很容易导致刚一开始崩溃,但如果已经用MMFF94优化过了,这种因素导致问题的可能性微乎其微。

student0618 发表于 Post on 1 hour ago
本帖最后由 student0618 于 2026-6-3 04:04 编辑

目的是什么?补参数?

要做非标准残基参数,只把残基拿出来封端算就行,不用算整条链。

参考Gromacs非标准氨基酸残基的蛋白MD建模:sobtop产生rtp文件
http://bbs.keinsci.com/forum.php ... 57309&fromuid=64740
(出处: 计算化学公社)


还要注意,短肽用AF未必准确。

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-6-3 05:50 , Processed in 0.325029 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list