qjy111 发表于 2023-6-29 12:21 不应当用sobtop处理含一堆分子的盒子的.pdb文件，sobtop处理的对象是单个分子 如果你是指靠sobtop来把pdb转成gro，根本没这个必要，grompp本来就能用pdb而不非得要求gro 水应当直接用现成的水模型，诸如SPC/E、OPC3，而不应当自己去搞它的拓扑文件 另外我不知道你用sobtop处理的是阴阳离子，还是对阴离子和阳离子分别处理的。必须分别处理，而且二者需分别优化，而且需要确保得到的是比较延展的构象。 |
sobereva 发表于 2023-6-28 00:35 好的好的，老师，我下次注意。我模拟的是离子液体水溶液的密度，因为采用相同的力场信息，所以模拟纯IL以及加水的都是采用RESP(0.5)计算出的电荷。在采用RESP(0.5)计算电荷中，我先是在# opt freq b3lyp/6-311++g(d) geom=connectivity下进行阴阳离子的几何优优化+振动分析并保留.fch，基于.fch文件保存为.gif文件。然后基于# B3LYP/TZVP em=GD3BJ opt→# B3LYP/def2TZVP em=GD3BJ geom=allcheck→# B3LYP/def2TZVP em=GD3BJ scrf=solvent=water geom=allcheck，进行计算然后取平均值（溶剂为水）。而拓扑文件是基于上述的.fch文件保存为.mol文件,采用sobtop例2的方法产生的拓扑文件。还有一个疑问是packmol产生的盒子.pdb文件，再用sobtop例2的方法产生的.gro文件，并没有进行取名称并进行分组编序，我试了gromac自带命令将.pdb文件转为.gro文件，可是依然没有名称，请问一下有什么解决办法嘛 |
|
不要用“sobtop的第二种方法”这种描述，根本不存在什么官方的第x种方法一说 提问时先交代自己算的是什么东西 纯物质密度算合理了再说加水的事 拓扑文件里原子顺序和结构文件对应就完了，不用管原子名 所有对结果可能产生影响的要素都交代清楚，诸如RESP(0.5)在不同级别、不同溶剂下计算结果明显不一样 认真看此文 在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚 http://sobereva.com/620（http://bbs.keinsci.com/thread-25787-1-1.html） |
pal 发表于 2023-6-27 10:58 是这个原因，但是看之前的论坛说只要你可以保证是对应上的，就可以忽略警告，然后我就忽略了，但是模拟出来的结果却误差大，所以我感觉还是有点问题。 |
qjy111 发表于 2023-6-27 10:38 第一个警告应该就是gro第二列的原子后面没有数字，跟itp里面的原子对应不起来 |
pal 发表于 2023-6-27 10:29 packmol生成mix.pdb文件，然后将mix.pdb文件用sobtop生成.gro文件 |
| 这个gro是怎么产生的，第一列全是1MOL，top文件肯定对应不起来 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-22 07:12 , Processed in 0.215656 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites