24小时热门版块排行榜    

查看: 1120  |  回复: 13

ouczhouyi

铜虫 (初入文坛)

[交流] 【求助】namd 密度控制 已有5人参与

本人是一新手,先正模仿一遍文献,做水在cnt中的md,文献中要求密度恒定,不知如何解决
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sjnyongle

金虫 (小有名气)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ghcacj(金币+2):谢谢 2010-11-16 14:01:46
密度恒定,就是要求在NPT系综下进行模拟,水在CNT中的模拟NAMD有个tutorial.你可以看一下!
2楼2010-11-15 17:45:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ouczhouyi

铜虫 (初入文坛)

引用回帖:
Originally posted by sjnyongle at 2010-11-15 17:45:32:
密度恒定,就是要求在NPT系综下进行模拟,水在CNT中的模拟NAMD有个tutorial.你可以看一下!

谢谢
3楼2010-11-16 21:33:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gzgzgz

铜虫 (小有名气)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ghcacj(金币+3):谢谢 2010-11-17 13:03:32
n=N/V,按定义应该是N和V是常数,那么密度为常数。所以真正的密度恒定应该是在NVT系综下实现的。然而模拟一般不可能一开始就达到想要的某一密度,所以先用NPT平衡调控到想要的密度平均值附近,然后转到NVT下模拟。除非你的“密度”不是体系的密度,那另说了。
引用回帖:
Originally posted by sjnyongle at 2010-11-15 17:45:32:
密度恒定,就是要求在NPT系综下进行模拟,水在CNT中的模拟NAMD有个tutorial.你可以看一下!

[ Last edited by gzgzgz on 2010-11-17 at 08:55 ]
4楼2010-11-17 08:51:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sjnyongle

金虫 (小有名气)

受教了,谢谢!
5楼2010-11-17 10:18:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

girljxj

新虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zh1987hs(金币+1):鼓励交流 2010-11-18 23:22:49
请问怎样才能在NPT条件下调控到目的密度的大概值?如果p固定,不是要在模拟过程中改变体积来达到压力的平衡吗?
6楼2010-11-18 18:16:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

席月调琴

新虫 (初入文坛)

还是密度问题——虔心求教

★
小木虫(金币+0.5):给个红包,谢谢回帖交流
NAMD模拟中,如果指定密度计算,怎么进行啊?比如计算体系的密度在0.001g/ml。多谢各位前辈
7楼2010-11-28 20:13:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ouczhouyi

铜虫 (初入文坛)

引用回帖:
Originally posted by gzgzgz at 2010-11-17 08:51:03:
n=N/V,按定义应该是N和V是常数,那么密度为常数。所以真正的密度恒定应该是在NVT系综下实现的。然而模拟一般不可能一开始就达到想要的某一密度,所以先用NPT平衡调控到想要的密度平均值附近,然后转到NVT下模拟。 ...

如果我做的是水在cnt中的MD,要控制水的密度,如何实现。再如果不控制密度,要做出密度与其他量的坐标图,如管子半径、温度,又该如何实现。问题有点多啊,麻烦您了,谢谢
8楼2010-12-01 20:52:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gzgzgz

铜虫 (小有名气)

★
小木虫(金币+0.5):给个红包,谢谢回帖交流
根据V=N/n,你可以算出你需要密度的盒子的体积吧?模拟程序一般会输出盒子大小改变的三维信息,自己可以算V=x*y*z.
引用回帖:
Originally posted by girljxj at 2010-11-18 18:16:17:
请问怎样才能在NPT条件下调控到目的密度的大概值?如果p固定,不是要在模拟过程中改变体积来达到压力的平衡吗?

9楼2010-12-02 01:13:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gzgzgz

铜虫 (小有名气)

★
小木虫(金币+0.5):给个红包,谢谢回帖交流
Av=6.02*1E-23
MW =  molecular weight, for example:  CH4 = 16 g/Mol
No. of molecules in 0.001 g = (0.001*Av/MW)
1mL = (0.01*1E10)^3 = 1E-24 Ang^3
n = (0.001*Av/MW) / (1E-24) = (1E21*Av/MW)   

所以,如果你有N个分子参与模拟,则要把你的体系体积控制在 N/n .
引用回帖:
Originally posted by 席月调琴 at 2010-11-28 20:13:24:
NAMD模拟中,如果指定密度计算,怎么进行啊?比如计算体系的密度在0.001g/ml。多谢各位前辈

[ Last edited by gzgzgz on 2010-12-2 at 01:29 ]
10楼2010-12-02 01:22:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ouczhouyi 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 WqtG0Aem8Ufl 2026-10-24 3/150 2026-10-25 20:24 by cLOojpOxkX0G
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 WqtG0Aem8Ufl 2026-10-24 3/150 2026-10-25 20:11 by cLOojpOxkX0G
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 WqtG0Aem8Ufl 2026-10-24 5/250 2026-10-25 20:07 by cLOojpOxkX0G
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 3xN9SualEUPt 2026-10-24 5/250 2026-10-25 20:01 by cLOojpOxkX0G
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 14:20 by XiYP04It3Evb
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 13:54 by XiYP04It3Evb
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 7xHjyB3bAiwg 2026-10-24 4/200 2026-10-25 13:15 by XiYP04It3Evb
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 7xHjyB3bAiwg 2026-10-24 6/300 2026-10-25 13:11 by XiYP04It3Evb
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-23 9/450 2026-10-25 13:02 by XiYP04It3Evb
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-25 12:59 by XiYP04It3Evb
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-25 12:51 by XiYP04It3Evb
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-25 12:49 by XiYP04It3Evb
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 11/550 2026-10-25 12:45 by XiYP04It3Evb
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 4/200 2026-10-25 12:36 by XiYP04It3Evb
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 GVGRDgNfv7P4 2026-10-23 4/200 2026-10-25 12:10 by XiYP04It3Evb
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-24 4/200 2026-10-25 01:22 by XcOYG53d3x4p
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 11/550 2026-10-24 23:35 by XcOYG53d3x4p
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-24 3/150 2026-10-24 18:23 by WqtG0Aem8Ufl
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 7xHjyB3bAiwg 2026-10-24 3/150 2026-10-24 15:41 by Equinoxhua
[考博] 郑州大学申博 +5 stars~ 2026-10-19 8/400 2026-10-22 15:06 by stars~
信息提示
请填处理意见