应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Gromacs/Amber/NAMD » 【求助】请教第二次 预平衡 问题
北京石油化工学院2026年研究生招生接收调剂公告
31/1返回列表
查看: 740  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

crispsoft

金虫 (小有名气)

[交流] 【求助】请教第二次 预平衡 问题

各位高手：
我最近按照《蛋白质结构预测实验指南》一书中的例子做一个蛋白的minimize 和平衡。
首先是固定CA碳进行第一次预平衡。这个没有问题。
问题出在第二次预平衡（全放开的结构优化过程）：
我安装书上的要求，将structure 后文件改为xxx.coor，按照.xsc文件内容修改盒子信息，并取消了CA碳的限制并取消minimize。
运行后提示：
charm++ FATAL ERROR：
FATAL ERROR: ABNORMAL EOF FOUND - buffer=*@@*

——请问这是为什么？
谢谢

附conf文件：
#############################################################
## JOB DESCRIPTION                                         ##
#############################################################
# Minimization and Equilibration of
# protein in a Water Box

#############################################################
## ADJUSTABLE PARAMETERS                                   ##
#############################################################
structure          ../2/eq1.coor                                  #书上说这步用coor文件取代原来的psf文件
coordinates        ../common/ionized.pdb
set temperature    310
set outputname     eq2
firsttimestep      0

#############################################################
## SIMULATION PARAMETERS                                   ##
#############################################################
# Input
paraTypeCharmm     on
parameters          ../common/par_all27_prot_lipid.inp
temperature         $temperature

# Force-Field Parameters
exclude             scaled1-4
1-4scaling          1.0
cutoff              12.0
switching           on
switchdist          10.0
pairlistdist        13.5

# Integrator Parameters
timestep            2.0  ;# 2fs/step
rigidBonds          all  ;# needed for 2fs steps
nonbondedFreq       2
fullElectFrequency  4  
stepspercycle       20

# Constant Temperature Control
langevin            on    ;# do langevin dynamics
langevinDamping     5     ;#  
langevinHydrogen    off    ;#

# Periodic Boundary Conditions  
cellBasisVector1    57.8367    0.   0.      #这里的信息我按照xsc文件改了
cellBasisVector2     0.   57.8367    0.
cellBasisVector3     0.    0   48.1973   
cellOrigin          15.01   5.22  -9.89
wrapAll             on

# PME (for full-system periodic electrostatics)
PME                 yes
PMEGridSpacing      1.0
#manual grid definition
PMEGridSizeX        72
PMEGridSizeY        72
PMEGridSizeZ        64  

# Constant Pressure Control (variable volume)
useGroupPressure      yes ;# needed for rigidBonds
useFlexibleCell       no
useConstantArea       no
langevinPiston        on
langevinPistonTarget  1.01325 ;#  in bar -> 1 atm
langevinPistonPeriod  200.0
langevinPistonDecay   100.0
langevinPistonTemp    $temperature
if {0} {
constraints        on
consref            ../common/ionized.pdb
conskfile          ../common/ionized.pdb
conskcol           B
constraintscaling  10.0
}
# Output
outputName          $outputname
restartname       eq2r
restartfreq         1000     ;#  
dcdfreq             1000
xstFreq             1000
outputEnergies      1000
outputPressure      1000

#############################################################
## EXTRA PARAMETERS                                        ##
#############################################################

#############################################################
## EXECUTION SCRIPT                                        ##
#############################################################
# Minimization
if {0} {
minimize            1000
reinitvels          $temperature
}
run 2500000 ;#

[ Last edited by crispsoft on 2011-3-20 at 16:36 ]
回复此楼

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2011-03-20 16:32:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crispsoft

金虫 (小有名气)
引用回帖:
Originally posted by crispsoft at 2011-03-20 16:32:56:
各位高手：
我最近按照《蛋白质结构预测实验指南》一书中的例子做一个蛋白的minimize 和平衡。
首先是固定CA碳进行第一次预平衡。这个没有问题。
问题出在第二次预平衡（全放开的结构优化过程）：
我安装书上 ...

