• 回答数

    2

  • 浏览数

    1776

  • 收藏数

    0

作者:素年未逝 发表于 2018-11-27 16:22:24
跳转到指定楼层
以前amber的in文件都是采用ntb=1来执行周期性盒子,但是从来没有注意过它的具体尺寸,现在想将amber移植到namd中MD,在设参数时碰到问题,就是红色部分的数值从哪可以获知:
TITLE
&cntrl
ntc=2, ntf=2
tol=0.0005,
nstlim=500,
ntpr=50,
ntwr=100,
ntwx=100,
dt=0.001,
tempi=300.,
cut=9.,
&end
&ewald        PME on # Use PME for electrostatic calculation
        # Orthogonal periodic box size
a=
62.23
,    cellBasisVector1 62.23 0 0
b=
62.23
,    cellBasisVector2 0 62.23 0
c=
62.23
,    cellBasisVector3 0 0 62.23
nfft1=
64
,      PMEGridSizeX 64
nfft2=
64
,      PMEGridSizeY 64
nfft3=
64
,     PMEGridSizeZ 64
ischrgd=
1
,     # NAMD doesn’t force neutralization of charge
分享:
回复

使用道具

该用户从未签到

新手上路

Rank: 1

积分
43
极客币
114
主题
16
帖子
47
注册时间
2018-10-19
在线时间
3 小时
性别
保密
发表于 2018-11-27 16:22:41 | 显示全部楼层
我猜是在你的输入坐标文件里面, .rst或者.inpcrd的最后一行-0.4101452 0.2378568 0.2090200 1.0447223 0.8526692 1.8107323
-0.9186951 0.1406647 0.5963050 0.0000000 0.0000000 0.0000000
65.067600060.396600059.6056000
90.000000090.000000090.0000000
我猜是紫色的那三个,但没用namd跑过amber的不完全清楚
回复

使用道具 举报

该用户从未签到

新手上路

Rank: 1

积分
44
极客币
115
主题
17
帖子
51
注册时间
2018-10-19
在线时间
2 小时
性别
保密
 楼主| 发表于 2018-11-27 16:22:50 | 显示全部楼层
回复

使用道具 举报

高级模式 评论
您需要登录后才可以回帖 登录 | 立即注册 微信登录