【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 2329  |  回复: 28
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 小范范1989 的 16 个金币 ,回帖就立即获得 1 个金币,每人有 1 次机会

小范范1989

木虫 (著名写手)


[交流] 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力场还缺少参数?具体应该如何做呢,有没有相应的帖子或者案例。谢谢

我想模拟的就是有机分子做成 的薄膜,不知道老师有什么好的建议,或者相对清晰的模拟过程。

谢谢老师的指点。

gromacs模拟咨询,谢谢指点
1.png


gromacs模拟咨询,谢谢指点-1
2.png
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

leafbrooks3

金虫 (正式写手)


★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小范范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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

星城辛未

木虫 (正式写手)



小范范1989(金币+1): 谢谢参与
祝早得解答
5楼2017-10-16 22:21:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leafbrooks3

金虫 (正式写手)



小范范1989(金币+1): 谢谢参与
说一下力场参数吧,不是所有的分子都有已有的力场,有些不常用的分子就没有,需要自己建参数。这个作者模拟的系统里就有这样的分子

发自小木虫Android客户端
11楼2017-10-17 10:52:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)


送红花一朵
引用回帖:
13楼: Originally posted by leafbrooks3 at 2017-10-17 11:08:52
1:周期性退火不是找能量最低点的,能量最低点是通过EM来找的。这里模拟的重点不是要得到0K时能量最小的点,对于大多数有机系统这样的能量最低点是很难找到的,而且也没有太大意义。模拟的重点是要找到实验温度和压 ...

非常感谢老师指点,谢谢老师。
自己初学gromacs,很多不明白的地方,多谢老师指点
19楼2017-10-17 14:39:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

青青白白123

新虫 (小有名气)



小范范1989(金币+1): 谢谢参与
引用回帖:
13楼: Originally posted by leafbrooks3 at 2017-10-17 11:08:52
1:周期性退火不是找能量最低点的,能量最低点是通过EM来找的。这里模拟的重点不是要得到0K时能量最小的点,对于大多数有机系统这样的能量最低点是很难找到的,而且也没有太大意义。模拟的重点是要找到实验温度和压 ...

老师,请问模拟退火的作用是什么呢?我要模拟膜的相变,看不同温度对膜的影响,怎么定义mdp文件来实现不同的温度呢?希望得到老师的解答!
26楼2018-01-12 14:27:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)


引用回帖:
26楼: Originally posted by 青青白白123 at 2018-01-12 14:27:41
老师,请问模拟退火的作用是什么呢?我要模拟膜的相变,看不同温度对膜的影响,怎么定义mdp文件来实现不同的温度呢?希望得到老师的解答!...

你好,我的理解是,退火模拟是一个温度变化的过程,是为了来寻找全局最小点,因为只是通过能量最小化,得到的是一个局部的极小点。
mdp中,你查阅一下,我这里
annealing = single
annealing_npoints = 2
annealing_time = 0 100 ;ps
annealing_temp = 0 298.15

具体含义,你查阅一下就明白了
27楼2018-01-13 16:12:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

青青白白123

新虫 (小有名气)


引用回帖:
27楼: Originally posted by 小范范1989 at 2018-01-13 16:12:35
你好,我的理解是,退火模拟是一个温度变化的过程,是为了来寻找全局最小点,因为只是通过能量最小化,得到的是一个局部的极小点。
mdp中,你查阅一下,我这里
annealing = single
annealing_npoints = 2
anne ...

那是通过看退火过程的不同温度对膜的影响吗?

发自小木虫Android客户端
28楼2018-01-13 22:55:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

随风zzd

银虫 (小有名气)



