计算化学公社

 找回密码 Forget password
 注册 Register
Views: 14733|回复 Reply: 9
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] gromacs拉伸C60小球与蛋白碰撞,距离蛋白质很近的地方小球会偏移是什么原因?

[复制链接 Copy URL]

9

帖子

0

威望

29

eV
积分
38

Level 2 能力者

各位大神们好,请问我想模拟拉伸C60小球与一个蛋白相互碰撞接触模拟,已经设置了沿一个方向拉伸C60小球,但是距离蛋白质很近的地方小球就会偏移,无法与蛋白接触是什么原因呢,谢谢大家

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

2#
发表于 Post on 2021-5-29 05:35:08 | 只看该作者 Only view this author
截图说明是什么情况,这样描述太抽象

模拟碰撞也不一定要用拉伸,直接在gro里设置初速度也可以

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示,求助帖必须在帖子标题明确体现出此帖内容是求助或提问,仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“gromacs拉伸配体与蛋白碰撞 ”改了,以后务必注意
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

9

帖子

0

威望

29

eV
积分
38

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2021-5-29 13:18:42 | 只看该作者 Only view this author
sobereva 发表于 2021-5-29 05:35
截图说明是什么情况,这样描述太抽象

模拟碰撞也不一定要用拉伸,直接在gro里设置初速度也可以

老师向上是y方向,已经把蛋白固定,设置了沿y方向牵引小球,但是蛋白与小球始终接触不上,一旦距离过近就会偏离y方向,学生有两个问题,1、为什么会接触不上2、明明已经设置了沿y方向,为什么小球还会沿别的方向偏移,谢谢老师

313

帖子

2

威望

3900

eV
积分
4253

Level 6 (一方通行)

4#
发表于 Post on 2021-5-29 13:43:39 | 只看该作者 Only view this author
正常吧,不可能碰撞上,理应会快接触的时候反弹。如果你的c60不带电,你提取c60 和蛋白之间的LJ作用应该就能看出来。gmx energy

9

帖子

0

威望

29

eV
积分
38

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2021-5-29 14:12:54 | 只看该作者 Only view this author
lyj714 发表于 2021-5-29 13:43
正常吧,不可能碰撞上,理应会快接触的时候反弹。如果你的c60不带电,你提取c60 和蛋白之间的LJ作用应该就 ...

之前看文献有的人是用lammps做的碰撞上了,我用的是gromacs就不知道为啥碰撞不上,还想请教一个问题我已经用energy提取了两者之间的时间与LJ势能的数据,用什么命令能求出球与蛋白之间的相互作用力与时间的数据呢?谢谢

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

6#
发表于 Post on 2021-5-30 01:32:36 | 只看该作者 Only view this author
我不知道你怎么定义碰撞。如果说肉眼看上去C60和蛋白质的原子都快挨到一起了,那明显不可能发生,要知道LJ势当中的互斥势部分是干什么的,当原子间距离小于平衡距离后就会感受到很强的斥力驱使它们分离。

而且你的体系也不是左右对称的,不可能始终恰好老老实实笔直地沿着Y方向运动,此体系中C60在X方向上感受到的外势都不是对称的。

头脑里有势函数的概念自然就会明白看到的现象
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

9

帖子

0

威望

29

eV
积分
38

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2021-5-30 14:32:28 | 只看该作者 Only view this author
sobereva 发表于 2021-5-30 01:32
我不知道你怎么定义碰撞。如果说肉眼看上去C60和蛋白质的原子都快挨到一起了,那明显不可能发生,要知道LJ ...

好的老师,我现在已经用gmx energy 求出来了他俩之间的相互作用能与时间的关系,那我想得出他俩之间的相互作用力该用什么命令呢,还是说把能量求导就可以得出相互作用力与时间的关系,谢谢老师

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

8#
发表于 Post on 2021-5-31 05:03:43 | 只看该作者 Only view this author
MrWang 发表于 2021-5-30 14:32
好的老师,我现在已经用gmx energy 求出来了他俩之间的相互作用能与时间的关系,那我想得出他俩之间的相 ...

没法简单得到,需要自己写程序对每一帧基于原子电荷和力场参数直接计算
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

313

帖子

2

威望

3900

eV
积分
4253

Level 6 (一方通行)

9#
发表于 Post on 2021-5-31 12:45:53 | 只看该作者 Only view this author
sobereva 发表于 2021-5-31 05:03
没法简单得到,需要自己写程序对每一帧基于原子电荷和力场参数直接计算

老师说的这个有文章理论参考么?我查过一些文献上的只看到过通过对相互作用能对距离求导数的那种方式

6万

帖子

99

威望

6万

eV
积分
125148

管理员

公社社长

10#
发表于 Post on 2021-6-1 04:16:39 | 只看该作者 Only view this author
lyj714 发表于 2021-5-31 12:45
老师说的这个有文章理论参考么?我查过一些文献上的只看到过通过对相互作用能对距离求导数的那种方式

自己编程,计算两个片段间每一对粒子的作用力然后加和
知道力场函数形式自然就知道怎么实现
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-23 01:21 , Processed in 0.204567 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list