计算化学公社

标题: 拉伸分子动力学问题 [打印本页]

作者
Author: 虔诚的人 时间: 2018-4-9 12:48
标题: 拉伸分子动力学问题
请问大家，用gmx做拉伸分子动力学，我想拉开我是分子，我实现了。但是移动的力应该是作用在质心上的。我想作用在分子边缘的几个原子上，可以吗？怎么实现？

就是力的作用点可以不以自己设定？

作者
Author: mol 时间: 2018-4-9 13:52
一个猜测，lz可以试试：可以多写几个pull_group，每个pull_group一个原子。

作者
Author: 虔诚的人 时间: 2018-4-9 15:09

mol 发表于 2018-4-9 13:52
一个猜测，lz可以试试：可以多写几个pull_group，每个pull_group一个原子。

“lz”是什么意思，好尴尬

我是这么理解的，可以把一个分子的几个原子写在index文件中，pull_group=index，这样好像可以。你觉得呢？

作者
Author: zsu007 时间: 2018-4-9 16:47
lz=楼主

作者
Author: sobereva 时间: 2018-4-9 18:15

虔诚的人发表于 2018-4-9 15:09
“lz”是什么意思，好尴尬
我是这么理解的，可以把一个分子的几个原子写在index文件中，pull_group ...

诸如pull-group1-name = left，则.ndx里定义的left组里若只包含某个局部区域的的原子，那么就只有这些原子会被拉

作者
Author: 虔诚的人 时间: 2018-4-9 20:14

zsu007 发表于 2018-4-9 16:47
lz=楼主

哈哈，好尴尬

作者
Author: 虔诚的人 时间: 2018-4-9 20:36

sobereva 发表于 2018-4-9 18:15
诸如pull-group1-name = left，则.ndx里定义的left组里若只包含某个局部区域的的原子，那么就只有这些 ...

谢谢sob老师

作者
Author: 青青青 时间: 2019-4-13 15:37
本帖最后由青青青于 2019-4-13 16:04 编辑

请问下图文献里的力常数1000 KJ mol−1 nm−2，是在mdp文件的pull_k1中设置的吗，意思是指在多肽和磷脂膜之间设置了恒力吧。谢谢 (, 下载次数 Times of downloads: 37)

作者
Author: 少年爱吃地三鲜 时间: 2019-4-13 15:53
请问你用的什么力场？

作者
Author: tjuptz 时间: 2019-8-13 15:45

青青青发表于 2019-4-13 15:37
请问下图文献里的力常数1000 KJ mol−1 nm−2，是在mdp文件的pull_k1中设置的 ...

这个是施加限制势做限制性动力学模拟

作者
Author: 青青青 时间: 2019-8-16 10:08

tjuptz 发表于 2019-8-13 15:45
这个是施加限制势做限制性动力学模拟

嗯嗯，谢谢

欢迎光临计算化学公社 (http://bbs.keinsci.com/)