小范范1989(金币+1): 谢谢参与
送红花一朵
请问你模拟出的薄膜能看出或是分析出成膜性吗,请求指点。
29楼2019-01-09 11:26:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2017-10-16 21:23   回复  
小范范1989(金币+1): 谢谢参与
发自小木虫Android客户端
dongmings3楼
2017-10-16 21:26   回复  
小范范1989(金币+1): 谢谢参与
发自小木虫IOS客户端
nono20094楼
2017-10-16 21:34   回复  
小范范1989(金币+1): 谢谢参与
发自小木虫Android客户端
朙天儿6楼
2017-10-16 23:07   回复  
小范范1989(金币+1): 谢谢参与
发自小木虫IOS客户端
gengxin607楼
2017-10-17 09:37   回复  
小范范1989(金币+1): 谢谢参与
skyish8楼
2017-10-17 09:51   回复  
小范范1989(金币+1): 谢谢参与
发自小木虫IOS客户端
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): 谢谢参与
2017-10-17 11:12   回复  
小范范1989(金币+1): 谢谢参与
k 发自小木虫IOS客户端
2017-10-17 11:12   回复  
小范范1989(金币+1): 谢谢参与
^_^ 发自小木虫IOS客户端
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   回复  
小范范1989(金币+1): 谢谢参与
发自小木虫Android客户端
2017-10-17 23:00   回复  
小范范1989(金币+1): 谢谢参与
zhgjiang25楼
2017-10-17 23:47   回复  
小范范1989(金币+1): 谢谢参与
发自小木虫IOS客户端
相关版块跳转 我要订阅楼主 小范范1989 的主题更新
提示: 如果您在30分钟内回复过其他散金贴,则可能无法领取此贴金币
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金和生小孩 +27 Ausy 2024-04-15 27/1350 2024-04-16 23:10 by 博csc分享
[论文投稿] The total environment格式 25+4 lyyyyds0 2024-04-13 8/400 2024-04-16 20:02 by 浅梦觉得
[考研] 345 求调剂 +4 南风轻语5 2024-04-12 6/300 2024-04-16 16:36 by chenweiwade
[考研] 295 +3 南山晤梅落 2024-04-12 8/400 2024-04-16 07:24 by waxx333
[考研] 一志愿北化085600 328求线上调剂 +8 Mrshark 2024-04-13 8/400 2024-04-15 22:32 by zzzhonghl
[考研] 调剂 +8 人间十三月@ 2024-04-10 10/500 2024-04-15 22:01 by zzzhonghl
[考研] 一志愿南京理工大学材料工程282求调剂 +13 abc16785 2024-04-11 13/650 2024-04-15 11:26 by coolfishwll
[考研] 280求调剂 +4 黑皮冰棒 2024-04-12 7/350 2024-04-14 23:06 by 永字号
[基金申请] 求问大神,教育部人文社科还有其他基金项目组成员是不是一定要有教授成员 +3 nanco123 2024-04-11 3/150 2024-04-13 11:56 by ljh9201
[教师之家] 表弟博士毕业,选30万安家费还是副教授??? +10 yingyu427 2024-04-11 12/600 2024-04-13 11:16 by karl_messi
[硕博家园] 生物炭采购 +3 锦鲤附体@ 2024-04-12 3/150 2024-04-12 22:10 by Danny614
[考研] 317求调剂 +4 xiuxii 2024-04-10 4/200 2024-04-12 06:55 by xlxu1103
[考研] 312求调剂 +16 dhdhjdhd 2024-04-10 20/1000 2024-04-11 21:34 by kaiyanjie
[考研] 290 生物与医药求调剂 +4 GPX4 2024-04-11 4/200 2024-04-11 21:28 by hunanzang
[考研] 085600 材料与化工 求调剂 总分 310 +22 a1845052693 2024-04-10 24/1200 2024-04-11 11:12 by houyaoxu
[考研] 求调剂 +5 para1527 2024-04-10 10/500 2024-04-10 23:21 by para1527
[硕博家园] 导师选择 +3 任晓霖 2024-04-10 3/150 2024-04-10 21:32 by 欢乐颂叶蓁
[考研] 求调剂 +3 LED绿光灯 2024-04-10 6/300 2024-04-10 17:22 by LED绿光灯
[考研] 290求调剂 +6 Z~N- 2024-04-10 7/350 2024-04-10 17:19 by Ghker
[考研] 材化专硕,一志愿中科大284求调剂 +16 玉上花 2024-04-10 17/850 2024-04-10 15:08 by qlm5820
信息提示
请填处理意见