数字完全不一样呀,应该是一样的嘛不是因为经过了换算嘛 |
luzujia 发表于 2024-11-5 16:00 单位不同把 |
钟老师,ztop在屏幕上显示的k和写入到top文件里的不一样,是两套不同的参数C和F? |
钟老师,能讲讲“该文件包含了所有的坐标和对应的量化能量,如果找到了就利用该文件中的数据进行拟合“ 这个拟合二面角的具体算法吗 |
刚读完钟老师的文章,还在想怎么引入二面角参数,就刷到这个帖子了 |
钟老师,您好,请问您在文中提到的“sob介绍的手动方法”是什么呢?有相关帖子么? |
本帖最后由 enthalpy 于 2023-2-27 21:42 编辑 很好的工具。请教以下,如果我已经得到QM以及MM(待拟合二面角的K设置为0)的能量,是不是也可用手动去调用脚本中函数fit_k_only和 fit_least_square 来拟合得到相关的参数? |
本帖最后由 Medivan 于 2022-5-28 04:25 编辑 请问老师,这种情况应该如何选择周期参数呢,period的1234排列组合都试过了,都不能把中间的峰压下去 |
ggdh 发表于 2022-4-20 19:15 明白了,谢谢老师! |
snljty 发表于 2022-4-20 17:42 因为最后二面角参数需要和RESP2(0.5)去匹配(共同决定二面角势能面) 所以拟合二面角参数需要用RESP2(0.5) |
钟老师好,我想请教一个问题,如果我想模拟某个分子在溶液中的动力学行为,那么我计算模拟使用的原子电荷时,可以考虑RESP2(0.5)。那我在去掉取代基,拟合二面角的时候,应该怎么考虑能量和原子电荷呢?之前和导师讨论的结果,是对二面角做柔性扫描,和对扫描的每个点做高精度单点计算的时候用真空下的结果即可。那就有个问题,拟合二面角我用的原子电荷应该是计算的真空下的原子电荷,还是诸如RESP2(0.5)等方法描述的溶液中的原子电荷呢?谢谢老师~ |
感谢大神分享!非常有用。 如果我想要调用的是OPLS-AA力场,是不是就不可以用ztop来模拟二面角参数了。。 |
ggdh 发表于 2021-11-6 09:24 谢谢钟老师,我来试试哈~ |
本帖最后由 ggdh 于 2021-11-6 09:37 编辑 好像不方便,因为tleap好像无法载入prmtop文件你可以用gromacs 加完水以后 再转成amber(用ztop就能直接转) |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-23 05:14 , Processed in 0.203879 second(s), 33 queries , Gzip On.