计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 分子模拟 (Molecular Modeling) › 用vmd显示sdf计算结果里面中心分子的每一个原子周围也有 ...

用vmd显示sdf计算结果里面中心分子的每一个原子周围也有溶剂分布

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

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

帖子模式

pg-zhd

发布时间: 2024-8-26 14:43

正文摘要:

我的gromacs版本是2021.6，vmd是1.9.3 我用以下代码做的sdf计算 gmx trjconv -s md.tpr -f md.xtc -n index.ndx -o md_cnt.xtc -boxcenter tric -ur compact -pbc none Select group for output Selected 0: 'Sy ...

回复 Reply

pg-zhd 发表于 Post on 2024-8-29 14:22:20

sobereva 发表于 2024-8-28 08:51
你算sdf时选的整个体系，包含了烯烃分子，当然在烯烃原子上也有等值面分布

跟归一化做法有关。默认情 ...

谢谢sob老师的指导，我去尝试一下

sobereva 发表于 Post on 2024-8-28 08:51:42

pg-zhd 发表于 2024-8-27 10:14
好的sob老师，我以后注意用词。
我以为isovalue的数值就是一定范围以外的分布。所以才那么说的。
我的 ...

你算sdf时选的整个体系，包含了烯烃分子，当然在烯烃原子上也有等值面分布

跟归一化做法有关。默认情况下spatial命令是将sdf除以了一个系数，以使得数值的数量级比较大。如果用-nodiv则不会做此处理

pg-zhd 发表于 Post on 2024-8-27 10:14:11

sobereva 发表于 2024-8-27 05:12
确认组选择得没问题。你的输入并没有直接体现“我设置的应该是一定范围以外的溶剂分布”

那叫命令，不叫 ...

好的sob老师，我以后注意用词。
我以为isovalue的数值就是一定范围以外的分布。所以才那么说的。
我的体系只有两种分子（1个烯烃，2000个溶剂）我的ndx文件是把md.gro文件里面的烯烃分子的所有原子改个名字，然后使用gmx make_ndx命令直接按q就保存了，在vmd里面我也看过了就是烯烃分子没有问题。
然后我在执行烯烃置于盒子中间，消除转动以及sdf命令的时候都是选择烯烃作为中心原子和另一个选项是整个体系（system）。
我还有一点很不理解的就是我的isovalue值特别大，我看别人最多就0.001或者0.01这个数量级的，可以麻烦老师解答一下吗？

sobereva 发表于 Post on 2024-8-27 05:12:33

确认组选择得没问题。你的输入并没有直接体现“我设置的应该是一定范围以外的溶剂分布”

那叫命令，不叫代码，看
计算化学中的一些常见不良写法和用词
http://sobereva.com/298（http://bbs.keinsci.com/thread-1358-1-1.html）

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