版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (571)
>考研 (67)
>导师招生 (46)
>虫友互识 (22)
>考博 (21)
>学术会议 (8)
>硕博家园 (3)
>休闲灌水 (2)
>材料工程 (1)
>博后之家 (1)
>金属 (1)
>基金申请 (1)
>教师之家 (1)
>健康生活 (1)
>文献求助 (1)
>药学 (1)

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

memgr

铜虫 (正式写手)

应助: 0 (幼儿园)
金币: 229.8
散金: 200
红花: 1
帖子: 330
在线: 101.4小时
虫号: 859962
注册: 2009-09-29
专业: 理论和计算化学

[交流] 【求助】计算结构大小已有1人参与

我想问一下怎么用gromacs分析聚集结构的大小，也就是如果我想用gromacs算一下我的表面活性剂聚集后每一簇所含的单体数目，这样的话怎么办呢

回复此楼

» 猜你喜欢

求调剂已经有14人回复
273求调剂已经有43人回复
307求调剂已经有13人回复
327求调剂已经有7人回复
338求调剂已经有6人回复
326分，一志愿沪9，求生物学调剂已经有3人回复
材料科学与工程320求调剂，080500 已经有9人回复
0703化学调剂 348分已经有10人回复
一志愿211，0703化学305分求调剂已经有15人回复
284求调剂已经有14人回复

1楼 2010-07-21 16:13:06

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

itismineok

捐助贵宾 (正式写手)

应助: 0 (幼儿园)
金币: 537.8
散金: 50
帖子: 372
在线: 97.5小时
虫号: 425200
注册: 2007-07-22
性别: GG
专业: 生物大分子结构与功能

★ ★ ★
雪狼乖乖(金币+3):谢谢希望能完整化 2010-07-28 21:08:00

解决这个问题分两步：
1. 首先使用top2psf.pl将gromacs的 .top file（或者） .tpr file转换成.psf文件

详细见http://www.gromacs.org/Developer_Zone/Programming_Guide/VMD

2. 然后将gromacs的.trr File(或者.trj File)和转换成的.psf文件输入到VMD中

编写tcl文件分析，举个例子：
set psf your.psf
set dcd your.dcd

mol load psf $psf dcd $dcd

set outfile [open cluster.dat w];

set nf [molinfo top get numframes]

# Select cluster.
#set sel [atomselect top "resname polar"]
# polar 表面活性剂的残基名

for {set i 1 } {$i < $nf } { incr i } {

# Select cluster number.
set cluster [atomselect top "resname polar and name only" frame $i]

set crd [$cluster get {x y z}]

#set clusters {}{}

foreach a $crd b $crd {
set diffr [veclength [vecsub $a $b]]

         #在这里根据表面活性剂分子之间的距离，生成clusters

}

   puts $outfile "$clusters"

   #set cluster_num {}
foreach a $clusters {
   #在这里统计每个cluster中分子的数量,

}

   puts $outfile "$cluster_num"

}

close $outfile

mol delete top

[ Last edited by itismineok on 2010-7-29 at 09:17 ]

赞一下(1人)