24小时热门版块排行榜    

查看: 889  |  回复: 8
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 coolrainbow 的 9 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

coolrainbow

木虫 (著名写手)

未来国家冻凉

[交流] 【求助】请问NAMD算自由能 已有2人参与

我照着教程里面那个Born模型的例子做了一下,结果总是出现:

ETITLE:      TS           BOND          ANGLE          DIHED          IMPRP               ELECT            VDW       BOUNDARY           MISC        KINETIC               TOTAL           TEMP      POTENTIAL         TOTAL3        TEMPAVG            PRESSURE      GPRESSURE         VOLUME       PRESSAVG      GPRESSAVG

ENERGY:       0         0.0000         0.0000         0.0000         0.0000         -50748.4978 9999999999.9999         0.0000         0.0000 9999999999.9999      9999999999.9999 9999999999.9999 9999999999.9999 9999999999.9999 9999999999.9999      9999999999.9999 9999999999.9999     27000.0000 9999999999.9999 9999999999.9999

我怀疑是两个Na原子重叠导致的能量无穷大造成的,但是我是按照教程里面做的啊,怎么会这样了,请大家提提意见好吗?
回复此楼

» 猜你喜欢

技术博客:http://hi.baidu.com/coolrainbow/blog
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

coolrainbow

木虫 (著名写手)

未来国家冻凉

引用回帖:
Originally posted by superdirac at 2010-04-18 17:40:13:
你拿什么算的自由能  ABF,FEP, or MCTI??

我用的是fep,按照A Tutorial on Alchemical Free Energy Perturbation Calculations in NAMD的第二个例子做的,可是怎么也做不好,总是能量无穷大,可是按照他的例子做的啊,诡异~~
技术博客:http://hi.baidu.com/coolrainbow/blog
3楼2010-04-18 18:42:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

superdirac

木虫 (正式写手)


coolrainbow(金币+1):谢谢回复 2010-04-18 18:40
lei0736(金币+1):谢谢 2010-04-18 22:23
你拿什么算的自由能  ABF,FEP, or MCTI??
我认为,酒一口一口喝,路一步一步走~步子迈大了,喀~容易扯着蛋
2楼2010-04-18 17:40:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superdirac

木虫 (正式写手)


lei0736(金币+1):谢谢 2010-04-18 22:23
namd 2.6 和 namd 2.7b2.对待 FEP的时候差别很大的!
 那你是用的那个??
我认为,酒一口一口喝,路一步一步走~步子迈大了,喀~容易扯着蛋
4楼2010-04-18 20:57:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

coolrainbow

木虫 (著名写手)

未来国家冻凉

引用回帖:
Originally posted by superdirac at 2010-04-18 20:57:14:
namd 2.6 和 namd 2.7b2.对待 FEP的时候差别很大的!
 那你是用的那个??

我用的是我自己编译的namd2.7b2,用intel+mpich2编译的
conf文件是:
##########################################################

fixedAtoms on
fixedAtomsFile ../sodium.fep
fixedAtomsCol B

#############################################################
## FEP  PARAMETERS                                         ##
#############################################################

source ../fep.tcl

alch on
alchtype fep
alchFile ../sodium.fep
alchCol B
alchOutFile alchemy.alchOutFile
alchOutFreq 10
alchEquilSteps 1000
set nSteps 11000
set dLambda 0.1

runFEP 0.0 1.0 $dLambda $nSteps

我的psf是用alchemify生成的,nnb部分
       1 NNB
       1
       0       1       1       1       1       1       1       1
       1       1       1       1       1       1       1       1
       1       1       1       1       1       1       1       1
       1       1       1       1       1       1       1       1
       1       1       1       1       1       1       1       1
       1       1       1       1       1       1       1       1
....

有什么问题吗?

[ Last edited by coolrainbow on 2010-4-19 at 07:56 ]
技术博客:http://hi.baidu.com/coolrainbow/blog
5楼2010-04-19 07:37:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 281求调剂(0805) +4 烟汐忆海 2026-03-16 11/550 2026-03-18 11:57 by djl2006
[考研] 278求调剂 +5 烟火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[考研] 302求调剂 +9 负心者当诛 2026-03-11 9/450 2026-03-17 17:13 by ruiyingmiao
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 8/400 2026-03-17 16:57 by ruiyingmiao
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
信息提示
请填处理意见