24小时热门版块排行榜    

查看: 780  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

小兵ai石头

铜虫 (初入文坛)

[求助] 在Linux上运行NAMD总是意外中断且不报错

我是分子模拟初学者。

在Linux上运行NAMD,总是会意外中断,查看.log文件也没发现报错,温度也正常;用VMD查看体系也没有爆开。
我用restart文件接着跑,过几个小时后又会断了,还是不报错。

另外,我跑的多个NAMD程序经常都在同一时间断开。由于是后台运行,所以应该可以排除断网的原因。

有没有虫友知道可能是什么原因导致的啊?可能是配置文件的问题吗(多个NAMD程序的配置文件相似)?
希望知道的虫友帮帮我吧!!

在这里把配置文件贴上:
##################################################
# MD SECTION
##################################################

numsteps             5000000
set temperature      298


# INPUT FILES
structure          ../../2-common/bcdsfla_wb.psf
coordinates        ../../2-common/bcdsfla_wb.pdb

firsttimestep      0
set inputname      emeq_s
binCoordinates     ../../3-emeqmd/sflabcd-cenC4/$inputname.restart.coor
binVelocities      ../../3-emeqmd/sflabcd-cenC4/$inputname.restart.vel  ;# remove the "temperature" entry if you use this!
extendedSystem           ../../3-emeqmd/sflabcd-cenC4/$inputname.restart.xsc
set outputname     abf_00
#temperature       $temperature


# FORCE FIELD
parameters          ../../2-common/fla_parm.inp
parameters          ../../2-common/aglc_parm.inp
paraTypeCharmm       on


# Force-Field Parameters
exclude              scaled1-4
1-4scaling           1.0
switching            on
switchdist           12.0                  
cutoff               14.0
pairlistdist         15.0                  
hgroupcutoff         2.8   


# Periodic Boundary Conditions
cellBasisVector1    35.0   0.0   0.0      
cellBasisVector2     0.0  35.0   0.0
cellBasisVector3     0.0    0   50.0
cellOrigin         0.132  -0.145  -2.452   
wrapAll             on                  
margin               3


# PME (for full-system periodic electrostatics)
PME                 yes
PMEGridSpacing      1.0


# Constant Temperature Control
langevin            on   
langevinDamping     1      ;# damping coefficient (gamma) of 1/ps
langevinTemp        $temperature
langevinHydrogen    off    ;# don't couple langevin bath to hydrogens


# Constant Pressure Control (variable volume)
useGroupPressure      yes                 ;# needed for rigidBonds

langevinPiston        on               
langevinPistonTarget  1.01325             ;# in bar -> 1 atm
langevinPistonPeriod  100.0               ;# fs!
langevinPistonDecay   50.0                ;# fs!
langevinPistonTemp    $temperature


# Integrator Parameters
timestep            2.0  
nonbondedFreq       1  
fullElectFrequency  2   
stepspercycle       10
splitpatch          hydrogen                   # how to assign atoms to patches


#Constraints and Restraints
rigidBonds          all  ;# needed for 2fs steps
rigidtolerance      0.000001
rigiditerations     400
rigiddieonerror     on

fixedAtoms           on
fixedAtomsFile       ../../3-emeqmd/sflabcd-cenC4/fixedatoms_s.pdb
fixedAtomsForces     off
fixedAtomsCol        B

constraints          on                                       
consexp              2                                          
consref              ../../3-emeqmd/sflabcd-cenC4/consref_s.pdb  
conskfile            ../../3-emeqmd/sflabcd-cenC4/consref_s.pdb   
conskcol             B                                       

selectConstraints    on              ;# restrain the positions of atoms to plane or a line in space
selectConstrX        on
selectConstrY        on
selectConstrZ        off

# OUTPUT FILES                                #!
outputname           abf_00_0
restartname          abf_00
dcdFile              abf_00.dcd

binaryoutput         no  
binaryrestart        yes

outputenergies       500
outputtiming         500
outputpressure       500
restartfreq          500
XSTFreq              500
dcdFreq              500

##################################################
# ABF SECTION   
##################################################
                     
colvars              on
colvarsConfig        distance.in
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小兵ai石头

铜虫 (初入文坛)

后来解决了,每次提交完任务之后都要退出提交界面,就不会停了~这个跟用的软件和系统有关~
3楼2014-10-28 11:39:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

小兵ai石头

铜虫 (初入文坛)

有没有可能是内存不足?
但是其他师兄在同一个节点上跑的就没有问题啊,只不过师兄跑的不是NAMD
2楼2014-06-04 10:38:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scutxl

新虫 (初入文坛)

楼主你好!
我最近也在学着用NAMD,也遇到和你差不多的问题,我是直接在自己计算机跑的,还没上服务器,
请问你说退出提交界面是什么意思,我直接用终端提交的,关了终端就会把NAMD停掉了。
4楼2016-04-10 16:01:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +4 a089 2026-08-07 4/200 2026-08-08 12:15 by arzu_hma
[基金申请] 关于代码变化问题,想知道的进来 +10 且听虎啸 2026-08-07 13/650 2026-08-08 12:09 by sakura12139
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] fileCode有新解读? +8 Tide man 2026-08-08 16/800 2026-08-08 09:59 by Tide man
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 08:59 by 6vVgjDL4CnGu
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 08:27 by 6vVgjDL4CnGu
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见