有同学指出
structure          ../2/eq1.coor不对，应该在bincoordinates 中指定。
我想知道除了这个命令外还需要其他命名吗？另外，如果有了这个命令，原来的structure还需要再指定吗？

谢谢
赞一下 回复此楼
2楼2011-03-21 12:33:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wally8962

木虫 (著名写手)
★ ★ ★ ★ ★ ★
crispsoft(金币+6): 谢谢帮助 2011-03-22 12:03:47
ghcacj(金币+6): 谢谢 2011-03-22 12:09:57
structure和coordinate不是一回事，你得搞清楚了。多看ug

# coordinates $ <$ coordinate PDB file $ >$
Acceptable Values: UNIX filename
Description: The PDB file containing initial position coordinate data. Note that path names can be either absolute or relative. Only one value may be specified.

# structure $ <$ PSF file $ >$
Acceptable Values: UNIX filename
Description: The X-PLOR format PSF file describing the molecular system to be simulated. Only one value may be specified.

bincoordinates $ <$ binary coordinate restart file $ >$
Acceptable Values: UNIX filename
Description: The binary restart file containing initial position coordinate data. A binary coordinate restart file is created as output from NAMD by activating the binaryrestart or binaryoutput options. Note that, in the current implementation at least, the bincoordinates option must be used in addition to the coordinates option, but the positions specified by coordinates will then be ignored.
赞一下(1人) 回复此楼
3楼2011-03-22 11:31:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 crispsoft 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 347材料专硕求调剂 +12 zj8215216 2026-04-06 12/600 2026-04-07 16:42 by 律香川照之
[考研] 化工学硕 285求调剂 +17 Wisjxn 2026-04-07 17/850 2026-04-07 15:55 by 乔哒哒哒
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
[考研] 385分 生物学（071000）求调剂 +12 qf626 2026-04-01 12/600 2026-04-06 17:36 by 52305043001
[考研] 复试调剂 +5 asdasdassda 2026-04-05 5/250 2026-04-06 09:32 by dongzh2009
[考研] 22408 总分320，一篇论文二作，两个国三，求调剂 +3 Leomulufu 2026-04-04 5/250 2026-04-05 19:04 by chongya
[考研] 288求调剂 一志愿哈工大 材料与化工 +13 洛神哥哥 2026-04-03 13/650 2026-04-05 17:27 by zzx2138
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-04-04 5/250 2026-04-05 14:06 by imissbao
[考研] 271分求调剂学校 +12 zph158488！ 2026-04-02 13/650 2026-04-05 10:13 by lqwchd
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
[考研] 296材料专硕求调剂 +21 202451007219 2026-04-02 22/1100 2026-04-04 21:48 by hemengdong
[考研] 357求调剂 +13 1050389037 2026-04-03 13/650 2026-04-03 22:27 by 无际的草原
[考研] 311求调剂 +11 勇敢的小吴 2026-04-02 11/550 2026-04-03 21:46 by qlm5820
[考研] 一志愿重庆大学085404，总分314分，求调剂 +4 zf83hn 2026-04-03 4/200 2026-04-03 21:25 by 啵啵啵0119
[考研] 材料科学与工程339求调剂 +12 hyz0119 2026-03-31 13/650 2026-04-03 18:33 by ls刘帅
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-04-02 6/300 2026-04-03 10:19 by 蓝云思雨
[考研] 求调剂 302分初试 0854 +5 伶可乐 2026-04-02 5/250 2026-04-02 17:53 by 笔落锦州
[考研] 321求调剂 一志愿 浙江工业大学生物医药 +5 嘿嘿HC 2026-04-01 6/300 2026-04-02 15:23 by sophie2180
[考研] 310分求调剂 +4 成功上岸wang 2026-04-01 4/200 2026-04-01 20:35 by liu823948201
[考研] 一志愿 南京航空航天大学 ，080500材料科学与工程学硕 +10 @taotao 2026-03-31 11/550 2026-04-01 09:43 by xiayizhi
信息提示
关闭
请填处理意见
关闭