牧生 发表于 2025-5-30 16:53 好的好的，感谢老师回复 |
oorangew 发表于 2025-5-30 16:22 用软件画好分子结构式以后，存为mol2格式，直接用sobtop调用openbabel就可以得到MMFF94电荷，不需要预先进行优化的。 http://sobereva.com/soft/Sobtop/
|
牧生 发表于 2023-5-30 22:45 老师你好，如果是碳氢氧为主的上千原子的木质素分子（结构也比较复杂）的优化有什么建议呢？分子大部分是苯环，木质素分子结构类似帖子中的贴图一，帖子链接：http://bbs.keinsci.com/thread-53626-1-1.html。后续是想用sobtop生成GAFF力场的拓扑，调用openbabel产生MMFF94电荷。 我目前的步骤是用MS建模之后用其的forcite中的compass II力场优化了一下之后得到pdb后，用openbabel转成mol2文件，就使用sobtop产生了MMFF94电荷的拓扑文件。 这样的步骤是否不得行？请老师指教  |
slxc920113 发表于 2023-7-7 15:33 理解了，谢谢！ |
本帖最后由 slxc920113 于 2023-7-7 15:35 编辑 somacule 发表于 2023-7-7 15:14 直接用MD的结构肯定是不会特别准，但是也不至于定性错误。因为需要平衡结构，你用最后一帧的结构引入了温度导致的振动微扰了，还需要做一下极小化（结构优化）才行。 |
slxc920113 发表于 2023-6-1 00:34 谢谢您的解答。还想请问一下，如果用MD之后的一帧的结构，不再优化，而是直接进行电荷的计算，是否可行呢 |
Geocmpt 发表于 2023-5-31 14:52 我算过大环分子和富勒烯结合的体系，400个原子，我是算的动的，就是得好几天，因为服务器比较多，所以等得起 |
slxc920113 发表于 2023-6-1 00:34 谢谢您的解答，很好。 |
| 你这个最开始的结构本身合理性就不高，直接优化很容易出问题。我之前帮你算过一次，流程是这样的：先用MMFF94电荷，GAFF力场做结构优化，然后换到XTB做结构优化，再换到ORCA用B97-3c做结构优化，最后用优化的结构计算AM1-BCC电荷（这一步基本没啥成本）。RESP电荷你就别指望了，我最近在写一个自动切片算法，未来有可能解决大体系的RESP电荷计算问题（仅限于中性体系，分片无法自动识别离子官能团）。 |
牧生 发表于 2023-5-30 22:45 谢谢您。我用的64核机子，试着没有直接用懒人脚本。首先自己手动xtb优化，然后ORCA继续优化，算不动，所以还是放弃了硬算的念头，转投了mmff94。 |
|
最好给出具体结构，或者如果需要保密，给出一个高度类似体系的结构 不同的分子往往有不同的最适合的处理方法，泛泛地说很难有针对性，有些特殊的技巧更是对于特定特征的分子才能用 |
本帖最后由 牧生 于 2023-5-30 22:54 编辑 Geocmpt 发表于 2023-5-30 16:00 我目前最多算过300原子的左右的，碳氢氧氮硫为主的化合物，睡一觉就算完了。 有个好一点的cpu，比如13900k，用orca计算resp电荷的懒人脚本，对于500多原子，如果硬造，也不是不行。有可能睡一觉就算完了，实在不行，睡两觉。 如果上千原子，可能选择余地就不大了，还得mm94。 别的一些建议：建模的时候，就尽量让其结构合理些。可以先用xtb优化以后，再用懒人脚本，大概率能节省一些时间。 |
GoldenBaby 发表于 2023-5-30 11:34 谢谢您的指导。还有一点，我这个是非均质，无定形的柔性有机大分子，可能无法像聚合物那样拆解“”精确计算电荷了。MMFF94电荷不依赖于具体构象，可能是个没办法的办法了。 |
|
大分子电荷肯定不能直接硬算的。 仔细阅读http://sobereva.com/soft/Sobtop/#ex5 并领会精神 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-14 08:48 , Processed in 0.231059 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites