24小时热门版块排行榜    

Znn3bq.jpeg
查看: 464  |  回复: 0

vallen

银虫 (正式写手)

[求助] Discover中怎样用BTCL编程计算两组分子之间的相互作用能

我现在有一段轨迹,想分析里面的两组分子之间作用能随时间的变化

网上有个例子是这样的
readFile returns the total energy of the ith step
# This loop will read all available MD records in the his file
for {$i = 1} {[set energy [readFile history frame = $i]] != ""} {incr i} {
# We can control which records will be calculated
if { $i < 3000 } {
# Output the step number to stdout
echo $i
# Write the calculation results to the output text file
puts $file "Step = $i $energy"
# Perform the nonbond energy calculation here,
# the subsets' names should be changed for a different system
EnergyAnalysis pv19,ureum -fid $file
# Output a newline
puts $file ""
} else {
break
}
}

这样的写法不知道对不对?
我不知道怎么从我的体系中选择两组分子  即  subset,用BTCL语句是怎么写的?
是不是要用subset define  命令啊?具体怎么写呢?
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 vallen 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 急需调剂 +5 绝不放弃22 2026-04-15 5/250 2026-04-15 23:11 by Tide man
[考研] 296求调剂 +12 汪!?! 2026-04-09 13/650 2026-04-15 20:01 by 学员JpLReM
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 化工学硕294分,求导师收留 +32 yzyzx 2026-04-12 36/1800 2026-04-14 17:45 by lhj2009
[考研] 290求调剂 +21 luoziheng 2026-04-10 23/1150 2026-04-14 15:49 by zs92450
[考研] 271求调剂 +35 2261744733 2026-04-11 41/2050 2026-04-14 15:36 by zs92450
[考研] 297求调剂 +23 ORCHID1 2026-04-10 26/1300 2026-04-14 13:52 by 陈皮皮
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] 302求调剂 +10 易!? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[考研] 一志愿西交机械专硕求调剂 +9 求上岸的小王 2026-04-10 9/450 2026-04-13 16:08 by jiangguiquan11
[考研] 材料复试求调剂 +24 xhhdjdjsjks 2026-04-09 24/1200 2026-04-13 15:49 by 幸免 ..
[考研] 生物学调剂 +11 小冉要努力 2026-04-10 13/650 2026-04-13 11:46 by 电化学及催化
[考研] +10 李多米lee. 2026-04-12 11/550 2026-04-12 22:58 by yuyin1233
[考研] 344 材料专业 求调剂211 无地域要求 +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
[考研] 277 数一104,学硕,求调剂 +21 瓶子PZ 2026-04-09 23/1150 2026-04-11 23:12 by labixiaoqiao
[考研] 一志愿厦大0856,306求调剂 +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
[考研] 282,求调剂 +12 jggshjkkm 2026-04-09 14/700 2026-04-11 09:39 by 猪会飞
[考研] 本科211 工科085400 280分求调剂 可跨专业 +11 LZH(等待调剂中 2026-04-10 11/550 2026-04-11 08:39 by zhq0425
[考研] 机械专368 有去处吗 +4 种大树 2026-04-10 4/200 2026-04-10 15:31 by jiajinhpu
信息提示
请填处理意见