24小时热门版块排行榜    

查看: 429  |  回复: 8
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yongbinz

金虫 (著名写手)

[交流] 【求助】求助:周期性边界模拟后没有氢键能部分?

我用MS模拟聚合物的聚集态,将3条高分子链放在一个周期性边界条件下进行分子模拟,分别经能量优化、几何优化和分子动力学优化后,数据列表中没有发现氢键能,请高手告知原因。是不是氢键的计算要另外设置呀?
Total enthalpy                  : 15274.343 kcal/mol
Total energy                    : 15295.459 kcal/mol
Temperature                     :  504.480 K
Pressure                        :   -0.004 GPa
Volume                          : 36924.219 A^3
Density                         :    0.979 g/cm^3

Contributions to total energy (kcal/mol):
  Valence energy (diag. terms)  : 10278.774
    Bond                        : 3657.819
    Angle                       : 3678.615
    Torsion                     : 2711.663
    Inversion                   :  230.678
  Valence energy (cross terms)  :   77.612
    Stretch-Stretch             :  225.128
    Stretch-Bond-Stretch        :  -65.384
    Stretch-Torsion-Stretch     :    0.000
    Separated-Stretch-Stretch   :    0.000
    Torsion-Stretch             :    0.000
    Bend-Bend                   :  -40.082
    Torsion-Bend-Bend           :  -42.050
    Bend-Torsion-Bend           :    0.000
  Non-bond energy               : 1549.700
    Hydrogen bond               :    0.000
    van der Waals               : 3169.431
    Electrostatic               : -1619.731

Cell parameters:     a:  33.37 A       b:  33.22 A        c:  33.31 A
                 alpha:  90.00 deg  beta:  90.00 deg  gamma:  90.00 deg
另外,我最后一步是用NPH进行的动力学优化,但优化后,温度怎么上升到了500K呢?我前一步NVT优化的温度可是300K呀/

我刚接触这个软件不久,请大家不吝赐教!本人不胜感谢!

[ Last edited by freshgirl on 2009-6-20 at 12:20 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongbinz

金虫 (著名写手)

我也想知道呀!
8楼2009-06-03 11:43:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

saitou

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yongbinz(金币+15,VIP+0):谢谢呀!非常感谢! 5-21 13:20
wuchenwf(金币+5,VIP+0):谢谢 5-21 21:34
1 不是没有氢键能,而是到了pcff和compass里氢键能不再作为单独的项列出,而这个条目是为了老的立场保留的。如果需要计算氢键能,可以使用其他办法,需要手动写perl脚本

2 NPH是不固定温度的吧,这和你之前用NVT平衡过是两码事。到了500k稳定那就是体系需要平衡到的状态,没有了热浴,温度保持不住也不奇怪
不抛弃,不放弃
2楼2009-05-21 13:01:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongbinz

金虫 (著名写手)

谢谢上面的热情回复!

引用回帖:
Originally posted by saitou at 2009-5-21 13:01:
1 不是没有氢键能,而是到了pcff和compass里氢键能不再作为单独的项列出,而这个条目是为了老的立场保留的。如果需要计算氢键能,可以使用其他办法,需要手动写perl脚本

2 NPH是不固定温度的吧,这和你之前用N ...

1您的意思是不是说,我的数据里的非键能部分就已经包括氢键能了,只是系统没有单独列出吗?
2虽然NPH不固定温度,但是经NPH模拟后的最终温度达500K,那么这个数据还能用吗?我的真实状态条件是300K的呀!
有时间的话,请您回复我,好吗?不胜感谢!
3楼2009-05-21 13:20:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saitou

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yongbinz(金币+15,VIP+0):谢谢回复! 5-22 00:10
qasd(金币+2,VIP+0):谢谢~ 6-3 22:10
1 正是如此,所以要把氢键能再剔出来我还是只会用笨办法,哈哈

2 如果要求温度稳定于300k的话 or 温度的影响重大,不可随意,那么建议考虑包含T的系综,NPT、NVT==,我没有做过固定H的模拟,结果如何解释一时间没有头绪
不抛弃,不放弃
4楼2009-05-22 00:03:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见