24小时热门版块排行榜    

查看: 1507  |  回复: 6

oudegeshen

银虫 (小有名气)

[求助] 在DPD Calculation设定Temperature=1.0,为什么呢?那么kBT=1.0意味着kB=1,T=1吗? 已有1人参与

在DPD Calculation设定Temperature=1.0,为什么呢?那么kBT=1.0意味着kB=1,T=1吗?可是KB不是波尔兹曼常数吗,怎么会等于1?温度怎么会等于1?
谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuntao118

新虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
oudegeshen: 金币+5, ★★★很有帮助 2015-05-09 11:13:46
abinitio: 金币+1, 鼓励交流 2015-05-11 08:30:54
不是这样的,温度在模拟式要无量刚化:模拟温度=模拟体系实际温度*波尔兹曼常数/体系选定的单位能量(e)来确定的,通常我们都是将kBT作为 单位能量了(1),比如实际温度298K,能量参数选为0.155kal/mol最终的模拟温度大约是3.8,因为这里面0.155(kal/mol)/kB=78.2K, 3.8=298/78.2.,个人认为,MD,DPD,MC中的初始温度方法都是一样的。
2楼2015-05-08 21:04:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oudegeshen

银虫 (小有名气)

引用回帖:
2楼: Originally posted by chuntao118 at 2015-05-08 21:04:03
不是这样的,温度在模拟式要无量刚化:模拟温度=模拟体系实际温度*波尔兹曼常数/体系选定的单位能量(e)来确定的,通常我们都是将kBT作为 单位能量了(1),比如实际温度298K,能量参数选为0.155kal/mol最终的模拟温 ...

O(∩_∩)O谢谢
模拟温度为什么默认为1.0呢?按你的意思298K时,T应该设定为3.8啊?但是T的范围只能在某一个范围(没记错的话好像是0.1-3.0)。您说能量参数选0.155kal/mol,请问为什么选这个值呢?

再次感!
3楼2015-05-09 11:13:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oudegeshen

银虫 (小有名气)

引用回帖:
2楼: Originally posted by chuntao118 at 2015-05-08 21:04:03
不是这样的,温度在模拟式要无量刚化:模拟温度=模拟体系实际温度*波尔兹曼常数/体系选定的单位能量(e)来确定的,通常我们都是将kBT作为 单位能量了(1),比如实际温度298K,能量参数选为0.155kal/mol最终的模拟温 ...

暂时给您5个金币,请您不要介意。
4楼2015-05-09 11:15:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuntao118

新虫 (著名写手)

【答案】应助回帖

★ ★ ★
oudegeshen: 金币+3, 有帮助 2015-05-11 12:52:09
引用回帖:
4楼: Originally posted by oudegeshen at 2015-05-09 11:15:11
暂时给您5个金币,请您不要介意。...

我只是给你举个例子,这是在水分子体系下选择的,不同体系有不同的选择,要看你模拟的体系是什么体系。体系变化就会跟着变化。
5楼2015-05-09 13:22:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oudegeshen

银虫 (小有名气)

引用回帖:
5楼: Originally posted by chuntao118 at 2015-05-09 13:22:17
我只是给你举个例子,这是在水分子体系下选择的,不同体系有不同的选择,要看你模拟的体系是什么体系。体系变化就会跟着变化。...

非常感谢您的指点。因为这几天我有期末考试,所以未能及时回复您。现在我还有一些疑问,如果您不介意的话还请指教。1.怎样知道我的体系的能量参数呢?2.我模拟时设定T=1.0(default),只是参考文献中别人这么用所以我也这么用了,但是不懂为什么。3.我的体系是药物分子&聚合物分子混合物,以及聚合物A&聚合物B混合物。O(∩_∩)O谢谢!
6楼2015-05-13 09:37:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuntao118

新虫 (著名写手)

引用回帖:
6楼: Originally posted by oudegeshen at 2015-05-13 09:37:00
非常感谢您的指点。因为这几天我有期末考试,所以未能及时回复您。现在我还有一些疑问,如果您不介意的话还请指教。1.怎样知道我的体系的能量参数呢?2.我模拟时设定T=1.0(default),只是参考文献中别人这么用所 ...

1.怎样知道我的体系的能量参数呢?
这个一般在势能函数中会体现出来,要是L-J势的话,前面的势阱EPSION就就选择单位能量单位的依据,这个值就像前面所举例子的0.155Kal/mol.
2.我模拟时设定T=1.0(default),只是参考文献中别人这么用所以我也这么用了,但是不懂为什么。
这个原因前面已经说过了。
3.我的体系是药物分子&聚合物分子混合物,以及聚合物A&聚合物B混合物?
对于这样的共混体系,选择主要依赖于不同分子之间的作用势的势阱来选择,你可以选择A中分子之间的作用势的势阱来作为单位能量,也可以考虑B中分子之间作用势的势阱来作为单位能量,但是不管怎么选择,一定只能以一个为为标准,其他的都通过选定的表示出来就行了。另外你可以追溯你做参考的文献,一定可以找到答案的,除非你所参考的文献的作者根本就不懂这个,只是套用别人的程序,别人的数据来糊弄人的,建议好好读读经典的教材。比如《The Art of Molecular Dynamics Simulation》、《 Introduction to Practice of Molecular Simulation 》后面这本书有专门一张是介绍DPD 的,但前面这本书李子类方法的经典。(个人观点),并且DPD就来源于MD.其实我不做DPD,只是看过DPD 的原创的文献,大概是92年的吧!
7楼2015-05-13 13:47:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 oudegeshen 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国社科又开始会评了,不知道这次命运如何 +3 雨打竹帘 2026-08-30 3/150 2026-08-30 09:43 by KNIGHT2011
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 06:46 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 06:11 by ZPa0EcMwuECS
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 06:00 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 13/650 2026-08-30 05:38 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 01:36 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 01:14 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:41 by ZPa0EcMwuECS
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 6/300 2026-08-29 18:10 by manplx
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见