| 查看: 2644 | 回复: 28 | |||
[交流]
gromacs模拟咨询，谢谢指点
|
|||
|
背景：别人一个gromacs模拟描述如下。（图一，图二） 我现在有以下几个问题不明白，望老师指点一二。 1：他的整个模拟过程应该是em，eq和prod。eq和prod采用了不同的热浴和压浴。但是，她整个过程没有进行周期性退火吧？只是单个的退火。这样能找到最小点吗？ 我理解的退火过程是，高温跨过势垒，然后退到0K的时候，才找到极小点，也就是周期性退火，多少个0K，就有多少个极小点，通过比较极小点，选出最小点。 但是他的这个过程，没有到0K。是不是对了是真实实验对比？或者是其他方面考虑？ 2：如果按照他这个过程模拟，最后选取结构，应该是对最后的那10ns进行了簇分析吧（我的猜测），然后选取了构型研究？或者是一些其他经常用的方法是什么？ 3：我看按照他说的，她整过过程下来30ns左右吧。一般模拟都这个时间范围左右？这个时间是不是偏长？ 4：最后一点，他的那20ns，其中前10ns和后10ns压浴和热浴的方式变了，这个在模拟过程中，是分开模拟前10ns，结束之后再进行后10ns，这样修改mdp文件。 还是说一个mdp文件，能够前10ns采用第一种方式，后10ns采用第二种方式？还是什么办法？ 5：MD模拟的开始，他说缺少参数，通过DFT计算得到了？这个是什么意思？我们对分子做出来的GAFF力场还缺少参数？具体应该如何做呢，有没有相应的帖子或者案例。谢谢 我想模拟的就是有机分子做成 的薄膜，不知道老师有什么好的建议，或者相对清晰的模拟过程。 谢谢老师的指点。  1.png  2.png |
回复此楼» 本帖已获得的红花（最新10朵）
随风zzd» 猜你喜欢
遇见不省心的家人很难过
已经有16人回复
-
退学或坚持读
已经有25人回复
-
博士延得我，科研能力直往上蹿
已经有4人回复
-
免疫学博士有名额，速联系
已经有14人回复
-
面上基金申报没有其他的参与者成吗
已经有4人回复
-
多组分精馏求助
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 抢金币啦！回帖就可以得到:
-
天津科技大学邓启良教授团队 招收2026年博士生
+1/80
-
结构动力学与结构健康监测方向欧盟玛丽居里全奖博士招聘
+1/78
-
非粮生物质能技术全国重点实验室合成生物学创新团队全球招聘博士/博士后
+1/74
-
中国石油大学（北京）国家级大人才团队博士招生2名：化学、材料、石油工程：油田化学
+2/61
-
香港科技大学计算物理及流体力学课题组招收全奖博士后及博士生（2026年9月入学）
+1/44
-
上海市“光探测材料与器件”工程技术研究中心（上海应用技术大学）招聘优秀研究人员
+1/28
-
山东科技大学招聘化学化工博士博士后
+1/28
-
教育部重点实验室和清华大学某国家重点实验室，联合培养硕生、博生，并长期招博士后
+1/27
-
考博求助
+1/14
-
中科院深圳先进院-免疫治疗方向-招收1名博士生(26年9月入学)
+1/13
-
西班牙巴塞罗那访学、博后、留学互动
+1/11
-
南京林业大学国家级青年人才团队招收2026年生物质转化和炭材料等方向博士生
+1/8
-
四川大学华西医院沈百荣教授课题组科研助理招聘启事
+1/8
-
M200-A摩擦磨损试验机
+1/7
-
山东第一医科大学第一附属医院招聘事业编制科研岗
+1/6
-
复旦大学化学系凡勇教授/张凡教授团队招聘博士后
+1/4
-
邵阳学院食品与化学工程学院硕士调剂
+1/3
-
英国兰卡斯特大学（Lancaster University）大模型、计算机视觉PhD招生
+1/3
-
双一流联合团队招聘团队青年人才与博后
+1/2
-
苏州大学医学院纳米生物医学方向招收申请考核制博士生1名
+1/1
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小范范1989: 金币+10, 非常感谢老师指点，谢谢老师 2017-10-17 14:25:50
smutao: 金币+3 2018-01-14 07:00:29
小范范1989: 金币+10, 非常感谢老师指点，谢谢老师 2017-10-17 14:25:50
smutao: 金币+3 2018-01-14 07:00:29
|
1：周期性退火不是找能量最低点的，能量最低点是通过EM来找的。这里模拟的重点不是要得到0K时能量最小的点，对于大多数有机系统这样的能量最低点是很难找到的，而且也没有太大意义。模拟的重点是要找到实验温度和压力下的平衡态结构。所以一般通过EM找到一个能量相对低的点，然后通过eq把系统加热到想要的温度，这个可以分几步，也可以一步到位。最后就是在NVT，或者NPT系统下，进行模拟找到系统的平衡结构。因为要加温，所以起始构型是不是能量最低不太重要，只要g构型差的不是太远，最后加温后的结构都会差不多。 2：这个看自己研究的需要，一般可能就从最后10ns里随机选一些构型作为初始构型进行下一步研究。 3：对于一般的系统来说，30ns不算长。就是一个初始的平衡过程。 4：是的，这个是分两步进行的，需要修改mdp文件 |
» 本帖已获得的红花（最新10朵）
13楼2017-10-17 11:08:52
5楼2017-10-16 22:21:40
11楼2017-10-17 10:52:02
19楼2017-10-17 14:39:18
26楼2018-01-12 14:27:41
27楼2018-01-13 16:12:35
28楼2018-01-13 22:55:45
29楼2019-01-09 11:26:27
简单回复
2017-10-16 21:23
回复
2017-10-16 21:26
回复
nono20094楼
2017-10-16 21:34
回复
小范范1989(金币+1): 谢谢参与
。 发自小木虫Android客户端
朙天儿6楼
2017-10-16 23:07
回复
2017-10-17 09:37
回复
小范范1989(金币+1): 谢谢参与
skyish8楼
2017-10-17 09:51
回复
小范范1989(金币+1): 谢谢参与
一 发自小木虫IOS客户端
一杯绿豆9楼
2017-10-17 10:45
回复
小范范1989(金币+1): 谢谢参与
^_^ 发自小木虫IOS客户端
岩拔t810楼
2017-10-17 10:47
回复
小范范1989(金币+1): 谢谢参与
M 发自小木虫IOS客户端
假大空12楼
2017-10-17 10:53
回复
小范范1989(金币+1): 谢谢参与
顶
syhorchid14楼
2017-10-17 11:11
回复
小范范1989(金币+1): 谢谢参与
诛仙青云15楼
2017-10-17 11:12
回复
小范范1989(金币+1): 谢谢参与
k 发自小木虫IOS客户端
迷丝忒陈16楼
2017-10-17 11:12
回复
小范范1989(金币+1): 谢谢参与
^_^ 发自小木虫IOS客户端
slytjiaofei17楼
2017-10-17 11:30
回复
小范范1989(金币+1): 谢谢参与
沙沙鱼18楼
2017-10-17 11:48
回复
小范范1989(金币+1): 谢谢参与
^_^ 发自小木虫IOS客户端
hongyingwu20楼
2017-10-17 15:03
回复
小范范1989(金币+1): 谢谢参与
qyj200621楼
2017-10-17 18:58
回复
小范范1989(金币+1): 谢谢参与
o 发自小木虫Android客户端
yuzhuce22楼
2017-10-17 20:11
回复
小范范1989(金币+1): 谢谢参与
祝福 发自小木虫IOS客户端
nicenicey23楼
2017-10-17 20:32
回复
sz_kebaoyuan24楼
2017-10-17 23:00
回复
小范范1989(金币+1): 谢谢参与
zhgjiang25楼
2017-10-17 23:47
回复