计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 分子模拟 (Molecular Modeling) › 求助：布洛芬100晶面（132个原子）如何生成GROMACS的输 ...

求助：布洛芬100晶面（132个原子）如何生成GROMACS的输入文件

查看数: 2502 | 评论数: 4 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

WuQs

发布时间: 2022-11-16 11:30

正文摘要:

本帖最后由 WuQs 于 2022-11-16 12:51 编辑 1.我有看过sob老师的生成拓扑文件的几种方法，想使用Gaussian+multiwfn+sobtop进行生成；但是132个原子用高斯计算较为困难，按照sob老师之前说的降低基组然后将em和opt ...

回复 Reply

WuQs 发表于 Post on 2022-11-18 14:25:23

sobereva 发表于 2022-11-17 02:17
只要单个分子的拓扑文件构建得得当，现实中能维持的簇结构自然就能维持，计算分子晶体时也能正常维持晶体 ...

阿里嘎多

sobereva 发表于 Post on 2022-11-17 02:17:04

WuQs 发表于 2022-11-16 21:51
好的，太着急了，谢谢老师；但是这个不算是复合物，因为是一个晶体结构，所以这四个分子是一类，都是布 ...

只要单个分子的拓扑文件构建得得当，现实中能维持的簇结构自然就能维持，计算分子晶体时也能正常维持晶体结构。

倘若你就挖一个簇出来单独在真空中跑MD，边缘原子还不做冻结，不严重变形才怪

WuQs 发表于 Post on 2022-11-16 21:51:53

本帖最后由 WuQs 于 2022-11-16 21:54 编辑

好的，太着急了，谢谢老师；但是这个不算是复合物，因为是一个晶体结构，所以这四个分子是一类，都是布洛芬分子，因为不清楚这样子计算是否会破坏晶体的结构，所以就用了一个单元去算

sobereva 发表于 Post on 2022-11-16 20:19:35

1 132原子的计算根本就没什么困难的，哪怕个人PC也不是个什么事。嫌DFT优化这步太耗时就用半经验或GFN-xTB凑合
要搞清楚，是每类分子各自建立一个拓扑文件，用sobtop等工具创建拓扑文件绝对不是对复合物而言的，sobtop主页的FAQ里都明确说了

再顺带一说，有些零基础的GROMACS初学者，居然试图用Sobtop产生整个模拟体系（含一大堆分子）的拓扑文件，然后直接用GROMACS跑，这种做法是极端野蛮粗暴不可理喻的。这么做可能在Sobtop里耗时很高，而且所有分子的拓扑信息都搅合在一起作为一个[ moleculetype ]出现时会特别难以管理和修改。

你对拓扑文件的基本理解有严重问题。先把基本原理想明白再去做，别赖程序

2 “GAFF力场好像并没有官方的力场可以下载”似乎你都不知道什么叫力场
但凡遇到什么问题，首先要想的是怎么正确地解决问题，而不是瞎试一些根本没用的白浪费时间

评分 Rate

参与人数 Participants 1	eV +3	收起理由 Reason
WuQs	+ 3	谢谢

查看全部评分 View all ratings

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register